Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 868,055 A→C 8.8% intergenic (+150/+27) PP_0750 → / ← mqo‑I hypothetical protein/malate:quinone oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947868,0550AC 8.8% 482.4 / 24.6 164intergenic (+150/+27)PP_0750/mqo‑Ihypothetical protein/malate:quinone oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (105/42);  new base C (12/3);  total (119/45)
Fisher's exact test for biased strand distribution p-value = 5.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 165 total aligned reads displayed.

TCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACAGCGCAGCATCGGTTGCCAGGGTCTTCTCGCGAGCCGGGAAG  >  NC_002947/867906‑868202
                                                                                                                                                     |                                                                                                                                                   
     aGCCTGTGCGGGCCTCTTGGGGGGGGCGCCGGCTCCCACTGGTGTCGCGCCGCACGTGCGGGAATTACCGAACGTGTGGGAGGGGGCATGCCCGCGATGAGGGCGGTACAGGCTGATCAAGCGCCAGGAACAAAAAACGCCCCGAGGGGc                                                                                                                                                <  1:4013961/150‑1 (MQ=255)
       ccAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCAGCGGGCATTACCGAACCAGTGGGAGCCGGCATGCCCCCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAAAGCCCCGGGGGGCGt                                                                                                                                              >  2:3331453/1‑150 (MQ=39)
       ccAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCCCACCTGCGGGCCTTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGCGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCgg                                                                                                                                              >  2:3331116/1‑149 (MQ=255)
        cAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACCGGAGTCGCGCCGCCCCAGCGGGCATTACCGAACCGGTGGGAGCGGGCCTGCCCGCGAAGCGGGCCGTACACGCGGATCAAGCGCCCGGCACAAAAAACGCCCCGGGGGGCggt                                                                                                                                             >  2:2020219/1‑148 (MQ=255)
          gggcgggCCTCGGCGCGGGCGCGCCCGCGCCCCCCGGTGTCGCGCCGCACCGGGGGGCATTCCCGAACCGGTGGGCGCGGGCAGGCCGGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGtttt                                                                                                                                           <  1:2715300/145‑1 (MQ=255)
           taccgccctcttcgcgggcatgcccgctcccacaggttcggtaatgcccgcaggtgcggcgcgactcctgtgggagcgggcgcgcccgcgaagaggGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACCCCCCGAGGGGCGttttt                                                                                                                                          >  1:3778858‑M2/97‑150 (MQ=255)
              cgccctcttcgcgggcatgcccgctcccacaggttcggtaatgcccgcaggtgcggcgcgactccggtgggagcgggcgcgcccgcgaagaggCCGGCACTGGCTGATCAAGCGCCAGGCACAAAAAAAGCCCCGAGGGGCGTTTTTGCg                                                                                                                                       >  1:721980‑M2/94‑150 (MQ=255)
               ggCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCGGTGGGCGCGGGCCTGCCCGCGAAGAGGGCGGTCCAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGt                                                                                                                                      >  2:2798565/1‑150 (MQ=255)
               ggCCTCTGCGCGGGCGCGCCCGCCCCCACCGGAGCCGCCCCGCACCTGCGGGCATCACCGACCCGGTGGGCGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCTCGGGGCGttttttgtg                                                                                                                                      <  1:1409795/150‑5 (MQ=255)
                gCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGtt                                                                                                                                     <  2:4407604/150‑1 (MQ=255)
                   TCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAAAGCCCCGAGGGGCGTTTTTGCGTTTCa                                                                                                                                  <  2:2642579‑M1/150‑1 (MQ=255)
                   TCTTCGCGGGCGCGCCCGCGCCCACAGGAGTCGCGCCGCCCCTGCGGGCATTAACGAACCTGGGGGCGCGGGGATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCACGCCCCAAAAAACCCCCGCGGGGCGTTTTTTCGGTTCa                                                                                                                                  <  1:354896‑M1/150‑1 (MQ=255)
                     ggaGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCAAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGa                                                                                                                                <  2:5085719‑M1/147‑1 (MQ=255)
                         cgGGCGCGCCCCCTCCCCCAGGAGGAGCGCCCCACCAGCGGGCATTAACGAACCGGTGGGGGCGGGCCTTACCCCGAACAGGGGGGGACACGCTGGTCAAGCGCCAGGCACAAAAAACCCCCCGCGGGGCGTTTTTTCGTTTGACAcccg                                                                                                                            >  2:1302795/1‑148 (MQ=255)
                              gcgcCCGCTCCCACCGGAGGCGCGCCGCACCTGCGGGCATTACCGACCCTGGGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTaa                                                                                                                       <  1:698368/150‑1 (MQ=255)
                               cgcCCGCTCCCCCCGGCGTCGCGCCGCCCCTGCGGGCATTACCGACCCTGGGGGCGCGGGCATGCCCGCGAAGAGGGCTGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAg                                                                                                                      <  1:3154150/150‑1 (MQ=255)
                                gcCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGc                                                                                                                     <  1:3604483/150‑1 (MQ=255)
                                  ccGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCCGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGt                                                                                                                   >  1:2029360/1‑150 (MQ=255)
                                   cgcgCCCCCAGGTGTCGCGCCGCACCTGCGGGCATTACCGAACCTGGGGGAGGGGCCAGGCCCGGGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTa                                                                                                                  <  1:1852506/150‑1 (MQ=255)
                                   cGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGGGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTa                                                                                                                  <  2:4311177/150‑1 (MQ=255)
                                    gctcccacaggttcggtaatgcccgcaggtgcggcgcgactcctgtgggagcgggcgcgcccgcgaagaggGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAg                                                                                                                 >  2:71662‑M2/72‑150 (MQ=255)
                                    gctcccacaggttcggtaatgcccgcaggtgcggcgcgactcctgtgggagcgggcgcgcccgcgaagaggGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAg                                                                                                                 >  2:3399450‑M2/72‑150 (MQ=255)
                                      tcccacaggttcggtaatgcccgcaggtgcggcgcgactcctgtgggagcgggcatgcccgcgaagaggGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGc                                                                                                                           <  2:40432‑M2/69‑1 (MQ=255)
                                      tcccacaggttcggtaatgcccgcaggtgcggcgcgactcctgtgggagcgggcatgcccgcgaagaggGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGc                                                                                                                           >  1:2594120‑M2/70‑138 (MQ=255)
                                      tcccacaggttcggtaatgcccgcaggtgcggcgcgactcctgtgggagcgggcatgcccgcgaagaggGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGc                                                                                                                           >  1:40432‑M2/70‑138 (MQ=255)
                                      tcccacaggttcggtaatgcccgcaggtgcggcgcgactcctgtgggagcgggcatgcccgcgaagaggGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTg                                                                                                               >  1:4772324‑M2/70‑150 (MQ=255)
                                      tCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCCCAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTg                                                                                                               >  1:213846/1‑150 (MQ=255)
                                        ccacaggttcggtaatgcccgcaggtgcggcgcgactcctgtgggagcgggcgcgcccgcgaagaggGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCt                                                                                                             >  1:4895612‑M2/68‑150 (MQ=255)
                                        ccACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCt                                                                                                             >  2:1315427/1‑150 (MQ=255)
                                         cacaggttcggtaatgcccgcaggtgcggcgcgactcctgtgggagcgggcgtgcccgcgaagaggGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTg                                                                                                           >  2:1707775‑M2/67‑150 (MQ=255)
                                         cacaggttcggtaatgcccgcaggtgcggcgcgactcctgtgggagcgggcgcgcccgcgaagaggCCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCtt                                                                                                            >  2:2086866‑M2/67‑150 (MQ=255)
                                         cacaggttcggtaatgcccgcaggtgcggcgcgactcctgtgggagcgggcgcgcccgcgaagaggCCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCtt                                                                                                            >  2:4111003‑M2/67‑150 (MQ=255)
                                         cacaggttcggtaatgcccgcaggtgcggcgcgactcctgtgggagcgggcatgcccgcgaagaggGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCtt                                                                                                            >  1:4133085‑M2/67‑150 (MQ=255)
                                         cacaGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCg                                                                                                                    <  2:668170/142‑1 (MQ=255)
                                         cacaGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCtt                                                                                                            >  1:3748367/1‑150 (MQ=255)
                                         cacaGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCtt                                                                                                            >  2:1113972/1‑150 (MQ=255)
                                          accggtgTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTg                                                                                                           <  1:1093887/144‑1 (MQ=255)
                                          acaGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTg                                                                                                           >  1:574792/1‑150 (MQ=255)
                                          acaGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTg                                                                                                           >  1:2409697/1‑150 (MQ=255)
                                          acaGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTg                                                                                                           >  2:4610667/1‑150 (MQ=255)
                                          acaGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTg                                                                                                           >  1:1632372/1‑150 (MQ=255)
                                          acaGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTg                                                                                                           >  1:387605/1‑150 (MQ=255)
                                          acaGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTg                                                                                                           >  1:2263393/1‑150 (MQ=255)
                                            aGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGgc                                                                                                         >  1:1244856/1‑150 (MQ=255)
                                              gttcggtaatgcccgcaggtgcggcgcgactcctgtgggagcgggcatgcccgcgaagaggGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGgccg                                                                                                       >  2:2348315‑M2/62‑150 (MQ=255)
                                              gAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGgccg                                                                                                       >  2:2292213/1‑150 (MQ=255)
                                              gAGTCGCGCCGCACCTGCGGGCAGTACCGAACCGGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAACAACGCCCCGCGGGGCGTATTTGCGGTTCCGACCGGCTTAAGAGTAGTGCTTGgccg                                                                                                       >  2:4863254/1‑150 (MQ=255)
                                                   gcgcCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTGCCGACCGGCTTAAGCGTAGTGCTTGGCCGCCggg                                                                                                  >  2:2698781/1‑150 (MQ=255)
                                                    cgccgcACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGgc                                                                                                 <  2:5059205/150‑1 (MQ=255)
                                                               cGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGcc                                                                                      >  2:3777166/1‑150 (MQ=255)
                                                               cGGCCATCACCGACCCTGGGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGcc                                                                                      <  1:2369282/150‑1 (MQ=255)
                                                                    aTTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTc                                                                                 >  2:4317559/1‑150 (MQ=255)
                                                                       aCCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTTCGTTTCAGACCGGCATAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAg                                                                              >  1:3664881/1‑150 (MQ=255)
                                                                        ccGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAgg                                                                             <  2:2905656/150‑1 (MQ=255)
                                                                        ccGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAgg                                                                             <  1:2457428/150‑1 (MQ=255)
                                                                              cTGTGGGCGCGGCCAGGCCGGCGAAGCGGGCGGTATAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCcgcg                                                                       <  1:493990/150‑1 (MQ=255)
                                                                               tgtgGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAATCGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGcgt                                                                      >  2:2882373/1‑150 (MQ=255)
                                                                               tgtgGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGcgt                                                                      >  2:2563152/1‑150 (MQ=255)
                                                                               tgtgGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGcgt                                                                      >  2:1858479/1‑150 (MQ=255)
                                                                                gggggCGCGGGCAGGCCCGCGAAGTGGGCGGTACAGGCTGAACACGCGCCAGGCACAAAACACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGcgtc                                                                     <  1:57467/148‑1 (MQ=255)
                                                                                  gggAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGcgtcgt                                                                   >  1:1650466/1‑150 (MQ=255)
                                                                                   ttcgCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGtt                                                                  >  2:2358911/4‑150 (MQ=255)
                                                                                   ggAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGtt                                                                  >  2:4142800/1‑150 (MQ=255)
                                                                                    gAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTa                                                                 >  2:3616288/1‑150 (MQ=255)
                                                                                     agcgggcgcgcccgcgaagaggGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTAt                                                                >  1:3189079‑M2/23‑150 (MQ=255)
                                                                                     aGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTAt                                                                >  1:3077567/1‑150 (MQ=255)
                                                                                          gCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCa                                                           <  1:834219/150‑1 (MQ=255)
                                                                                          gCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCa                                                           >  1:3396084/1‑150 (MQ=255)
                                                                                                ccgcgaagaggGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATc                                                     >  1:1250229‑M2/12‑150 (MQ=255)
                                                                                                ccgcgaagaggGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGACAGGTCCAGGGCCGCGTCGTTATCGGCACTGATc                                                     >  1:3699469‑M2/12‑150 (MQ=255)
                                                                                                  gcgaagaggGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCtt                                                   >  1:896267‑M2/10‑150 (MQ=255)
                                                                                                     aagaggGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTtgtg                                                <  2:1632372‑M2/144‑1 (MQ=255)
                                                                                                     aagaggGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTtgtg                                                >  1:3341166‑M2/7‑150 (MQ=255)
                                                                                                          gGCGGTACAGGCTGAACAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCATGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACa                                           >  1:4147732‑M2/2‑150 (MQ=255)
                                                                                                               tACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACAGCgca                                      >  2:4579859/1‑150 (MQ=255)
                                                                                                               tACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACAGCgca                                      >  2:4579910/1‑150 (MQ=255)
                                                                                                                  aGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACAGCgcagca                                   >  2:4945117/1‑150 (MQ=255)
                                                                                                                  aGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACAGCgcagca                                   >  1:1013047/1‑150 (MQ=255)
                                                                                                                    gCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACAGCGCAGCATc                                 >  2:3607685/1‑150 (MQ=255)
                                                                                                                    gCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACAGCGCAGCATc                                 <  2:4068778/150‑1 (MQ=255)
                                                                                                                     cTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACAGCgca                                      >  1:4226955/1‑144 (MQ=255)
                                                                                                                     cTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACAGCgca                                      <  2:4226955/144‑1 (MQ=255)
                                                                                                                        aTCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACAGCGCAGCATCGGtt                             >  2:2019051/1‑150 (MQ=255)
                                                                                                                         tCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACAGCGCAGCATCGGTTg                            >  2:4741524/1‑150 (MQ=255)
                                                                                                                           aaGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACAGCGCAGCATCGGTTGcc                          >  1:3397674/1‑150 (MQ=255)
                                                                                                                            aGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACAGCGCAGCATCGGTTGCCa                         >  1:74970/1‑150 (MQ=255)
                                                                                                                            aGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGGGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACAGCGCAGCATCGGTTGCCa                         >  2:314956/1‑150 (MQ=255)
                                                                                                                             gcgcCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACAGCGCAGCATCGGTTGCCAg                        >  2:453009/1‑150 (MQ=255)
                                                                                                                              cgcCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACAGCGCAGCATCGGTTGCCAgg                       <  1:2798579/150‑1 (MQ=255)
                                                                                                                               gcCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACAGCGCAGCATCGGTTGCCAggg                      >  1:615331/1‑150 (MQ=255)
                                                                                                                                      acaAAAAACGCCCCGAGGGGCGTGTTTGCGTGTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACAGCGCAGCAGCGGTTGCCAGGGTCTtctc               >  2:4726925/1‑150 (MQ=255)
                                                                                                                                       caAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACAGCGCAGCATCGGTTGCCAGGGTCTTCTcg              >  2:2164927/1‑150 (MQ=255)
                                                                                                                                        caaaaaCGCCCCTCGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACAGCGCAGCATCGGTTGCCCGGGTCTTCTcgc             >  1:4546138/2‑150 (MQ=255)
                                                                                                                                         aaaaaCGCCCCTCGGGGCGTTTTTGCGTTTCAGACCGGCATAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCCTCGTTATCGGCACAGATCTTGTGGTACAGCGCAGCAACGGTTTCCAGGGTCTTCTcgcg            >  1:451894/1‑150 (MQ=255)
                                                                                                                                          acaaCGCCCCGCGGGGCGATTTTGCGTTTCAGCCCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACAGCGCCGCTTCGGTTGCCAGGGTCTTCTCGCGa           <  1:2549297/148‑1 (MQ=255)
                                                                                                                                          aaaaCGCCCCTCGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGGAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACAGCGCAGCATCGGTTGCCAGGGTCTTCTCGCGa           >  1:1357129/1‑150 (MQ=255)
                                                                                                                                             aCGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACAGCGCAGCATCGGTTGCCAGGGTCTTCTCGCGAGcc        <  1:1014247/150‑1 (MQ=255)
                                                                                                                                                 cccGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACAGCGCAGCATCGGTTGCCAGGGTCTTCTCGCGAGCCGGGa    >  2:1943164/1‑150 (MQ=255)
                                                                                                                                                 cccGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACAGCGCAGCATCGGTTGCCAGGGTCTTCTCGCGAGCCGGGa    >  1:381210/1‑150 (MQ=255)
                                                                                                                                                 cccGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACAGCGCAGCATCGGTTGCCAGGGTCTTCTCGCGAGCCGGGa    >  1:1079672/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACAGCGCAGCATCGGTTGCCAGGGTCTTCTCGCGAGCCGGGAAG  >  NC_002947/867906‑868202

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: