Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,236,982 C→T 10.7% intergenic (+15/+23) ybaK → / ← PP_1078 Cys‑tRNAPro deacylase/ABC transporter ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,236,9820CT10.7% 490.5 / 38.9 177intergenic (+15/+23)ybaK/PP_1078Cys‑tRNAPro deacylase/ABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (84/74);  new base T (5/14);  total (89/88)
Fisher's exact test for biased strand distribution p-value = 3.06e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

Only 100 of 190 total aligned reads displayed.

TTTATTGATGACTCGGCGCAGAATTTTGCAACTATCCATGTCAGCGCCGGCAGGCGTGGCCTGGAAGTGGAGTTGGCCGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGTGGCAACTGTGCCGCAACGCTACCCGGTACCGGGAATAGCACGAGGTGATCGGCCTTCACTGCAATCTCCACATTCTCACCGATC  >  NC_002947/1236833‑1237131
                                                                                                                                                      |                                                                                                                                                     
 ttATTGATGACTCGGCGCAGAATTTTGCAACTATCCATGTCAGCGCCGGCAGGCGTGGCCTGGAAGTGGAGTTGGCCGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCa                                                                                                                                                      >  2:4739378/1‑150 (MQ=255)
  tATTGATGACTCGGCGCAGAATTTTGCAACTATCCATGTCAGCGCCGGCAGGCGTGGCCTGGAAGTGGAGTTGGCCGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCaa                                                                                                                                                     >  1:684410/1‑150 (MQ=255)
   aTTGATGACTCGGCGCAGAATTTTGCAACTATCCATGTCAGCGCCGGCAGGCGTGGCCTGGAAGTGGAGTTGGCCGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAg                                                                                                                                                    >  2:3100482/1‑150 (MQ=255)
   aTTGATGACTCGGCGCAGAATTTTGCAACTATCCATGTCAGCGCCGGCAGGCGTGGCCTGGAAGTGGAGTTGGCCGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAg                                                                                                                                                    >  2:3592433/1‑150 (MQ=255)
       atgaCTCGGCGCAGAATTTTGCAACTATCCATGTCAGCGCCGGCAGGCGTGGCCTGGAAGTGGAGTTGGCCGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGc                                                                                                                                                <  1:2197574/150‑1 (MQ=255)
       atgaCTCGGCGCAGAATTTTGCAACTATCCATGTCAGCGCCGGCAGGCGTGGCCTGGAAGTGGAGTTGGCCGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGc                                                                                                                                                >  2:480998/1‑150 (MQ=255)
           cTCGGCGCAGAATTTTGCAACTATCCATGTCAGCGCCGGCAGGCGTGGCCTGGAAGTGGAGTTGGCCGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCctc                                                                                                                                            >  2:4566254/1‑150 (MQ=255)
            tCGGCGCAGAATTTTGCAACTATCCATGTCAGCGCCGGCAGGCGTGGCGTGGAAGTGGAGTTGGCCGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCctct                                                                                                                                           >  2:4127007/1‑150 (MQ=255)
            tCGGCGCAGAATTTTGCAACTATCCATGTCAGCGCCGGCAGGCGTGGCCTGGAAGTGGAGTTGGCCGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCTTGGCGGCCccct                                                                                                                                           <  2:122344/150‑1 (MQ=255)
            tCGGCGCAGAATTTTGCAACTATCCATGTCAGCGCCGGCAGGCGTGGCCTGGAAGTGGAGTTGGCCGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCctct                                                                                                                                           <  1:827004/150‑1 (MQ=255)
            tCGGCGCAGAATTTTGCAACTATCCATGTCAGCGCCGGCAGGCGTGGCCTGGAAGTGGAGTTGGCCGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCctct                                                                                                                                           >  1:2128062/1‑150 (MQ=255)
            tCGGCGCAGAATTTTGCAACTATCCATGTCAGCGCCGGCAGGCGTGGCCTGGAAGTGGAGTTGGCCGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGAGGGCCGCTTGGCGGCCccct                                                                                                                                           <  2:5093980/150‑1 (MQ=255)
            tCGGCGCAGAATTTTGCAACTATCCATGTCAGCGCCGGCAGGCGTGGCCTGGAAGTGGAGTTGGCCGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGAGGGCCGCTTGGCGGCCccct                                                                                                                                           <  1:5244357/150‑1 (MQ=255)
                 gcAGAATTTTGCAACTATCCATGTCAGCGCCGGCAGGCGTGGCCTGGAAGTGGAGTTGGCCGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTc                                                                                                                                      >  2:5137243/1‑150 (MQ=255)
                 gcAGAATTTTGCAACTATCCATGTCAGCGCCGGCAGGCGTGGCCTGGAAGTGGAGTTGGCCGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTc                                                                                                                                      >  1:3747433/1‑150 (MQ=255)
                    gAATTTTGCAACTATCCATGTCAGCGCCGGCAGGCGTGGCCTGGAAGTGGAGTTGGCCGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACCAGAGGGCCGCTTGGCGGCCCCCTgttctca                                                                                                                                    <  1:1653100/150‑7 (MQ=38)
                     caTTTTGCAACTCGCCATGTCAGCGCCGGCAGGCGTGGCCTGGAAGTGGAGTTGGCCGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGaa                                                                                                                                  <  2:4680326/149‑1 (MQ=255)
                     aaTTTTGCAACTATCCATGTCAGCGCCGGCAGGCGTGGCCTGGAAGTGGAGTTGGCCGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGaa                                                                                                                                  <  2:5343765/150‑1 (MQ=255)
                     aaTTTTGCAACTATCCATGTCAGCGCCGGCAGGCGTGGCCTGGAAGTGGAGTTGGCCGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGaa                                                                                                                                  <  2:2013480/150‑1 (MQ=255)
                     aaTTTTGCAACTATCCATGTCAGCGCCGGCAGGCGTGGCCTGGAAGTGGAGTTGGCCGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGaa                                                                                                                                  <  2:684410/150‑1 (MQ=255)
                     aaTTTTGCAACTATCCATGTCAGCGCCGGCAGGCGTGGCCTGGAAGTGGAGTTGGCCGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGaa                                                                                                                                  <  1:4421369/150‑1 (MQ=255)
                     aaTTTTGCAACTATCCATGTCAGCGCCGGCAGGCGTGGCCTGGAAGTGGAGTTGGCCGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGAGGGCCGCTTGGCGGCCccct                                                                                                                                           <  2:4634833/141‑1 (MQ=255)
                     aaTTTTGCAACTATCCATGTCAGCGCCGGCAGGCGTGGCCTGGAAGTGGAGTTGGCCGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGAGGGCCGCTTGGCGGCCccct                                                                                                                                           >  1:4634892/1‑141 (MQ=255)
                       ttttGCAACTATCCATGTCAGCGCCGGCAGGCGTGGCCTGGAAGTGGAGTTGGCCGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGaagt                                                                                                                                >  2:3221876/1‑150 (MQ=255)
                       ttttGCAACTATCCATGTCAGCGCCGGCAGGCGTGGCCTGGAAGTGGAGTTGGCCGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGaagt                                                                                                                                >  2:494558/1‑150 (MQ=255)
                          tGCAACTATCCATGTCAGCGCCGGCAGGCGTGGCCTGGAAGTGGAGTTGGCCGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGaagtcaa                                                                                                                             <  2:255687/150‑1 (MQ=255)
                          tGCAACTATCCATGTCAGCGCCGGCAGGCGTGGCCTGGAAGTGGAGTTGGCCGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGaagtcaa                                                                                                                             <  2:862614/150‑1 (MQ=255)
                             aaCTATCCATGTCAGCGCCGGCAGGCGTGGCCTGGAAGTGGAGTTGGCCGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGaagtcaagtc                                                                                                                          >  1:437316/1‑150 (MQ=255)
                             aaCTATCCATGTCAGCGCCGGCAGGCGTGGCCGGGAAGTGGAGTTGGCCGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGaagtcaagtc                                                                                                                          >  1:1481342/1‑150 (MQ=255)
                               cTATCCATGTCAGCGCCGGCAGGCGTGGCCTGGAAGTGGAGTTGGCCGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGc                                                                                                                        >  1:2259818/1‑150 (MQ=255)
                                tATCCATGTCAGCGCCGGCAGGCGTGGCCTGGAAGTGGAGTTGGCCGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGcc                                                                                                                       <  1:4279784/150‑1 (MQ=255)
                                 aTCCATGTCAGCGCCGGCAGGCGTGGCCTGGAAGTGGAGTTGGCCGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGt                                                                                                                                        <  2:3818779/132‑1 (MQ=255)
                                 aTCCATGTCAGCGCCGGCAGGCGTGGCCTGGAAGTGGAGTTGGCCGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGt                                                                                                                                        >  1:3818779/1‑132 (MQ=255)
                                  tCCATGTCAGCGCCGGCAGGCGTGGCCTGGAAGTGGAGTTGGCCGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGa                                                                                                                     >  1:3157997/1‑150 (MQ=255)
                                    cATGTCAGCGCCGGCAGGCGTGGCCTGGAAGTGGAGTTGGCCGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAAc                                                                                                                   >  1:1302137/1‑150 (MQ=255)
                                    cATGTCAGCGCCGGCAGGCGTGGCATGGAAGTGGAGTTGGCCGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAAc                                                                                                                   >  1:502840/1‑150 (MQ=255)
                                      tGTCACAGCCGGAAGTCGTGGCATGTAAGTGGAGTTGGCCGCAACGGTCCTCGCGGAGTATAAGCAGGGCACGTTTGCCGGGATTGGTCGATGCTGAGAACAGGGGGCCGCCACGCGGCCCTCTGGTTCAGAAGTCAAGTCGCCGAACTg                                                                                                                 >  2:4150050/1‑150 (MQ=255)
                                        tCAGCGCCGGCAGGCGTGGCCTGGAAGTGGAGTTGGCCGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTa                                                                                                               <  1:1922965/150‑1 (MQ=255)
                                          aGCGCCGGCAGGCGTGGCCTGGAAGTGGAGTTGGCCGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATc                                                                                                             <  2:1826125/150‑1 (MQ=255)
                                           gcgcCGGCAGGCGTGGCCTGGAAGTGGAGTTGGCCGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCg                                                                                                            >  2:5181731/1‑150 (MQ=255)
                                             gcCGGCAGGCGTGGCCTGGAAGTGGAGTTGGCCGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGac                                                                                                          <  2:3459436/150‑1 (MQ=255)
                                                   aGGCGTGGCCTGGAAGTGGAGTTGGCCGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGacgcacgc                                                                                                    <  2:437316/150‑1 (MQ=255)
                                                       gTGGCCTGGAAGTGGAGTTGGCCGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCAtt                                                                                                <  1:5094974/150‑1 (MQ=255)
                                                       gTGGCCTGGAAGTGGAGTTGGCCGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCAtt                                                                                                <  1:5094875/150‑1 (MQ=255)
                                                        tGGCCTGGAAGTGGAGTTGGCCGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATtc                                                                                               <  1:4495328/150‑1 (MQ=255)
                                                         ggCCTGGAAGTGGAGTTGGCCGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATtct                                                                                              >  2:2234554/1‑150 (MQ=255)
                                                             tGGAAGTGGAGTTGGCCGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTg                                                                                          >  2:2101880/1‑150 (MQ=255)
                                                                aaGTGGAGTTGGCCGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCg                                                                                       >  1:4962277/1‑150 (MQ=255)
                                                                  gTGGAGTTGGCCGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGTg                                                                                     <  2:2333118/150‑1 (MQ=255)
                                                                      aGTTGGCCGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGTGGCaa                                                                                 <  2:3793701/150‑1 (MQ=255)
                                                                        ttGGCCGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGTGGCAACt                                                                               <  1:3681789/150‑1 (MQ=255)
                                                                        ttGGCCGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGTGGCAACt                                                                               >  2:852806/1‑150 (MQ=255)
                                                                        ttGGCCGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGTGGCAACt                                                                               >  2:1735218/1‑150 (MQ=255)
                                                                        ttGGCCGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGTGGCAACt                                                                               >  2:4206142/1‑150 (MQ=255)
                                                                        ttGGCCGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGTGGCAACt                                                                               <  1:3681924/150‑1 (MQ=255)
                                                                          ggCCGCAGCGGTGCTGGCGGAGTATTCGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGTGGCAACtgt                                                                             <  2:2903324/150‑1 (MQ=255)
                                                                          ggCCGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGTGGCAACtgt                                                                             <  1:4291855/150‑1 (MQ=255)
                                                                          ggCCGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGTGGCAACtgt                                                                             >  2:3433599/1‑150 (MQ=255)
                                                                            ccGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGTGGCAACtgtg                                                                            >  2:4455401/1‑149 (MQ=255)
                                                                            ccGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGTGGCAACTGTGc                                                                           <  1:2166657/150‑1 (MQ=255)
                                                                            ccGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGTGGCAACTGTGc                                                                           <  1:108531/150‑1 (MQ=255)
                                                                            ccGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGTGGCAACTGTGc                                                                           >  1:4099202/1‑150 (MQ=255)
                                                                            ccGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGTGGCAACTGTGc                                                                           <  1:5201881/150‑1 (MQ=255)
                                                                             cGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTg                                                                                          <  1:3146591/134‑1 (MQ=255)
                                                                             cGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTg                                                                                          >  2:3146591/1‑134 (MQ=255)
                                                                             cGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGTGGCAACTGTGcc                                                                          >  1:1780305/1‑150 (MQ=255)
                                                                                 gCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGTGGCAACTGTGCCGCaa                                                                      <  1:355719/150‑1 (MQ=255)
                                                                                  cGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGTGGCAACTGTGCCGCAAc                                                                     <  2:1397936/150‑1 (MQ=255)
                                                                                    gTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGGCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGTGGCAACTGTGCCGCAACGc                                                                   >  1:4707876/1‑150 (MQ=255)
                                                                                       cTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGAGGGCCGCTTGGCGGcccc                                                                                                                                             <  2:466730/73‑2 (MQ=38)
                                                                                        tGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGTGGCAACTGTGCCGCAACGCTAcc                                                               >  1:2500965/1‑150 (MQ=255)
                                                                                        tGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGTGGCAACTGTGCCGCAACGCTAcc                                                               >  2:2574415/1‑150 (MQ=255)
                                                                                         ggCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGTGGCAACTGTGCCGCAACGCTAccc                                                              <  2:2743472/150‑1 (MQ=255)
                                                                                         ggCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGTGGCAACTGTGCCGCAACGCTAccc                                                              <  1:4686009/150‑1 (MQ=255)
                                                                                            ggAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGTGGCAACTGTGCCGCAACGCTACCCGGt                                                           >  1:610621/1‑150 (MQ=255)
                                                                                              aGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGTGGCAACTGTGCCGCAACGCTACCCGGTAc                                                         <  2:456021/150‑1 (MQ=255)
                                                                                                 ataCGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGt                                                                                      >  2:3727642/1‑118 (MQ=255)
                                                                                                 ataCGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGTGGCAACTGTGCCGCAACGCTACCCGGTACCgg                                                      <  1:3069194/150‑1 (MQ=255)
                                                                                                 ataCGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGTGGCAACTGTGCCGCAACGCTACCCGGTACCgg                                                      >  2:1440549/1‑150 (MQ=255)
                                                                                                 ataCGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGTGGCAACTGTGCCGCAACGCTACCCGGTACCgg                                                      <  1:4857480/150‑1 (MQ=255)
                                                                                                    cGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGTGGCAACTGTGCCGCAACGCTACCCGGTACCGGGaa                                                   >  2:842140/1‑150 (MQ=255)
                                                                                                      cAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGTGGCAACTGTGCCGCAACGCTACCCGGTACCGGGAATa                                                 <  2:2689897/150‑1 (MQ=255)
                                                                                                         ggCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGTGGCAACTGTGCCGCAACGCTACCCGGTACCGGGAATAGCa                                              <  1:3100482/150‑1 (MQ=255)
                                                                                                                    cGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGTGGCAACTGTGCCGCAACGCTACCCGGTACCGGGAATAGCACGAGGTGATCg                                   >  1:1385614/1‑150 (MQ=255)
                                                                                                                    cGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTCGCGTGGCAACTGTGCCGCAACGCTACCCGGTACCGGGAATAGCACGAGGTGATCg                                   <  2:1327392/150‑1 (MQ=255)
                                                                                                                       gATTGGTCGAGGCTGAGAACGGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGTGGCAACTGTGCCGCAACGCTACCCGGTACCGGGAATAGCACGAGGTGATCGGcc                                >  2:4829766/1‑150 (MQ=255)
                                                                                                                       gATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGTGGCAACTGTGCCGCAACGCTACCCGGTACCGGGAATAGCACGAGGTGATCGGcc                                >  2:1506160/1‑150 (MQ=255)
                                                                                                                          tGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGTGGCAACTGTGCCGCAACGCTACCCGGTACCGGGAATAGCACGAGGTGATCGGCCTTc                             >  2:588767/1‑150 (MQ=255)
                                                                                                                           ggTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGTGGCAACTGTGCCGCAACGCTACCCGGTACCGGGAATAGCACGAGGTGATCGGCCTTCa                            <  2:3635374/150‑1 (MQ=255)
                                                                                                                                     gagaACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGTGGCAACTGTGCCGCAACGCTACCCGGTACCGGGAATAGCACGAGGTGATCGGCCTTCACTGCAAtctc                  >  2:4438688/1‑150 (MQ=255)
                                                                                                                                        aaCAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGTGGCAACTGTGCCGCAACGCTACCCGGTACCGGGAATAGCACGAGGTGATCGGCCTTCACTGCAATCTCcac               >  2:2791124/1‑150 (MQ=255)
                                                                                                                                        aaCAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGTGGCAACTGTGCCGCAACGCTACCCGGTACCGGGAATAGCACGAGGTGATCGGCCTTCACTGCAATCTCcac               <  1:4127019/150‑1 (MQ=255)
                                                                                                                                         aCAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGTGGCAACTGTGCCGCAACGCTACCCGGTACCGGGAATAGCACGAGGTGATCGGCCTTCACTGCAATCTCcaca              <  1:480998/150‑1 (MQ=255)
                                                                                                                                           cAGAGGGCCGCTTGGCGGCCCCCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGTGGCAACTGTGCCGCAACGCTACCCGGTACCGGGAATAGCACGAGGTGATCGGCCTTCACTGCAATCTCCACAt             >  1:4671404/1‑150 (MQ=255)
                                                                                                                                               gggCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGTGGCAACTGTGCCGCAACGCTACCCGGTACCGGGAATAGCACGAGGTGATCGGCCTTCACTGCAATCTCCACATtctc         >  1:2067021/1‑150 (MQ=255)
                                                                                                                                               gggCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGTGGCAACTGTGCCGCAACGCTACCCGGTACCGGGAATAGCACGAGGTGATCGGCCTTCACTGCAATCTCCACATtctc         >  2:3235094/1‑150 (MQ=255)
                                                                                                                                                 gccgccAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGTGGCAACTGTGCCGCAACGCTACCCGGTACCGGGAATAGCACGAGGTGATCGGCCTTCACTGCAATCTCCACATTCTCAc       >  1:1268417/1‑150 (MQ=255)
                                                                                                                                                   cgccAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGTGGCAACTGTGCCGCAACGCTACCCGGTACCGGGAATAGCACGAGGTGATCGGCCTTCACTGCAATCTCCACATTCTCACCg     >  1:218157/1‑150 (MQ=255)
                                                                                                                                                    gccAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGTGGCAACTGTGCCGCAACGCTACCCGGTACCGGGAATAGCACGAGGTGATCGGCCTTCACTGCAATCTCCACATtct          <  1:4522240/144‑1 (MQ=255)
                                                                                                                                                      cAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGTGGCAACTGTGCCGCAACGCTACCCGGTACCGGGAATAGCACGAGGTGATCGGCCTTCACTGCAATCTCCACATTCTCACCGAtc  <  1:5137243/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
TTTATTGATGACTCGGCGCAGAATTTTGCAACTATCCATGTCAGCGCCGGCAGGCGTGGCCTGGAAGTGGAGTTGGCCGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGTGGCAACTGTGCCGCAACGCTACCCGGTACCGGGAATAGCACGAGGTGATCGGCCTTCACTGCAATCTCCACATTCTCACCGATC  >  NC_002947/1236833‑1237131

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: