Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,093,040 G→A A215T (GCT→ACT)  hisC → histidinol‑phosphate aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,093,0400GA98.3% 181.5 / NA 60A215T (GCT→ACT) hisChistidinol‑phosphate aminotransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (31/27);  minor base C (0/1);  total (31/29)
Fisher's exact test for biased strand distribution p-value = 4.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTGCTGGAGTTAACCCGCGGTAAGGCGATTGTGGTTGCCGATGAAGCCTATATCGAGTTTTGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAGCTTTTGCTCTGG‑CGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGTGATCGCCCCCTACCCGCTCTCGACGCCGGTTGCCGACATTGCGGCCCAGGCGTTAAGCCCACAGGGAATC  >  NC_000913/2092907‑2093183
                                                                                                                                     |                                                                                                                                                
ccTGCTGGAGTTAACCCGCGGTAAGGCGATTGTGGTTGCCGATGAAGCCTATATCGAGTTTTGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGG‑Cgg                                                                                                                                  >  1:940560/1‑149 (MQ=255)
 ctgctgGAGTTAACCCGCGGTAAGGCGATTGTGGTTGCCGATGAAGCCTATATCGAGTTTTGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGG‑Cggg                                                                                                                                 <  2:784305/149‑1 (MQ=255)
     tgGAGTTAACCCGCGGTAAGGCGATTGTGGTTGCCGATGAAGCCTATATCGAGTTTTGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGG‑CGGGGCtt                                                                                                                             >  1:965624/1‑149 (MQ=255)
       gAGTTAACCCGCGGTAAGGCGATTGTGGTTGCCGATGAAGCCTATATCGAGTTTTGCCCGCAGGCATCGCTGGCCGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTATCGACAACTTTTCCTCTGG‑CGGGGCCTCg                                                                                                                           >  1:720869/1‑149 (MQ=255)
       gAGTTAACCCGCGGTAAGGCGAGTGTGGTTGCCGATGAAGCCTATATAGAGTTTTGCCCGCAGGCATCGCGGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTCTTTTACGCACCCTGCCGGAAACTTTGGCTCTGG‑GGGGGCTTCg                                                                                                                           >  2:470541/1‑149 (MQ=255)
         gTTAACCCGCCGTAACGCGATTGTGGTTGCCGATGAAGCCTATAACGAGTTTTGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGG‑CGGGGCTTCGtt                                                                                                                         <  1:632680/149‑1 (MQ=255)
             aCCCGCGGTAAGGCGATTGTGGTTGACGATGAAGCCTATATCGAGTTATGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTTGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGG‑CGGGGCTTCGTTGCgg                                                                                                                     >  1:612901/1‑149 (MQ=255)
               ccGCGGTAAGGCGATTGTGGTTGCCGATGAAGCCTATATCGAGTTTTGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGG‑CGGGGCTTCGTTGCGGAt                                                                                                                   >  2:1012095/1‑149 (MQ=255)
                      aaGGCGATTGTGGTTGCCGATGAAGCCTATATCGAGTTTTGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGG‑CGGGGCTTCGTTGCGGATTTACGCt                                                                                                            >  1:219635/1‑149 (MQ=255)
                      aaGGCGATTGTGGTTGCCGATGAAGCCTATATCGAGTTTTGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGG‑CGGGGCTTCGTTGCGGATTTACGCt                                                                                                            <  1:628074/149‑1 (MQ=255)
                       aGGCGATTGGGGTTGCCGATGAAGCCTATATCGAGTTTTGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGG‑CGGGGCTTCGTTGCGGATTTACGCTg                                                                                                           <  1:160870/149‑1 (MQ=255)
                         gCGATTGGGGTTGCCGATGAAGCCTATAACGCGTTTGGCCCGCGGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGG‑CGGGGCTTCGTTGCGGATTTACGCTGGc                                                                                                         >  2:461941/1‑149 (MQ=255)
                           gATTGTGCTTGCCGATGATGCCTATATCGAGTTTTGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGG‑CGGGGCTTCGTTGCGGATTTACGCTGGCaa                                                                                                       <  1:16428/149‑1 (MQ=255)
                                tgGTTGCCGATGAAGCCTACACCGAGTTTTGCCCGCAGGCATCGCTGCCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGG‑CGGGGCTTCGTTGCGGATTTACGCTGGCAAACgaa                                                                                                  <  2:259573/149‑1 (MQ=255)
                                  gTTGCCGATGAAGCCTATATCGAGTTTTGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGG‑CGGGGCTTCGTTGCGGATTTACGCTGGCAAACgaaga                                                                                                >  1:1026657/1‑149 (MQ=255)
                                   ttGCCGATGAAGCCTATATCGAGTTTTGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGG‑CGGGGCTTCGTTGCGGATTTACGCTGGCAAACgaagaa                                                                                               <  1:500153/149‑1 (MQ=255)
                                        gATGAAGCCTATATCGAG‑TTTGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGG‑CGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGtcatc                                                                                         <  1:341541/149‑1 (MQ=255)
                                               ccTATATCGAGTTTTGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGG‑CGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACctg                                                                                   <  1:671129/149‑1 (MQ=255)
                                                  atatCGAGTTTTGCCCGCAGGCATCTCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGG‑CGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACctgctg                                                                                <  1:1012095/149‑1 (MQ=255)
                                                  atatCGAGTTTTGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGGACTCTGGGTTG‑TGCTGATTTAAGCTGGCGGACGAAGAAGTCATCAACctgctg                                                                                >  1:790518/1‑149 (MQ=255)
                                                  atatCGAGTTTTGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGG‑CGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACctgctg                                                                                <  2:925666/149‑1 (MQ=255)
                                                   catCGAGTTTTGCCCGCGGGCACCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGG‑CGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCtga                                                                               <  2:592477/148‑1 (MQ=255)
                                                    atCGAGTTTTGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGG‑CGGGGCTTCGTTGCGGATTTACGCTGGCAAACGCAGGAGTCATCAACCTGCcgat                                                                              >  1:39229/1‑149 (MQ=255)
                                                        aGTTTTGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGG‑CGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGaaa                                                                          <  1:503219/149‑1 (MQ=255)
                                                               cccGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGG‑CGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGGCATCAACCTGCTGATGAAAGTGATCg                                                                   >  1:693347/1‑149 (MQ=255)
                                                                ccGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGG‑CGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGTGATCGc                                                                  >  2:789631/1‑149 (MQ=255)
                                                                 cGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGG‑CGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGTGATCGcc                                                                 >  2:779321/1‑149 (MQ=255)
                                                                  gCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGG‑CGGGGCTTCGTTGCGGAGTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGTGATCGccc                                                                >  1:846693/1‑149 (MQ=255)
                                                                  gCAGGCATCGCTGGCTGGCGGGCTGGCCGAACATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGG‑CGGGGCTTCGTTGCGGAGTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGTGATCGccc                                                                <  2:486280/149‑1 (MQ=255)
                                                                   cAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGG‑CGGGGCTTCGTTGCGGACTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGTGATCGcccc                                                               <  1:752651/149‑1 (MQ=255)
                                                                      gCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGG‑CGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGGTGAAAGTGATCGCCGCCTa                                                            >  2:215562/1‑149 (MQ=255)
                                                                        aTCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGG‑CGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGTGATCGCCCCCTAcc                                                          <  1:395931/149‑1 (MQ=255)
                                                                                 tggctggctggcGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGG‑CGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGGGATCGCCCCCTACCCGCTCCCGa                                                 >  2:692469/1‑149 (MQ=255)
                                                                                 tggctggccggcGGAATATCCGCACCTGGCTAATTTACGCACACTGTCGAAAACTTTTGCTCTGG‑CGGGGCTTCGTTGCGTATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGTGATCGCCCCCTACCCGCTCTCGa                                                 <  2:612709/149‑1 (MQ=255)
                                                                                     tggctggcGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGG‑CGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGTGATCGCCCCCTACCCGCTCTCGACGcc                                             >  2:241745/1‑149 (MQ=255)
                                                                                       gctggcGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGG‑CGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGTGATCGCCCCCTACCCGCTCTCGACGCCgg                                           >  1:510469/1‑149 (MQ=255)
                                                                                         tggcTGAATCTCCGCTCCTGCCTATTTTCCGCACCCTGGCGAAAACTTTCGCTCTGG‑CGGGGCTTCGTTCCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGTGATCGCCCCCTACCCGCTCTCGACGCCGGtt                                         <  1:129343/149‑1 (MQ=255)
                                                                                          ggcGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGT‑CGGGGCTTCGTTGCGGATTTACGCTGGCGAACGAAGAAGTCATCAACCTGCCGATGCAAGTGCTCGCCCCCTACCCGCTCTCGCCGCCTGTTg                                        >  2:3415/1‑149 (MQ=255)
                                                                                           gccgAATATCCGCACCTGGCTATCTTACGCACACTGTCGAAACCTTTTGCTCTGG‑CGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGTGATCGCCCCCTACCCGCTCTCGACGCCGGTTGc                                       <  2:184542/146‑1 (MQ=255)
                                                                                           gcGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGG‑CGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGTGATCGCCCCCTACCCGCTCTCGACGCCGGTTGc                                       >  1:918239/1‑149 (MQ=255)
                                                                                            cGGAATATCCGCACCTGGCTATTTTACGCCCACTGTCGAAAACTTTTGCTCTGG‑CGGGGCTTCGTTGCGGATTTACGCTGGCAAACGACGAAGTCAGCAACCTGCTGCTGAAAGTGATCGCCCACTACCCGCTCTCGCCGCCGGTCGcc                                      >  2:308983/1‑149 (MQ=255)
                                                                                              gAACATCCGCACCTGTCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGG‑CGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGCCATCAACCTGCTGCTGAAAGTGATCGCCCCCTACCCGCTCTCGACGCCGGTTGCCGa                                    <  2:254109/149‑1 (MQ=255)
                                                                                                atatCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGG‑CGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGTGATCGCCCCCTACCCGCTCTCGACGCCGGTTGCCGACa                                  >  1:1091268/1‑149 (MQ=255)
                                                                                                atatCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGG‑CGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGTGATCGCCCCCTACCCGCTCTCGACGCCGGTTGCCGACa                                  <  1:461941/149‑1 (MQ=255)
                                                                                                   tCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGG‑CGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGTGATCGCCCCCTACCCGCTCTCGACGCCGGTTGCCGACATTg                               <  2:588806/149‑1 (MQ=255)
                                                                                                     cGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGG‑CGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGTGATCGCCCCCTACCCGCTCTCGACGCCGGTTGCCGACATTGCg                             <  1:599851/149‑1 (MQ=255)
                                                                                                     cGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGG‑CGGGGCTTCGTTGCGGATGTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGTGATCGCCCCCTACCCGCTCTCGACGCCGGTTGCCGACATTGCg                             <  2:790518/149‑1 (MQ=255)
                                                                                                      gCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGG‑CGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGTGATCGCCCCCTCCCCGCACTCTACGCCGGTTGCCGACATTGCgg                            >  2:226885/1‑149 (MQ=255)
                                                                                                          cTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGG‑CGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGTGATCGCCCCCTACCCGCTCTCGACGCCGGTTGCCGACATTGCGGCCCa                        <  2:510469/149‑1 (MQ=255)
                                                                                                            ggCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGG‑CGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGTGATCGCCCCCTACCCGCTCTCGACGCCGGTTGCCGACATTGCGGCCCAgg                      <  2:792827/149‑1 (MQ=255)
                                                                                                             gCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGG‑CGGGGCTTCGTTGCGGATGTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGTGATCGCCCCCTACCCGCTCTCGACGCCGGTTGCCGACATTGCGGCCCAGGc                     >  1:310984/1‑149 (MQ=255)
                                                                                                                aTTTTACGCACACTGTCGAAAACTTTTGCTCTGG‑CGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGTGATCGCCCCCTACCCGCTCTCGACGCCGGTTGCCGACATTGCGGCCCAGGCGtt                  >  1:541904/1‑149 (MQ=255)
                                                                                                                aTTTTACGCACACTGTCGAAAACTTTTGCTCTGG‑CGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGTGATCGCCCCCTACCCGCTCTCGACGCCGGTTGCCGACATTGCGGCCCAGGCGtt                  >  1:577314/1‑149 (MQ=255)
                                                                                                                 ttttACGCACACTGTCGAAAACTTTTGCTCTGG‑CGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGTGATCGCCCCCTACCCGCTCTCGACGCCGGTTGCCGACATTGCGGCCCAGGCGTTa                 >  1:714600/1‑149 (MQ=255)
                                                                                                                 ttttACGCACACTGTCGAAAACTTTTGCTCTGG‑CGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGTGATCGCCCCCTACCCGCTCTCGACGCCGGTTGCCGACATTGCGGCCCAGGCGTTa                 <  1:271853/149‑1 (MQ=255)
                                                                                                                    tACGCACACTGTCGAAAACTTGTGCTCTGG‑CGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGTGATCGCCCCCTACCCGCTCTCGACGCCGGTTGCCGACATTGCGGCCCAGGCGTTAAGc              <  1:214858/149‑1 (MQ=255)
                                                                                                                         acagTGTCTATATCTTTGGCTCTGG‑CGGGGCTTCGTTTCGGATTTACGCTGTTTAACAAAGAAGTCATCAATCTGCTTTTGAGCGTGGTCGCCCCCTGCCCGCTCTCGACGTCGGTTGCCGACATTGCGGCTCTGTCGTGAAGCCcaca         <  2:467519/149‑1 (MQ=255)
                                                                                                                             tGTCGAAAACTTTTGCTCTGG‑CGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGTGATCGCCCCCTACCCGCTCTCGCCGCCGGTTGCCGACATTGCGGCCCAGGCGTTAAGCCCACAGGGa     >  1:26320/1‑149 (MQ=255)
                                                                                                                               tCGAAAACTTTTGCTCTGG‑CGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGTGATCGCCCCCTACCCGCTCTCGACGCCGGTTGCCGACATTGCGGCCCAGGCGTTAAGCCCACAGGGAAt   >  1:949834/1‑149 (MQ=255)
                                                                                                                               tCGAAAACTTTTGCTCTGG‑CGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGTGATCGCCCCCTACCCGCTCTCGACGCCGGTTGCCGACATTGCGGCCCAGGCGTTAAGCCCACAGGGAAt   >  1:693567/1‑149 (MQ=255)
                                                                                                                                caaaaaCTTTTTCCCGGG‑CGGGCCTTCGTTGCGGATTTACGCTGGCAACCGAAGAAGTCATCAACCTGCTGATGAAAGTGATCGCCCCCTACCCGCTCTCGACGCCGGTTGCCGACATTGCGGCCCAGGCGTTAAGCCCACAGGGAAtc  <  1:71164/147‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                                
CCTGCTGGAGTTAACCCGCGGTAAGGCGATTGTGGTTGCCGATGAAGCCTATATCGAGTTTTGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAGCTTTTGCTCTGG‑CGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGTGATCGCCCCCTACCCGCTCTCGACGCCGGTTGCCGACATTGCGGCCCAGGCGTTAAGCCCACAGGGAATC  >  NC_000913/2092907‑2093183

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑