Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,106,285 T→C F156L (TTT→CTT)  ymdB → 2'‑O‑acetyl‑ADP‑ribose deacetylase, regulator of RNase III activity

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,106,2850TC96.8% 100.9 / ‑4.3 31F156L (TTT→CTT) ymdB2'‑O‑acetyl‑ADP‑ribose deacetylase, regulator of RNase III activity
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (18/12);  total (19/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

GTGGCTTTTCCTGCAATCAGTACTGGGGTTTATGGTTACCCTCGTGCGGCAGCGGCTGAAATCGCAGTAAAAACCGTTTCAGAATTTATTACCCGTCACGCTTTACCCGAACAGGTATACTTTGTCTGTTATGATGAAGAAAACGCCCACCTCTACGAAAGACTCCTTACCCAACAAGGAGATGAATGATTTGCCCCGGCTGGCGAGCGCGGTGCTGCCACTGTGTTCGCAACATCCCGGTCAGTGTG  >  NC_000913/1106165‑1106412
                                                                                                                        |                                                                                                                               
gTGGCTTTTCCTGCAATCAGTACTGGGGTTTATGGTTACCCTCGTGCGGCAGCGGCTGAAATCGCAGTAAAAACCGTTTCAGAATTTATTAGCCGTCACGGTTTGCCCGAACAGGTAGACTTCGTCTGTTatg                                                                                                                     >  2:193063/1‑133 (MQ=255)
gTGGCTTTTCCTGCAATCAGTACTGGGGTTTATGGTTACCCTCGTGCGGCAGCGGCTGAAATCGCAGTAAAAACCGTTTCAGAATTTATTACCCGTCACGCTTTACCCGAACAGGTATACCTTGTCTGTTatg                                                                                                                     >  1:20110/1‑133 (MQ=255)
gTGGCTTTTCCTGCAATCAGTACTGGGGTTTATGGTTACCCTCGTGCGGCAGCGGCTGAAATCGCAGTAAAAACCGTTTCAGAATTTATTACCCGTCACGCTTTACCCGAACAGGTATACCTTGTCTGTTatg                                                                                                                     >  1:451491/1‑133 (MQ=255)
gTGGCTTTTCCTGCAATCAGTACTGGGGTTTATGGTTACCCTCGTGCGGCAGCGGCTGAAATCGCAGTAAAAACCGTTTCAGAATTTATTACCCGTCACGCTCTACCCGAGCAGATAGACCTTGACTGTTatg                                                                                                                     >  1:165212/1‑133 (MQ=255)
       ttCCTGCAATCAGTACTGGGGTTTATGGTTACCCTCGTGCGGCAGCGGCTGAAATCGCAGTAAAAACCGTTTCAGAATTTATTACCCGTCACGCTTTACCCGAACAGGTATACCTTGTCTGTTATGATgaaga                                                                                                              <  2:20110/133‑1 (MQ=255)
           tGCAATCAGTACTGGGGTTTATGGTTACCCTCGTGCGGCAGCGGCTGAAATCGCAGTAAAAACCGTTTCAGAATTTATTACCCGTCACGCTTTACCCGAACAGGTATACCTTGTCTGTTATGATGAACAAAAc                                                                                                          >  2:130368/1‑133 (MQ=255)
                           gTTTATGGTTACCCTAGTGCGGCACCGGATGACATCGCAGTCAAAACCGTTTCAGAATTTATTACCCGTCACGCTTTACCCGAACAGGTATACCTTGTCTGTTATGATGAAGAAAACGCCCACCTCTACGaaa                                                                                          <  2:165212/133‑1 (MQ=255)
                              tATGGTTACCCTCGTGCGGCAGCGGCTGAAATCGCAGTAAAAACCGTTTCAGAATTTATTACCCGTCACGCTTTACCCGAACAGGTATACCTTGTCTGTTATGATGAAGAAAACGCCCACCTCTACGAAAGAc                                                                                       >  1:153727/1‑133 (MQ=255)
                              tATGGTTACCCTCGTGCGGCAGCGGCTGAAATCGCAGTAAAAACCGTTTCAGAATTTATTACCCGTCACGCTTTACCCGAACAGGTATACCTTGTCTGTTATGATGAAGAAAACGCCCACCTCTACGAAAGAc                                                                                       >  2:425854/1‑133 (MQ=255)
                              tATGGTTACCCTCGTGCGGCAGCGGCTGAAATCGCAGTAAAAACCGTTTCAGAATTTATTACCCGTCACGCTTTACCCGAACAGGTATACCTTGTCTGTTATCATGAAGATAACTCCTACCTCTACGAAAGAc                                                                                       >  1:289130/1‑133 (MQ=255)
                               aTGGTTACCCTCGTGCGGCAGCGGCTGAAATCGCAGTAAAAACCGTTTCAGAATTTATTACCCGTCACGCTTTACCCGAACAGGTATACCTTGTCTGTTATGATGAAGAAAACGCCCACCTCTACGAAAGACt                                                                                      <  2:398021/133‑1 (MQ=255)
                                      cccTCGTGCGGCAGCGGCTGAAATCGCAGTAAAAACCGTTTCAGAATTTATTACCCGTCACGCTTTACCCGAACAGGTATACCTTGTCTGTTATGATGAAGAAAACGCCCACCTCTACGAAAGACTCCTTAcc                                                                               <  1:471198/133‑1 (MQ=255)
                                      cccTCGTGCCGCAGCAGCTGAAATCGCAGTAAGAACCGTTTCAGAATTTATTACCCGTCACGCTTTACCCGAACAGGTATACCTTGTCTGTTATGATGAAGAAAACGCCCACCTCTACGAAAGACTCCTTAcc                                                                               <  2:289130/133‑1 (MQ=255)
                                                 cAGCGGCTGAAATCGCAGTAAAAACCGTTTCAGAATTTATTACCCGTCACGCTTTACCCGAACAGGTATACCTTGTCTGTTATGATGAAGAAAACGCCCACCTCTACGAAAGACTCCTTACCCAACAAGgaga                                                                    >  1:575431/1‑133 (MQ=255)
                                                    cGGCTGAAATCGCAGTAAAAACCGTTTCAGAATTTATTACCCGTCACGCTTTACCCGAACAGGTATACCTTGTCTGTTATGATGAAGAAAACGCCCACCTCTACGAAAGACTCCTTACCCAACAAGGAGatga                                                                 >  2:509789/1‑133 (MQ=255)
                                                                             ttCAGAATTTATTACCCGTCACGCTTTACCCGAACAGGTATACCTTGTCTGTTATGATGAAGAAAACGCCCACCTCTACGAAAGACTCCTTACCCAACAAGGAGATGAATGATTTGCCCCAGCTGGCGAgcgc                                        <  2:451491/133‑1 (MQ=255)
                                                                                     ttattaCCCGTCACGCTTTACCCGAACAGGTATACCTTGTCTGTTATGATGAAGAAAACGCCCACCTCTACGAAAGACTCCTTACCCAACAAGGAGATGAATGATTTGCCCCGGCTGGCGAGCGCGGtgctgc                                >  1:484547/1‑133 (MQ=255)
                                                                                          aCCCGTCACGCTTTACCCGAACAGGTATACCTTGTCTGTTATGATGAAGAAAACGCCCACCTCTACGAAAGACTCCTTACCCAACAAGGAGATGAATGATTTGCGCCGGCTGGCGAGCGCGGAGCTGCCACtg                           >  2:224332/1‑133 (MQ=255)
                                                                                                cACGCTTTACCCGAACAGGTATACCTTGTCTGTTATGATGAAGAAAATGCCCACCTCTACGAAAGACTCCTTACCCAACAAGGAGATGAATGATTTGCCCCGGCTGGCGAGCGCGGTGCTGCCACTGTGTTCg                     >  2:518798/1‑133 (MQ=255)
                                                                                                cACGCTTTACCCGAACAGGTATACCTTGTCTGTTATGATGAAGAAAACGCCCACCTCTACGAAAGACTCCTTACCCAACAAGGAGATGAATGATTTGCCCCGGCTGGCGAGCGCGGTGCTGCCACTGTGTTCg                     >  2:142262/1‑133 (MQ=255)
                                                                                                cACGCTTTACCCGAACAGGTATACCTTGTCTGTTATGATGAAGAAAACGCCCACCTCTACGAAAGACTCCTTACCCAACAAGGAGATGAATGATTTGCCCCGACTGGCGAGCGCGGTGCTGCCACTGTGTTCg                     >  1:211712/1‑133 (MQ=255)
                                                                                                  cGCTTTACCCGAACAGGTATACCTTGTCTGTTATGATGAAGAAAACGCCCACCTCTACGAAAGACTCCTTACCCAACAAGGAGATGAATGATTTGCCCCGGCTGGCGAGCGCGGTGCTGCCACTGTGTTCGCa                   <  2:538633/133‑1 (MQ=255)
                                                                                                  cGCTTTACCCGAACAGGTATACCTTGTCTGTTATGATGAAGAAAACGCCCACCTCTACGAAAGACTCCTTACCCAACAAGGAGATGAATGATTTGCCCCGGCTGGCGAGCGCGGTGCTGCCACTGTGTTCGCa                   <  2:575431/133‑1 (MQ=255)
                                                                                                           cGAACAGGTATACCTTGTCTGTTATGATTAAGTAAACGCCCGCCTCTACGAAAGACTCCTTACCCAACAAGGAGATGAATGATTTGCCCCGGCTGGCGAGCGCGGTGCTGCCACTGTGTTCGCAACATCCTgg          <  2:206837/133‑1 (MQ=255)
                                                                                                            gAACAGGTATACCTTGTCTGTTATGATGAAGAAAACGCCCACCTCTACGAAAGACTCCTTACCCAACAAGGAGATGAATGATTTGCCCCGGCTGGCGAGCGCGGTGCTGCCACTGTGTTCGCAACATCCCGGt         <  1:509789/133‑1 (MQ=255)
                                                                                                             aaCAGGTATACCTTGTCTGTTATGATGAAGAAAACGCCCACCTCTACGAAAGACTCCTTACCCAACAAGGAGATGAATGATTTGCCCCGGCTGGCGAGCGCGGTGCTGCCACTGTGTTCGCAACATCCCGGTc        >  2:587826/1‑133 (MQ=255)
                                                                                                             aaCAGGTATACCTTGTCTGTTATGATGAAGAAAACGCCCACCTCTACGAAAGACTCCTTACCCAACAAGGAGATGAATGATTTGCCCCGGCTGGCGAGCGCGGTGCTGCCACTGTGTTCGCAACATCCCGGTc        >  2:149500/1‑133 (MQ=255)
                                                                                                             aaCAGGTATACCTTGTCTGTTATGATGAAGAAAACGCCCACCTCTACGAAAGACTCCTTACCCAACAAGGAGATGAATGATTTGCCCCGGCTGGCGAGCGCGGTGCTGCCACTGTGTTCGCAACATCCCGGTc        >  1:334568/1‑133 (MQ=255)
                                                                                                             aaCAGGTATACCTTGTCTGTTATGATGAAGAAAACGCCCACCTCTACGAAAGACTCCTTACCCAACAAGGAGATGAATGATTTGCCCCGGCTGGCGAGCGCGCTGCTGCCCCTGTGTACGCAACCTCCCGGTc        >  1:206837/1‑133 (MQ=255)
                                                                                                                 ggTATACCTTGTCTGTTATGATGAAGAAAACGCCCACCTCTACGAAAGACTCCTTACCCAACAAGGAGATGAATGATTTGCCCCGGCTGGCGAGCGCGGTGCTGCCACTGTGTTCGCAACATCCCGGTCAgtg    >  1:53639/1‑133 (MQ=255)
                                                                                                                  gTATACCTTGTCTGTTATGATGAAGAAAACGCCCACCTCTACGAAAGACTCCTTACCCAACAAGGAGATGAATGATTTGCCCCGGCTGGCGAGCGCGGTGCTGCCACTGTGTTCGCAACATCCCGGTCAgtgt   <  1:142262/133‑1 (MQ=255)
                                                                                                                   tataCCTTGTCTGTTATGATGAAGAAAACGCCCACCTCTACGAAAGACTCCTTACCCAACAAGGAGATGAATGATTTGCCCCGGCTGGCGAGCGCGGTGCTGCCACTGTGTTCGCAACATCCCGGTCAGtgtg  <  2:463465/133‑1 (MQ=255)
                                                                                                                        |                                                                                                                               
GTGGCTTTTCCTGCAATCAGTACTGGGGTTTATGGTTACCCTCGTGCGGCAGCGGCTGAAATCGCAGTAAAAACCGTTTCAGAATTTATTACCCGTCACGCTTTACCCGAACAGGTATACTTTGTCTGTTATGATGAAGAAAACGCCCACCTCTACGAAAGACTCCTTACCCAACAAGGAGATGAATGATTTGCCCCGGCTGGCGAGCGCGGTGCTGCCACTGTGTTCGCAACATCCCGGTCAGTGTG  >  NC_000913/1106165‑1106412

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