Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,085,635 C→T 5.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,085,6350CT 5.4% 266.0 / 3.9 93intergenic (+301/‑294)PP_1864/PP_1865hypothetical protein/transposase
Reads supporting (aligned to +/- strand):  ref base C (44/44);  new base T (3/2);  total (47/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.31e-01

GGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTACAGTCCGTTT  >  NC_002947/2085492‑2085781
                                                                                                                                               |                                                                                                                                                    
gggCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCaa                                                                                                                                                >  2:1449569/1‑150 (MQ=35)
 ggCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCaaa                                                                                                                                               <  1:1778255/150‑1 (MQ=35)
 ggCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCaaa                                                                                                                                               <  1:1777794/150‑1 (MQ=35)
  gCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAg                                                                                                                                              >  2:947183/1‑150 (MQ=35)
    aTGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTa                                                                                                                                            >  1:1846754/1‑150 (MQ=35)
          gCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTCGTTACTTTGGTCGTGGCCAAAGTAAccaaa                                                                                                                                      <  1:2135989/150‑4 (MQ=11)
            ccGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACc‑‑cgccg                                                                                                                                  >  1:2158383/1‑150 (MQ=35)
            ccGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACc‑‑cgccg                                                                                                                                  >  2:1998325/1‑150 (MQ=35)
             cGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCCCGACCAAAGTAACC‑‑CGCCgg                                                                                                                                 >  2:972454/1‑150 (MQ=25)
              gTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGa                                                                                                                                >  2:149411/1‑150 (MQ=35)
                 aaGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAgg                                                                                                                             <  2:46560/150‑1 (MQ=35)
                  aGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAggg                                                                                                                            >  1:1812912/1‑150 (MQ=35)
                         ggCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAg                                                                                                                     <  1:2341879/150‑1 (MQ=35)
                         ggCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGAAGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGaaa                                                                                                                      <  2:2361384/149‑1 (MQ=25)
                         ggCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGAAGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGaaa                                                                                                                      >  1:2361384/1‑149 (MQ=25)
                         ggCCCCAGGATGGGGACTGAAGGGAGGTCCAAAAGGAGTTAAGTCTGGAAGTAGCGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAg                                                                                                                     <  2:1841721/150‑1 (MQ=11)
                           ccccAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGt                                                                                                                   <  2:1201342/150‑1 (MQ=35)
                           ccccAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGt                                                                                                                   <  1:1233554/150‑1 (MQ=35)
                           ccccAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGt                                                                                                                   >  1:1418047/1‑150 (MQ=35)
                              cAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGTCCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGaa                                                                                                                >  2:789082/1‑150 (MQ=25)
                                ggATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCACCGAACAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATa                                                                                                              <  1:789082/150‑1 (MQ=255)
                                 gATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAg                                                                                                             <  2:108912/150‑1 (MQ=35)
                                 gATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAg                                                                                                             >  1:2267577/1‑150 (MQ=35)
                                  aTGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAgg                                                                                                            >  1:617117/1‑150 (MQ=35)
                                    ggggCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGgcg                                                                                                          >  1:1134974/1‑150 (MQ=35)
                                     gggCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGgcgc                                                                                                         <  1:1449569/150‑1 (MQ=35)
                                       gCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCt                                                                                                       <  2:661690/150‑1 (MQ=35)
                                         cTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGc                                                                                                     >  2:1136105/1‑150 (MQ=35)
                                          tGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGaa                                                                                                                >  2:292825/1‑138 (MQ=37)
                                          tGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCAGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGaa                                                                                                                <  1:292825/138‑1 (MQ=25)
                                              gcgcGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTACCC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGt                                                                                                <  2:1846754/150‑1 (MQ=255)
                                              gcgcGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCAGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACAAAAGTAAAC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGt                                                                                                <  2:1163885/150‑1 (MQ=255)
                                                 cgGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGtgaa                                                                                             <  1:1330237/150‑1 (MQ=39)
                                                    tcctccCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCATGACCAAAGTAACC‑‑CGCCGGAAGGGCGGaa                                                                                                                       <  1:511582/121‑1 (MQ=37)
                                                    tcctccCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCATGACCAAAGTAACC‑‑CGCCGGAAGGGCGGaa                                                                                                                       >  2:511582/1‑121 (MQ=37)
                                                    tcctccCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGtgaatga                                                                                          >  1:1605240/1‑150 (MQ=38)
                                                      ctccCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTgaatgaat                                                                                        >  2:457294/1‑150 (MQ=38)
                                                      ctccCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTgaatgaat                                                                                        <  2:2267577/150‑1 (MQ=38)
                                                        cccGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGt                                                                                      <  2:2084042/150‑1 (MQ=38)
                                                         ccGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTc                                                                                     >  1:17408/1‑150 (MQ=38)
                                                            ggAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTAc                                                                                  <  2:17408/150‑1 (MQ=39)
                                                                   gCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGc                                                                                                                           <  2:546467/102‑1 (MQ=38)
                                                                   gCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGc                                                                                                                           >  1:546467/1‑102 (MQ=38)
                                                                    cccGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGAtgtg                                                                          >  1:2345174/1‑150 (MQ=40)
                                                                     ccGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGa                                                                         >  2:558855/1‑150 (MQ=40)
                                                                         aGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAgc                                                                     <  2:1134974/150‑1 (MQ=255)
                                                                          gTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGTGGAAAGGTGAATAGGCGCCTGCATCGTgaatgaat                                                                                        <  2:2001623/130‑1 (MQ=255)
                                                                          gTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGTGGAAAGGTGAATAGGCGCCTGCATCGTgaatgaat                                                                                        >  1:2001623/1‑130 (MQ=255)
                                                                           tGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGTCGTGGCCAAAGTaa                                                                                                                                           >  1:346087/1‑80 (MQ=14)
                                                                           tGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAgcgc                                                                   <  2:1354230/150‑1 (MQ=255)
                                                                            gAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGTGGGAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGca                                                                  <  2:2158383/150‑1 (MQ=255)
                                                                            gAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGca                                                                  <  2:617117/150‑1 (MQ=255)
                                                                             aGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTATTC‑‑CTCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGAATCTTTAATGAATGTCTACTTGATGTGAAAACGcaa                                                                 >  1:531689/1‑149 (MQ=14)
                                                                               gggACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACaa                                                               >  2:2284162/1‑150 (MQ=40)
                                                                                    cccGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGa                                                                                                                 <  1:2149962/95‑1 (MQ=39)
                                                                                    cccGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGa                                                                                                                 >  2:2149962/1‑95 (MQ=39)
                                                                                    cccGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGtgaatga                                                                                          >  2:138902/1‑118 (MQ=255)
                                                                                    cccGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGtgaatga                                                                                          <  1:138902/118‑1 (MQ=255)
                                                                                    cccGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCaa                                                          >  1:252043/1‑150 (MQ=255)
                                                                                    cccGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCaa                                                          >  1:896269/1‑150 (MQ=255)
                                                                                     ccGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCaaa                                                         <  1:516640/150‑1 (MQ=255)
                                                                                      cGGCGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATCTCTACTTGATGTGAAAGCGCACAACCCAAAg                                                        <  2:772447/150‑1 (MQ=255)
                                                                                       ggAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAgg                                                       <  1:947183/150‑1 (MQ=255)
                                                                                       ggAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAgg                                                       >  2:1059720/1‑150 (MQ=255)
                                                                                             tagcgtagGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCcaaa                                                 >  1:1443962/1‑150 (MQ=255)
                                                                                              agcgtagGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTACCC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCcaaaa                                                <  2:438620/150‑1 (MQ=255)
                                                                                                      ggCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAGCC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCcaaaagctcaaaa                                        <  1:1136105/150‑1 (MQ=255)
                                                                                                      ggCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCcaaaagctcaaaa                                        >  2:1184715/1‑150 (MQ=255)
                                                                                                        ccGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCcaaaagctcaaaagc                                      >  1:1395109/1‑150 (MQ=255)
                                                                                                         cGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCcaaaagctcaaaagct                                     <  1:1184715/150‑1 (MQ=255)
                                                                                                           gATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCttt                                   <  1:637579/150‑1 (MQ=255)
                                                                                                                 gggAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGcc                             >  1:1252708/1‑150 (MQ=255)
                                                                                                                  ggAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCg                            >  2:2378144/1‑150 (MQ=255)
                                                                                                                   gAGCCGGCGGTTTTGGTTACTTGGCCCACGACCAGAGTCCCC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCgg                           <  2:1812912/150‑1 (MQ=17)
                                                                                                                    aGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGt                          >  1:1310989/1‑150 (MQ=255)
                                                                                                                       cGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGACTAGGAGCAATTATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCt                       <  1:147064/150‑1 (MQ=17)
                                                                                                                         aCGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCcaaaagctcaaaagct                                     >  1:1393120/1‑134 (MQ=255)
                                                                                                                         aCGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCcaaaagctcaaaagct                                     <  2:1393120/134‑1 (MQ=255)
                                                                                                                           ggCTTTGGTTACTTTGGCCACGACCAAAGTCCCC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATc                   <  1:1704563/150‑1 (MQ=255)
                                                                                                                           ggCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATc                   <  1:423563/150‑1 (MQ=255)
                                                                                                                             cTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCgg                 <  2:1605240/150‑1 (MQ=255)
                                                                                                                             cTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCgg                 <  1:219968/150‑1 (MQ=255)
                                                                                                                                  gTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTa            <  1:2284162/150‑1 (MQ=255)
                                                                                                                                   ttACTTTGGTCGTGGCCAAAGTAACCAAAGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGt             >  2:668423/1‑150 (MQ=18)
                                                                                                                                   ttACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTgaatgaat                                                                                        >  2:1214022/1‑73 (MQ=35)
                                                                                                                                   ttACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTgaatgaat                                                                                        <  1:1214022/73‑1 (MQ=35)
                                                                                                                                   ttACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAACGGCCCAAAAGCTCAAAAGCTTTATCGCCGCTTCTGATCGGGCGTAc           >  2:388757/1‑150 (MQ=255)
                                                                                                                                   ttACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTAc           >  2:453596/1‑150 (MQ=255)
                                                                                                                                   ttACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCGCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTAc           >  2:1309736/1‑150 (MQ=255)
                                                                                                                                    tACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTACa          <  1:457294/150‑1 (MQ=255)
                                                                                                                                       tttGGCCACGACCAAAGTAACC‑‑AGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTACAGTc       >  2:868280/1‑150 (MQ=255)
                                                                                                                                            ccACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTACAGTCCGttt  >  1:177467/1‑150 (MQ=255)
                                                                                                                                            ccACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTACAGTCCGttt  >  2:1195491/1‑150 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
GGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTACAGTCCGTTT  >  NC_002947/2085492‑2085781

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