Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,087,557 C→T 6.4% intergenic (+301/‑294) 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,5570CT 6.4% 314.6 / 8.7 111intergenic (+301/‑294)PP_1864/PP_1865hypothetical protein/transposase
Reads supporting (aligned to +/- strand):  ref base C (79/24);  new base T (4/3);  total (83/28)
Fisher's exact test for biased strand distribution p-value = 3.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 112 total aligned reads displayed.

GTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTACAGTCCGTTTGCG  >  NC_002947/2087409‑2087706
                                                                                                                                                     |                                                                                                                                                       
gTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACg                                                                                                                                                        >  1:275987/1‑150 (MQ=34)
gTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACg                                                                                                                                                        >  1:275899/1‑150 (MQ=34)
  gACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGAc                                                                                                                                                      >  2:1442226/1‑150 (MQ=34)
   aCGGGCATCCATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGTCGGCTTTGGTTACTTTGGCCACGAcc                                                                                                                                                     >  1:807429/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:646511/1‑150 (MQ=34)
                  cGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCgg                                                                                                                                    >  1:2909552/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                                                                                                                    >  2:784860/1‑150 (MQ=34)
                                  ccAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGa                                                                                                                    >  1:1735518/1‑150 (MQ=34)
                                  ccAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGa                                                                                                                    >  2:2590571/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‑‑CGCCGGAAGGGCGGACAGGTGAATAGGCGcc                                                                                                           >  2:3811246/1‑150 (MQ=21)
                                           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)
                                                  aGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCg                                                                                                    >  1:2525230/1‑150 (MQ=37)
                                                   gggcgGTCCCCCGGGGAGCACGCGGGGTGTGGGGGTCCCCGGTGGGGTCGCTTGGGGGCGGTATGGGGGGAGCGGCGGGCTTTGGTTACTTTGCCAACGACAAATTTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGt                                                                                                   <  2:3531626/148‑1 (MQ=2)
                                                    cgcgGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGCCAGGTGAATAGGCGCCTGCATCGtg                                                                                                  >  1:4005610/1‑150 (MQ=255)
                                                     gcgGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGtga                                                                                                 >  2:3986829/1‑150 (MQ=37)
                                                      gggTCCGCCGGGGTGCACCCCGGGTGTGGGGGGACCCGGGTGGGTTGGTTCGGGGCGGGTTGTGGGGTGCCGCCGGGTTTGTTTTGTTGGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGtgaa                                                                                                <  2:3538484/149‑1 (MQ=2)
                                                         tcctccGGGGAGCAAGCCGGGAGTGAGGGGACCCCGGAGCGTAGCGTTGGGGCCGGATGAGGGGAGCCGACGGCTTTGTTTACTTTGGCCACGACAAAAGTACCC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGtgaatga                                                                                             <  2:2239208/150‑1 (MQ=11)
                                                            tccCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGaatgaatg                                                                                          >  2:886031/1‑150 (MQ=37)
                                                            tccCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGaatgaatg                                                                                          >  1:530895/1‑150 (MQ=37)
                                                            tccCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGaatgaatg                                                                                          >  2:4126470/1‑150 (MQ=37)
                                                             ccggggAGCACGCCGGTCGTGTGGGGACCCCGGAGCTTCGCTTTGGGGCGGTATGCGGGGAGCCGCGGGCTTGGGTTAGTTTGGCCACGACAAATTTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGt                                                                                         <  2:1692777/147‑1 (MQ=2)
                                                              ccGGGAGCAAGTCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGACAGGTGAATAGGCGCCTGCATCGTGAATGAATGTc                                                                                        >  1:3955630/1‑150 (MQ=255)
                                                              ccGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTc                                                                                        >  1:4516723/1‑150 (MQ=37)
                                                              ccGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTc                                                                                        >  2:255014/1‑150 (MQ=37)
                                                              ccGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTc                                                                                        >  2:1514378/1‑150 (MQ=37)
                                                              ccGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTc                                                                                        >  2:1514361/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)
                                                                          ccGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGAAACC‑‑CGCCGGAAGGGCGGAAACGTGAATAGTCGCCAGCATCGTCAATGAAAGACAACATGATGTGa                                                                            >  1:3531044/1‑150 (MQ=14)
                                                                                   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)
                                                                                                                        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:875716/1‑150 (MQ=255)
                                                                                                                                       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:3273068/1‑150 (MQ=255)
                                                                                                                                                    aCGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTACAGTCCGTTTGc   >  2:1655127/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                       
GTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTACAGTCCGTTTGCG  >  NC_002947/2087409‑2087706

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: