Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,094,683 C→T intergenic (‑411/+40) rhaD ← / ← rhaA rhamnulose‑1‑phosphate aldolase/L‑rhamnose isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,094,6830CT92.3% 32.5 / NA 13intergenic (‑411/+40)rhaD/rhaArhamnulose‑1‑phosphate aldolase/L‑rhamnose isomerase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (6/6);  minor base A (1/0);  total (7/6)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTT  >  NC_000913/4094575‑4094817
                                                                                                            |                                                                                                                                      
aCAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCGACAACTATGGCCTACCAGTAGGCCGGATAAAGCGCTCGCGCCGCATcc                                                                                                         >  2:213964/1‑140 (MQ=11)
                     gCCGGATAAGGCACTTGTGCCGTATCCGGCATTCATTGCTTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGa                                                                                    <  2:312383/140‑1 (MQ=14)
                                  cTTGTGCCGCATCCGGCAATCAATGCCTGATGCGATGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAAtttctt                                                                       <  1:596662/140‑1 (MQ=255)
                                        ccGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCGCTCGCACCGCATCCGGCATTGTTCACCCGCGACGACAAAAAATTTATTTCTCAt                                                                 >  2:378686/1‑140 (MQ=14)
                                                           ccTGATGCGACGCTGACGCGTCGTAGCAGTCCTACAATTATTGCCGACCTGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAg                                              <  1:378686/140‑1 (MQ=25)
                                                               aTGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCAt                                          >  1:123735/1‑140 (MQ=255)
                                                                tGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCAtt                                         <  2:123735/140‑1 (MQ=255)
                                                                         gTCGCGTCTTATCAGTCCTACAGCTATTGCCGACCTGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACctg                                <  1:339419/140‑1 (MQ=255)
                                                                                          cTACAACTATTGCCGACCTGTAAGCCGGACAAGGCGCACGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGt               <  1:213964/140‑1 (MQ=39)
                                                                                             cAACTATTGCCGACCTGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATCCGCGACATGCTGGCGTATCGTGACGTTgg            >  1:299595/1‑140 (MQ=38)
                                                                                                  attGCCGACCTGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATa       >  1:164953/2‑140 (MQ=255)
                                                                                                     gCCGACCTGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGGCTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACAt    >  2:427743/1‑140 (MQ=255)
                                                                                                       cGACCTGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACAttt  >  2:739325/1‑140 (MQ=255)
                                                                                                            |                                                                                                                                      
ACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTT  >  NC_000913/4094575‑4094817

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: