Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,053,758 T→A 6.6% W16R (TGG→AGG)  exbB → biopolymer transport protein ExbB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,053,7580TA 6.6% 381.2 / 11.4 120W16R (TGG→AGG) exbBbiopolymer transport protein ExbB
Reads supporting (aligned to +/- strand):  ref base T (67/45);  new base A (5/3);  total (72/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

Only 100 of 121 total aligned reads displayed.

TTTTGTTAACAAGAATTACTTGCAATAATGACCGCCAAATTGTTCTCGGCCTACATCGAGGCCTTGGAGAACGCTCACTTTTTTCTTCATCAGGTCCGGCCAGCATGACACGTACTCAACCTTCCGCTTCGCCAACCCCGTCGCGCGCCTGGCGCGCCATCGCCGCGCTGACACTCAGCCTGGTGCTGGCCCCGGTGGCCATGGCCGATGAGCCAACCACCAACGCCGCCACACCTGCTGCTGCCACTGCCCCGGCCACACCTGCTGCTGCCCCGGCGGATGCTCCGGCGCCTGAT  >  NC_002947/6053609‑6053904
                                                                                                                                                     |                                                                                                                                                  
ttttGTTAACAAGAATTACTTGCAATAATGACCGCCAAATTGTTCTCGGCCTACATCGAGGCCTTGGAGAACGCTCACTTTTTTCTTCATCAGGTCTGGCCAGCATGACACGTACTCAACCTTCCGCTTCGCCAAACCCGTCGCGCCCCt                                                                                                                                                    >  2:1379507/1‑150 (MQ=255)
        aCAAGAATTACTTGCAATAATGACCGCCAAATTGTTCTCGGCCTACATCGAGGCCTTGGAGAACGCTCACTTTTTTCTTCATCAGGTCCGGCCAGCATGACACGTACTCAACCTTCCGCTTCGCCAACCCCGTCGCGCGCCTGGCGCGcc                                                                                                                                            >  1:396289/1‑150 (MQ=255)
        aCAAGAATTACTTGCAATAATGACCGCCAAATTGTTCTCGGCCTACATCGAGGCCTTGGAGAACGCTCACTTTTTTCTTCATCAGGTCCGGCCAGCATGACACGTACTCAACCTTCCGCTTCGCCAACCCCGTCGCGCGCCTGGCGCGcc                                                                                                                                            >  1:75604/1‑150 (MQ=255)
         cAAGAATTAGTTTCAATAATGACCGCCAAATTGTTCTCGGCCTACATCGAGGCCTTGGAGAACGCTCACTTTTTTCTTCATCAGGTCCGGCCAGCATGACACGTACTCAACCTTCCGCTTCGCCAACCCCGTCGCGCGCCAGgcgcgcga                                                                                                                                           <  1:1447795/150‑3 (MQ=255)
            gAATTACTTGCAATAATGACCGCCAAATTGTTCTCGGCCTACATCGAGGCCTTGGAGAACGCTCACTTTTTTCTTCATCAGGTCCGGCCAGCATGACACGTACTCAACCTTCCGCTTCGCCAACCCCGTCGCGCGCCAGgcgcgcgacgg                                                                                                                                        <  1:1030611/150‑6 (MQ=255)
             aaTTACTTGCAATAATGACCGCCAAATTGTTCTCGGCCTACATCGAGGCCTTGGAGAACGCTCACTTTTTTCTTCATCAGGTCCGGCCAGCATGACACGTACTCAACCTTCCGCTTCGCCAACCCCGTCGCGCGCCTGGCGCGCCATcgc                                                                                                                                       >  2:951844/1‑150 (MQ=255)
              aTTACTTGCAATAATGACCGCCAAATTGTTCTCGGCCTACATCGAGGCCTTGGAGAACGCTCACTTTTTTCTTCATCAGGTCCGGCCAGCATGACACGTACTCAACCTTCCGCTTCGCCAACCCCGTGGCGCGCCAGgcgcgcgacgggg                                                                                                                                      <  2:322063/150‑8 (MQ=255)
              aTTACTTGCAATAATGACCGCCAAATTGTTCTCGGCCTACATCGAGGCCTTGGAGAACGCTCACTTTTTTCTTCATCAGGTCCGGCCAGCATGACACGTACTCAACCTTCCGCTTCGCCAACCCCGTCGCGCGCCTGGCGCGCCATcgcc                                                                                                                                      <  2:562504/150‑1 (MQ=255)
                tACTTGCAATAATGACCGCCAAATTGTTCTCGGCCTACATCGAGGCCTTGGAGAACGCTCACTTTTTTCTTCATCAGGTCCGGCCAGCATGACACGTACTCAACCTTCCGCTTCGCCAACCCCGTCGCGCGCCTGGCGCGCCATcgccgc                                                                                                                                    <  2:1764657/150‑1 (MQ=255)
                 aCTTGCAATAATGACCGCCAAATTGTTCTCGGCCTACATCGAGGCCTTGGAGAACGCTCACTTTTTTCTTCATCAGGTCCGGCCAGCATGACACGTACTCAACCTTCCGCTTCGCCAACCCCGTCGCGCGCCTGGCGCGCCATCGCcgcg                                                                                                                                   <  1:1694948/150‑1 (MQ=255)
                       aataatGACCGCCAAATTGTTCTCGGCCTACATCGAGGCCTTGGAGAACGCTCACTTTTTTCTTCATCAGGTCCGGCCAGCATGACACGTACTCAACCTTCCGCTTCGCCAACCCCGTCGCGCGCCTGGCGCGCCATCGCCGCGCTGaca                                                                                                                             >  1:319629/1‑150 (MQ=255)
                          aatGACCGCCAAATTGTTCTCGGCCTACATCGAGGCCTTGGAGAACGCTCACTTTTTTCTTCATCAGGTCCGGCCAGCATGACACGTACTCAACCTTCCGCTTCGCCAACCCCGTCGCGCGCCTGGCGCGCCATCGCCGCGCTGACACTc                                                                                                                          <  1:951844/150‑1 (MQ=255)
                                cGCCAAATTGTTCTCGGCCTACATCGAGGCCTTGGAGCACGCTCACTTTTTTCTTCATCAGGTCCGGCCAGCATGACACGTACTCAACCTTCCGCTTCGCCAACCCCGTCGCGCGCCTGGCGCGCCATCGCCGCGCTGACACTCAGCCTg                                                                                                                    >  1:962025/1‑150 (MQ=255)
                                 gCCAAATTGTTCTCGGCCTACATCGAGGCCTTGGAGAACGCTCACTTTTTTCTTCATCAGGTCCGGCCAGCATGACACGTACTCAACCTTCCGCTTCGCCAACCCCGTCGCGCGCCTGGCGCGCCATCGCCGCGCTGACACTCAGCCTgg                                                                                                                   >  1:37343/1‑150 (MQ=255)
                                  ccAAATTGTTCTCGGCCTACATCGAGGCCTTGGAGAACGCTCACTTTTTTCTTCATCTATCCCGGCCAGCATGACACGTACTCAACCTTCCGCTTCGCCAACCCCGTCGCGCGCCTGGCGCGCCATCGCCGCGCTGACACTCAGCCTGGt                                                                                                                  <  1:1601718/150‑1 (MQ=255)
                                  ccAAATTGTTCTCGGCCTACATCGAGGCCTTGGAGAACGCTCACTTTTTTCTTCATCAGGTCCGGCCAGCATGACACGTACTCAACCTTCCGCTTCGCCAACCCCGTCGCGCGCCTGGCGCGCCATCGCCGCGCTGACACTCAGCCTGGt                                                                                                                  <  2:849324/150‑1 (MQ=255)
                                   cAAATTGTTCTCGGCCTACATCGAGGCCTTGGAGAACGCTCACTTTTTTCTTCATCAGGTCCGGCCAGCATGACACGTACTCAACCTTCCGCTTCGCTAACCCCGTCGCGCGCCTGGCGCGCCATCGCCGCGCTGACACTCAGCCTGGTg                                                                                                                 >  2:1874343/1‑150 (MQ=255)
                                      aTTGTTCTCGGCCTACATCGAGGCCTTGGAGAACGCTCACTTTTTTCTTCATCAGGTCCGGCCAGCATGACACGTACTCAACCTTCCGCTTCGCCAACCCCGTCGCGCGCCTGGCGCGCCATCGCCGCGCTGACACTCAGCCTGGTGCTg                                                                                                              >  2:1081661/1‑150 (MQ=255)
                                      aTTGTCATCGGCGAACTTCGAGGCCATGTAGTAAGCGCACTGTTTTCTTCATCAGGTCCGGCCAGCATGACACGTACTCAACCTTCCGCTTCGCCAACCCCGTCGCGCGCCTGGCGCGCCATCGCCGCGCTGACACTCAGCCTGGTGCTg                                                                                                              >  2:1353106/1‑150 (MQ=255)
                                             tcGGCCTACATCGAGGCCTTGGAGAACGCTCACTTTTTTCTTCATCAGGTCCGGCCAGCATGACACGTACTCAACCTTCCGCTTCGCCAACCCCGTCGCGCGCCTGGCGCGCCATCGCCGCGCTGACACTCAGCCTGGTGCTGGCCCCgg                                                                                                       >  2:956043/1‑150 (MQ=255)
                                             tcGGCCTACATCGAGGCCTTGGAGAACGCTCACTTTTTTCTTCATCAGGTCCGGCCAGCATGACACGTACTCAACCTTCCGCTTCGCCAACCCCGTCGCGCGCCTGGCGCGCCATCGCCGCGCTGACACTCAGCCTGGTGCTGGCCCCgg                                                                                                       >  1:2189920/1‑150 (MQ=255)
                                              cGGCCTACATCGAGGCCTTGGAGAACGCTCACTTTTTTCTTCATCAGGTCCGGCCAGCATGACACGTACTCAACCTTCCGCTTCGCCAACCCCGTCGCGCGCCTGGCGCGCCATCGCCGCGCTGACACTCAGCCTGGTGCTGGCCCCGGt                                                                                                      <  2:75604/150‑1 (MQ=255)
                                                   tACATCGAGGCCTTGGAGAACGCTCACTTTTTTCTTCATCAGGTCCGGCCAGCATGACACGTACTCAACCTTCCGCTTCGCCAACCCCGTCGCGCGCCTGGCGCGCCATCGCCGCGCTGACACTCAGCCTGGTGCTGGCCCCGGTGGCCa                                                                                                 <  1:1081661/150‑1 (MQ=255)
                                                     cagcgtggCCTTGGAGAACGCTCACTTTTTTCTTCATCAGGTCCGGCCAGCATGACACGTACTCAACCTTCCGCTTCGCCAACCCCGTCGCGCGCCTGGCGCGCCATCGCCCCGCTGCCACTCAGCCTGGTGCTGGCCCCGGTCGCCCTg                                                                                               >  2:408720/7‑150 (MQ=255)
                                                      aTCGAGGCCTTGGAGAACGCTCACTTTTTTCTTCATCAGGTCCGGCCAGCATGACACGTACTCAACCTTCCGCTTCGCCAACCCCGTCGCGCGCCTGGCGCGCCATCGCCGCGCTGACACTCAGCCTGGTGCTGGCCCCGGTGGCCATgg                                                                                              >  2:331026/1‑150 (MQ=255)
                                                      aTCGAGGCCTTGGAGAACGCTCACTTTTTTCTTCATCAGGTCCGGCCAGCATGACACGTACTCAACCTTCCGCTTCGCCAACCCCGTCGCGCGCCTGGCGCGCCATCGCCGCGCTGACACTCAGCCTGGTGCTGGCCCCGGTGGCCATgg                                                                                              <  2:396289/150‑1 (MQ=255)
                                                      aTCGAGGCCTTGGAGAACGCTCACTTTTTTCTTCATCAGGTCCGGCCAGCATGACACGTACTCAACCTTCCGCTTCGCCAACCCCGTCGCGCGCCTGGCGCGCCATCGCCGCGCTGACACTCAGCCTGGTGCTGGCCACCGGGGCCagca                                                                                              <  2:319629/150‑4 (MQ=255)
                                                        cGAGGCCTTGGAGAACGCTCACTTTTTTCTTCACCAGGTCCGGCCAGCATGACACGTACTCAACCTTCCGCTTCGCTAACCCCGTCGCGCGCCTGTCGCGCCATCGCCGCGCTGACACTCCGCCTGGTGCTGGCCCCGGTGGCCATGGcc                                                                                            <  1:1379507/150‑1 (MQ=255)
                                                              cTTGGAGAACGCTCACTTTTTTCTTCATCAGGTCCGGCCAGCATGACACGTACTCAACCTTCCGCTTCGCCAACCCCGTCGCGCGCCTGGCGCGcc                                                                                                                                            >  1:1946591/1‑96 (MQ=255)
                                                                      aCGCTCACTTTTTTCTTCATCAGGTCCGGCCAGCATGACACGTACTCAACCTTCCGCTTCGCCAACCCCGTCGCGCGCCTGGCGCGCCATCGCCGCGCTGACACTCAGCCTGGTGCTGGCCCCGGTGGCCATGGCCGATGAGCCAaccac                                                                              <  1:956043/150‑1 (MQ=255)
                                                                        gCTCACTTTTTTCTTCATCAGGTCCGGCCAGCATGACACGTACTCAACCTTCCGCTTCGCCAACCCCGTCGCGCGCCTGGCGCGCCATCGCCGCGCTGACACTCAGCCTGGTGCTGGCCCCGGTGGCCATGGCCGa                                                                                          <  1:469313/136‑1 (MQ=255)
                                                                        gCTCACTTTTTTCTTCATCAGGTCCGGCCAGCATGACACGTACTCAACCTTCCGCTTCGCCAACCCCGTCGCGCGCCTGGCGCGCCATCGCCGCGCTGACACTCAGCCTGGTGCTGGCCCCGGTGGCCATGGCCGa                                                                                          >  2:469313/1‑136 (MQ=255)
                                                                        gCTCACTTTTTTCTTCATCAGGTCCGGCCAGCATGACACGTACTCAACCTTCCGCTTCGCCAACCCCGTCGCGCGCCTGGCGCGCCATCGCCGCGCTGACACTCAGCCTGGTGCTGGCCCCGGTGGCCATGGCCGATGAGCCAAccacca                                                                            >  1:1177608/1‑150 (MQ=255)
                                                                        gCTCACTTTTTTCTTCATCAGGTCCGGCCAGCATGACACGTACTCAACCTTCCGCTTCGCCAACCCCGTCGCGCGCCTGGCGCGCCATCGCCGCGCTGACACTCAGCCTGGTGCTGGCCCCGGTGGCCATGGCCGATGAGCCAAccacca                                                                            >  1:253369/1‑150 (MQ=255)
                                                                          tCACTTTTTTCTTCATCAGGTCCGGCCAGCATGACACGTACTCAACCTTCCGCTTCGCCAACCCCGTCGCGCGCCTGGCGCGCCATCGCCGCGCTGACACTCAGCCTGGTGCTgg                                                                                                             >  2:1006725/1‑115 (MQ=255)
                                                                          tCACTTTTTTCTTCATCAGGTCCGGCCAGCATGACACGTACTCAACCTTCCGCTTCGCCAACCCCGTCGCGCGCCTGGCGCGCCATCGCCGCGCTGACACTCAGCCTGGTGCTgg                                                                                                             <  1:1006725/115‑1 (MQ=255)
                                                                            aCTTTTTTCTTCATCAGGTCCGGCCAGCATGACACGTACTCAACCTTCCGCTTCGCCAACCCCGTCGCGCGCCTGGCGCGCCATCGCCGCGCTGACACTCAGCCTGGTGCTGGCCCCGGTGGCCATGGCCGATGAGCCAACCACCAAcgc                                                                        >  2:399084/1‑150 (MQ=255)
                                                                              ttttttCTTCATCAGGTCCGGCCAGCATGACACGTACTCAACCTTCCGCTTCGCCAACCCCGTCGCGCGCCTGGCGCGCCATCGCCGCGCTGACACTCAGCCTGGTGCTGGCCCCGGTGGCCa                                                                                                 >  1:55608/1‑123 (MQ=255)
                                                                              ttttttCTTCATCAGGTCCGGCCAGCATGACACGTACTCAACCTTCCGCTTCGCCAACCCCGTCGCGCGCCTGGCGCGCCATCGCCGCGCTGACACTCAGCCTGGTGCTGGCCCCGGTGGCCATGGCCGATGAGc                                                                                     >  2:1390730/1‑135 (MQ=255)
                                                                              ttttttCTTCATCAGGTCCGGCCAGCATGACACGTACTCAACCTTCCGCTTCGCCAACCCCGTCGCGCGCCTGGCGCGCCATCGCCGCGCTGACACTCAGCCTGGTGCTGGCCCCGGTGGCCATGGCCGATGAGCCAACCACCAAcgccg                                                                      >  1:363842/1‑150 (MQ=255)
                                                                               tttttCTTCATCAGGTCCGGCCAGCATGACACGTACTCAACCTTCCGCTTCGCCAACCCCGTCGCGCGCCTGGCGCGCCATCGCCGCGCTGACACTCAGCCTGGTGCTGGCCCCGGTGGCCATGGCCGATGAGCCAACCACCAAcgccgc                                                                     <  2:2213142/150‑1 (MQ=255)
                                                                                   tcttcATCAGGTCCGGCCAGCATGACACGTACTCAACCTTCCGCTTCGCCAACCCCGTCGCGCGCCTGGCGCGCCATCGCCGCGCTGACACTCAGCCTGGTGCTGGCCCCGGTGGCCATGGCCGATGAGCCAACCACCAACGCCGCcaca                                                                 >  2:1492733/1‑150 (MQ=255)
                                                                                     ttcATCAGGTCCGGCCAGCATGACACGTACTCAACCTTCCGCTTCGCCAACCCCGTCGCGCGCCTGGCGCGcc                                                                                                                                            <  1:511893/73‑1 (MQ=255)
                                                                                     ttcATCAGGTCCGGCCAGCATGACACGTACTCAACCTTCCGCTTCGCCAACCCCGTCGCGCGCCTGGCGCGcc                                                                                                                                            >  2:511893/1‑73 (MQ=255)
                                                                                        atcaGGTCCGGCCAGCATGACACGTACTCAACCTTCCGCTTCGCCAACCCCGTCGCGCGCCTGGCGCGCCATCGCCGCGCTGACACTCAGCCTGGTGCTGGCCCCGGTGGCCATGGCCGATGAGCCAACCACCAACGCCGCCACACctgc                                                            >  1:1972575/1‑150 (MQ=255)
                                                                                         tcaGGTCCGGCCAGCATGACACGTACTCAACCTTCCGCTTCGCCAACCCCGTCGCGCGCCTGGCGCGCCATCGCCGCGCTGACACTCAGCCTGGTGCTGCCCCCGGTGGCCATGGCCGATGAGCCAACCACCAACGCCGCCACACctgct                                                           >  1:2331141/1‑150 (MQ=255)
                                                                                           aGGTCCGGCCAGCATGACACGTACTCAACCTTCCGCTTCGCCAACCCCGTCGCGCGCCTGGCGCGCCATCGCCGCGCTGACACTCAGCCTGGTGCTGGCCCCGGTGGCCATGGCCGATGAGCCAACCACCAACGCCGCCACACctgctgc                                                         >  1:1357629/1‑150 (MQ=255)
                                                                                            ggTCCGGCCAGCATGACACGTACTCAACCTTCCGCTTCGCCAACCCCGTCGCGCGCCTGGCGCGCCATCGCCGCGCTGACACTCAGCCTGGTGCTGGCCCCGGTGGCCATGGCCGATGAGCCAACCACCAACGCCGCCACACctgctgct                                                        >  2:1033951/1‑150 (MQ=255)
                                                                                             gTCCGGCCAGCATGACACGTACTCAACCTTCCGCTTCGCCAACCCCGTCGCGCGCCTGGCGCGCCATCGCCGCGCTGACACTCAGCCTGGTGCTGGCCCCGGTGGCCATGGCCGATGAGCCAACCACCAACGCCGCCACACctgctgctg                                                       >  1:755172/1‑150 (MQ=255)
                                                                                             gTCCGGCCAGCATGACACGTACTCAACCTTCCGCTTCGCCAACCCCGTCGCGCGCCTGGCGCGCCATCGCCGCGCTGACACTCAGCCTGGTGCTGGCCCCGGTGGCCATGGCCGATGAGCCAACCACCAACGCCGCCACACctgctgctg                                                       >  1:1517235/1‑150 (MQ=255)
                                                                                             gTCCGGCCAGCATGACACGTACTCAACCTTCCGCTTCGCCAACCCCGTCGCGCGCCTGGCGCGCCATCGCCGCGCTGACACTCAGCCTGGTGCTGGCCCCGGTGGCCATGCACGATGAGCCAACCACCAACGCCGCCACACctgctgctg                                                       >  1:1231253/1‑150 (MQ=255)
                                                                                              tCCGGCCAGCATGACACGTACTCAACCTTCCGCTTCGCCAACCCCGTCGCGCGCCTGGCGCGCCATCGCCGCGCTGACACTCACCCTGGTGCTGGCCCCGGTGGCCATGGCCGATGAGCCAACCACCAACGCCGCCACACCtgctgctgc                                                      <  2:1142595/150‑1 (MQ=255)
                                                                                                cGGCCAGCATGACACGTACTCAACCTTCCGCTTCGCCAACCCCGTCGCGCGCCTGGCGCGCCATCGCCGCGCTGACACTCAGCCTg                                                                                                                    >  1:566199/1‑86 (MQ=255)
                                                                                                cGGCCAGCATGACACGTACTCAACCTTCCGCTTCGCCAACCCCGTCGCGCGCCTGGCGCGCCATCGCCGCGCTGACACTCAGCCTg                                                                                                                    <  2:566199/86‑1 (MQ=255)
                                                                                                        aTGACACGTACTCAACCTTCCGCTTCGCCAACCCCGTCGCGCGCCTGGCGCGCCATCGCCGCGCTGACACTCAGCCTGGTGCTGGCCCCGGTGGCCATGGCCGATGAGCCAACCACCAACGCCGCCACCACTACTGCTGCCACTGCCCCg                                            >  2:1641028/1‑150 (MQ=255)
                                                                                                        aTGACACGTACTCAACCTTCCGCTTCGCCAACCCCGTCGCGCGCCTGGCGCGCCATCGCCGCGCTGACACTCAGCCTGGTGCTGGCCCCGGTGGCCATGGCCGATGAGCCAACCACCAACGCCGCCACACCTGCTGCTGCCACTGCCCCg                                            >  2:1294791/1‑150 (MQ=255)
                                                                                                          gACACGTACTCAACCTTCCGCTTCGCCAACCCCGTCGCGCGCCTGGCGCGCCATCGCCGCGCTGACACTCAGCCTGGTGCTGGCCCCGGTGGCCATGGCCGATGAGCCAACCACCAACGCCGCCACACCTGCTGCTGCCACTGCCCCGGc                                          <  2:1842889/150‑1 (MQ=255)
                                                                                                           acacGTACTCAACCTTCCGCTTCGCCAACCCCGTCGCGCGCCTGGCGCGCCATCGCCGCGCTGACACTCAGCCTGGTGCTGGCCCCGGTGGCCATGGCCGATGAGCCAACCACCAACGCCGCCACACCTGCTGCTGCCACTGCCCCGGcc                                         <  2:1177608/150‑1 (MQ=255)
                                                                                                            cacGTACTCAACCTTCCGCTTCGCCAACCCCGTCGCGCGCCTGGCGCGCCATCGCCGCGCTGACACTCAGCCTGGTGCTGGCCCCGGTGGCCATGGCCGATGAGCCAACCACCAACGCCGCCACACCTGCTGCTGCCACTGCCCCGGCca                                        <  1:1052947/150‑1 (MQ=255)
                                                                                                             acGTACTCAACCTTCCGCTTCGCCAACCCCGTCGCGCGCCTGGCGCGCCATCGCCGCGCTGACACTCAGCCTGGTGCTGGCCCCGGTGGCCATGGCCGATGAGCCAACCACCAACGCCGCCACACCTGCTGCTGCCACTGCCCCGGCcac                                       >  2:621599/1‑150 (MQ=255)
                                                                                                             acGTACTCAACCTTCCGCTTCGCCAACCCCGTCGCGCGCCTGGCGCGCCATCGCCGCGCTGACACTCAGCCTGGTGCTGGCCCCGGTGGCCATGGCCGATGAGCCAACCACCAACGCCGCCACACCTGCTGCTGCCACTGCCCCGGCcac                                       >  1:1796634/1‑150 (MQ=255)
                                                                                                               gTACTCAACCTTCCGCTTCGCCAACCCCGTCGCGCGCCTGGCGCGCCATCGCCGCGCTGACACTCAGCCTGGTGCTGGCCCCGGTGGCCATGGCCGATGAGCCAACCACCAACGCCGCCACACCTGCTGCTGCCACTGCCCCGGCCacac                                     <  1:1874343/150‑1 (MQ=255)
                                                                                                                  cTCAACCTTCCGCTTCGCCAACCCCGTCGCGCGCCTGGCGCGCCATCGCCGCGCTGACACTCAGCCTGGTGCTGGCCCCGGTGGCCATGGCCGATGAGCCAACCACCAACGCCGCCACACCTGCTGCTGCCACTGCCCCGGCCACACctg                                  <  1:1559045/150‑1 (MQ=255)
                                                                                                                   tCAACCTTCCGCTTCGCCAACCCCGTCGCGCGCCTGGCGCGCCATCGCCGCGCTGACACTCAGCCTGGTGCTGGCCCCGGTGGCCATGGCCGATGAGCCAACCACCAACGCCGCca                                                                   <  2:1298186/116‑1 (MQ=255)
                                                                                                                   tCAACCTTCCGCTTCGCCAACCCCGTCGCGCGCCTGGCGCGCCATCGCCGCGCTGACACTCAGCCTGGTGCTGGCCCCGGTGGCCATGGCCGATGAGCCAACCACCAACGCCGCca                                                                   >  1:1298186/1‑116 (MQ=255)
                                                                                                                   tCAACCTTCCGCTTCGCCAACCCCGTCGCGCGCCTGGCGCGCCATCGCCGCGCTGACACTCAGCCTGGTGCTGGCCCCGGTGGCCATGGCCGATGAGCCAACCACCAACGCCGCCACACCTGCTGCTGCCACTGCCCCGGCCACACctgc                                 >  2:1183285/1‑150 (MQ=255)
                                                                                                                     aaCCTTCCGCTTCGCCAACCCCGTCGCGCGCCTGGCGCGCCATCGCCGCGCTGACACTCAGCCTGGTGCTGGGCCCGGTGGCCATGGCCGATGAGCCAACCACCAACGCCGCCACACCTGCTGCTGCCACTGCCCCGGCCACACctgctg                               >  2:277299/1‑150 (MQ=255)
                                                                                                                     aaCCTTCCGCTTCGCCAACCCCGTCGCGCGCCTGGCGCGCCATCGCCGCGCTGACACTCAGCCTGGTGCTGGCCCCGGTGGCCATGGc                                                                                             >  1:1319582/1‑88 (MQ=255)
                                                                                                                     aaCCTTCCGCTTCGCCAACCCCGTCGCGCGCCTGGCGCGCCATCGCCGCGCTGACACTCAGCCTGGTGCTGGCCCCGGTGGCCATGGCCGATGAGCCAACCACCAACGCCGCCACACCTGCTGCTGCCACTGCCCCGGCCACACctgctg                               >  2:250880/1‑150 (MQ=255)
                                                                                                                     aaCCTTCCGCTTCGCCAACCCCGTCGCGCGCCTGGCGCGCCATCGCCGCGCTGACACTCAGCCTGGTGCTGGCCCCGGTGGCCATGGCCGATGAGCCAACCACCAACCCCGCCACACCTGCTGCTGCCACTGCCCCGGCCACAcctcctc                               >  2:989775/1‑149 (MQ=255)
                                                                                                                       ccTTCCGCTTCGCCAACCCCGTCGCGCGCCTGGCGCGCCATCGCCGCGCTGACACTCAGCCTGGTGCTGGCCCCGGTGGCCATGGCCGATGAGCCAACCACCAACGCCGCCACACCTGCTGCTGCCACTGCCCCGGCCACACctgctgct                             <  2:1067227/150‑1 (MQ=255)
                                                                                                                        cTTCCGCTTCGCCAACCCCGTCGCGCGCCTGGCGCGCCATCGCCGCGCTGACACTCAGCCTGGTGCTGGCCCCGGTGGCCATGGCCGATGAGCCAACCACCAACGCCGCCACACCTGCTGCTGCCACTGCCCCGGCCACACctgctgctg                            >  1:2207234/1‑150 (MQ=255)
                                                                                                                        cTTCCGCTTCGCCAACCCCGTCGCGCGCCTGGCGCGCCATCGCCGCGCTGACACTCAGCCTGGTGCTGGCCCCGGTGGCCATGGCCGATGAGCCAACCACCAACGCCGCCACACCTGCTGCTGCCACTGCCCCGGCCACACct                                   >  1:1004277/1‑143 (MQ=255)
                                                                                                                        cTTCCGCTTCGCCAACCCCGTCGCGCGCCTGGCGCGCCATCGCCGCGCTGACACTCAGCCTGGTGCTGGCCCCGGTGGCCATGGCCGATGAGCCAACCACCAACGCCGCCACACCTGCTGCTGCCACTGCCCCGGCCACACct                                   <  2:1004277/143‑1 (MQ=255)
                                                                                                                        cTTCCGCGTCGCCAACCCCGTCGCGCGCCTGGCGCGCCATCGCCGCGCTGACACTCAGCCTGGTGCTGGCCCCGGTGGCCATGGCCGATGAGCCAACCACCAACGCCGCCACACCTGCTGCTGCCACTGCCCCGGCCACACctgctgctg                            >  1:2289490/1‑150 (MQ=255)
                                                                                                                         ttCCGCTTCGCCAACCCCGTCGCGCGCCTGGCGCGCCATCGCCTCGCTGACACTCAGCCTGGTGCTGGCCCCGGTGGCCATGGCCGAt                                                                                         <  2:208665/88‑1 (MQ=255)
                                                                                                                         ttCCGCTTCGCCAACCCCGTCGCGCGCCTGGCGCGCCATCGCCTCGCTGACACTCAGCCTGGTGCTGGCCCCGGTGGCCATGGCCGAt                                                                                         >  1:208665/1‑88 (MQ=255)
                                                                                                                         ttCCGCTTCGCCAACCCCGTCGCGCGCCTGGCGCGCCATCGCCGCGCTGACACTCAGCCTGGTGCTGGCCCCGGTGGCCATGGCCGATGAGCCAACCACCAACGCCGCCACACCtgctgctgc                                                      >  1:1166903/1‑123 (MQ=255)
                                                                                                                         ttCCGCTTCGCCAACCCCGTCGCGCGCCTGGCGCGCCATCGCCGCGCTGACACTCAGCCTGGTGCTGGCCCCGGTGGCCATGGCCGATGAGCCAACCACCAACGCCGCCACACCtgctgctgc                                                      <  2:1166903/123‑1 (MQ=255)
                                                                                                                         ttCCGCTTCGCCAACCCCGTCGCGCGCCTGGCGCGCCATCGCCGCGCTGACACTCAGCCTGGTGCTGGCCCCGGTGGCCATGGCCGATGAGCCAACCACCAACGCCGCCACACCTGCTGCTGCCACTGCCCCGGCCACACCtgctgctgc                           <  1:823629/150‑1 (MQ=255)
                                                                                                                         ttCCGCTTCGCCAACCCCGTCGCGCGCCTGGCGCGCCATCGCCGCGCTGACACTCAGCCTGGTGCTGGCCCCGGTGGCCATGGCCGATGAGCCAACCACCAACGCCGCCACACCTGCTGCTGCCACTGCCCCGGCCACACCtgctgctgc                           <  2:2207234/150‑1 (MQ=255)
                                                                                                                           ccGCTTCGCCAACCCCGTCGCGCGCCTGGCGCGCCATCGCCGCGCTGACACTTAGCCTGGTGCTGGCCCCGGTGGCCATGGCCGATGAGCCAACCACCAACGCCGCCACACCTGCTGCTGCCACTGCCCCGGCCACACCTGCTGCTGccc                         >  1:1266863/1‑150 (MQ=255)
                                                                                                                           ccGCTTCGCCAACCCCGTCGCGCGCCTGGCGCGCCATCGCCGCGCTGACACTCAGCCTGGTGCTGGCCCCGGTGGCCATGGCCGATGAGCCAACCACCAACGCCGCCACACCTGCTGCTGCCACTGCCCCGGCCACACCTGCTGCTGtcc                         <  1:1033951/150‑1 (MQ=255)
                                                                                                                             gCTTCGCCAACCCCGTCGCGCGCCTGGCGCGCCATCGCCGCGCTGACACTCAGCCTGGTGCTGGCCCCGGTGGCCATGGCCGATGAGCCAACCACCAACGCCGCCacac                                                                >  2:266490/1‑109 (MQ=255)
                                                                                                                             gCTTCGCCAACCCCGTCGCGCGCCTGGCGCGCCATCGCCGCGCTGACACTCAGCCTGGTGCTGGCCCCGGTGGCCATGGCCGATGAGCCAACCACCAACGCCGCCacac                                                                <  1:266490/109‑1 (MQ=255)
                                                                                                                                   ccAACCCCGTCGCGCGCCTGGCGCGCCATCGCCGCGCTGACACTCAGCCTGGTGCTGGCCCCGGTGGCCATGGCCGATGAGCCAACCACCAACGCCGCCACACctg                                                             <  2:1277473/106‑1 (MQ=255)
                                                                                                                                   ccAACCCCGTCGCGCGCCTGGCGCGCCATCGCCGCGCTGACACTCAGCCTGGTGCTGGCCCCGGTGGCCATGGCCGATGAGCCAACCACCAACGCCGCCACACctg                                                             >  1:1277473/1‑106 (MQ=255)
                                                                                                                                   ccAACCCCGTCGCGCGCCTGGCGCGCCATCGCCGCGCTGACACTCAGCCTGGTGCTGGCCCCGGTGGCCATGGCCGATGAGCCAACCACCAACGCCGCCACACCTGCTGCTGCCACTGCCCCGGCCACACCTGCTGCTGCCCCGGCGGAt                 >  1:1517570/1‑150 (MQ=255)
                                                                                                                                    cAACCCCGTCGCGCGCCTGGCGCGCCATCGCCgcgc                                                                                                                                  >  2:2192273/1‑36 (MQ=255)
                                                                                                                                    cAACCCCGTCGCGCGCCTGGCGCGCCATCGCCGCGCTGACACTCAGCCTGGTGCTGGCCCCGGTGGCCATGGCCGATGAGCCAACCACCAACGCCGCCACACCTGCTGCTGCCACTGCCCCGGCCACACCTGCTGCTGCCCCGGCGGATg                >  2:1890039/1‑150 (MQ=255)
                                                                                                                                        cccGTCGCGCGCCTGGCGCGCCATCGCCGCGCTGACACTCAGCCTGGTGCTGGCCCCGGTGGCCATGGCCGATGAGCCAACCACCAACGCCGCCACACctg                                                             <  1:1993502/101‑1 (MQ=255)
                                                                                                                                        cccGTCGCGCGCCTGGCGCGCCATCGCCGCGCTGACACTCAGCCTGGTGCTGGCCCCGGTGGCCATGGCCGATGAGCCAACCACCAACGCCGCCACACctg                                                             >  2:1993502/1‑101 (MQ=255)
                                                                                                                                        cccGTCGCGCGCCTGGCGCGCCATCGCCGCGCTGACACTCAGCCTGGTGCTGGCCCCGGTGGCCATGGCCGATGAGCCAACCACCAACGCCGCCACACCTGCTGCTGCCACTGCCCCGGCCACACCTGCTGCTGCCCCGGCGGATGCTcc            >  1:995705/1‑150 (MQ=255)
                                                                                                                                          gatggcgcgcCAGGCGCGCCATCGCCGCGCTGACACTCAGCCTGGTGCTGGCCCCGGTGGCCATGGCCGATGAGCCAACCACCAACGCCGCCACACCTGCTGCTGCCACTGCCCCGGCCACACCTGCTGCTGCCCCGGCGGATGCTCCgg          >  2:244938/5‑150 (MQ=255)
                                                                                                                                          gatggcgcgcCAGACGCGCGATCGCCGCGCTGACACTCAGCCTGGTGCTGGCCCCGGTGGCCATGGCCGATGAGCCAACCACCAACGCCGCCACACCTGCTGCTGCCACTGCCCCGGCCACACCTGCTGCTGCCCCGGCGGATGCTCCgg          >  2:69547/5‑150 (MQ=255)
                                                                                                                                           atggcgcgcCAGGCGCGCGATCGCCGCGCTGACACTCAGCCTGGTGCTGGCCCCGGTGGCCATGGCCGATGAGCCAACCACCAACGCCGCCACACCTGCTGCTGCCACTGCCCCGGCCACACCTGCTGCTGCCCCGGCGGATGCTCCGgc         >  2:1476955/4‑150 (MQ=255)
                                                                                                                                           atggcgcgcCAGGCGCGCCATCGCCGCGCTGACACTCAGCCTGGTGCTGGCCCCGGTGGCCATGGCCGATGAGCCAACCACCAACGCCGCCACACCTGCTGCTGCCACTGCCCCGGCCACACCTGCTGCTGCCCCGGCGGATGCTCCGgc         >  1:1512300/4‑150 (MQ=255)
                                                                                                                                                gcgcCTGCCGCGCCATCGCCGCGCTGACACTCAGCCTGGTGCTGGCCCCGGTGGCCATGGCCGATGAGCCAACCACCAACGCAGCCACACCTGCTCCTGCCACTGCCCCGGCCACACCTGCTACTGCCCCGGCGGATGCTCCGGCGCCTg    >  2:244347/1‑150 (MQ=255)
                                                                                                                                                  gcCTGGCGCGCCATCGCCGCGCTGACACTCAGCCTGGTGCTGGCCCCGGTGGCCATGGCCGATGAGCCAACCACCAACGCCGCCACACCTGCTGCTGCCACTGCCCCGGCCACACCTGCTGCTGCCCCGGCGGATGCTCCGGCGCCTGAt  <  2:2189920/150‑1 (MQ=255)
                                                                                                                                                  gcCTGGCGCGCCATCGCCGCGCTGACACTCAGCCTGGTGCTGGCCCCGGTGGCCATGGCCGATGAGCCAACCACCAACGCCGCCACACCTGCTGCTGCCACTGCCCCGGCCACACCTGCTGCTGCCCCGGCGGATGCTCCGGCGCCTGAt  <  1:250880/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
TTTTGTTAACAAGAATTACTTGCAATAATGACCGCCAAATTGTTCTCGGCCTACATCGAGGCCTTGGAGAACGCTCACTTTTTTCTTCATCAGGTCCGGCCAGCATGACACGTACTCAACCTTCCGCTTCGCCAACCCCGTCGCGCGCCTGGCGCGCCATCGCCGCGCTGACACTCAGCCTGGTGCTGGCCCCGGTGGCCATGGCCGATGAGCCAACCACCAACGCCGCCACACCTGCTGCTGCCACTGCCCCGGCCACACCTGCTGCTGCCCCGGCGGATGCTCCGGCGCCTGAT  >  NC_002947/6053609‑6053904

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