Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,409,410:1 +A 8.9% coding (136/732 nt) rlmB → 23S rRNA 2'‑O‑ribose G2251 methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,409,4101.A 8.9% 668.0 / 66.6 179coding (136/732 nt)rlmB23S rRNA 2'‑O‑ribose G2251 methyltransferase
Reads supporting (aligned to +/- strand):  ref base . (77/86);  new base A (11/5);  total (88/91)
Fisher's exact test for biased strand distribution p-value = 1.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.17e-01

Only 100 of 186 total aligned reads displayed.

CAACGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGATCGCTTCGCTCGATCAACC  >  NC_000913/4409261‑4409558
                                                                                                                                                      |                                                                                                                                                    
caacGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGtctc                                                                                                                                                       >  1:971576/1‑150 (MQ=255)
  acGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTGCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAg                                                                                                                                                    >  1:1606266/1‑150 (MQ=255)
   cGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAgg                                                                                                                                                   >  2:587998/1‑150 (MQ=255)
      ttACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGTCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCg                                                                                                                                                <  1:2207216/149‑1 (MQ=255)
      gTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCg                                                                                                                                                >  1:1702026/1‑150 (MQ=255)
       tACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTGAGACTCAAGGGCGt                                                                                                                                               <  1:1510647/150‑1 (MQ=255)
         cATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTgg                                                                                                                                             <  2:2335746/150‑1 (MQ=255)
         cATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTgg                                                                                                                                             <  1:697934/150‑1 (MQ=255)
         cATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTgg                                                                                                                                             >  2:1199813/1‑150 (MQ=255)
         cATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTgg                                                                                                                                             <  1:372576/150‑1 (MQ=255)
          aTTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGt                                                                                                                                            <  1:1949072/150‑1 (MQ=255)
           ttAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTGAGACTCAAGGGCGTgaat                                                                                                                                           <  2:118274/150‑4 (MQ=255)
            tAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTa                                                                                                                                          >  2:1052412/1‑150 (MQ=255)
            tAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTGAGACTCAAGGGCGTgaatc                                                                                                                                          <  2:1491501/150‑5 (MQ=255)
             aaTGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTAt                                                                                                                                         >  2:1731904/1‑150 (MQ=255)
             aaTGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATc                                                                                                                                        <  1:1874261/150‑1 (MQ=255)
              aTGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATc                                                                                                                                        <  2:615731/150‑1 (MQ=255)
                gAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCa                                                                                                                                      >  1:1102169/1‑150 (MQ=255)
                  gCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGt                                                                                                                                    >  2:1832001/1‑150 (MQ=255)
                    gAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTg                                                                                                                                  <  2:1666334/150‑1 (MQ=255)
                       aTGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCa                                                                                                                               <  1:2153654/150‑1 (MQ=255)
                       aTGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCa                                                                                                                               >  1:749524/1‑150 (MQ=255)
                          aTTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAAc                                                                                                                            >  2:308478/1‑150 (MQ=255)
                            ttACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCg                                                                                                                          <  1:1732175/150‑1 (MQ=255)
                             tACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGc                                                                                                                         >  2:1974752/1‑150 (MQ=255)
                              aCGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGcc                                                                                                                        <  2:2487711/150‑1 (MQ=255)
                               cGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGTAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCa                                                                                                                       >  1:258653/1‑150 (MQ=255)
                                ggCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCaa                                                                                                                      <  2:1903544/150‑1 (MQ=255)
                                 gCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAat                                                                                                                     <  1:1956382/150‑1 (MQ=255)
                                  cATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAata                                                                                                                    >  2:1698158/1‑150 (MQ=255)
                                  cATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAata                                                                                                                    >  1:2360404/1‑150 (MQ=255)
                                      cACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATAtctc                                                                                                                >  1:1297726/1‑150 (MQ=255)
                                         gCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTcgac                                                                                                             <  2:1343909/150‑1 (MQ=255)
                                              gCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGaa                                                                                                        >  2:1206909/1‑150 (MQ=255)
                                              gCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGaa                                                                                                        >  2:920933/1‑150 (MQ=255)
                                                aGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGaaaa                                                                                                      >  2:1502916/1‑150 (MQ=255)
                                                 ggCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAg                                                                                                     <  1:2147982/150‑1 (MQ=255)
                                                   cccTCCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCg                                                                                                   >  2:1537872/1‑150 (MQ=255)
                                                    ccTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGa                                                                                                  >  2:1887338/1‑150 (MQ=255)
                                                       gctgGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACgg                                                                                               >  1:1067680/1‑150 (MQ=255)
                                                       gctgGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACgg                                                                                               <  2:2137110/150‑1 (MQ=255)
                                                                cgcCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCa                                                                                      <  2:971576/150‑1 (MQ=255)
                                                                  ccccTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATc                                                                                    >  1:466601/1‑150 (MQ=255)
                                                                  ccccTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATc                                                                                    >  1:189708/1‑150 (MQ=255)
                                                                      tGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAggg                                                                                >  1:1480168/1‑150 (MQ=255)
                                                                              ttCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCg                                                                        <  1:164882/150‑1 (MQ=255)
                                                                                cAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGcc                                                                      >  1:1644252/1‑150 (MQ=255)
                                                                                 agcaaGTATTTATTTTAAAATGAAATGAAGATAAACTTATGTTACACCTGATTAAAGAAATTAAATCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGccc                                                                     >  2:1800564/4‑150 (MQ=255)
                                                                                        cTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGaa                                                              >  1:1615412/1‑150 (MQ=255)
                                                                                         tttATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAg                                                             <  2:1962859/150‑1 (MQ=255)
                                                                                          ttATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGc                                                            <  2:2436845/150‑1 (MQ=255)
                                                                                           tATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGcc                                                           >  2:396864/1‑150 (MQ=255)
                                                                                           tATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGcc                                                           >  2:787489/1‑150 (MQ=255)
                                                                                                 aaaaGGCAGGGAAGATAAACGTCTGTTACCGCTGATTCTCGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACg                                                     <  2:1730573/150‑1 (MQ=255)
                                                                                                   aaGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTc                                                   <  2:35584/150‑1 (MQ=255)
                                                                                                    aGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCa                                                  <  1:1042687/150‑1 (MQ=255)
                                                                                                     ggCCGTGAAGATTAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAg                                                 <  2:1723064/150‑1 (MQ=255)
                                                                                                     ggCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAg                                                 <  2:1846157/150‑1 (MQ=255)
                                                                                                       ccGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGACCGCGTGAAGCCAGGACGTCAGTa                                               >  1:733499/1‑150 (MQ=255)
                                                                                                          tGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCa                                            >  1:1826618/1‑150 (MQ=255)
                                                                                                           gAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAg                                           <  1:1800477/150‑1 (MQ=255)
                                                                                                           gAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAg                                           <  1:550488/150‑1 (MQ=255)
                                                                                                            aaGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAgg                                          <  1:547076/150‑1 (MQ=255)
                                                                                                                 aaaCGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAAc                                     <  1:280674/150‑1 (MQ=255)
                                                                                                                 aaaCGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAAc                                     >  1:755041/1‑150 (MQ=255)
                                                                                                                 aaaCGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGGCAGTATCAGGAAAAc                                     <  1:253408/150‑1 (MQ=255)
                                                                                                                       cTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTg                               <  2:206466/150‑1 (MQ=255)
                                                                                                                         gTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGcc                             <  2:1480168/150‑1 (MQ=255)
                                                                                                                           tACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGACAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCgg                           >  1:1391306/1‑150 (MQ=255)
                                                                                                                            aCCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGa                          <  1:397156/150‑1 (MQ=255)
                                                                                                                             ccGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGAt                         <  1:1199813/150‑1 (MQ=255)
                                                                                                                             ccGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGAt                         >  2:960317/1‑150 (MQ=255)
                                                                                                                              cGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCTTTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATc                        <  1:606945/150‑1 (MQ=255)
                                                                                                                              cGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATc                        <  2:1771633/150‑1 (MQ=255)
                                                                                                                               gCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCt                       <  1:2250289/150‑1 (MQ=255)
                                                                                                                               gCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCt                       <  2:162969/150‑1 (MQ=255)
                                                                                                                                gataaccACGCCCTGAGACTCAAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTg                      >  2:1378609/7‑150 (MQ=255)
                                                                                                                                  taaccACGCCCTGAGACTCAAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGAt                    >  2:1024507/5‑150 (MQ=255)
                                                                                                                                    ttCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGAtcg                  >  2:741598/1‑150 (MQ=255)
                                                                                                                                    ttCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGAtcg                  <  2:2390041/150‑1 (MQ=255)
                                                                                                                                    accACGCCCTGAGACTCAAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGAtcg                  >  2:488702/3‑150 (MQ=255)
                                                                                                                                    accACGCCCTGAGACTCAAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGAtcg                  >  2:825316/3‑150 (MQ=255)
                                                                                                                                    accACGCCCTGAGACTCAAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGAtcg                  >  2:2236209/3‑150 (MQ=255)
                                                                                                                                     tCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGAtcgc                 >  2:1448987/1‑150 (MQ=255)
                                                                                                                                     ccACGCCCTGAGACTCAAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGAtcgc                 >  1:889920/2‑150 (MQ=255)
                                                                                                                                      cACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGAtcgct                >  1:2458332/1‑150 (MQ=255)
                                                                                                                                      cACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGAtcgct                >  1:2005955/1‑150 (MQ=255)
                                                                                                                                        cGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGAtcgcttc              <  1:308478/150‑1 (MQ=255)
                                                                                                                                         gCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGAtcgcttcg             >  1:1749451/1‑150 (MQ=255)
                                                                                                                                           ccTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGAtcgcttcgct           <  1:1022222/150‑1 (MQ=255)
                                                                                                                                           ccTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCACCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGAtcgcttcgct           >  1:1349391/1‑150 (MQ=255)
                                                                                                                                           ccTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCACCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGAtcgcttcgct           >  1:1340518/1‑150 (MQ=255)
                                                                                                                                               gAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGATCGCTTCGCTCGAt       <  1:2374184/150‑1 (MQ=255)
                                                                                                                                               gAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGATCGCTTCGCTCGAt       <  1:1284149/150‑1 (MQ=255)
                                                                                                                                               gAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGATCGCTTCGCTCGAt       <  1:949772/150‑1 (MQ=255)
                                                                                                                                                aGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGATCGCTTCGCTCGATc      <  1:1832001/150‑1 (MQ=255)
                                                                                                                                                   ctcAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGATCGCTTCGCTCGATCAAc   >  2:1461703/1‑150 (MQ=255)
                                                                                                                                                    tcAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGATCGCTTCGCTCGATCAAcc  <  1:1729401/150‑1 (MQ=255)
                                                                                                                                                    tcAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGATCGCTTCGCTCGATCAAcc  <  2:1067680/150‑1 (MQ=255)
                                                                                                                                                    tcAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGATCGCTTCGCTCGATCAAcc  <  1:787489/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
CAACGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGATCGCTTCGCTCGATCAACC  >  NC_000913/4409261‑4409558

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