Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,111,259 Δ1 bp coding (718/900 nt) rfbD ← dTDP‑4‑dehydrorhamnose reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,111,2590G.93.3% 56.2 / ‑1.2 15coding (718/900 nt)rfbDdTDP‑4‑dehydrorhamnose reductase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base . (10/4);  total (11/4)
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.63e-01

TCATCATTCCTTTTAATTCATCTTGCTCCACCATCACGAACAAGATGCAAAAACTATTAAATTGCTGTAGTCGTAAATAATTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGGGAATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGATCGTTAATA  >  NC_000913/2111021‑2111443
                                                                                                                                                                                                                                              |                                                                                                                                                                                        
tcatcaTTCCTTTTAATTCATCTTGCTCCACCATCACGAACAAGATGCAAAAACTATTAAATTGCTGTAGTCGTAAATAATTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTg                                                                                                                                                                             >  2:80435/1‑251 (MQ=255)
             tAATTCATCTTGCTCCACCATCACGAACAAGATCCAAAAACTATTAAATTGTTGTAGTCGTAAATAATTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑ACTGCCTGCTTTGCGCGCCTCTTCa                                                                                                                                                                 >  2:223618/1‑250 (MQ=255)
                                  cACGAACAAGATGCAAAAACTATTAAATTGCTGTAGTCGTAAATAATTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCg                                                                                                                                            <  2:172635/250‑1 (MQ=255)
                                                         tAAATTGCTGTAGTCGTAAATAATTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCt                                                                                                                    >  1:54694/1‑251 (MQ=255)
                                                         tAAATTGCTGTAGTCGTAAATAATTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCt                                                                                                                    >  2:52454/1‑251 (MQ=255)
                                                         tAAATTGCTGTAGTCGTAAAGAATTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGc                                                                                                                     <  1:91098/250‑1 (MQ=255)
                                                           aaTTGTTGGAGTAGTATATACTTCATTGAGCATCCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGAGCACAGTTCTGCTGAAATTTTTCTGTCTTAAGGCGAGCGTTATGTGGACGACGAGATGGTGTAGGATAGGCTGTTGTTGGTAATGCGTTTAGCTTGTTGAGTG‑AAAGGGAATGCCTGCTTTGCGCGCCTATTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGATAc                                                                                                                  >  2:237629/1‑251 (MQ=255)
                                                                     gTCGTAAATAATTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTCGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTAc                                                                                                        >  1:35507/1‑251 (MQ=255)
                                                                            ataatTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGCAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTg                                                                                                 >  2:60240/1‑251 (MQ=255)
                                                                                                  cACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGc                                                                           <  2:3234/251‑1 (MQ=255)
                                                                                                                                         ctgctgAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCATTTc                                    >  1:13675/1‑251 (MQ=255)
                                                                                                                                                                ggCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACt             >  2:143816/1‑251 (MQ=255)
                                                                                                                                                                 gCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTg            >  2:18664/1‑251 (MQ=255)
                                                                                                                                                                   gagaGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGAt          <  1:174493/251‑1 (MQ=255)
                                                                                                                                                                            gtgGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCCGCACCTGTTGGCGCACCAAACTGATCGTtaata  >  1:44256/1‑250 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                        
TCATCATTCCTTTTAATTCATCTTGCTCCACCATCACGAACAAGATGCAAAAACTATTAAATTGCTGTAGTCGTAAATAATTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGGGAATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGATCGTTAATA  >  NC_000913/2111021‑2111443

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