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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,409,4101.A 5.6% 608.0 / 29.2 144coding (136/732 nt)rlmB23S rRNA mG2251 2'‑O‑ribose methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base . (69/67);  new base A (4/4);  total (73/71)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.97e-01

Only 100 of 154 total aligned reads displayed.

ACGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGATCGCTTCGCTCGATCAACC  >  NC_000913/4409263‑4409558
                                                                                                                                                    |                                                                                                                                                    
acGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑a                                                                                                                                                     <  3:882554/149‑1 (MQ=255)
acGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑a                                                                                                                                                     <  2:482108/149‑1 (MQ=255)
  gAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑Agg                                                                                                                                                   <  1:1055915/149‑1 (MQ=255)
  gAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑Agg                                                                                                                                                   <  3:452425/149‑1 (MQ=255)
  gAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑Agg                                                                                                                                                   >  1:799937/1‑149 (MQ=255)
       cATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTg                                                                                                                                              <  4:1499233/149‑1 (MQ=255)
        aTTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTgg                                                                                                                                             <  3:1070689/149‑1 (MQ=255)
        aTTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTgg                                                                                                                                             <  4:353728/149‑1 (MQ=255)
         ttAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGt                                                                                                                                            <  3:1237263/149‑1 (MQ=255)
             tGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATc                                                                                                                                        <  2:184410/149‑1 (MQ=255)
               aGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCa                                                                                                                                      <  3:259610/149‑1 (MQ=255)
                gCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAg                                                                                                                                     >  4:1134874/1‑149 (MQ=255)
                gCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAg                                                                                                                                     >  2:1468501/1‑149 (MQ=255)
                gCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTAACCGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTAAGTCTC‑AGGGCGTGGTTATCCAg                                                                                                                                     >  2:994981/1‑149 (MQ=255)
                    aaTGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTgg                                                                                                                                 >  1:1164604/1‑149 (MQ=255)
                               gCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCaa                                                                                                                      <  3:976076/149‑1 (MQ=255)
                                 aTCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑Agg                                                                                                                                                   <  4:693242/118‑1 (MQ=255)
                                 aTCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑Agg                                                                                                                                                   >  3:693242/1‑118 (MQ=255)
                                    cACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATAtct                                                                                                                 >  2:1153489/1‑149 (MQ=255)
                                     aCGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATAtctc                                                                                                                >  2:588741/1‑149 (MQ=255)
                                     aCGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATAtctc                                                                                                                <  1:243909/149‑1 (MQ=255)
                                        cAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTcgac                                                                                                             <  4:48816/149‑1 (MQ=255)
                                        cAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTcgac                                                                                                             <  2:1563478/149‑1 (MQ=255)
                                             cAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGaa                                                                                                        >  4:393858/1‑149 (MQ=255)
                                              aGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑a                                                                                                                                                     >  2:157428/1‑103 (MQ=255)
                                               ggCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGaaaa                                                                                                      <  3:1134874/149‑1 (MQ=255)
                                                gCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAg                                                                                                     >  2:1234612/1‑149 (MQ=255)
                                                 cccTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGc                                                                                                    <  3:275967/149‑1 (MQ=255)
                                                  ccTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATAtc                                                                                                                  <  1:550572/134‑1 (MQ=255)
                                                  ccTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCg                                                                                                   >  4:1559867/1‑149 (MQ=255)
                                                   ctgctgGAGCGCGCCCCTGAACGTTTTCAGGACGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGa                                                                                                  <  4:1279173/149‑1 (MQ=255)
                                                    tgctgGAGCGCGCCCCCGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGAc                                                                                                 <  3:361908/149‑1 (MQ=255)
                                                       tgGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACgg                                                                                               >  4:1123405/1‑148 (MQ=255)
                                                            cgcgcCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGt                                                                                         <  1:1114008/149‑1 (MQ=255)
                                                              cgcCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGc                                                                                       >  3:736724/1‑149 (MQ=255)
                                                               gcCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCa                                                                                      >  1:174447/1‑149 (MQ=255)
                                                                 cccTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGc                                                                                            <  4:1284298/141‑1 (MQ=255)
                                                                 cccTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATc                                                                                    >  2:684912/1‑149 (MQ=255)
                                                                 cccTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATc                                                                                    >  3:148501/1‑149 (MQ=255)
                                                                  ccTGAACGTGTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGTCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCa                                                                                   <  4:964081/149‑1 (MQ=255)
                                                                   cTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCa                                                                                   >  2:440245/1‑148 (MQ=255)
                                                                   cTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCa                                                                                   <  1:440245/148‑1 (MQ=255)
                                                                    tGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAgg                                                                                 >  4:553841/1‑149 (MQ=255)
                                                                     gAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGc                                                                                       <  1:1535046/142‑1 (MQ=255)
                                                                        cGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCAt                                                                             <  2:911317/149‑1 (MQ=255)
                                                                          ttttCAGGAAGTCTTTACTTGAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGCCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTa                                                                           <  4:736724/149‑1 (MQ=255)
                                                                                 gAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAat                                                                                                                     <  1:1325902/100‑1 (MQ=255)
                                                                                 gAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCcg                                                                    >  2:1485645/1‑149 (MQ=255)
                                                                                  aaGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCGCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCcgc                                                                   >  1:1040832/1‑149 (MQ=255)
                                                                                      cTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGa                                                               >  2:271369/1‑149 (MQ=255)
                                                                                       tttATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGaaaa                                                                                                      >  3:1104904/1‑109 (MQ=255)
                                                                                       tttATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGaaaa                                                                                                      <  4:1104904/109‑1 (MQ=255)
                                                                                       tttATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGaa                                                              >  1:1089395/1‑149 (MQ=255)
                                                                                       tttATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGaa                                                              >  1:758641/1‑149 (MQ=255)
                                                                                       tttATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGaa                                                              <  3:371067/149‑1 (MQ=255)
                                                                                       tttATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGaa                                                              <  2:592855/149‑1 (MQ=255)
                                                                                         tATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGc                                                            <  2:291789/149‑1 (MQ=255)
                                                                                           ttttAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCa                                                          <  1:1153489/149‑1 (MQ=255)
                                                                                            tttAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAg                                                         <  2:207231/149‑1 (MQ=255)
                                                                                               aaaaGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGAc                                                      <  4:664674/149‑1 (MQ=255)
                                                                                                  aGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTc                                                   <  4:752976/149‑1 (MQ=255)
                                                                                                  aGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTc                                                   >  3:1292834/1‑149 (MQ=255)
                                                                                                  aGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTc                                                   >  1:1566965/1‑149 (MQ=255)
                                                                                                   ggCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGACa                                                  <  1:1351904/149‑1 (MQ=255)
                                                                                                    gCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAg                                                 <  3:553841/149‑1 (MQ=255)
                                                                                                       gTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGaaaa                                                                                                      <  1:446505/93‑1 (MQ=255)
                                                                                                       gTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTAt                                              >  4:559136/1‑149 (MQ=255)
                                                                                                         gAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCa                                            <  4:523131/149‑1 (MQ=255)
                                                                                                          aaGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAg                                           <  3:266928/149‑1 (MQ=255)
                                                                                                           aGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAgg                                          <  3:1515949/149‑1 (MQ=255)
                                                                                                              tAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGaaa                                       >  1:883595/1‑149 (MQ=255)
                                                                                                               aaaCGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGaaaa                                      <  2:747774/149‑1 (MQ=255)
                                                                                                               aaaCGTCTGTTACCGCTGATTCACGCCCTTGAGGCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGt                                                                                              >  4:430255/1‑93 (MQ=255)
                                                                                                                aaCGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAAc                                     <  1:671388/149‑1 (MQ=255)
                                                                                                                    tCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCa                                            <  3:1491004/138‑1 (MQ=255)
                                                                                                                    tCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCa                                            >  4:1491004/1‑138 (MQ=255)
                                                                                                                    tCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATc                                 >  2:953271/1‑149 (MQ=255)
                                                                                                                       gTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGc                              >  3:499012/1‑149 (MQ=255)
                                                                                                                           ccGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGt                                                                 >  3:403860/1‑110 (MQ=255)
                                                                                                                           ccGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGa                          <  1:994981/149‑1 (MQ=255)
                                                                                                                            cGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGAt                         >  2:1287185/1‑149 (MQ=255)
                                                                                                                             gCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGTAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATc                        <  3:63252/149‑1 (MQ=255)
                                                                                                                             gCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATc                        <  3:566285/149‑1 (MQ=255)
                                                                                                                             gCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATc                        <  1:904337/149‑1 (MQ=255)
                                                                                                                             gCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATc                        >  4:1226787/1‑149 (MQ=255)
                                                                                                                                 aTTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGAt                    <  1:30543/149‑1 (MQ=255)
                                                                                                                                  ctCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGAtc                   <  3:138323/148‑1 (MQ=255)
                                                                                                                                   ccACGCCC‑TGAGACTCAAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATCATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGaaaa                                      <  3:631069/128‑1 (MQ=255)
                                                                                                                                   ccACGCCC‑TGAGACTCAAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATCATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGaaaa                                      >  4:631069/2‑129 (MQ=255)
                                                                                                                                    cACGCCC‑TGAGACTCAAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGAtcg                  >  4:603578/1‑148 (MQ=255)
                                                                                                                                    cACGCCC‑TGAGACTCAAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGAtc                   <  2:671142/147‑1 (MQ=255)
                                                                                                                                    cACGCCC‑TGAGACTCAAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGAtc                   >  1:671142/1‑147 (MQ=255)
                                                                                                                                     aCGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGAtcgct                >  2:414829/1‑149 (MQ=255)
                                                                                                                                       gCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGaa                                                              <  2:370860/101‑1 (MQ=255)
                                                                                                                                         ccTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGAtcgcttcgc            >  2:480390/1‑149 (MQ=255)
                                                                                                                                          cTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGAtcgcttcgct           >  4:1065523/1‑149 (MQ=255)
                                                                                                                                             gAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGATCGCTTCGCTCGa        <  2:799937/149‑1 (MQ=255)
                                                                                                                                                  tc‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGATCGCTTCGCTCGATCAAc   <  4:147335/149‑1 (MQ=255)
                                                                                                                                                  tc‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGATCGCTTCGCTCGATCAAc   <  1:513445/149‑1 (MQ=255)
                                                                                                                                                   c‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGATCGCTTCGCTCGATCAAcc  >  2:1100743/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ACGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGATCGCTTCGCTCGATCAACC  >  NC_000913/4409263‑4409558

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑