Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,171,699 C→A 6.2% intergenic (‑551/‑230) PP_5483 ← / → PP_1921 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,171,6990CA 6.2% 1127.8 / 42.6 335intergenic (‑551/‑230)PP_5483/PP_1921hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (146/168);  new base A (9/12);  total (155/180)
Fisher's exact test for biased strand distribution p-value = 8.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 339 total aligned reads displayed.

CTCCCTCGATCAGAAAGTTACTCTCAGCAACCTATTGCAGAGGGTTTAAGGATCGCCGTTGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGGCGCGGGTGCAGATTCTGTCGAGCGAAGCAGGAGAGAGGAGGGTTGGTGTTATTTTGCTACTATGACTATGCGCACGACGGGCGGT  >  NC_002947/2171550‑2171848
                                                                                                                                                     |                                                                                                                                                     
ctcCCTCGATCAGAAAGTTACTCTCAGCAACCTATTGCAGAGGGTTTAAGGATCGCCGTTGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGc                                                                                                                                                       >  2:3274782/1‑150 (MQ=255)
 tcCCTCGATCAGAAAGTTACTCTCAGCAACCTATTGCAGAGGGTTTAAGGATCGCCGTTGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCt                                                                                                                                                      >  2:1994214/1‑150 (MQ=255)
     tCGATCAGAAAGTTACTCTCAGCAACCTATTGCAGAGGGTTTAAGGATCGCCGTTGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTcta                                                                                                                                                  <  2:5342102/150‑3 (MQ=255)
     tCGATCAGAAAGTTACTCTCAGCAACCTATTGCAGAGGGTTTAAGGATCGCCGTTGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCa                                                                                                                                                  >  1:8034296/1‑150 (MQ=255)
     tCGATCAGAAAGTTACTCTCAGCAACCTATTGCAGAGGGTTTAAGGATCGCCGTTGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCa                                                                                                                                                  <  2:3127610/150‑1 (MQ=255)
      cGATCAGAAAGTTACTCGCAGCAACCTATTGCAGAGGGTTTAAGGATCGCCGTTGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAg                                                                                                                                                 <  1:797037/150‑1 (MQ=255)
        aTCAGAAAGTTACTCTCAGCAACCTATTGCAGAGGGTTTAAGGATCGCCGTTGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGt                                                                                                                                               <  1:1852823/150‑1 (MQ=255)
         tCAGAAAGTTACTCTCAGCAACCTATTGCAGAGGGTTTAAGGATCGCCGTTGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGtt                                                                                                                                              <  1:3642516/150‑1 (MQ=255)
         tCAGAAAGTTACTCTCAGCAACCTATTGCAGAGGGTTTAAGGATCGCCGTTGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGAAGCCAGGtt                                                                                                                                              <  1:2034129/150‑1 (MQ=255)
          cAGAAAGTTACTCTCAGCAACCTATTGCAGAGGGTTTAAGGATCGCCGTTGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTg                                                                                                                                             <  1:7506553/150‑1 (MQ=255)
            gAAAGTTACTCTCAGCAACCTATTGCAGAGGGTTTAAGGATCGCCGTTGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGAAGCCAGGTTgat                                                                                                                                           <  2:602911/150‑3 (MQ=255)
               aGTTACTCTCAGCAACCTATTGCAGAGGGTTTAAGGATCGCCGTTGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGAAGCCAGGTTgattga                                                                                                                                        <  2:3941066/150‑6 (MQ=255)
                   aCTCTCAGCAACCTATTGCAGAGGGTTTAAGGATCGCCGTTGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCgta                                                                                                                                    <  1:6383701/150‑1 (MQ=255)
                    ctctCAGCAACCTATTGCAGAGGGTTTAAGGATCGCCGTTGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCgtag                                                                                                                                   <  2:7644948/150‑1 (MQ=255)
                     tctcAGCAACCTATTGCAGAGGGTTTAAGGATCGCCGTTGGCCATTGTGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCgtagg                                                                                                                                  >  1:3782267/1‑150 (MQ=255)
                     tctcAGCAACCTATTGCAGAGGGTTTAAGGATCGCCGTTGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCgtagg                                                                                                                                  >  1:6662749/1‑150 (MQ=255)
                     tctcAGCAACCTATTGCAGAGGGTTTAAGGATCGCCGTTGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCgtagg                                                                                                                                  <  1:6131095/150‑1 (MQ=255)
                     tctcAGCAACCTAGTGCAGAGGGTTTAAGGATCGCCGTGGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCgtagg                                                                                                                                  <  1:7682700/150‑1 (MQ=255)
                       tcAGCAACCTATTGCAGAGGGTTTAAGGATCGCCGTTGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCgtaggga                                                                                                                                <  2:6470225/150‑1 (MQ=255)
                         aGCAACCTATTGCAGAGGGTTTAAGGATCGCCGTTGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCgtagggagt                                                                                                                              >  2:841752/1‑150 (MQ=255)
                             aCCTATTGCAGAGGGTTTAAGGATCGCCGTTGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCgtagggagtaggg                                                                                                                          <  2:7963509/150‑1 (MQ=255)
                              ccTATTGCAGAGGGTTTAAGGATCGCCGTTGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCgtagggagtaggga                                                                                                                         >  1:3883526/1‑150 (MQ=255)
                              ccTATTGCAGAGGGTTTAAGGATCGCCGTTGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCgtagggagtaggga                                                                                                                         >  1:477079/1‑150 (MQ=255)
                              ccTATTGCAGAGGGTTTAAGGATCGCCGTTGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCgtagggagtaggga                                                                                                                         >  1:7317769/1‑150 (MQ=255)
                               cTATTGCAGAGGGTTTAAGGATCGCCGTTGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCgtagggagtagggag                                                                                                                        <  2:4290953/150‑1 (MQ=255)
                                tATTGCAGAGGGTTTAAGGATCGCCGTTGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCgtagggagtagggagt                                                                                                                       >  1:4479322/1‑150 (MQ=255)
                                tATTGCAGAGGGTTTAAGGATCGCCGTTGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCgtagggagtagggagt                                                                                                                       >  1:2569442/1‑150 (MQ=255)
                                 aTTGCAGAGGGTTTAAGGATCGCCGTTGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCgtagggagtagggagta                                                                                                                      <  2:2258403/150‑1 (MQ=255)
                                  ttGCAGAGGGTTTAAGGATCGCCGTTGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCgtagggagtagggagtag                                                                                                                     <  2:7443491/150‑1 (MQ=255)
                                     cAGAGGGTTTAAGGATCGCCGTTGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCgtagggagtagggagtaggga                                                                                                                  <  2:3202554/150‑1 (MQ=255)
                                     cAGAGGGTTTAAGGATCGCCGTTGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCgtagggagtagggagtaggga                                                                                                                  >  2:2862672/1‑150 (MQ=255)
                                          gttttAAGGTTCGCCTTTGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGATTCCAGGTTGCAGGTCgtagggagtagggagtagggagtatg                                                                                                             <  1:5676534/148‑3 (MQ=255)
                                          ggTTTAAGGATCGCCGATGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCgtagggagtagggagtagggagtagg                                                                                                             <  1:3576006/150‑1 (MQ=255)
                                          ggTTTAAGGAGCGCCGTGGGCCATTAAGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCgtagggagtagggagtagggagtagg                                                                                                             <  1:3808076/150‑1 (MQ=255)
                                               aaGGATCGCCGTTGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCGTAgggagtagggagtagggagtagggagta                                                                                                        >  2:4938124/1‑150 (MQ=255)
                                               aaGGATCGCCGTTGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCGTAgggagtagggagtagggagtagggagta                                                                                                        <  2:910501/150‑1 (MQ=255)
                                                  gATCGCCGTTGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCGTAGGGagtagggagtagggagtagggagtaggg                                                                                                     >  1:8720427/1‑150 (MQ=255)
                                                  gATCGCCGTTGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCGTAGGGagtagggagtagggagtagggagtaggg                                                                                                     <  2:2542309/150‑1 (MQ=255)
                                                       ccGTTGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAgg                                                                                                >  1:8444161/1‑150 (MQ=255)
                                                       ccGTTGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAgg                                                                                                <  2:7296356/150‑1 (MQ=255)
                                                        cGTTGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAggg                                                                                               >  2:1676537/1‑150 (MQ=255)
                                                        cGTTGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAggg                                                                                               <  1:2137266/150‑1 (MQ=255)
                                                              ccATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTc                                                                                         >  2:4541514/1‑150 (MQ=255)
                                                                 ttATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGgc                                                                                      >  2:8722785/1‑150 (MQ=255)
                                                                 ttATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGgc                                                                                      >  2:3561449/1‑150 (MQ=255)
                                                                   aTGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGgcgc                                                                                    <  1:5771965/150‑1 (MQ=255)
                                                                     gtgAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGGCGCgg                                                                                  <  2:7846749/150‑1 (MQ=255)
                                                                     gtgAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGGCGCgg                                                                                  >  1:2805156/1‑150 (MQ=255)
                                                                     gtgAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCGGAGGGAGTCGGGAGTAGGGAGTCGGGAGTCGGGGCAGGGAGAGTCGGCGCgg                                                                                  >  2:5737579/1‑150 (MQ=255)
                                                                       gAAACGACCCGCTCCTGGAGCGGGATGATCGGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGGCGCGGGt                                                                                >  1:2243885/1‑150 (MQ=255)
                                                                           catcccGCTCCAGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCGTAGGGAGTAGGGGGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGGCGCGGGTGCAg                                                                            >  1:6316163/4‑150 (MQ=255)
                                                                            gACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGGCGCGGGTGCAGa                                                                           >  1:3696984/1‑150 (MQ=255)
                                                                            gACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGGCGCGGGTGCAGa                                                                           <  1:153608/150‑1 (MQ=255)
                                                                             tcccGCTCCAGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGGCGCGGGTGCAGAt                                                                          >  1:5009301/2‑150 (MQ=255)
                                                                             tcccGCTCCAGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGGCGCGGGTGCAGAt                                                                          >  1:297527/2‑150 (MQ=255)
                                                                             aCCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGGCGCGGGTGCAGAt                                                                          <  2:3626948/150‑1 (MQ=255)
                                                                              cccGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGGCGCGGGTGCAGAtt                                                                         <  2:1259988/150‑1 (MQ=255)
                                                                              cccGCTCCAGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCGGAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGGCGCGGGTGCAGAtt                                                                         >  1:5524179/1‑150 (MQ=255)
                                                                               ccGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGGCGCGGGTGCAGATTc                                                                        <  1:4520168/150‑1 (MQ=255)
                                                                                cGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGGCGCGGGTGCAGATTCt                                                                       >  2:3784123/1‑150 (MQ=255)
                                                                                  cTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGGCGCGGGTGCAGATTCTGt                                                                     <  1:4932781/150‑1 (MQ=255)
                                                                                   tCCTGGAGCGGGATGATCCGAGGGCTGTTGGTGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGGCGCGGGTGCAGATTCTGTc                                                                    <  1:3315959/150‑1 (MQ=255)
                                                                                   tCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGGCGCGGGTGCAGATTCTGTc                                                                    <  1:978005/150‑1 (MQ=255)
                                                                                   tCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGGCGCGGGTGCAGATTCTGTc                                                                    >  2:699610/1‑150 (MQ=255)
                                                                                     cTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGGCGCGGGTGCAGATTCTGTCGa                                                                  <  2:3302263/150‑1 (MQ=255)
                                                                                     cTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGGCGCGGGTGCAGATTCTGTCGa                                                                  <  2:6637694/150‑1 (MQ=255)
                                                                                      gggAGCGGGATGCTCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGGCGCGGGTGCAGATTCTGTCGAg                                                                 <  2:7223810/149‑1 (MQ=255)
                                                                                      gggAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGGCGCGGGTGCAGATTCTGTCGAg                                                                 <  1:839829/149‑1 (MQ=255)
                                                                                          gCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGGCGCGGGTGCAGATTCTGTCGAGCGaa                                                             <  1:7406770/150‑1 (MQ=255)
                                                                                              gatgatCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGGCGCGGGTGCAGATTCTGTCGAGCGAAGCAg                                                         >  1:1121546/1‑150 (MQ=255)
                                                                                                   tCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGGCGCGGGTGCAGATTCTGTCGAGCGAAGCAGgagag                                                    >  1:1241314/1‑150 (MQ=255)
                                                                                                    ccGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGGCGCGGGTGCAGATTCTGTCGAGCGAAGCAGgagaga                                                   >  1:8341034/1‑150 (MQ=255)
                                                                                                      gAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGGCGCGGGTGCAGATTCTGTCGAGCGAAGCAGGAGAgagg                                                 >  2:3143123/1‑150 (MQ=255)
                                                                                                       aGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGGCGCGGGTGCAGATTCTGTCGAGCGAAGCAGGAGAgagga                                                <  1:5257790/150‑1 (MQ=255)
                                                                                                       aGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGGCGCGGGTGCAGATTCTGTCGAGCGAAGCAGGAGAgagga                                                <  2:5080114/150‑1 (MQ=255)
                                                                                                          gCTGTTGGTGTTTCGCTGTTTGTGTGATGATCATTCTACCTGGCTTACTTGTTGCAGGTCGTAGGGATTAGGGAGTTGGGAGTTGGTAGTATGTGCAGGGTGATTCGGCGCGTGTGCAGATTCTGTCGAGCGAAGCAGGTGAgaggatgt                                             <  2:1234232/150‑4 (MQ=255)
                                                                                                          gCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGGCGCGGGTGCAGATTCTGTCGAGCGAAGCAGGAGAGAGGAggg                                             <  1:3305454/150‑1 (MQ=255)
                                                                                                            tGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGGCGCGGGTGCAGATTCTGTCGAGCGAAGCAGGAGAGAGGAGGGtt                                           <  1:1769294/150‑1 (MQ=255)
                                                                                                               tGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCGTAGGGAGTAGGGTGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGGCGCGGGTGCAGATTCTGTCGAGCGAAGCAGGAGAGAGGAGGGTTGgt                                        <  2:6316163/150‑1 (MQ=255)
                                                                                                                        gCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGGCGCGGGTGCAGATTCTGTCGAGCGAAGCAGGAGAGAGGAGGGTTGGTGTTATTTTg                               >  1:5847230/1‑150 (MQ=255)
                                                                                                                        gCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGGCGCGGGTGCAGATTCTGTCGAGCGAAGCAGGAGAGAGGAGGGTTGGTGTTATTTTg                               <  2:1783283/150‑1 (MQ=255)
                                                                                                                         cTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGGCGCGGGTGCAGATTCTGTCGAGCGAAGCAGGAGAGAGGAGGGTTGGTGTTATTTTGc                              <  1:5925986/150‑1 (MQ=255)
                                                                                                                          tGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGGCGCGGGTGCAGATTCTGTCGAGCGAAGCAGGAGAGAGGAGGGTTGGTGTTATTTTGct                             <  1:5440719/150‑1 (MQ=255)
                                                                                                                          tGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGGCGCGGGTGCAGATTCTGTCGAGCGAAGCAGGAGAGAGGAGGGTTGGTGTTATTTTGct                             >  1:4686622/1‑150 (MQ=255)
                                                                                                                              tGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGGCGCGGGTGCAGATTCTGTCGAGCGAAGCAGGAGAGAGGAGGGTTGGTGTTATTTTGctacta                         <  1:3993620/150‑1 (MQ=255)
                                                                                                                               gAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGGCGCGGGTGCAGATTCTGTCGAGCGAAGCAGGAGAGAGGAGGGTTGGTGTTATTTTGCTactat                        <  2:6522674/150‑1 (MQ=255)
                                                                                                                                     aGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGGCGCGGGTGCAGATTCTGTCGAGCGAAGCAGGAGAGAGGAGGGTTGGTGTTATTTTGCTactatgactat                  >  2:4433765/1‑150 (MQ=255)
                                                                                                                                       cTCAATCAACCTGGCTTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGGCGCGGGTGCAGATTCTGTCGAGCGAAGCAGGAGAGAGGAGGGTTGGTGTTATTTTGCTACTATGACTATgc                >  2:402525/1‑150 (MQ=255)
                                                                                                                                        tcaatcaaCCTGGCTTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGGCGCGGGTGCAGATTCTGTCGAGCGAAGCAGGAGAGAGGAGGGTTGGTGTTATTTTGCTACTATGACTATgcg               >  1:1480221/1‑150 (MQ=255)
                                                                                                                                         caatcaaCCTGGCTTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGGCGCGGGTGCAGATTCTGTCGAGCGAAGCAGGAGAGAGGAGGGTTGGTGTTATTTTGCTACTATGACTATgcgc              <  2:7524413/150‑1 (MQ=255)
                                                                                                                                          ctgcaaCCTGGAAGCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGGCGCGGGTGCAGATTCTGTCGAGCGAAGCAGGAGAGAGGAGGGTTGGTGTTATTTTGCTACTATGACTATGCGCa             >  2:2336544/4‑150 (MQ=255)
                                                                                                                                          ctgcaaCCTGGAAGCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGGCGCGGGTGCAGATTCTGTCGAGCGAAGCAGGAGAGAGGAGGGTTGGTGTTATTTTGCTACTATGACTATGCGCa             >  2:8256353/4‑150 (MQ=255)
                                                                                                                                          ctgcaaCCTGGAAGCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGGCGCGGGTGCAGATTCTGTCGAGCGAAGCAGGAGAGAGGAGGGTTGGTGTTATTTTGCTACTATGACTATGCGCa             >  2:1392850/4‑150 (MQ=255)
                                                                                                                                          aatcaaCCTGGCTTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGGCGCGGGTGCAGATTCTGTCGAGCGAAGCAGGAGAGAGGAGGGTTGGTGTTATTTTGCTACTATGACTATGCGCa             >  2:2060262/1‑150 (MQ=255)
                                                                                                                                          aagccaCCTGGCTTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTCGGGGCAGGGAGAGTCGGCGCGGGTGCAGATTCTGTCGAGCGAAGCAGGAGAGAGGAGGGTTGGTGTTATTTTGCTACTATGACTATGCGCa             >  2:8363298/6‑150 (MQ=255)
                                                                                                                                           tgcaaCCTGGAAGCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGGCGCGGGTGCAGATTCTGTCGAGCGAAGCAGGAGAGAGGAGGGTTGGTGTTATTTTGCTACTATGACTATGCGCac            >  1:3145334/3‑150 (MQ=255)
                                                                                                                                           tgcaaCCTGGAAGCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGGCGCGGGTGCAGATTCTGTCGAGCGAAGCAGGAGAGAGGAGGGTTGGTGTTATTTTGCTACTATGACTATGCGCac            >  1:3001173/3‑150 (MQ=255)
                                                                                                                                            tcaaCCTGGCTTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGGCGCGGGTGCAGATTCTGTCGAGCGAAGCAGGAGAGAGGAGGGTTGGTGTTATTTTGCTACTATGACTATGCGCacg           >  1:4324578/1‑150 (MQ=255)
                                                                                                                                                    gCTTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGGCGCGGGTGCAGATTCTGTCGAGCGAAGCAGGAGAGAGGAGGGTTGGTGTTATTTTGCTACTATGACTATGCGCACGACGGGCgg   <  1:6987980/150‑1 (MQ=255)
                                                                                                                                                     cTTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGGCGCGGGTGCAGATTCTGTCGAGCGAAGCAGGAGAGAGGAGGGTTGGTGTTATTTTGCTACTATGACTATGCGCACGACGGGCGGt  <  1:6515411/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTCCCTCGATCAGAAAGTTACTCTCAGCAACCTATTGCAGAGGGTTTAAGGATCGCCGTTGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGCTTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGGCGCGGGTGCAGATTCTGTCGAGCGAAGCAGGAGAGAGGAGGGTTGGTGTTATTTTGCTACTATGACTATGCGCACGACGGGCGGT  >  NC_002947/2171550‑2171848

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