Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,942,835 G→A 6.4% P55P (CCG→CCA)  PP_5210 → alcohol dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,942,8350GA 6.4% 181.2 / 2.1 62P55P (CCG→CCA) PP_5210alcohol dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (34/24);  new base A (2/2);  total (36/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.97e-01

CCTGGGTCCGGCGCGCAACCTGGTGCTGGAAGAGGTGGCCAGCCCACTGCCGAAGAAGAACGAGATCCTGCT‑GGACGTGCAAGCTGCCGGGGTCAACTTCCCCGACACCTTGATCATCGAAGGCAAGTACCAGTTCCAGCCGCCGCTGCCGTTCTCTCCAGGTGGCGAGGCAGCGGGCGTGGTGGCCGCCGTCGGCGAAAAGGCCGGCGCGTTCAAGGTCGGCGACCGGGTCATGGCGCTCACCGGCTGGGGTGCGTTCGCCGAGCAAGTGGCGGTGCCGTTCTATAAC  >  NC_002947/5942694‑5942982
                                                                                                                                              |                                                                                                                                                   
ccTGGGTCCGGCGCGCAACCTGGTGCTGGAAGAGGTGGCCAGCCCACTGCCGAAGAAGAACGAGATCCTGCT‑GGACGTGCAAGCTGCCGGGGTCAACTTCCCCGACACCTTGATCATCGAAGGCAAGTACCAGTTCCAGCCGCCGCTGcc                                                                                                                                             >  1:1123169/1‑150 (MQ=255)
ccTGGGTCCGGCGCGCAACCTGGTGCTGGAAGAGGTGGCCAGCCCACTGCCGAAGAAGAACGAGATCCTGCT‑GGACGTGCAAGCTGCCGGGGTCAACTTCCCCGACACCTTGATCATCGAAGGCAAGTACCAGTTCCAGCCGCCGCTGcc                                                                                                                                             >  1:1626952/1‑150 (MQ=255)
         ggCGCGCAACCTGGTGCTGGAAGAGGTGGCCAGCCCACTGCACTGGAGGTCAGACGT‑GTGCTCTTCCGATCAAGCTGCCGGGGTCAACTTCCCCGACACCTTGATCATCGAAGGCAAGTACCAGTTCCAGCCGCCGCTGCCGTTCTCTcc                                                                                                                                    <  1:717151/150‑1 (MQ=255)
             cgcAACCTGGTGCTGGAAGAGGTGGCCAGCCCACTGCCGAAGAAGAACGAGATCCTGCT‑GGACGTGCAAGCTGCCGGGGTCAACTTCCCCGACACCTTGATCATCGAAGGCAAGTACCAGTTCCAGCCGCCGCTGCCGTTCTCTCCAGGt                                                                                                                                <  2:466022/150‑1 (MQ=255)
                    tGGTGCTGGAAGAGGTGGCCAGCCCACTGCCGAAGAAGAACGAGATCCTGCT‑GGACGTGCAAGCTGCCGGGGTCAACTTCCCCGACACCTTGATCATCGAAGGCAAGTACCAGTTCCAGCCGCCGCTGCCGTTCTCTCCAGGTGGCGAgg                                                                                                                         >  1:263925/1‑150 (MQ=255)
                    tGGTGCTGGAAGAGGTGGCCAGCCCACTGCCGAAGAAGAACGAGATCCTGCT‑GGACGTGCAAGCTGCCGGGGTCAACTTCCCCGACACCTTGATCATCGAAGGCAAGTACCAGTTCCAGCCGCCGCTGCCGTTCTCTCCAGGTGGCGAgg                                                                                                                         <  2:341047/150‑1 (MQ=255)
                    tGGTGCTGGAAGAGGTGGCCAGCCCACTGCCGAAGAAGAACGAGATCCTGCT‑GGACGTGCAAGCTGCCGGGGTCAACTTCCCCGACACCTTGATCATCGAAGGCAAGTACCAGTTCCAGCCGCCGCTGCCGTTCTCTCCAGGTGGCGAgg                                                                                                                         <  1:1818235/150‑1 (MQ=255)
                          tGGAAGAGGTGGCCAGCCCACTGCCGAAGAAGAACGAGATCCTGCT‑GGACGTGCAAGCTGCCGGGGTCAACTTCCCCGACACCTTGATCATCGAAGGCAAGTACCAGTTCCAGCCGCCGCTGCCGTTCTCTCCAGGTGGCGAGGCAGCgg                                                                                                                   >  2:1211246/1‑150 (MQ=255)
                          tGGAAGAGGTGGCCAGCCCACTGCCGAAGAAGAACGAGATCCTGCT‑GGACGTGCAAGCTGCCGGGGTCAACTTCCCCGACACCTTGATCATCGAAGGCAAGTACCAGTTCCAGCCGCCGCTGCCGTCCTCTCCAGGTGGCGAGGCCGCgg                                                                                                                   >  1:294941/1‑150 (MQ=255)
                           ggAAGAGGTGGCCAGCCCACTGCCGAAGAAGAACGAGATCCTGCT‑GGACGTGCAAGCTGCCGGGGTCAACTTCCCCGACACCTTGATCATCGAAGGCAAGTACCAGTTCCAGCCGCCGCTGCCGTTCTCTCCAGGTGGCGAGGCAGCggg                                                                                                                  <  2:1548972/150‑1 (MQ=255)
                            gAAGAGGTGGCCAGCCCACTGCCGAAGAAGAACGAGATCCTGCT‑GGACGTGCAAGCTGCCGGGGTCAACTTCCCCGACACCTTGATCATCGAAGGCAAGTACCAGTTCCAGCCGCCGCTGCCGTTCTCTCCAGGTGGCGAGGCAGCGGGc                                                                                                                 >  2:1917156/1‑150 (MQ=255)
                                 ggTGGCCAGCCCACTGCCGAAGAAGAACGAGATCCTGCT‑GGACGTGCAAGCTGCCGGGGTCAACTTCCCCGACACCTTGATCATCGAAGGCAAGTACCAGTTCCAGCCGCCGCTGCCGTTCTCTCCAGGTGGCGAGGCAGCGGGCgtggt                                                                                                            >  1:1789435/1‑150 (MQ=255)
                                        aGCCCACTGCCGAAGAAGAACGAGATCCTGCT‑GGACGTGCAAGCTGCCGGGGTCAACTTCCCCGACACCTTGATCATCGAAGGCAAGTACCAGTTCCAGCCGCCGCTGCCGTTCTCTCCAGGTGGCGAGGCAGCGGGCGTGGTGgccgcc                                                                                                     >  1:1932402/1‑150 (MQ=255)
                                           gcACTTCGGAAGAAGAACGAGATCCTGCT‑GGACGTGCAAGCTGCCGGGGTCAACTTCCCCGACACCTTGATCATCGAAGGCAAGTACCAGTTCCAGCCGCCGCTGCCGTTCTCTCCAGGTGGCGAGGCAGCGGGCGTGGTGGCCGCCGTc                                                                                                  >  2:1530531/2‑150 (MQ=255)
                                            cACTGCCGAAGAAGAACGAGATCCTGCT‑GGACGTGCAAGCTGCCGGGGTCAACTTCCCCGACACCTTGATCATCGAAGGCAAGTACCAGTTCCAGCCGCCGCTGCCGTTCTCTCCAGGTGGCGAGGCAGCGGGCGTGGTGGCCGCCGTCg                                                                                                 >  2:340200/1‑150 (MQ=255)
                                             aCTGCCGAAGAAGAACGAGATCCTGCT‑GGACGTGCAAGCTGCCGGGGTCAACTTCCCCGACACCTTGATCATCGAAGGCAAGTACCAGTTCCAGCCACCGCTGCCGTTCTCTCCAGGTGGCGAGGCAGCGGGCGTGGTGGCCGCCGTCgg                                                                                                >  1:269648/1‑150 (MQ=255)
                                             aCTGCCGAAGAAGAACGAGATCCTGCT‑GGACGTGCAAGCTGCCGGGGTCAACTTCCCCGACACCTTGATCATCGAAGGCAAGTACCAGTTCCAGCCACCGCTGCCGTTCTCTCCAGGTGGCGAGGCAGCGGGCGTGGTGGCCGCCGTCgg                                                                                                >  1:269688/1‑150 (MQ=255)
                                              cTGCCGAAGAAGAACGAGATCCTGCT‑GGACGTGCAAGCTGCCGGGGTCAACTTCCCCGACACCTTGATCATCGAAGGCAAGTACCAGTTCCAGCCGCCGCTGCCGTTCTCTCCAGGTGGCGAGGCAGCGGGCGTGGTGGCCGCCGTCGGc                                                                                               <  1:1356070/150‑1 (MQ=255)
                                                   gaagaagaaCGAGATCCTGCT‑GGACGTGCAAGCTGCCGGGGTCAACTTCCCCGACACCTTGATCATCGAAGGCAAGTACCAGTTCCAGCCGCCGCTGCCGTTCTCTCCAGGTGGCGAGGCAGCGGGCGTGGTGGCCGCCGTCGGCGaaaa                                                                                          >  2:386151/1‑150 (MQ=255)
                                                    aagaagaaCGAGATCCTGCT‑GGACGTGCAAGCTGCCGGGGTCAACTTCCCCGACACCTTGATCATCGAAGGCAAGTACCAGTTCCAGCCGCCGCTGCCGTTCTCTCCAGGTGGCGAGGCAGCGGGCGTGGTGGCCGCCGTCGGCGAAAAg                                                                                         >  1:686379/1‑150 (MQ=255)
                                                         gaaCGAGATCCTGCT‑GGACGTGCAAGCTGCCGGGGTCAACTTCCCCGACACCTTGATCATCGAAGGCAAGTACCAGTTCCAGCCGCCGCTGCCGTTCTCTCCAGGTGGCGAGGCAGCGGGCGTGGTGGCCGCCGTCGGCGAAAAGGCCgg                                                                                    >  1:1234553/1‑150 (MQ=255)
                                                                    tgct‑gGACGTGCAAGCTGCCGGGGTCAACTTCCCCGACACCTTGATCATCGAAGGCAAGTACCAGTTCCAACCGCCGCTGCCGTTCTCTCCAGGTGGCGAGGCAGCGGGCGTGGTGGCCGCCGTCGGCGAAAAGGCCGGCGCGTTCAAgg                                                                         >  2:880578/1‑150 (MQ=255)
                                                                      ct‑gGACGTGCAAGCTGCCGGGGTCAACTTCCCCGACACCTTGATCATCGAAGGCAAGTACCAGTTCCAGCCGCCGCTGCCGTTCTCTCCAGGGGGCGAGGCAGCGGGCGTGGTGGCCGCCGTCGGCGAAAAGGCCGGCGCGTTCAAGGTc                                                                       <  2:263925/150‑1 (MQ=255)
                                                                       t‑gGACGTGCAAGCTGCCGGGGTCAACTTCCCCGACACCTTGATCATCGAAGGCAAGTACCAGTTCCAGCCGCCGCTGCCGTTCTCTCCAGGACGCCATGCAGCGGGCGTGGTGGCCGCCGTCGGCGAAAAGGCCGGCGCGTTCAAGGTCg                                                                      <  2:1626952/150‑1 (MQ=255)
                                                                             gTGCAAGCTGCCGGGGTCAACTTCCCCGACACCTTGATCATCGAAGGCAAGTACCAGTTCCAGCCGCCGCTGCCGTTCTCTCCAGGTGGCGAGGCAGCGGGCGTGGTGGCCGCCGTCGGCGAAAAGGCCGGCGCGTTc                                                                             <  2:701998/138‑1 (MQ=255)
                                                                             gTGCAAGCTGCCGGGGTCAACTTCCCCGACACCTTGATCATCGAAGGCAAGTACCAGTTCCAGCCGCCGCTGCCGTTCTCTCCAGGTGGCGAGGCAGCGGGCGTGGTGGCCGCCGTCGGCGAAAAGGCCGGCGCGTTc                                                                             >  1:701998/1‑138 (MQ=255)
                                                                             gTGCAAGCTGCCGGGGTCAACTTCCCCGACACCTTGATCATCGAAGGCAAGTACCAGTTCCAGCCGCCGCTGCCGTTCTCTCCAGGTGGCGAGGCAGCGGGCGTGGTGGCCGCCGTCGGCGAAAAGGCCGGCGCGTTCAAGGTCGGCGAc                                                                 >  1:1405486/1‑150 (MQ=255)
                                                                              tGCAAGCTGCCGGGGTCAACTTCCCCGACACCTTGATCATCGAAGGCAAGTACCAGTTCCAGCCGCCGCTGCCGTTCTCTCCAGGTGGCGAGGCAGCGGGCGTGGTGGCCGCCGTCGGCGAAAAGGCCGGCGCGTTCAAGGTCGGCGAcc                                                                >  1:625572/1‑150 (MQ=255)
                                                                                     tGCCGGGGTCAACTTCCCCGACACCTTGATCATCGAAGGCAAGTACCAGTTCCAGCCGCCGCTGCCGTTCTCTCCAGGTGGCGAGGCAGCGGGCGTGGTGGCCGCCGTCGGCGAAAAGGCCGGCGCGTTCAAGGTCGGCGACCGGGTCAt                                                         >  2:238084/1‑150 (MQ=255)
                                                                                       ccGGGGTCAACTTCCCCGACACCTTGATCATCGAAGGCAAGTACCAGTTCCAGCCGCCGCTGCCGTTCTCTCCAGGTGGCGAGGCAGCGGGCGTGGTGGCCGCCGTCGGCGAAAAGGCCGGCGCGTTCAAGGTCGGCGACCGGGTCATgg                                                       <  2:1405486/150‑1 (MQ=255)
                                                                                             tCAACTTCCCCGACACCTTGATCATCGAAGGCAAGTACCAGTTCCAGCCGCCGCTGCCGTTCTCTCCAGGTGGCGAGGCAGCGGGCGTGGTGGCCGCCGTCGGCGAAAAGGCCGGCGCGTTCAAGGTCGGCGACCGGGTCATGTCGCTCa                                                 >  2:541651/1‑150 (MQ=255)
                                                                                              agttgaccccGACACCTTGATCATCGAAGGCAAGTACCAGTTCCAGCCGCCGCTGCCGTTCTCTCCAGGTGGCGAGGCAGCGGGCGTGGTGGCCGCCGTCGGCGAAAAGGCCGGCGCGTTCAAGGTCGGCGACCGGGTCATGGCGCTCAc                                                >  2:1370547/7‑150 (MQ=255)
                                                                                                aCTTCCCCGACACCTTGATCATCGAAGGCAAGTACCAGTTCCAGCCGCCGCTGCCGTTCTCTCCAGGTGGCGAGGCAGCGGGCGTGGTGGCCGCCGTCGGCGAAAAg                                                                                         <  1:803479/107‑1 (MQ=255)
                                                                                                aCTTCCCCGACACCTTGATCATCGAAGGCAAGTACCAGTTCCAGCCGCCGCTGCCGTTCTCTCCAGGTGGCGAGGCAGCGGGCGTGGTGGCCGCCGTCGGCGAAAAg                                                                                         >  2:803479/1‑107 (MQ=255)
                                                                                                aCTTCCCCGACACCTTGATCATCGAAGGCAAGTACCAGTTCCAGCCGCCGCTGCCGTTCTCTCCAGGTGGCGAGGCAGCGGGCGTGGTGGCCGCCGTCGGCGAAAAGGCCGGCGCGTTCAAGGTCGGCGACCg                                                               <  1:1162435/133‑1 (MQ=255)
                                                                                                aCTTCCCCGACACCTTGATCATCGAAGGCAAGTACCAGTTCCAGCCGCCGCTGCCGTTCTCTCCAGGTGGCGAGGCAGCGGGCGTGGTGGCCGCCGTCGGCGAAAAGGCCGGCGCGTTCAAGGTCGGCGACCg                                                               >  2:1162435/1‑133 (MQ=255)
                                                                                                     ctcGACACCTTGATCATCGAAGGCAAGTACCAGTTCCAGCCGCCGCTGCCGTTCTCTCCAGGTGGCGAGGCAGCGGGCGTGGTGGCCGCCGTCGGCGAAAAGGCCGGCGCGTTCAAGGTCGGCGACCGGGTCATGGCGCTCACCGGCTgg                                         >  2:378030/3‑150 (MQ=255)
                                                                                                        gACACCTTGATCATCGAAGGCAAGTACCAGTTCCAGCCGCCGCTGCCGTTCTCTCCAGGTGGCGAGGCAGCGGGCGTGGTGGCCGCCGTCGGCGAAAAGGCCGGCGCGTTCAAGGTCGGCGACCGGGTCATGGCGCTCa                                                 >  1:1642213/1‑139 (MQ=255)
                                                                                                        gACACCTTGATCATCGAAGGCAAGTACCAGTTCCAGCCGCCGCTGCCGTTCTCTCCAGGTGGCGAGGCAGCGGGCGTGGTGGCCGCCGTCGGCGAAAAGGCCGGCGCGTTCAAGGTCGGCGACCGGGTCATGGCGCTCa                                                 <  2:1642213/139‑1 (MQ=255)
                                                                                                          cacCTTGATCATCGAAGGCAAGTACCAGTTCCAGCCGCCGCTGCCGTTCTCTCCAGGTGGCGAGGCAGCGGGCGTGGTGGCCGCCGTCGGCGAAAAGGCCGGCGCGTTCAAGGTCGGCGACCGGGTCATGGCGCTCACCGGCTGGGGTGc                                    <  2:1789435/150‑1 (MQ=255)
                                                                                                              ttGATCATCGAAGGCAAGTACCAGTTCCAGCCGCCGCTGCCGTTCTCTCCAGGTGGCGAGGCAGCGGGCGTGGTGGCCGCCGTCGGCGAAAAGGCCGGCGCGTTCAAGGTCGGCGACCGGGTCATGGCGCTCACCGGCTGGGGTGCGTTc                                >  1:1936337/1‑150 (MQ=255)
                                                                                                              ttGATCATCGAAGGCAAGTACCAGTTCCAGCCGCCGCTGCCGTTCTCTCCAGGTGGCGAGGCAGCGGGCGTGGTGGCCGCCGTCGGCGAAAAGGCCGGCGCGTTCAAGGTCGGCGACCGGGTCATGGCGCTCACCGGCTGGGGTGCGTTc                                >  2:715522/1‑150 (MQ=255)
                                                                                                                gATCATCGAAGGCAAGTACCAGTTCCAGCCGCCGCTGCCGTTCTCTCCAGGTGGCGAGGCAGCGGGCGTGGTGGCCGCCGTCGGCGAAAAGGCCGGCGCGTTCAAGGTCGGCGACCGGGTCATGGCGCTCACCGGCTGGGGTGCGTTCGc                              <  1:715522/150‑1 (MQ=255)
                                                                                                                gATCATCGAAGGCAAGTACCAGTTCCAGCCACCGCTGCCGTTCTCTCCAGGTGGCGAGGCAGCGGGCGTGGTGGCCGCCGTCGGCGAAAAGGCCGGCGCGTTCAAGGTCGGCGACCGGGTCATGGCGCTCACCGGCTGGGGTGCGTTCGc                              <  2:269688/150‑1 (MQ=255)
                                                                                                                gATCATCGAAGGCAAGTACCAGTTCCAGCCACCGCTGCCGTTCTCTCCAGGTGGCGAGGCAGCGGGCGTGGTGGCCGCCGTCGGCGAAAAGGCCGGCGCGTTCAAGGTCGGCGACCGGGTCATGGCGCTCACCGGCTGGGGTGCGTTCGc                              <  2:269648/150‑1 (MQ=255)
                                                                                                                 atcatcGAAGGCAAGTACCAGTTCCAGCCGCCGCTGCCGTTCTCTCCAGGTGGCGAGGCAGCGGGCGTGGTGGCCGCCGTCGGCGAAAAGGCCGGCGCGTTCAAGGTCGGCGACCGGGTCATGGCGCTCACCGGCTGGGGTGCGTTCGcc                             <  2:686379/150‑1 (MQ=255)
                                                                                                                   catcGAAGGCAAGTACCAGTTCCAGCCGCCGCTGCCGTTCTCTCCAGGTGGCGAGGCAGCGGGCGTGGTGGCCGCCGTCGGCGAAAAGGCCGGCGCGTTCAAGGTCGGCGACCGGGTCATGGCGCTCACCGGCTGGGGTGCGTTCGCCGa                           <  1:1371653/150‑1 (MQ=255)
                                                                                                                    atcGAAGGCAAGTACCAGTTCCAGCCGCCGCTGCCGTTCTCTCCAGGTGGCGAGGCAGCGGGCGTGGTGGCCGCCGTCGGCGAAAAGGCCGGCGCGTTCa                                                                            >  2:208607/1‑100 (MQ=255)
                                                                                                                    atcGAAGGCAAGTACCAGTTCCAGCCGCCGCTGCCGTTCTCTCCAGGTGGCGAGGCAGCGGGCGTGGTGGCCGCCGTCGGCGAAAAGGCCGGCGCGTTCa                                                                            <  1:208607/100‑1 (MQ=255)
                                                                                                                     tcGAAGGCAAGTACCAGTTCCAGCCGCCGCTGCCGTTCTCTCCAGGTGGCGAGGCAGCGGGCGTGGTGGCCGCCGTCGGCGAAAAGGCCGGCGCGTTCAAGGTCGGCGACCGGGTCATGGCGCTCACCGGCTGGGGTGCGTTCGcc                             <  1:1512074/146‑1 (MQ=255)
                                                                                                                     tcGAAGGCAAGTACCAGTTCCAGCCGCCGCTGCCGTTCTCTCCAGGTGGCGAGGCAGCGGGCGTGGTGGCCGCCGTCGGCGAAAAGGCCGGCGCGTTCAAGGTCGGCGACCGGGTCATGGCGCTCACCGGCTGGGGTGCGTTCGcc                             >  2:1512074/1‑146 (MQ=255)
                                                                                                                     tcGAAGGCAAGTACCAGTTCCAGCCGCCGCTGCCGTTCTCTCCAGGTGGCGAGGCAGCGGGCGTGGTGGCCGCCGTCGGCGAAAAGGCCGGCGCGTTCAAGGTCGGCGACCGGGTCATGGCGCTCACCGGCTGGGGTGCGTTCGCCGAGc                         >  2:1620586/1‑150 (MQ=255)
                                                                                                                       gAAGGCAAGTACCAGTTCCAGCCGCCGCTGCCGTTCTCTCCAGGTGGCGAGGCAGCGGGCGTGGTGGCCGCCGTCGGCGAAAAGGCCGGCGCGTTCAAGGTCGGCGACCGGGTCATGGCGCTCACCGGCTGGGGTGCGTTCGCCGAGCaa                       <  2:1936337/150‑1 (MQ=255)
                                                                                                                       gAAGGCAAGTACCAGTTCCAGCCGCCGCTGCCGTTCTCTCCAGGTGGCGAGGCAGCGGGCGTGGTGGCCGCCGTCGGCGAAAAGGCCGGCGCGTTCAAGGTCGGCGACCGGGTCATGGCGCTCACCGGCTGGGGTGCGTTCGCCGAGCaa                       >  1:806639/1‑150 (MQ=255)
                                                                                                                         aGGCAAGTACCAGTTCCAGCCGCCGCTGCCGTTCTCTCCAGGTGGCGAGGCAGCGGGCGTGGTGGCCGCCGTCGGCGAAAAGGCCGGCGCGTTCAAGGTCGGCGACCGGGTCATGGCGCTCACCGGCTGGGGTGCGTTCGCCGAGCAAGt                     >  2:250889/1‑150 (MQ=255)
                                                                                                                             aaGTACCAGTTCCAGCCGCCGCTGCCGTTCTCTCCAGGTGGCGAGGCAGCGGGCGTGGTGGCCGCCGTCGGCGAAAAGGCCGGCGCGTTCAAGGTCGGCGACCGGGTCATGGCGCTCACCGGCTGGGGTGCGTTCGCCGAGCAAGTGGCg                 <  1:340200/150‑1 (MQ=255)
                                                                                                                                tACCAGTTCCAGCCGCCGCTGCCGTTCTCTCCAGGTGGCGAGGCAGCGGGCGTGGTGGCCGCCGTCGGCGAAAAGGCCGGCGCGTTCAAGGTCGGCGACCGGGTCATGGCGCTCACCGGCTGGGGTGCGTTCGCCGAGCAAGTGGCGGTg              >  2:1083260/1‑150 (MQ=255)
                                                                                                                                     gTTCCAGCCGCCGCTGCCGTTCTCTCCAGGTGGCGAGGCAGCGGGCGTGGTGGCCGCCGTCGGCGAAAAGGCCGGCGCGTTCAAGGTCGGCGACCGGGTCATGGCGCTCACCGGCTGGGGTGCGTTCGCCGAGCAAGTGGCGGTGCCGtt         <  1:378030/150‑1 (MQ=255)
                                                                                                                                      ttCCAGCCGCCGCTGCCGTTCTCTCCAGGTGGCGAGGCAGCGGGCGTGGTGGCCGCCGTCGGCGAAAAGGCCGGCGCGTTCAAGGTCGGCGACCGGGTCATGGCGCTCACCGGCTGGGGTGCGTTCGCCGAGCAAGTGGCGGTGCCGTTc        >  1:436645/1‑150 (MQ=255)
                                                                                                                                         cAGCCGCCGCTGCCGTTCTCTCCAGGTGGCGAGGCAGCGGGCGTGGTGGCCGCCGTCGGCGAAAAGGCCGGCGCGTTCAAGGTCGGCGACCGGGTCATGGCGCTCACCGGCTGGGGTGCGTTCGCCGAGCAAGTGGCGGTGCCGTTCtat     <  1:1370547/150‑1 (MQ=255)
                                                                                                                                          aGCCGCCGCTGCCGTTCTCTCCAGGTGTCGAGGCAGCGGGCGTGGTGGCCGCCGTCGGCGAAAAGGCCGGCGCGTTCAAGGTCGGCGACCGGGTCATGGCGCTCACCGGCTGGGGTGCGTTCGCCGATCAAGTGGCGGTGCCGTTCtata    <  2:740977/150‑1 (MQ=255)
                                                                                                                                          aGCCGCCGCTGCCGTTCTCTCCAGGTGGCGAGGCAGCGGGCGTGGTGGCCGCCGTCGGCGAAAAGGCCGGCGCGTTCAAGGTCGGCGCCTGGGTCATGGCGCTCACCCGGTGGGGTGCGTTCGCCGAGCAAGTGGCGGTGTCGTTCtata    >  1:786072/1‑150 (MQ=255)
                                                                                                                                           gccgccGCTGCCGTTCTCTCCAGGTGGCGAGGCAGCGGGCGTGGTGGCCGCCGTCGGCGAAAAGGCCGGCGCGTTCAAGGTCGGCGACCGGGTCATGGCGCTCACCGGCTGGGGTGCGTTCGCCGAGCAAGTGGCGGTGCCGTTCTATaa   <  1:238084/150‑1 (MQ=255)
                                                                                                                                            ccgccgCTGCCGTTCTCTCCAGGTGGCGAGGCAGCGGGCGTGGTGGCCGCCGTCGGCGAAAAGGCCGGCGCGTTCAAGGTCGGCGACCGGGTCATGGCGCTCACCGGCTGGGGTGCGTTCGCCGAGCAAGTGGCGGTGCCGTTCTATAAc  <  1:541651/150‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                   
CCTGGGTCCGGCGCGCAACCTGGTGCTGGAAGAGGTGGCCAGCCCACTGCCGAAGAAGAACGAGATCCTGCT‑GGACGTGCAAGCTGCCGGGGTCAACTTCCCCGACACCTTGATCATCGAAGGCAAGTACCAGTTCCAGCCGCCGCTGCCGTTCTCTCCAGGTGGCGAGGCAGCGGGCGTGGTGGCCGCCGTCGGCGAAAAGGCCGGCGCGTTCAAGGTCGGCGACCGGGTCATGGCGCTCACCGGCTGGGGTGCGTTCGCCGAGCAAGTGGCGGTGCCGTTCTATAAC  >  NC_002947/5942694‑5942982

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