Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 170,841 A→C 5.5% intergenic (+22/‑548) PP_0159 → / → PP_16SA family III CoA‑transferase/16S ribosomal RNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947170,8410AC 5.5% 506.7 / 14.4 187intergenic (+22/‑548)PP_0159/PP_16SAfamily III CoA‑transferase/16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (91/80);  new base C (5/7);  total (96/91)
Fisher's exact test for biased strand distribution p-value = 5.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

Only 100 of 189 total aligned reads displayed.

CGGAGACGCCAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAAGTAGTTAAGGC  >  NC_002947/170692‑170987
                                                                                                                                                     |                                                                                                                                                  
cGGAGACGCCAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCa                                                                                                                                                    >  1:3443574/1‑150 (MQ=255)
 ggAGACGCCAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAg                                                                                                                                                   <  1:1772079/150‑1 (MQ=255)
  gagaCGCCAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAgg                                                                                                                                                  >  1:5202446/1‑150 (MQ=255)
  gagaCGCCAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGAGGAAAAGCGGGGTGGCCGTCAgg                                                                                                                                                  >  1:4837253/1‑150 (MQ=255)
   agaCGCCAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGc                                                                                                                                                 <  1:4512861/150‑1 (MQ=255)
       gCCAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTcc                                                                                                                                             >  2:1290513/1‑150 (MQ=255)
        ccAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCTGACGGCCTccc                                                                                                                                            <  2:1054572/150‑1 (MQ=255)
        ccAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTccc                                                                                                                                            <  2:2597356/150‑1 (MQ=255)
        ccAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTccc                                                                                                                                            <  1:305728/150‑1 (MQ=255)
         cAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCTGACGGCCTcccc                                                                                                                                           <  2:6495228/150‑1 (MQ=255)
          aGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCg                                                                                                                                          >  1:5679091/1‑150 (MQ=255)
            tGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCTGACGGCCTCCCCGCt                                                                                                                                        <  1:4330494/150‑1 (MQ=255)
             ggtgTATCGCCGGGCGCCACCAGGGCGGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCGGGATGCTGGCGAGGTACAGCGTCTGCGTGATGGTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCCCCCCGCtt                                                                                                                                       <  2:4018953/147‑1 (MQ=255)
             ggAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCAGGATGCTGGCGAGGTACAGCTTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGGGGGGAGGCCGTCAGGCCTCCCCGCtt                                                                                                                                       <  1:3930861/150‑1 (MQ=255)
              gAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGATGTGTTGGGGCTGGATGCTGGTGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCTGACGGCCTCCCCGCttt                                                                                                                                      <  1:1009962/150‑1 (MQ=255)
                gTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGt                                                                                                                                    >  1:3366908/1‑150 (MQ=255)
                 tATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGtt                                                                                                                                   >  1:6643918/1‑150 (MQ=255)
                  aTCGCCGGGCGCCCCCAGTGCGGGGCGAGCATACCGACGTGGTGCTGGGGGAAGTGTTGGGGCGGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGCGGCCGTCAGGCCTCCCCGCTTTGttt                                                                                                                                  <  1:5663262/150‑1 (MQ=255)
                  aTCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGttt                                                                                                                                  <  2:4693494/150‑1 (MQ=255)
                   tCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGtttt                                                                                                                                 >  1:3810136/1‑150 (MQ=255)
                          gcgcCCCCACTGCTGGGCGAGCATACCGACGTGGTGCTGGGGGACGTGTTGGGGCTGGTTGCTGGCGAGGTACAGCGTCTGCGTGAGGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAgcgc                                                                                                                          <  1:5017810/150‑1 (MQ=39)
                          gcgcCACCACTGCTGGGCGAGCATACCGACGTGGTGGTGGGGGAAGTGTTGGGGCTGGATGCGGGCGAGGTACCGCGTCTGCGTGATGCGGGGGTGCTCTTGATCAAAAGCGGGGCGGCCGACAGGCCGCCCCGCTTTGTTTTTCAgcgc                                                                                                                          <  1:1538450/150‑1 (MQ=255)
                           cgcCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCt                                                                                                                         <  2:5029622/150‑1 (MQ=255)
                                 cacTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCaaa                                                                                                                   <  2:4556596/150‑1 (MQ=255)
                                 cacTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCACCCCTTTGTTTTTAAGCGCTATCaaa                                                                                                                   >  1:5585946/1‑150 (MQ=255)
                                 cacTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGAACAGCGTCTGCTTGATGCTGGGGTGCTCTAGATCAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCaaa                                                                                                                   <  2:1916178/150‑1 (MQ=255)
                                  acTGCTGGGCGAGCATACCGACGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAg                                                                                                                  <  2:6810504/150‑1 (MQ=255)
                                    tgctgGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAATGt                                                                                                                >  2:3453014/1‑150 (MQ=255)
                                    tgctgGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGt                                                                                                                >  1:4369304/1‑150 (MQ=255)
                                      ctgGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCCATCAAAGGTaa                                                                                                              >  2:4724287/1‑150 (MQ=255)
                                         ggCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTg                                                                                                           <  2:3639113/150‑1 (MQ=255)
                                         ggCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTg                                                                                                           >  1:3912251/1‑150 (MQ=255)
                                                aTACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTaa                                                                                                    >  2:6038736/1‑150 (MQ=255)
                                                aTACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCACCCCGCGTTGTTTTTAAGCGCAATAAAAGGTAATTGAAATTaa                                                                                                    >  1:3708099/1‑150 (MQ=255)
                                                 tccggatgtggtggtgGGGTATGTGTGGGGGCTGTTTGGTGGCGTGTTACAGGGTCTGCGTGTTGCGGGGTTGCTATATATTAAATGGGGGGTGGCGGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTaaa                                                                                                   <  1:5271985/143‑1 (MQ=255)
                                                 tACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTaaa                                                                                                   <  2:6404446/150‑1 (MQ=255)
                                                 tACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCGGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTaaa                                                                                                   <  2:1662804/150‑1 (MQ=255)
                                                    cGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTAAGCGCTATCAAAGGTAATTGAAATTAAAGGt                                                                                                >  1:6461981/1‑150 (MQ=255)
                                                    cGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCCCGTTGTTTTTCCGCGCTATCAAAGGTAATTTAAATTAAAGGt                                                                                                >  2:5820261/1‑150 (MQ=255)
                                                     gAAGTGGTGCTGGGGGACGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGtt                                                                                               <  2:5645442/150‑1 (MQ=255)
                                                     gAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGACGCCCTCATGCCACCCCCCATTGTTTTTAAGCCCAAAAAAAGGTAATataaaataaaattt                                                                                               >  1:2071858/1‑146 (MQ=255)
                                                     catgtggtgCTGGGGGATGTGTTGGGGCTGTATGCTGGCGTGTTACAGCGTCTGCTTGTTGCGGGGGTGCTATAGATTAAATGGGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGtt                                                                                               <  1:3122809/147‑1 (MQ=255)
                                                        gtggtgCTGGGGGATGTGTTGGGGCGGGATGCTGGCGAGGTACCGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGAc                                                                                            <  2:2071858/150‑1 (MQ=255)
                                                        gtggtgCTGGGGGATGGGTTGGGGCTGGATGCTGGCTAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCTGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGAc                                                                                            <  2:2869131/150‑1 (MQ=255)
                                                            tgCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGccc                                                                                        <  2:1567312/150‑1 (MQ=255)
                                                              cTGGGGGACGTGTTGGGGCTGGATGCTGGCGAGGTACCGGGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCCGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCctc                                                                                      <  2:1401088/150‑1 (MQ=255)
                                                                ggggtATGTGTTGGGGCTGGATGGTGGCGCGGTACTGCGTCGGGGTGTTGGTGGGGTGCTCTAGATGAAACGCGGGGCGGCCGCCCGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCtctc                                                                                    <  1:1375860/150‑1 (MQ=255)
                                                                gggggTAGTGTTGGGGCTGGATGCGGGCGTGGTACAGCGTCTGCGTGATGGTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACTCCCtctc                                                                                    <  2:5291556/150‑1 (MQ=39)
                                                                gggggATGTGTGGGGGGTGGATGCTGGCGAGGTACAGCGTCTGCGTGTTGCGGGGGTGCTCTAGATGAAAAGGGGGGCGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCtctc                                                                                    <  1:4797546/150‑1 (MQ=255)
                                                                gggggAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATTAAAAGCGGGGCGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCtctc                                                                                    <  1:6415654/150‑1 (MQ=255)
                                                                gggggAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCtctc                                                                                    <  2:5669417/150‑1 (MQ=255)
                                                                   ggAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCATTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGaa                                                                                 >  1:3288120/1‑150 (MQ=255)
                                                                   ggAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTACTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGaa                                                                                 >  2:374452/1‑150 (MQ=255)
                                                                    tttgtgtTGGGGGTGTATGGTGGCGTGTTACCGCTTGTGGTTTTTGCTGGGGTGATCTAGATTAAAAGGGGGGTGGCCGTCAGGCCCCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAAt                                                                                <  1:295239/147‑1 (MQ=255)
                                                                     aaGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATc                                                                               <  2:2249289/150‑1 (MQ=255)
                                                                     aaGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATc                                                                               <  1:3974031/150‑1 (MQ=255)
                                                                      cgtgtGGGGGCGGTCTGCTGGCGCGGCACCCCGGCGGCGTCAGGCGGGGGGCCCCTAGCGCAACCGCGGGGCCGCCGTCAGGCGGCCCCGCTTTGTTTGTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATcc                                                                              <  1:5798870/149‑1 (MQ=255)
                                                                      aGTGTCGGGGCTGCATGCTGGCGAGGGACAGCGTCGGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTACAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCACTCGAATcc                                                                              <  2:3256920/150‑1 (MQ=39)
                                                                       gtgtTGGGGCTGGTTGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATccc                                                                             <  2:41321/150‑1 (MQ=255)
                                                                        tgtTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTTTCGAATcccc                                                                            <  2:4902175/150‑1 (MQ=255)
                                                                         gtTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCt                                                                           >  2:5539195/1‑150 (MQ=255)
                                                                          ttGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCtt                                                                          >  2:3491348/1‑150 (MQ=255)
                                                                            ggggCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTtat                                                                        >  1:2796945/1‑150 (MQ=255)
                                                                                 tGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCACTCGAATCCCCTTATAATgc                                                                   >  1:1635254/1‑150 (MQ=255)
                                                                                   gATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGGGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATAAAAGGTAATTGTAATTAAAGGTTGTCGCCCTCTCGAATCCCCTTATAATgcgc                                                                 >  1:4153527/1‑150 (MQ=255)
                                                                                     tGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCGTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGccc                                                               >  2:3799648/1‑150 (MQ=255)
                                                                                      gCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGcccc                                                              <  2:3904918/150‑1 (MQ=255)
                                                                                        tGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTGAGGCCTCCCCGGTTTGTTTTTGAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCGTATAATGCGCCCCAc                                                            >  2:1359161/1‑150 (MQ=255)
                                                                                        tGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCAc                                                            >  2:5516461/1‑150 (MQ=255)
                                                                                             aGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCa                                                       <  2:778209/150‑1 (MQ=255)
                                                                                             aGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCGTTGTTTTTAAGCGCTATCAAAGGTAATTGAAATTAAAGGTTTAAGCCCTCTCGAATACCCGTATAATGCGCCCCACGTCCa                                                       >  1:3264335/1‑150 (MQ=255)
                                                                                               gTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGCGGCCGGCAGGCCGCCCCCCGTTGTTTTTCAGCGCTATAAAAGGTAATTGAAATTAAAGGTTGACGCCCACTCGCATCCCCATATAAATCGCCCCACATCCAGc                                                     >  1:2004296/1‑150 (MQ=37)
                                                                                                  cAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGGGGCCGGCAGGGCTCCCCCCGTTGTTTTTAAGCGCGATAAAAGGTAATTGAAATTAAAGGTTTAAGCCCACAAGAATACCCCTATAATTCGCCCCCCATACAGCGAc                                                  >  1:2451179/1‑150 (MQ=255)
                                                                                                     cGTCGGGGTGTTGCTGGGGTGCACTAGATCAAACGGGGGGAGGCCGTCTGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATc                                               <  1:211972/150‑1 (MQ=25)
                                                                                                        cTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGa                                            <  2:3366908/150‑1 (MQ=255)
                                                                                                         tGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGaa                                           <  2:1141137/150‑1 (MQ=255)
                                                                                                                 gCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAAc                                   <  1:4543484/150‑1 (MQ=255)
                                                                                                                 gCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCCCTTTGTTTTTAAGCGCTATAAAATGTAATTGAAATTAAATGTTGAAGCCCTCTCGAATACCCTTATAAATCGCCCCACATCCAGCCACAACCGAACCAAAAAc                                   >  1:6660227/1‑150 (MQ=255)
                                                                                                                   tGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTc                                 <  1:4911995/150‑1 (MQ=255)
                                                                                                                   tGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACATCCAGCGACATCGGAACGAAAAACTc                                 >  2:3231815/1‑150 (MQ=255)
                                                                                                                    ggggTGCTCTAGATGAAAAGCGGGGAGGCCGTCCGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTTAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTcc                                <  1:4619411/150‑1 (MQ=255)
                                                                                                                      ggTGCTCTAGATGAAAAGCGGGGCGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCtt                              <  1:531661/150‑1 (MQ=255)
                                                                                                                      ggTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCtt                              <  1:2183036/150‑1 (MQ=255)
                                                                                                                        tGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTga                            <  2:1250093/150‑1 (MQ=255)
                                                                                                                          ctctAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTgaga                          <  2:5110278/150‑1 (MQ=255)
                                                                                                                           tctAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTgagct                         >  2:4493377/1‑148 (MQ=255)
                                                                                                                           tctAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAAATAAAGGTTGACGCCCTCTCGAATCCCCTTATAATTCGCCCCACTTCCAGCGACATCGGAACGAAAAACACCTAGAGAt                         >  2:4944515/1‑150 (MQ=255)
                                                                                                                             tAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACCTCGGAACGAAAAACTCCTTGCGATTc                       >  2:3815769/1‑150 (MQ=255)
                                                                                                                               gATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCaa                     <  2:1735472/150‑1 (MQ=255)
                                                                                                                               gATGAAAAGCGGGGAGGCCGGCAGGCCTCCCCCCGTTGTTTTTCCGCGCGATCAAAGGTAATTGAAATTAAAGGTTGACGCCCGCTCGAAACCCCTTATAATGCGCCCCACATCCAGCGACCACGGACCGAAAAACACCGTGCGATTcca                     >  2:1487050/1‑148 (MQ=21)
                                                                                                                                     aaGCGGGGAGGCCTGACGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAAACCCCGTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAAATCCATGAGATTCAATGAGtt               >  1:1010973/1‑150 (MQ=18)
                                                                                                                                     aaGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGtt               >  2:2694587/1‑150 (MQ=255)
                                                                                                                                       gCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTaa             >  1:1695946/1‑150 (MQ=255)
                                                                                                                                       gCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTaa             <  2:3049938/150‑1 (MQ=255)
                                                                                                                                         ggggAGGCCTGACGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAagt           >  2:6319052/1‑150 (MQ=255)
                                                                                                                                               gCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAAGTAGTTaa     <  1:928980/150‑1 (MQ=255)
                                                                                                                                               gCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAAGTAGTTaa     >  1:676268/1‑150 (MQ=255)
                                                                                                                                                ccttCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAAGTAGTTAAg    <  1:1290513/147‑1 (MQ=255)
                                                                                                                                                 cGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAAGTATTTAAgg   >  2:4051346/1‑150 (MQ=255)
                                                                                                                                                 cGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAAGTAGTTAAgg   >  2:5224587/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
CGGAGACGCCAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAAGTAGTTAAGGC  >  NC_002947/170692‑170987

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