Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,955,264 G→C 5.7% intergenic (+51/‑142) PP_1752 → / → PP_1753 hypothetical protein/DMT superfamily permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,955,2640GC 5.7% 670.2 / 19.4 202intergenic (+51/‑142)PP_1752/PP_1753hypothetical protein/DMT superfamily permease
Reads supporting (aligned to +/- strand):  ref base G (83/107);  new base C (5/7);  total (88/114)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 212 total aligned reads displayed.

CGCTGGTGGAGGCCAATGACAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTCAACAACCCCCACCACGGGCTTGATCGAGGATTTTATGCACAC  >  NC_002947/1955116‑1955413
                                                                                                                                                    |                                                                                                                                                     
cGCTGGTGGAGGCCAATGACAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAgg                                                                                                                                                      <  1:4392000/150‑1 (MQ=255)
 gCTGGTGGAGGCCAATGACAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAggg                                                                                                                                                     >  2:4230188/1‑150 (MQ=255)
 gCTGGTGGAGGCCAATGACAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAggg                                                                                                                                                     >  1:212121/1‑150 (MQ=255)
  cTGGTGGAGGCCAATGACAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGt                                                                                                                                                    <  2:2711999/150‑1 (MQ=255)
    ggtggAGGCCAATGACAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCa                                                                                                                                                  <  2:3727836/150‑1 (MQ=255)
         tggCCAATGACAAGCGCGGGGCGGGGGTTATGGTGGTGCGCACAATCAACCAGGCGCTTGTGGCGGGGTTTGCCAAGCAAAATAAATAGCCGAAGGTGTGGTTGCCTGTGCTGACCCTTTCGCGGTTAAACCCGCGAAGCGGCCAGaaca                                                                                                                                             <  1:3202135/149‑1 (MQ=16)
         cggCCAATGACACGCGGTGGGCGGGGGTGATGGTGGTGCGCACAATCAACCAGGCGGTTGTGGGGGGGTTTGCCAAGCAAAACAAATAGCCGATGGTGTGGTTGCCTGTGCTGACCCTTCCGCGGGTAAACCCGCGAAGCGGCCAGaaca                                                                                                                                             <  1:2408827/149‑1 (MQ=12)
         aGGCCAATGACAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGcaca                                                                                                                                             >  2:1936309/1‑150 (MQ=255)
         aGGCCAATGACAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGcaca                                                                                                                                             >  1:3094491/1‑150 (MQ=255)
          ggCCATTGACAAGCGGTGGGCGGGGGTGATGGTGGTGCGCACAATCAACCAGGCGCTTGAGGGGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTGCTGACCCTTTCGCGGTTAAACCCGCGAAGCGGCCAGAACAg                                                                                                                                            <  1:3969113/150‑1 (MQ=21)
          ggCCAATGACAAGCGCTGGGCGCGGGTGATGGTGGTGCGCAAAATCAACCTGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTGCTGACCCTTTCGCGGTTAAACCCGCGAAGCGGCCAGAACAg                                                                                                                                            <  1:96264/150‑1 (MQ=25)
          ggCCAATGACAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAg                                                                                                                                            >  1:2949659/1‑150 (MQ=255)
            ccAATGACAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGCGGCGGCGTTGGCCACGCACAAGAAAGAGCCGACGGTGAGGTGGCCTGTGCTGACCCTTTCGCGGGTAAACCCGCGAAGCGGCCAGAACAGGc                                                                                                                                          <  2:968598/150‑1 (MQ=21)
             cAATGACAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCa                                                                                                                                         >  1:2243833/1‑150 (MQ=255)
             cAATGACAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCa                                                                                                                                         <  2:823931/150‑1 (MQ=255)
             cAATGACAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCa                                                                                                                                         <  2:3092973/150‑1 (MQ=255)
             cAATGACAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCa                                                                                                                                         >  2:727262/1‑150 (MQ=255)
             cAATGACAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCa                                                                                                                                         >  2:420586/1‑150 (MQ=255)
              aaTGACAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCaa                                                                                                                                        <  1:1500953/150‑1 (MQ=255)
              aaTGACAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCaa                                                                                                                                        >  2:2309010/1‑150 (MQ=255)
              aaTGACAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCaa                                                                                                                                        <  1:3022397/150‑1 (MQ=255)
                tGACAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAaca                                                                                                                                      <  1:2994877/150‑1 (MQ=255)
                   cAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCt                                                                                                                                   <  2:3741561/150‑1 (MQ=255)
                      gcgcTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCa                                                                                                                                <  2:29194/150‑1 (MQ=255)
                       cgcTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCAc                                                                                                                               >  2:2180427/1‑150 (MQ=255)
                          tGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCa                                                                                                                            >  1:1237734/1‑150 (MQ=255)
                           gggCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAg                                                                                                                           <  1:2340214/150‑1 (MQ=255)
                              cgcgGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCa                                                                                                                        >  1:2285889/1‑150 (MQ=255)
                                cgGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATg                                                                                                                      <  1:2465371/150‑1 (MQ=255)
                                        ggtggtGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCaaaa                                                                                                              <  2:4371086/150‑1 (MQ=255)
                                        ggtggtGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCaaaa                                                                                                              <  2:2491196/150‑1 (MQ=255)
                                         gtggtgCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAg                                                                                                             <  2:3109528/150‑1 (MQ=255)
                                           ggtgCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGaa                                                                                                           <  2:2157407/150‑1 (MQ=255)
                                           ggtgCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGaa                                                                                                           >  1:4030719/1‑150 (MQ=255)
                                            gtgCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGTACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAAt                                                                                                          <  1:373442/150‑1 (MQ=255)
                                            gtgCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAAt                                                                                                          <  2:2570055/150‑1 (MQ=255)
                                            gtgCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAAt                                                                                                          <  1:126432/150‑1 (MQ=255)
                                            gtgCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAAt                                                                                                          <  2:536971/150‑1 (MQ=255)
                                            gtgCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAAt                                                                                                          >  2:1623603/1‑150 (MQ=255)
                                                 cacaATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCt                                                                                                     >  2:1371454/1‑150 (MQ=255)
                                                    aaTCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCtgat                                                                                                  <  2:2906365/150‑1 (MQ=255)
                                                    aaTCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTAACCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCtgat                                                                                                  >  1:611872/1‑150 (MQ=255)
                                                    aaTCAACCAGGCGCTTGAGGCGGCGTTTGCAAAGCACAATAAATAGCCTAAGTTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCtgat                                                                                                  <  1:2014824/150‑1 (MQ=255)
                                                      tCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCtgatga                                                                                                >  1:328630/1‑150 (MQ=255)
                                                        aaCCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGaa                                                                                                           <  2:638635/137‑1 (MQ=255)
                                                        aaCCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGaa                                                                                                           >  1:638635/1‑137 (MQ=255)
                                                        aaCCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAAt                                                                                              <  2:1348470/150‑1 (MQ=255)
                                                        aaCCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAAt                                                                                              >  1:737359/1‑150 (MQ=255)
                                                          ccAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTc                                                                                            >  2:3419809/1‑150 (MQ=255)
                                                           cAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCt                                                                                           >  1:2558219/1‑150 (MQ=255)
                                                            aGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCtt                                                                                          <  2:1841403/150‑1 (MQ=255)
                                                             ggCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCttt                                                                                         <  2:1966076/150‑1 (MQ=255)
                                                             ggCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCttt                                                                                         <  2:1966150/150‑1 (MQ=255)
                                                             ggCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTAGCGGGTTTAACCCCGAAACGGTCAGCACAGGCAACACCACCACCCAACCAAGCCCCAAAAGAAATACCAGATGaaatatat                                                                                         >  1:1285033/1‑144 (MQ=255)
                                                              gcgcTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCtttt                                                                                        <  2:4136562/150‑1 (MQ=255)
                                                                  ttGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAAAGACCTGATGAATTCTTTTCTCg                                                                                    >  2:4316001/1‑150 (MQ=255)
                                                                  ttGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTAACCGCGAAAGGGTCAGCACAGGCAACACCACCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTATTTTCTCg                                                                                    >  1:493356/1‑150 (MQ=255)
                                                                   tGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCgg                                                                                   <  2:3453924/150‑1 (MQ=255)
                                                                   tGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCgg                                                                                   <  2:4404950/150‑1 (MQ=255)
                                                                   tGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCgg                                                                                   <  2:2252802/150‑1 (MQ=255)
                                                                       gcggcgTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGAACTGATGAATTATTTTCTCGGCAct                                                                               >  1:831676/1‑150 (MQ=255)
                                                                       gcggcgTTAGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTGTACCCCCGAAACGGTCCGCACCGGCAACACCCCCACCCAGCCAAGCCCCAAAAGAATGACCAGAAGAAGTCTTTTATCGGCAct                                                                               >  1:1361136/1‑150 (MQ=255)
                                                                             tttGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCAc                                                                         <  1:1609503/150‑1 (MQ=255)
                                                                             tttGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCAc                                                                         <  1:1024005/150‑1 (MQ=255)
                                                                                  cAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGtct                                                                    <  1:3651567/150‑1 (MQ=255)
                                                                                  cAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGtct                                                                    >  2:1632449/1‑150 (MQ=255)
                                                                                   aaGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGtctc                                                                   >  2:3269449/1‑150 (MQ=255)
                                                                                     gcccaAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCCCAGGCAACACCTCCACCCAGCCATGCCCCAAAAGGATGGGCTGATGAATTCGTTTGGCGGCACTCGCCCCCGTCTCTg                                                                 >  1:4191439/4‑150 (MQ=255)
                                                                                      cacaAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGt                                                                >  2:3171341/1‑150 (MQ=255)
                                                                                         aaGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGtcct                                                             >  2:4453016/1‑150 (MQ=255)
                                                                                          aGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGtcctc                                                            <  1:3048077/150‑1 (MQ=255)
                                                                                                 gCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAAcccc                                                     <  2:2539611/150‑1 (MQ=255)
                                                                                                 gCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAAcccc                                                     <  1:278402/150‑1 (MQ=255)
                                                                                                  ccGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGACTGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCg                                                    <  1:4458915/150‑1 (MQ=255)
                                                                                                         tGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGc                                                                                                                                          >  1:4396131/1‑57 (MQ=255)
                                                                                                          gAGGTTGGCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTc                                            <  2:2913707/150‑1 (MQ=255)
                                                                                                          gAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTc                                            >  2:1590327/1‑150 (MQ=255)
                                                                                                            ggTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGTTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTcaa                                          <  1:1204105/150‑1 (MQ=255)
                                                                                                             gTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTcaac                                         >  2:2227963/1‑150 (MQ=255)
                                                                                                                             gCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTCAACAACCCCCACCACGGGc                         <  1:3845518/150‑1 (MQ=255)
                                                                                                                             gCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTCAACAACCCCCACCACGGGc                         <  1:3457112/150‑1 (MQ=255)
                                                                                                                              tttCGCGGGTAAACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTCAACAACCCCCACCACGGGCt                        >  2:4246582/2‑150 (MQ=255)
                                                                                                                               ttCGCGGGTAAACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTCAACAACCCCCACCACGGGCtt                       >  1:2825288/1‑150 (MQ=255)
                                                                                                                                 cgcgGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTCAACAACCCCCACCACGGGCTTGa                     <  2:3901380/150‑1 (MQ=255)
                                                                                                                                  gcgGGTTTACCCGCGAAAGGGTCAGCACAGGCAGCACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTCAACAACCCCCACCACGGGCTTGAt                    >  1:3804297/1‑150 (MQ=255)
                                                                                                                                  gcgGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTCAACAACCCCCACCACGGGCTTGAt                    >  2:2330606/1‑150 (MQ=255)
                                                                                                                                  gcgGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCAAGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGACCTCCAACCCCCTAGAATTCAACAACCCCCACCACGGGCTTGAt                    >  2:466137/1‑150 (MQ=255)
                                                                                                                                  gcgGGTGTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTCAACAACCCCCACCACGGGCTTGAt                    <  2:1119567/150‑1 (MQ=255)
                                                                                                                                  gcgGGTAAACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTCAACAACCCCCACCACGGGCTTGAt                    <  1:4246582/150‑1 (MQ=255)
                                                                                                                                    gggTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTCAACAACCCCCACCACGGGCTTGATCg                  >  1:3914657/1‑150 (MQ=255)
                                                                                                                                     ggTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTCAACAACCCCCACCACGGGCTTGATCGa                 <  2:1003220/150‑1 (MQ=255)
                                                                                                                                     ggTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTCAACAACCCCCACCACGGGCTTGATCGa                 <  1:4429551/150‑1 (MQ=255)
                                                                                                                                     ggTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTCAACAACCCCCACCACGGGCTTGATCGa                 <  1:4230188/150‑1 (MQ=255)
                                                                                                                                      gTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTCAACAACCCCCACCACGGGCTTGATCGAg                <  2:527407/150‑1 (MQ=255)
                                                                                                                                         tACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTCAACAACCCCCACCACGGGCTTGATCGAGGAt             >  2:575148/1‑150 (MQ=255)
                                                                                                                                            ccGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTCAACAACCCCCACCACGGGCTTGATCGAGGAtttt          >  2:2364065/1‑150 (MQ=255)
                                                                                                                                              gcgAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTCAACAACCCCCACCACGGGCTTGATCGAGGATTTTAt        <  1:2180427/150‑1 (MQ=255)
                                                                                                                                              gcgAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTCAACAACCCCCACCACGGGCTTGATCGAGGATTTTAt        <  1:124319/150‑1 (MQ=255)
                                                                                                                                                 aaaGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTCAACAACCCCCACCACGGGCTTGATCGAGGATTTTATGca     >  2:2276836/1‑150 (MQ=255)
                                                                                                                                                    gggTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGGAGAATTCAACAACCCCCACCACGGGCTTGATCGAGGATTTTATGCacac  <  2:2768043/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CGCTGGTGGAGGCCAATGACAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTCAACAACCCCCACCACGGGCTTGATCGAGGATTTTATGCACAC  >  NC_002947/1955116‑1955413

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: