Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,409,407 T→A 5.5% S45T (TCT→ACT)  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,4070TA 5.5% 590.5 / 15.3 175S45T (TCT→ACT) rlmB23S rRNA mG2251 2'‑O‑ribose methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base T (83/82);  new base A (8/2);  total (91/84)
Fisher's exact test for biased strand distribution p-value = 1.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.26e-01

Only 100 of 182 total aligned reads displayed.

AACAACGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGATCGCTTCGCTCGATCA  >  NC_000913/4409259‑4409555
                                                                                                                                                    |                                                                                                                                                     
  caacGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGtct                                                                                                                                                     >  2:1699531/1‑149 (MQ=255)
   aacGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGtctc                                                                                                                                                    >  2:389230/1‑149 (MQ=255)
   aacGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGtctc                                                                                                                                                    >  3:1339493/1‑149 (MQ=255)
    acGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑a                                                                                                                                                  <  3:17490/149‑1 (MQ=255)
           cATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTg                                                                                                                                           <  2:1099025/149‑1 (MQ=255)
           cATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTg                                                                                                                                           <  3:863180/149‑1 (MQ=255)
           cATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTg                                                                                                                                           <  1:1705151/149‑1 (MQ=255)
            aTTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTgg                                                                                                                                          >  2:589416/1‑149 (MQ=255)
               aaTGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTa                                                                                                                                       >  1:199967/1‑149 (MQ=255)
                 tGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATc                                                                                                                                     <  1:1226958/149‑1 (MQ=255)
                  gAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATcc                                                                                                                                    >  1:431244/1‑149 (MQ=255)
                     cGAAATGATTTACGGCATCCACGCAGTTCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGt                                                                                                                                 >  3:1530344/1‑149 (MQ=255)
                      gAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGtt                                                                                                                                <  4:1295402/149‑1 (MQ=255)
                      gAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGtt                                                                                                                                <  1:598307/149‑1 (MQ=255)
                        aaTGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGCACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTgg                                                                                                                              <  2:133377/149‑1 (MQ=255)
                        aaTGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTgg                                                                                                                              >  2:771354/1‑149 (MQ=255)
                        aaTGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAAGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTgg                                                                                                                              >  3:996051/1‑149 (MQ=255)
                           gATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGTAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCaa                                                                                                                           >  2:866788/1‑149 (MQ=255)
                           gATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCaa                                                                                                                           >  4:1704243/1‑149 (MQ=255)
                                 cGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGcc                                                                                                                     <  2:199967/149‑1 (MQ=255)
                                   gCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCaa                                                                                                                   <  2:1585738/149‑1 (MQ=255)
                                     aTCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAata                                                                                                                 >  4:909310/1‑149 (MQ=255)
                                        cACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTGTTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATAtct                                                                                                              >  1:1129326/1‑149 (MQ=255)
                                        cACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGtt                                                                                                                                <  3:360700/131‑1 (MQ=255)
                                         aCGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAGGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATAtctc                                                                                                             >  2:528620/1‑149 (MQ=255)
                                         aCGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATAtctc                                                                                                             >  4:384443/1‑149 (MQ=255)
                                         aCGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACATATTCAGGAAGTCTTAATGCTAAACGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTAAGTCTC‑AGGGCGTGGTTCTCCAGTTAGCAAACCGCCAATAactc                                                                                                             >  2:1370021/1‑149 (MQ=255)
                                           gCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTcga                                                                                                           <  4:566334/149‑1 (MQ=255)
                                              gTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTcgacga                                                                                                        <  4:343893/149‑1 (MQ=255)
                                                   ggCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGaaaa                                                                                                   <  3:922984/149‑1 (MQ=255)
                                                        tgctgGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGAc                                                                                              <  3:1235410/149‑1 (MQ=255)
                                                              aGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGcc                                                                                        <  1:1061197/149‑1 (MQ=255)
                                                                cgcgcCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGt                                                                                      >  2:1301729/1‑149 (MQ=255)
                                                                   gcCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCa                                                                                   >  4:195252/1‑149 (MQ=255)
                                                                        tGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAgg                                                                              <  1:689388/149‑1 (MQ=255)
                                                                          aaCGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGc                                                                            <  4:1339493/149‑1 (MQ=255)
                                                                                   aGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGcc                                                                   >  4:683210/1‑149 (MQ=255)
                                                                                     gAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCcg                                                                 >  2:1515440/1‑149 (MQ=255)
                                                                                     gAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCcg                                                                 >  1:1527051/1‑149 (MQ=255)
                                                                                       aGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCcgcg                                                               >  1:14108/1‑149 (MQ=255)
                                                                                        gTCTTTATTTTAAAAGGCCGTGTAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGt                                                              <  1:866788/149‑1 (MQ=255)
                                                                                           tttATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGaa                                                           <  2:839973/149‑1 (MQ=255)
                                                                                           tttATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGaa                                                           >  3:549531/1‑149 (MQ=255)
                                                                                           tttATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGaa                                                           >  2:1349328/1‑149 (MQ=255)
                                                                                            ttATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAg                                                          >  1:303289/1‑149 (MQ=255)
                                                                                               ttttAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCa                                                       >  3:891781/1‑149 (MQ=255)
                                                                                               ttttAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCa                                                       >  2:674820/1‑149 (MQ=255)
                                                                                               ttttAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCa                                                       >  1:366545/1‑149 (MQ=255)
                                                                                                 ttAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAgg                                                     >  1:36923/1‑149 (MQ=255)
                                                                                                 ttAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAgg                                                     <  4:580112/149‑1 (MQ=255)
                                                                                                      aGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTc                                                <  1:867672/149‑1 (MQ=255)
                                                                                                       ggCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTcgac                                                                                                          >  3:1245309/1‑90 (MQ=255)
                                                                                                       ggCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTcgac                                                                                                          <  4:1245309/90‑1 (MQ=255)
                                                                                                       ggCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCa                                               <  2:366545/149‑1 (MQ=255)
                                                                                                        gCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAg                                              <  3:195252/149‑1 (MQ=255)
                                                                                                            tGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATc                                          <  2:303289/149‑1 (MQ=255)
                                                                                                             gAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCa                                         <  3:1372366/149‑1 (MQ=255)
                                                                                                             gAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCa                                         <  2:538467/149‑1 (MQ=255)
                                                                                                                   aaaCGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGaaaa                                   >  2:501339/1‑149 (MQ=255)
                                                                                                                     aCGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCa                                               <  2:1690112/135‑1 (MQ=255)
                                                                                                                     aCGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCa                                               >  1:1690112/1‑135 (MQ=255)
                                                                                                                      cGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGa                                <  4:662869/149‑1 (MQ=255)
                                                                                                                       gTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGACGTGCATCAGGGCATTAACGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGAt                               >  3:972425/1‑149 (MQ=255)
                                                                                                                        tCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCGAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATc                              <  2:1079461/149‑1 (MQ=255)
                                                                                                                        tCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATc                              <  1:589416/149‑1 (MQ=255)
                                                                                                                        tCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATc                              <  1:835452/149‑1 (MQ=255)
                                                                                                                        tCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATc                              >  3:375357/1‑149 (MQ=255)
                                                                                                                            ttACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGcc                          >  4:1411122/1‑149 (MQ=255)
                                                                                                                                cGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAggg                                                                             <  4:1197058/94‑1 (MQ=255)
                                                                                                                                cGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGAt                      <  1:51147/149‑1 (MQ=255)
                                                                                                                                cGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGAt                      >  1:1346973/1‑149 (MQ=255)
                                                                                                                                 gCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGacc                     <  2:1396736/149‑3 (MQ=255)
                                                                                                                                 gCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATc                     <  2:120241/149‑1 (MQ=255)
                                                                                                                                  gataaccACGCCC‑TGAGACTCAAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCAATATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAAc                                  >  3:5880/7‑135 (MQ=255)
                                                                                                                                      ttCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGAtc                >  1:665067/1‑149 (MQ=255)
                                                                                                                                       tCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGAtcg               >  3:1403155/1‑149 (MQ=255)
                                                                                                                                       ccACGCCC‑TGAGACTCAAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGAtcg               >  4:948759/2‑149 (MQ=255)
                                                                                                                                        cACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGAc                                                                                              >  4:1563182/1‑69 (MQ=255)
                                                                                                                                        cACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGAc                                                                                              <  3:1563182/69‑1 (MQ=255)
                                                                                                                                        cACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGAt                               <  1:655578/132‑1 (MQ=255)
                                                                                                                                        cACGCCC‑TGAGACTCAAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGAtcgc              >  3:848688/1‑149 (MQ=255)
                                                                                                                                         aCGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGGGTGAAGCCAg                                                      <  3:1045022/108‑1 (MQ=255)
                                                                                                                                         aCGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAg                                                      >  4:1045022/1‑108 (MQ=255)
                                                                                                                                         aCGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATc                     <  2:569343/141‑1 (MQ=255)
                                                                                                                                         aCGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGAtcgct             >  1:557489/1‑149 (MQ=255)
                                                                                                                                         aCGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGAtcgct             >  4:1370174/1‑149 (MQ=255)
                                                                                                                                         aCGCCC‑TGAGACTCAAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGAtcgct             >  4:1066543/1‑149 (MQ=255)
                                                                                                                                         aCGCCC‑TGAGACTCAAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGAtcgct             >  3:1595216/1‑149 (MQ=255)
                                                                                                                                           gCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGAtcgcttc           <  4:1609487/149‑1 (MQ=255)
                                                                                                                                             ccTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGAtcgcttcgc         >  2:566527/1‑149 (MQ=255)
                                                                                                                                             ccTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGAtcgcttcgc         >  4:1093780/1‑149 (MQ=255)
                                                                                                                                              cTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCt                             <  2:1577112/128‑1 (MQ=255)
                                                                                                                                              cTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCt                             >  1:1577112/1‑128 (MQ=255)
                                                                                                                                              cTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGAtcgcttcgct        >  3:293628/1‑149 (MQ=255)
                                                                                                                                              cTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGAtcgcttcgct        <  3:749658/149‑1 (MQ=255)
                                                                                                                                               ttGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCa                                         >  3:1135544/1‑115 (MQ=255)
                                                                                                                                               ttGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCa                                         <  4:1135544/115‑1 (MQ=255)
                                                                                                                                               ttGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGATCGCTTCGCTc       >  4:249432/1‑149 (MQ=255)
                                                                                                                                                 gAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGATCGCTTCGCTCGa     <  3:118384/149‑1 (MQ=255)
                                                                                                                                                 gAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGATCGCTTCGCTCGa     <  2:813209/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AACAACGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGATCGCTTCGCTCGATCA  >  NC_000913/4409259‑4409555

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