Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 549,670 G→A 9.3% P114P (CCG→CCA tufB → elongation factor Tu‑B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947549,6700GA 9.3% 521.2 / 32.8 182P114P (CCG→CCAtufBelongation factor Tu‑B
Reads supporting (aligned to +/- strand):  ref base G (71/94);  new base A (5/12);  total (76/106)
Fisher's exact test for biased strand distribution p-value = 3.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.93e-01

Only 100 of 185 total aligned reads displayed.

ACCGCTCACGTCGAGTACAACTCGAACATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCC  >  NC_002947/549521‑549819
                                                                                                                                                     |                                                                                                                                                     
     tCACGTCGAGTACAACTCGAACATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAAc                                                                                                                                                  >  2:4802643/1‑150 (MQ=34)
          tCGAGTACAACTCGAACATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTg                                                                                                                                             >  1:2678890/1‑150 (MQ=34)
           cGAGTACAACTCGAACATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGa                                                                                                                                            <  1:3389877/150‑1 (MQ=34)
           cGAGTACAACTCGAACATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGa                                                                                                                                            >  1:3218112/1‑150 (MQ=34)
           cGAGTACAACTCGAACATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGa                                                                                                                                            <  1:4143559/150‑1 (MQ=34)
           cGAGTACAACTCGAACATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGa                                                                                                                                            <  2:667593/150‑1 (MQ=34)
           cGAGTACAACTCGAACATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGa                                                                                                                                            <  2:647421/150‑1 (MQ=34)
           cGAGTACAACTCGAACATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGa                                                                                                                                            <  2:5297885/150‑1 (MQ=34)
              gttgtaCTCGAACATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGca                                                                                                                                         >  2:1192803/6‑150 (MQ=21)
                aCAACTCGAACATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGcaca                                                                                                                                       >  2:1767351/1‑150 (MQ=34)
                  aaCTCGAACATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATc                                                                                                                                     <  2:824375/150‑1 (MQ=34)
                  aaCTCGAACATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATc                                                                                                                                     <  1:5411982/150‑1 (MQ=34)
                  aaCTCGAACATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATc                                                                                                                                     <  2:57429/150‑1 (MQ=34)
                  aaCTCGAACATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATc                                                                                                                                     >  2:591708/1‑150 (MQ=34)
                     tCGAACATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCctt                                                                                                                                  <  2:704934/150‑2 (MQ=21)
                     tCGAACATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCctg                                                                                                                                  >  1:5101423/1‑150 (MQ=34)
                       gAACATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCctgct                                                                                                                                <  2:1529482/150‑1 (MQ=35)
                       gAACATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCctgct                                                                                                                                >  1:354782/1‑150 (MQ=35)
                       gAACATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCctgct                                                                                                                                >  1:354795/1‑150 (MQ=35)
                               gTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTc                                                                                                                        >  1:2301326/1‑150 (MQ=34)
                               gTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTc                                                                                                                        >  1:315105/1‑150 (MQ=34)
                                  aCTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAgg                                                                                                                     >  2:1206205/1‑150 (MQ=34)
                                   cTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGt                                                                                                                    >  2:2672153/1‑150 (MQ=34)
                                    tACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTa                                                                                                                   <  1:5097753/150‑1 (MQ=34)
                                         tCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGt                                                                                                              >  1:3638470/1‑150 (MQ=34)
                                         tCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGt                                                                                                              >  1:4962985/1‑150 (MQ=34)
                                         tCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGt                                                                                                              >  1:1337406/1‑150 (MQ=34)
                                            cGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTcc                                                                                                           >  1:3577833/1‑150 (MQ=34)
                                                gACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTAc                                                                                                       >  1:2359316/1‑150 (MQ=34)
                                                  cTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACAt                                                                                                     <  1:2297035/150‑1 (MQ=34)
                                                    gCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCg                                                                                                   >  1:1188997/1‑150 (MQ=34)
                                                    gCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCg                                                                                                   <  1:3285081/150‑1 (MQ=34)
                                                    gCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCg                                                                                                   >  1:1188918/1‑150 (MQ=34)
                                                    gCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCg                                                                                                   <  1:2539827/150‑1 (MQ=34)
                                                      ccAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCTCGTCAGGTAGGCGTTCCGTACATCGTg                                                                                                 <  2:574357/150‑1 (MQ=21)
                                                        aGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGt                                                                                               <  2:3990747/150‑1 (MQ=34)
                                                            cACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTc                                                                                           <  2:1121109/150‑1 (MQ=34)
                                                              cGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCt                                                                                         <  1:327763/150‑1 (MQ=34)
                                                               gCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTg                                                                                        >  2:3960075/1‑150 (MQ=34)
                                                                    cTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACaa                                                                                   <  2:5067016/150‑1 (MQ=34)
                                                                      aCGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAgg                                                                                 <  2:5101423/150‑1 (MQ=34)
                                                                        gTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGcc                                                                               >  2:5194281/1‑150 (MQ=35)
                                                                          gaagaaCATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGa                                                                             >  2:2084506/1‑150 (MQ=37)
                                                                            agaaCATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGAcc                                                                           <  1:2687798/150‑1 (MQ=37)
                                                                                cATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGt                                                                       <  2:1448833/150‑1 (MQ=37)
                                                                                 aTGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTa                                                                      <  2:3050406/150‑1 (MQ=37)
                                                                                 aTGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTa                                                                      <  2:3944874/150‑1 (MQ=37)
                                                                                    aTCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAgac                                                                   <  1:2398082/150‑1 (MQ=37)
                                                                                     tCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAgacg                                                                  >  1:4907852/1‑150 (MQ=37)
                                                                                          ggTGCTGCCCAGATGGACGGCGCGAACCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGACTTCCTGAACAAGTCCGACCTGGTAGACGACGCt                                                             >  2:2539921/1‑150 (MQ=255)
                                                                                          ggTGCTGCACAGATGGACGGCGCGATCCTGGTTGTTTCGGCAGCTGACGGCCCGATGCCACAGACCCGTGAGCACATCCTGCTCGCTCGTCAGGTTGGCGTTCCTGCGATCGTCGTGTTCCTGAACAAGTGCGATCAGGTTGACGATGCt                                                             >  1:3063394/1‑150 (MQ=2)
                                                                                               tgcACAGATGGACGGCGCGATCCTGGTTGTTTCGGCAGCTGACGGCCCGATGCCACAGACCCGTGAGCACATCCTGCTCGCTCGTCAGGTTGGCGTTCCTGCGATCGTCGTGTTCCTGAACAAGTGCGATCAGGTTGACGATGCTGAact                                                        >  1:570552/1‑150 (MQ=2)
                                                                                                 cccAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAgctgc                                                      <  1:995926/150‑1 (MQ=37)
                                                                                                 cccAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAgctgc                                                      >  2:4097632/1‑150 (MQ=37)
                                                                                                  ccAGATGGACGGCGCTATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAgctgct                                                     >  1:3158865/1‑150 (MQ=37)
                                                                                                  ccAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAgctgct                                                     >  1:747200/1‑150 (MQ=37)
                                                                                                  ccAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAgctgct                                                     >  2:3459174/1‑150 (MQ=37)
                                                                                                   cAGATGGACGGCGCGATCCTGGTTGTTTCGGCAGCTGACGGCCCGATGCCACAGACCCGTGAGCACATCCTGCTCGCTCGTCAGGTTGGCGTTCCTGCGATCGTCGTGTTCCTGAACAAGTGCGATCAGGTTGACGATGCTGAacttctc                                                    <  1:2741506/150‑2 (MQ=2)
                                                                                                    aGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTgg                                                   <  1:2561002/150‑1 (MQ=37)
                                                                                                      aTGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGaa                                                 <  1:2237055/150‑1 (MQ=38)
                                                                                                       tGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACACGGCCGACCTGGTAGACGACGCTGAGCTGCTGGacg                                                >  1:2828477/1‑148 (MQ=37)
                                                                                                       tGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAAc                                                <  1:2766772/150‑1 (MQ=38)
                                                                                                       tGGACGGCGCGATCCTGGTTGTTTCGGCAGCTGACGGCCCGATGCCACAGACCCGTGAGCACATCCTGCTCGCTCGTCAGGTTGGCGTTCCTGCGATCGTCGTGTTCCTGAACAAGTGCGATCAGGTTGACGATGCTGAACTTCTCGagc                                                >  1:850276/1‑148 (MQ=2)
                                                                                                           cGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGt                                            >  2:3686878/1‑150 (MQ=39)
                                                                                                           cGGCGCGATCCTGGTTGTTTCGGCAGCTGACGGCCCGATGCCACAGACCCGTGAGCACATCCTGCTCGCTCGTCAGGTTGGCGTTCCTGCGATCGTCGTGTTCCTGAACAAGTGCGATCAGGTTGACGATGCTGAACTTCTCGAACTCGt                                            <  1:5284870/150‑1 (MQ=255)
                                                                                                           cGGCGCGATCCTGGTTGTTTCGGCAGCTGACGGCCCGATGCCACAGACCCGTGAGCACATCCTGCTCGCTCGTCAGGTTGGCGTTCCTGCGATCGTCGTGTTCCTGAACAAGTGCGATCAGGTTGACGATGCTGAACTTCTCGAACTCGt                                            <  1:2930077/150‑1 (MQ=255)
                                                                                                                cgATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCgaga                                       <  2:4579730/150‑1 (MQ=39)
                                                                                                                 gATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGAt                                      <  2:2230726/150‑1 (MQ=39)
                                                                                                                 gATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGAt                                      <  2:2230723/150‑1 (MQ=39)
                                                                                                                 gATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGAt                                      <  2:1052779/150‑1 (MQ=39)
                                                                                                                   tCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATgg                                    <  1:3319624/150‑1 (MQ=39)
                                                                                                                   tCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATgg                                    <  2:1998407/150‑1 (MQ=39)
                                                                                                                   tCCTGGTTGTTTCGGCAGCTGACGGCCCGATGCCACAGACCCGTGAGCACATCCTGCTCGCTCGTCAGGTTGGCGTTCCTGCGATCGTCGTGTTCCTGAACAAGTGCGATCAGGTTGACGATGCTGAACTTCTCGAACTCGTAgaactcg                                    <  1:5369862/150‑6 (MQ=255)
                                                                                                                     cTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGaa                                  <  1:4004716/150‑1 (MQ=39)
                                                                                                                     cTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGaa                                  <  2:354782/150‑1 (MQ=39)
                                                                                                                        gTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGtt                               <  1:3229547/150‑1 (MQ=39)
                                                                                                                         tttGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTTGTCGAGATGGAAGTTc                              >  2:2737200/1‑150 (MQ=255)
                                                                                                                         tttGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTc                              <  1:1767351/150‑1 (MQ=39)
                                                                                                                            gCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTcgcg                           <  2:3253137/150‑1 (MQ=39)
                                                                                                                             cTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCg                                                                                                          <  2:5212988/70‑1 (MQ=255)
                                                                                                                             cTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCg                                                                                                          >  1:5212988/1‑70 (MQ=255)
                                                                                                                             cTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGa                          <  2:2085160/150‑1 (MQ=39)
                                                                                                                             cTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGa                          >  2:1282437/1‑150 (MQ=39)
                                                                                                                              tCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGAc                         >  1:4986312/1‑150 (MQ=39)
                                                                                                                               cGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGAcc                        >  1:1626241/1‑150 (MQ=39)
                                                                                                                               cGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGAcc                        >  2:2143582/1‑150 (MQ=39)
                                                                                                                                ggCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACct                       >  1:1934870/1‑150 (MQ=39)
                                                                                                                                ggCAGCTGACGGCCCGATGCCACAGACCCGTGAGCACATCCTGCTCGCTCGTCAGGTTGGCGTTCCTGCGATCGTCGTGTTCCTGAACAAGTGCGATCAGGTTgacgatg                                                               <  2:2965963/110‑3 (MQ=255)
                                                                                                                                ggCAGCTGACGGCCCGATGCCACAGACCCGTGAGCACATCCTGCTCGCTCGTCAGGTTGGCGTTCCTGCGATCGTCGTGTTCCTGAACAAGTGCGATCAGGTTgacgatg                                                               <  2:2965964/110‑3 (MQ=255)
                                                                                                                                  ccgccgATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACctgc                     >  1:2142348/1‑150 (MQ=39)
                                                                                                                                    gccgATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACctgctg                   <  1:1638368/150‑1 (MQ=39)
                                                                                                                                      cgATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTc                 >  1:2884663/1‑150 (MQ=39)
                                                                                                                                       gATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTcc                <  1:2452379/150‑1 (MQ=39)
                                                                                                                                        aTGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCAGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCa               <  1:591917/150‑1 (MQ=255)
                                                                                                                                          ggcccGATGCCACAGACCCGTGAGCACATCCTGCTCGCTCGTCAGGTTGGCGTTCCTGCGATCGTCGTGTTCCTGAACAAGTGCGATCAGGTTGACGATGCTGAACTTCTCGAGCTCGTAGAACTCGAAGTTCGCGAACTGCTTTcgaag             >  2:5028980/4‑146 (MQ=255)
                                                                                                                                              cGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGt                                                                       <  2:4953369/88‑1 (MQ=255)
                                                                                                                                               gATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGa        <  2:2301326/150‑1 (MQ=39)
                                                                                                                                                aTGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGAc       <  2:2850856/150‑1 (MQ=39)
                                                                                                                                                    cGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTcc   >  1:1487984/1‑150 (MQ=38)
                                                                                                                                                     gCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTccc  >  2:1557656/1‑150 (MQ=38)
                                                                                                                                                     |                                                                                                                                                     
ACCGCTCACGTCGAGTACAACTCGAACATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCC  >  NC_002947/549521‑549819

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