Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,504,076 T→A 6.0% intergenic (+34/‑229) rpsI → / → petA 30S ribosomal protein S9/ubiquinol‑cytochrome c reductase iron‑sulfur subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,504,0760TA 6.0% 819.2 / 23.9 257intergenic (+34/‑229)rpsI/petA30S ribosomal protein S9/ubiquinol‑cytochrome c reductase iron‑sulfur subunit
Reads supporting (aligned to +/- strand):  ref base T (119/120);  new base A (6/10);  total (126/131)
Fisher's exact test for biased strand distribution p-value = 4.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 274 total aligned reads displayed.

CCCTGCGTGGCGCTCTGCGTCGTGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAGTGTGTGGATTTGTTTCTTACCATGACG  >  NC_002947/1503927‑1504225
                                                                                                                                                     |                                                                                                                                                      
   tGCGTGGCGCTCTGCGTCGTGCTGGCTACGTCACCCGCGACGCTCGTGTTGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCTGTCGTAAAAAGCCCGGTTGCTTgg                                                                                                                                                     <  2:5345635/150‑1 (MQ=255)
   tGCGTGGCGCTCTGCGTCGTGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTgg                                                                                                                                                     <  2:6238547/150‑1 (MQ=255)
        ggCGCTCTGCGTCGTGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCg                                                                                                                                                <  1:7352181/150‑1 (MQ=255)
         gcgcTCTGCGTCGTGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCgg                                                                                                                                               <  2:6634762/150‑1 (MQ=255)
         gcgcTCTGCGTCGTGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCgg                                                                                                                                               <  2:1293380/150‑1 (MQ=255)
         gcgcTCTGCGTCGTGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCAAGGAACCgg                                                                                                                                               <  2:87690/150‑1 (MQ=255)
         gcgcTCTGCGTCGTGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCAAGGAACCgg                                                                                                                                               <  1:584088/150‑1 (MQ=255)
          cgcTCTGCGTCGTGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCggg                                                                                                                                              <  2:8060923/150‑1 (MQ=255)
          cgcTCTGCGTCGTGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCggg                                                                                                                                              <  2:6210312/150‑1 (MQ=255)
          cgcTCTGCGTCGTGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCggg                                                                                                                                              <  1:3798904/150‑1 (MQ=255)
          cgcTCTGCGTCGTGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCAAGGAACCggg                                                                                                                                              <  1:1440138/150‑1 (MQ=255)
              ctGCGTCGTGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCttt                                                                                                                                         <  2:3190548/150‑1 (MQ=255)
                gCGTCGTGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCttttt                                                                                                                                       <  2:649780/150‑1 (MQ=255)
                gCGTCGTGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCttttt                                                                                                                                       <  2:4007461/150‑1 (MQ=255)
                gCGTCGTGCTGGCTACGTCACACGCGACGCTCGTGAAGTTGAGCGGAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCttttt                                                                                                                                       >  1:8055697/1‑150 (MQ=255)
                  gtcgtGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAACGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTGAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCttttttg                                                                                                                                     >  2:7596309/1‑149 (MQ=255)
                  gtcgtGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCttttttc                                                                                                                                     <  1:4926577/150‑2 (MQ=255)
                     gtGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTGGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTa                                                                                                                                  <  1:3406775/150‑1 (MQ=255)
                      tGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAg                                                                                                                                 >  2:3564447/1‑150 (MQ=255)
                        cTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGcc                                                                                                                               >  2:1681922/1‑150 (MQ=255)
                         tGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAATGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCt                                                                                                                              >  2:1759635/1‑150 (MQ=255)
                         tGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCt                                                                                                                              >  2:8889962/1‑150 (MQ=255)
                         tGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCt                                                                                                                              >  2:1111247/1‑150 (MQ=255)
                          ggCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCtt                                                                                                                             <  1:8171131/150‑1 (MQ=255)
                          ggCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCtt                                                                                                                             <  2:8620852/150‑1 (MQ=255)
                            cTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGa                                                                                                                           >  2:4659988/1‑150 (MQ=255)
                            cTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGa                                                                                                                           <  1:1813431/150‑1 (MQ=255)
                             tACGTCACCCGCGACGCTCGTGAAGTTGAGCGTGAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGaa                                                                                                                          >  1:2202325/1‑150 (MQ=255)
                             tACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGaa                                                                                                                          >  1:9090931/1‑150 (MQ=255)
                             tACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGaa                                                                                                                          >  1:2743862/1‑150 (MQ=255)
                                  cACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATCCCGCTTCGGCAGTCGAAAAAAGCCCGGGTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAAc                                                                                                                     <  1:4082351/150‑1 (MQ=255)
                                  cACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAAc                                                                                                                     <  1:986143/150‑1 (MQ=255)
                                    cccGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCt                                                                                                                   <  2:506152/150‑1 (MQ=255)
                                     ccGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTg                                                                                                                  >  1:6730442/1‑150 (MQ=255)
                                     ccGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTg                                                                                                                  >  1:7185896/1‑150 (MQ=255)
                                      cgcgACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGt                                                                                                                 <  2:4217728/150‑1 (MQ=255)
                                        cgACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTca                                                                                                               <  2:7415749/150‑1 (MQ=255)
                                         gACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTcag                                                                                                              <  1:888160/150‑1 (MQ=255)
                                              tCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGcagtc                                                                                                         >  2:2693008/1‑150 (MQ=255)
                                               cGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGcagtca                                                                                                        <  1:5162681/150‑1 (MQ=255)
                                                  gAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTg                                                                                                     <  1:3645031/150‑1 (MQ=255)
                                                   aaGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTgg                                                                                                    >  2:6088516/1‑150 (MQ=255)
                                                      ttGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCg                                                                                                 >  1:8400675/1‑150 (MQ=255)
                                                          gCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAATGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTTGCCTTGAAGTATCCTGTCAGCAGTCAGTGGCCGCTTg                                                                                             >  1:7648214/1‑150 (MQ=255)
                                                          gCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTg                                                                                             <  2:7070366/150‑1 (MQ=255)
                                                          gCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTg                                                                                             >  1:8164790/1‑150 (MQ=255)
                                                            gTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTc                                                                                           <  1:1142182/150‑1 (MQ=255)
                                                               aGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCa                                                                                        <  2:6647831/150‑1 (MQ=255)
                                                               aGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAAACTGTCAGCAGTCAGTGGCCGCGTGTCGCa                                                                                        >  1:6249584/1‑150 (MQ=255)
                                                                 aaaGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATa                                                                                      >  2:5603830/1‑150 (MQ=255)
                                                                         tCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGa                                                                              >  1:5343584/1‑150 (MQ=255)
                                                                            gCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCa                                                                           >  2:8378811/1‑150 (MQ=255)
                                                                             cGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCaa                                                                          <  2:5843750/150‑1 (MQ=255)
                                                                              gTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAg                                                                         <  2:6469779/150‑1 (MQ=255)
                                                                              gTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAg                                                                         <  2:3239568/150‑1 (MQ=255)
                                                                              gTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAg                                                                         <  2:6725952/150‑1 (MQ=255)
                                                                                  aGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTATCCGTGAACTAACCGGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCaaa                                                                     >  1:8699063/1‑150 (MQ=255)
                                                                                     gcgTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCg                                                                  >  2:5457525/1‑150 (MQ=255)
                                                                                        tAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTgg                                                               <  1:6132927/150‑1 (MQ=255)
                                                                                         aaGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTggg                                                              <  2:558958/150‑1 (MQ=255)
                                                                                          aGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGt                                                             >  2:685028/1‑150 (MQ=255)
                                                                                            cGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTg                                                           <  1:4819385/150‑1 (MQ=255)
                                                                                             gTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAATTGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTTATGGTTGTACTAAAGTGTCAGGAGTGAGTGGTGGCTTGTCGGATAGAGGTAGATCAATAAAATCGTGGGttt                                                           >  2:4276694/1‑149 (MQ=255)
                                                                                                cTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGc                                                       >  2:8336787/1‑150 (MQ=255)
                                                                                                   aGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCt                                                    >  1:316195/1‑150 (MQ=255)
                                                                                                     ctgcccAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTg                                                  >  2:2313672/5‑150 (MQ=255)
                                                                                                        tCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAAc                                                                                                                     >  1:4769716/1‑80 (MQ=255)
                                                                                                        tCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAAc                                                                                                                     <  2:4769716/80‑1 (MQ=255)
                                                                                                        tCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCg                                               >  2:2893644/1‑150 (MQ=255)
                                                                                                        tCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCg                                               >  1:5089051/1‑150 (MQ=255)
                                                                                                        ccgaaGCGGTATTAACGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCg                                               >  2:9158315/4‑150 (MQ=255)
                                                                                                          cAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGTACTAAACTGTCAGCAGTCAGTGGCCGCGTTTCGCATAGAAGTAGATAAAGCAAAGCGTGGGTTGTCGCCCTTGACggg                                             >  1:7000899/1‑150 (MQ=255)
                                                                                                          cAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCggg                                             <  1:1575226/150‑1 (MQ=255)
                                                                                                          cAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCggg                                             <  2:7763222/150‑1 (MQ=255)
                                                                                                          cAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCggg                                             >  1:5749321/1‑150 (MQ=255)
                                                                                                          cAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCggg                                             <  1:325365/150‑1 (MQ=255)
                                                                                                          cAAGCGTTAAGACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTTCCggg                                             <  2:8168787/150‑1 (MQ=255)
                                                                                                           aaGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCgggg                                            >  1:6500369/1‑150 (MQ=255)
                                                                                                           aaGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCGTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCgggg                                            >  1:2300020/1‑150 (MQ=255)
                                                                                                            aGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGt                                           <  2:5199138/150‑1 (MQ=255)
                                                                                                               gTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAAt                                        <  1:59797/150‑1 (MQ=255)
                                                                                                               gTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCATTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAAt                                        <  2:9108420/150‑1 (MQ=255)
                                                                                                                 tAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATc                                       <  2:7237111/150‑1 (MQ=255)
                                                                                                                     aCCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCtt                                  >  1:7819110/1‑150 (MQ=255)
                                                                                                                         cTTCGGCAGTCGAAAAAAGCCCGGGTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCa                              <  1:1059385/150‑1 (MQ=255)
                                                                                                                               gctaaaaaaaaaGCCCGGTTCCAAGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGGCGTAGATCAAGCAAAGCGGGGGGTGTCGCCCTTGCCGGGGTAATCACCATGTCAgtgtgt                        >  2:1935244/7‑150 (MQ=38)
                                                                                                                               gctaaaaaaaaaGCCCGGTTCCAAGGAACCGGGCTTATTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAgtgtgt                        <  1:5678113/144‑1 (MQ=39)
                                                                                                                                aGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAGtgtgtg                       <  2:1462453/150‑1 (MQ=255)
                                                                                                                                  tCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAGTGTGTGGa                     >  2:7982110/1‑150 (MQ=255)
                                                                                                                                   cGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAGTGTGTGGAt                    <  1:4177684/150‑1 (MQ=255)
                                                                                                                                       aaatGCCCGGTTCCTAGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGGATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCGCTTGCGGGGGTAATCAGCTTGTGAGTGTGTGGATTTGt                >  2:8245717/1‑150 (MQ=255)
                                                                                                                                          aGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAGTGTGTGGATTTGTTTc             >  1:4389190/1‑150 (MQ=255)
                                                                                                                                            cccGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAGTGTGTGGATTTGTTTCtt           >  2:5765027/1‑150 (MQ=255)
                                                                                                                                            cccGGTTCCAAGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAGTGTGTGGATTTGTTTCtt           >  2:7960235/1‑150 (MQ=255)
                                                                                                                                             ccGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAGTGTGTGGATTTGTTTCTTa          >  2:261134/1‑150 (MQ=255)
                                                                                                                                             ccGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAGTGTGTGGATTTGTTTCTTa          <  1:7441016/150‑1 (MQ=255)
                                                                                                                                                 ttCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAGTGTGTGGATTTGTTTCTTACCAt      <  2:6567914/150‑1 (MQ=255)
                                                                                                                                                   ccTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAGTGTGTGGATTTGTTTCTTACCATGa    <  1:5702218/150‑1 (MQ=255)
                                                                                                                                                     ttGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAGTGTGTGGATTTGTTTCTTACCATGACg  >  1:500522/1‑150 (MQ=255)
                                                                                                                                                     ttGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAGTGTGTGGATTTGTTTCTTACCATGACg  >  1:2874879/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
CCCTGCGTGGCGCTCTGCGTCGTGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAGTGTGTGGATTTGTTTCTTACCATGACG  >  NC_002947/1503927‑1504225

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