Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,087,554 C→T 6.7% intergenic (+298/‑297) 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,5540CT 6.7% 286.4 / 8.8 104intergenic (+298/‑297)PP_1864/PP_1865hypothetical protein/transposase
Reads supporting (aligned to +/- strand):  ref base C (73/24);  new base T (4/3);  total (77/27)
Fisher's exact test for biased strand distribution p-value = 3.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 106 total aligned reads displayed.

GTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTACAGTC  >  NC_002947/2087409‑2087698
                                                                                                                                                  |                                                                                                                                                  
gTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACg                                                                                                                                                >  1:275987/1‑150 (MQ=34)
gTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACg                                                                                                                                                >  1:275904/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)
           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:2909552/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:2590571/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)
                                            gCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGGGACTAGGCGccg                                                                                                  >  2:502719/1‑149 (MQ=21)
                                             ccTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTg                                                                                                 >  1:3593890/1‑150 (MQ=35)
                                              cTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGc                                                                                                >  2:130343/1‑150 (MQ=35)
                                                 ctgcggGGTCCCCCGGGGAGCAAGCCCGGTGTGTGGGGCCCCCGGCGCGTCGGGTTGGGGCCGGATGCGGGGAGCCGAGGGCTTTGGTTCCTTGGGCCACGCCCAATGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATc                                                                                             <  2:4009047/148‑1 (MQ=11)
                                                 cAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTTGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGc                                                                                                                                                    <  2:1961141/97‑1 (MQ=37)
                                                 cAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGc                                                                                                                                                    >  1:1961141/1‑97 (MQ=37)
                                                  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                                                                                  >  2:4126470/1‑150 (MQ=37)
                                                            tccCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGaatgaatg                                                                                  >  2:886031/1‑150 (MQ=37)
                                                            tccCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGaatgaatg                                                                                  >  1:530895/1‑150 (MQ=37)
                                                             ccggggAGCACGCCGGTCGTGTGGGGACCCCGGAGCTTCGCTTTGGGGCGGTATGCGGGGAGCCGCGGGCTTGGGTTAGTTTGGCCACGACAAATTTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGt                                                                                 <  2:1692777/147‑1 (MQ=2)
                                                              ccGGGAGCAAGTCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGACAGGTGAATAGGCGCCTGCATCGTGAATGAATGTc                                                                                >  1:3955630/1‑150 (MQ=255)
                                                              ccGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTc                                                                                >  2:1514378/1‑150 (MQ=37)
                                                              ccGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTc                                                                                >  2:255014/1‑150 (MQ=37)
                                                              ccGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTc                                                                                >  2:1514361/1‑150 (MQ=37)
                                                              ccGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTc                                                                                >  1:4516723/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:326432/1‑150 (MQ=39)
                                                                          ccGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGa                                                                    >  2:4132192/1‑150 (MQ=39)
                                                                          ccGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGAAACC‑‑CGCCGGAAGGGCGGAAACGTGAATAGTCGCCAGCATCGTCAATGAAAGACAACATGATGTGa                                                                    >  1:3531044/1‑150 (MQ=14)
                                                                                 gtggggACCCCGGAGCGTAGCTTTGGGGCCGGTTGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGca                                                             <  2:1464478/148‑1 (MQ=255)
                                                                                   ggggACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATATGCGCCTGGATCGTGAATGAATGTCTACTTGATTTGAAATCGcaca                                                           >  1:1706359/1‑150 (MQ=255)
                                                                                   ggggACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGcaca                                                           >  2:897721/1‑150 (MQ=39)
                                                                                   ggggACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGcaca                                                           >  2:897730/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)
                                                                                           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)
                                                                                                                                     ttGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGc          <  2:2182649/150‑1 (MQ=255)
                                                                                                                                       ggTTACTTTGGTCGTGGCCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGt        >  1:875711/1‑150 (MQ=255)
                                                                                                                                       ggTTACTTTGGTCGTGGCCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGt        >  1:875716/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)
                                                                                                                                                  |                                                                                                                                                  
GTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTACAGTC  >  NC_002947/2087409‑2087698

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: