Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,087,559 A→G 6.3% intergenic (+303/‑292) PP_1864 → / → PP_1865 hypothetical protein/transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,087,5590AG 6.3% 316.9 / 8.5 111intergenic (+303/‑292)PP_1864/PP_1865hypothetical protein/transposase
Reads supporting (aligned to +/- strand):  ref base A (78/24);  new base G (4/3);  total (82/29)
Fisher's exact test for biased strand distribution p-value = 3.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

Only 100 of 112 total aligned reads displayed.

GACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTACAGTCCGTTTGCGAT  >  NC_002947/2087411‑2087708
                                                                                                                                                     |                                                                                                                                                       
gACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGAc                                                                                                                                                        >  2:1442226/1‑150 (MQ=34)
 aCGGGCATCCATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGTCGGCTTTGGTTACTTTGGCCACGAcc                                                                                                                                                       >  1:807429/1‑150 (MQ=255)
 aCGGGCATCCATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGGCGGCTTTGGTTACTTTGGCCACGAcc                                                                                                                                                       >  1:807464/1‑150 (MQ=255)
      cttGGAGGCCCGACAGGCGGGGGGCCCCTGGTGGGGGCCGGCTGGGGGGCCCCCCCGGGACCACGCCCGGTGTGTGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGTCGTGGCCAAAGt                                                                                                                                                  <  1:3783927/148‑1 (MQ=2)
        tGGATGCCCGTCAAGCGCTGGGCCCCAGGTTGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGTCGTGGCCAAAGTaa                                                                                                                                                <  1:4026118/150‑1 (MQ=11)
         ggATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAAc                                                                                                                                               >  2:3582448/1‑150 (MQ=34)
         ggATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAAc                                                                                                                                               >  1:262297/1‑150 (MQ=34)
          gATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGTAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAAcc                                                                                                                                              <  2:1558593/150‑1 (MQ=21)
             gCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACc‑‑cgc                                                                                                                                         >  1:1921401/1‑150 (MQ=34)
                cGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCgg                                                                                                                                      >  1:646511/1‑150 (MQ=34)
                cGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGCAGCAACCCCGGAGTGAGGGGACCCCGCAGCGAAGCGTAGGGGCCGGATCAGGGGAGCCGACGGCTTGGGTTACTTTGGTCGTGGCCAAAGTAAccaaagcc                                                                                                                                        <  2:1099166/150‑7 (MQ=11)
                    aaGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAgg                                                                                                                                  >  2:2276441/1‑150 (MQ=34)
                          tGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTCGCGAAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAACGGCGGaa                                                                                                                            >  1:2110978/1‑150 (MQ=255)
                            ggCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGAGGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTACCC‑‑CGCCGGAAGGGCGGAAAg                                                                                                                          >  2:1238858/1‑150 (MQ=21)
                                ccAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGa                                                                                                                      >  1:1735518/1‑150 (MQ=34)
                                ccAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGa                                                                                                                      >  2:784860/1‑150 (MQ=34)
                                 cAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGaa                                                                                                                     >  1:1610864/1‑150 (MQ=34)
                                  atcctGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCCAAGTAACC‑‑CGCCGGAAGGGCGGACAGGTGAAt                                                                                                                    >  2:3822133/5‑150 (MQ=255)
                                    gggggggCCTGCCGCGCGTCCCCCCCTTGCGCTCGCCGGGGGTGGGGGGCCCGCGGCGGGTGGCGTTGGGGGGGGTTGGGGGGCGCCGAGGTCTTTGGTTTTTTGGGCCACGACAAAATTACCC‑‑CGCCGGAAGGGCGGAAAGGTGAATAg                                                                                                                  <  2:3162251/147‑1 (MQ=2)
                                     ttGGGGCCGGCCGCGCGGTCCTCCCGGGAGCACGCCGGGTGTGTGGGGCCCCCGGTGCGTTGCGTTGGGGCCGGTTGCGGGGCGCCGAGGGCTTTGGTTACTTTGGCCACGACAAATGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAgg                                                                                                                 <  2:356171/149‑1 (MQ=11)
                                        gggCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGgcgc                                                                                                              >  2:3384998/1‑150 (MQ=34)
                                         ggCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGcc                                                                                                             >  1:1426821/1‑150 (MQ=34)
                                           ccTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTg                                                                                                           >  1:3593890/1‑150 (MQ=35)
                                            cTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGc                                                                                                          >  2:130343/1‑150 (MQ=35)
                                               ctgcggGGTCCCCCGGGGAGCAAGCCCGGTGTGTGGGGCCCCCGGCGCGTCGGGTTGGGGCCGGATGCGGGGAGCCGAGGGCTTTGGTTCCTTGGGCCACGCCCAATGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATc                                                                                                       <  2:4009047/148‑1 (MQ=11)
                                                aGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCg                                                                                                      >  1:2525230/1‑150 (MQ=37)
                                                 gggcgGTCCCCCGGGGAGCACGCGGGGTGTGGGGGTCCCCGGTGGGGTCGCTTGGGGGCGGTATGGGGGGAGCGGCGGGCTTTGGTTACTTTGCCAACGACAAATTTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGt                                                                                                     <  2:3531626/148‑1 (MQ=2)
                                                  cgcgGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGCCAGGTGAATAGGCGCCTGCATCGtg                                                                                                    >  1:4005610/1‑150 (MQ=255)
                                                  cgcgGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGtg                                                                                                    >  2:437336/1‑150 (MQ=37)
                                                   gcgGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGtga                                                                                                   >  2:3986829/1‑150 (MQ=37)
                                                    gggTCCGCCGGGGTGCACCCCGGGTGTGGGGGGACCCGGGTGGGTTGGTTCGGGGCGGGTTGTGGGGTGCCGCCGGGTTTGTTTTGTTGGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGtgaa                                                                                                  <  2:3538484/149‑1 (MQ=2)
                                                     ggTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGtgaat                                                                                                 >  2:2557710/1‑150 (MQ=37)
                                                       tcctccGGGGAGCAAGCCGGGAGTGAGGGGACCCCGGAGCGTAGCGTTGGGGCCGGATGAGGGGAGCCGACGGCTTTGTTTACTTTGGCCACGACAAAAGTACCC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGtgaatga                                                                                               <  2:2239208/150‑1 (MQ=11)
                                                       tcctccCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGtgaatga                                                                                               >  2:3662782/1‑150 (MQ=37)
                                                          tccCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGaatgaatg                                                                                            >  1:530895/1‑150 (MQ=37)
                                                          tccCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGaatgaatg                                                                                            >  2:4126470/1‑150 (MQ=37)
                                                          tccCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGaatgaatg                                                                                            >  2:886031/1‑150 (MQ=37)
                                                           ccggggAGCACGCCGGTCGTGTGGGGACCCCGGAGCTTCGCTTTGGGGCGGTATGCGGGGAGCCGCGGGCTTGGGTTAGTTTGGCCACGACAAATTTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGt                                                                                           <  2:1692777/147‑1 (MQ=2)
                                                            ccGGGAGCAAGTCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGACAGGTGAATAGGCGCCTGCATCGTGAATGAATGTc                                                                                          >  1:3955630/1‑150 (MQ=255)
                                                            ccGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTc                                                                                          >  2:255014/1‑150 (MQ=37)
                                                            ccGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTc                                                                                          >  2:1514361/1‑150 (MQ=37)
                                                                  gCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTg                                                                                    >  2:2118819/1‑150 (MQ=37)
                                                                    aaGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTAATTGAt                                                                                  >  1:1523556/1‑150 (MQ=37)
                                                                     aGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGAtg                                                                                 >  2:2844959/1‑150 (MQ=37)
                                                                       cccGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGAtgtg                                                                               >  2:3603607/1‑150 (MQ=38)
                                                                        ccGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGa                                                                              >  2:4132192/1‑150 (MQ=39)
                                                                        ccGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGa                                                                              >  2:326432/1‑150 (MQ=39)
                                                                               gtggggACCCCGGAGCGTAGCTTTGGGGCCGGTTGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGca                                                                       <  2:1464478/148‑1 (MQ=255)
                                                                                 ggggACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATATGCGCCTGGATCGTGAATGAATGTCTACTTGATTTGAAATCGcaca                                                                     >  1:1706359/1‑150 (MQ=255)
                                                                                 ggggACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGcaca                                                                     >  2:897730/1‑150 (MQ=39)
                                                                                 ggggACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGcaca                                                                     >  2:897721/1‑150 (MQ=39)
                                                                                 ggggACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGcaca                                                                     >  2:1269932/1‑150 (MQ=39)
                                                                                  gggACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACaa                                                                    >  1:1461345/1‑150 (MQ=39)
                                                                                  gggACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACGTGATGTGAAAGCGCACaa                                                                    >  1:1559801/1‑150 (MQ=255)
                                                                                   ggACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACTGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACAAGATGTGAAAGCTCACAAc                                                                   >  2:298759/1‑150 (MQ=255)
                                                                                       cccGGTGCGTCGCTTTGGGGCCGGATGTGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACAAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCaa                                                               <  1:3704704/150‑1 (MQ=14)
                                                                                       cccGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCaa                                                               >  2:2653020/1‑150 (MQ=39)
                                                                                       cccGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTAATTTATGTTAAAGCGCACAACCCaa                                                               >  1:4157692/1‑150 (MQ=255)
                                                                                        ccGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGCAAGGTGAATAGGCGCCGGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCaaa                                                              >  1:2996811/1‑150 (MQ=255)
                                                                                        ccGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCaaa                                                              >  1:463591/1‑150 (MQ=39)
                                                                                        ccGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCaaa                                                              >  2:2555391/1‑150 (MQ=39)
                                                                                         cGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGACTGTCTACTTGATGTGAAAGCGCACAACCCAAAg                                                             >  2:3820656/1‑150 (MQ=255)
                                                                                                 tagcgtagGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCcaaa                                                      >  2:3824923/1‑150 (MQ=255)
                                                                                                  tggtttggggCGGGTTGTGGGGAGCGGAGGGATTTGTTTATTTTGGCCAGTACAAAATTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCcaaaa                                                     <  1:3302572/144‑1 (MQ=14)
                                                                                                  agcgtagGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCcaaaa                                                     >  2:3557288/1‑150 (MQ=255)
                                                                                                   gcgtcgGGGCCGGATGCGGGGAGCCCCCGGCTTTGGTTCCTGTGGCCACGCCCAACGTACCC‑‑CGCCGGAAGGGCGGAAAGGGGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCcaaaag                                                    <  2:2023025/150‑1 (MQ=14)
                                                                                                    cgtagGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGAGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCcaaaagc                                                   >  1:3874102/1‑150 (MQ=255)
                                                                                                         gggCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGACAGGTGAATAGGCGCCTGCATCGTGTATGAATTTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCcaaaagctcaaa                                              >  1:2839369/1‑150 (MQ=255)
                                                                                                          ggCCGGTTGTGGGGAGCCGACGGCTTTGTTTATTTTGGCCACTACAAATGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCcaaaagctcaaaa                                             <  2:1738995/150‑1 (MQ=14)
                                                                                                             cGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCcaaaagctcaaaagct                                          >  2:2504249/1‑150 (MQ=255)
                                                                                                              ggATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCtt                                         >  2:1363014/1‑150 (MQ=255)
                                                                                                               gttgtggggAGCGGCGGGTTTTGTTTACTTGGGCCACGACCAATTTACCC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCttt                                        <  1:1189493/145‑1 (MQ=14)
                                                                                                                 tggggggCGCGGCCGGCTTGGTTTACTTGGGCCACGACAAAATAACCC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAg                                      <  2:232989/147‑1 (MQ=14)
                                                                                                                 tgcggggCGCCGACGGCTTTGTTTAATTTGCCCCCCACCAATGAAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAg                                      <  2:232996/147‑1 (MQ=14)
                                                                                                                    ggggAGCCGACCGCTTTGGTAACTATGCCCACTACAACAGTCACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAACGGCCCAAAAGCTCAAAAGCTTTAgcgc                                   <  1:3365877/150‑1 (MQ=17)
                                                                                                                      ggAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGGAAGGTGAATAGGCGCCTGCATCGGGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAACCTCAAAAGCTTTAGCGCCg                                 >  1:496372/1‑150 (MQ=255)
                                                                                                                      ggAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCg                                 >  1:3638293/1‑150 (MQ=255)
                                                                                                                          ccGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTGTAGCGCCGGTTc                             >  1:2979502/1‑150 (MQ=255)
                                                                                                                             ccGGTTTTGTTTCCTTTGGCCACTAACAATGTAACC‑‑CGCCGTAAGTGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATTTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGa                          <  2:1577145/149‑1 (MQ=17)
                                                                                                                             aCGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGTATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGa                          >  1:1858797/1‑150 (MQ=255)
                                                                                                                                 ggcgggTTACTTTGGTCGTGGCCAAAGTAACCAAAGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATc                        >  2:3315810/5‑150 (MQ=17)
                                                                                                                                     ggTTACTTTGGTCGTGGCCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGt                  >  1:875711/1‑150 (MQ=255)
                                                                                                                                      gtgccTTGGGCCACGACCAAAGTACCC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTa                 <  2:3233570/146‑1 (MQ=255)
                                                                                                                                      gTTACTTTGGCCACGACCAACGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTa                 <  2:807429/150‑1 (MQ=255)
                                                                                                                                      gTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGTGCCGGTTCTGATCGGGCGTa                 <  2:807464/150‑1 (MQ=255)
                                                                                                                                       ttACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGAAGTGAAAGCGCACAACCCAAAGGCCCAAAAGCACAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTAc                >  1:1953419/1‑150 (MQ=255)
                                                                                                                                        tACTTTGGCCCCGACAAATTTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTACa               <  1:2390132/150‑1 (MQ=255)
                                                                                                                                        tACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCgc                                                                                                     <  2:631992/64‑2 (MQ=25)
                                                                                                                                        tACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCgc                                                                                                     >  1:631992/1‑63 (MQ=25)
                                                                                                                                         aCTTTGGTCGTGGCCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTACAg              >  1:1168547/1‑150 (MQ=255)
                                                                                                                                           tttGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGGAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGGTCTGATCGGGCGGACAGTc            >  1:3442724/1‑150 (MQ=255)
                                                                                                                                                 cACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGCAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTACAGTCCGTTTg      >  1:311531/1‑150 (MQ=255)
                                                                                                                                                 cACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTACAGTCCGTTTg      >  2:2820932/1‑150 (MQ=255)
                                                                                                                                                 cACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTACAGTCCGTTTg      >  2:2016790/1‑150 (MQ=255)
                                                                                                                                                  aCGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTACAGTCCGTTTGc     <  1:2962218/150‑1 (MQ=255)
                                                                                                                                                  aCGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTACAGTCCGTTTGc     >  2:3119397/1‑150 (MQ=255)
                                                                                                                                                  aCGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTACAGTCCGTTTGc     >  2:1655127/1‑150 (MQ=255)
                                                                                                                                                    gACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTACAGTCCGTTTGCGa   <  1:3384998/150‑1 (MQ=255)
                                                                                                                                                    gACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTACAGTCCGTTTGCGa   >  2:633759/1‑150 (MQ=255)
                                                                                                                                                     aCCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCGGCCTCGGGAATGAATGTCTACTTGCTGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCGTTAGCGCCGGTTCTGAGCGGGCGGACAGTCCGGGTGCGag  >  2:178992/1‑149 (MQ=18)
                                                                                                                                                     |                                                                                                                                                       
GACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTACAGTCCGTTTGCGAT  >  NC_002947/2087411‑2087708

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: