Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,335,123 C→T intergenic (+294/‑25) cysB → / → ymiA DNA‑binding transcriptional dual regulator CysB/uncharacterized protein YmiA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,335,1230CT83.3% 22.9 / ‑2.8 12intergenic (+294/‑25)cysB/ymiADNA‑binding transcriptional dual regulator CysB/uncharacterized protein YmiA
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (5/5);  total (6/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

TATATTCTCTTTTGTTACTTATTTCTATACCGTGCACAATATTTTTATCCCTCATATTTATAGGGTAAATTACCTGTAAAGTAGTGTCTAATTACCCGCCTCGCAAACTCCAGTTTTCAAATATCTATTTCAGGTCAAAAGATTGAATTATTCACTGGAGACGATTCGTAAATTCGCTGGATTTTTGACTAAAGTTTCTTTAGGATTTGTCTCATCGTATGATTAATTACACCAATCGTTCAACGTCTTGCTGATAAGTGATGGCGATTATATGAGGTTAGCAATGCCTTCTGGAAATCAGGAACCCCGCCGCGATCCTGAATTAAAACGTAAAGCCTGGCTGGCGGTTTTTCTTGGTTCTGCACTTTTCTGGGTGGTTGTCGCACTGCTGATTTGGAAAGTGTGGGGATAACTATGGTCGGTCAGGAGCAACTGGAGTCGTCACCATTATGCCAGCATAGTGACAATGAAACAGAAACGAAACGGGAATGTT  >  NC_000913/1334877‑1335369
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      
aatatTTTCTTTTGTTACTTATTTCTATACCGTGCACAAAATTTTTATCCCTCATATTTAAAGGGTAAATTTCCTGTAAAATAGTGTCTAATTACCCGGCTCCCAAACTCCAGTTTTCAAATATCTAGTTCAGGTCAAAAGATTTAATTAATCACTGGAGACGATTCGTAAATTCGCTGTATTTTTGACTAAAGTTGCTTTAGGATTTGTCTCATCGTATGATTAATTACACCAATCGTTCAACGTTTTGc                                                                                                                                                                                                                                                    <  1:51093/250‑1 (MQ=255)
            tGTTACTTATTTCTATACCGTGCACAATATTTTTATCCCTCATATTTATAGGGTAAATTACCTGTAAAGTAGTGTCTAATTACCAGCCTCGCAAACTCCAGTTTTCAAATATCTATTTCAGGTCAAAAGATTGAATTATTCACTGGAGACGATTCGTCAATTCGCTGGATTTTTGACTACAGTTTCTTTAGGATTTGTCTCATCGTATGATTAATTACACCAATCGTTCAATGTTTTGCTGATATGTGATg                                                                                                                                                                                                                                        >  2:474685/1‑251 (MQ=255)
                   tATTTCTATACCGTGCACAATATTTTTATCCCTCATATTTATAGGGTAAATTTCCTGTAAAGTAGTGTCTAATTACCCGCCTCGCAAACTCCAGTTTTCAAATATCTATTTCAGGTCAAAAGATTGAATTATTCACTGGAGACGATTCGTAAATTCGCTGGATTTTTGACTAAAGTTTCTTTAGGATTTGTCTCATCGTATGATTAATTACACCAATCGTTCAACGTTTTGCTGATAAGTGATGGCGAtt                                                                                                                                                                                                                                  <  1:180978/250‑1 (MQ=255)
                      ttCTATACCGTGCACAATATTTTTATCCCTCATATTTATAGGGTAAATTACCTGTAAAGTAGTGTCTAATTACCCGCCTCGCAAACTCCAGTTTTCAAATATCTATTTCAGGTCAAAAGATTGAATTATTCACTGGAGACGATTCGTAAATTCGCTGGATTTTTGACTAAAGTTTCTTTAGGATTTGTCTCATCGTATGATTAATTACACCAATCGTTCAACGTTTTGCTGATAAGTGATGGCGATTatat                                                                                                                                                                                                                              >  1:345964/1‑251 (MQ=255)
                             ccGTGCACAATATTTTTATACCTCATATTTATAGGGTAAATTACCTGTAAAGTAGAGTCTAATTACCCGCCTCGCAAACTCCAGTTTTCAAAAATCTATTTCAGGTCAAAAGATTGCATTCTTCACTGGAGACGATTCGTAAATTCGCTGGATTTTTGAATAAAGTTTCTTTAGGATTTGTCTCATCGTATGATTATTTACACACATCGTTCACCGTTTTGCTGATAAGTGTTGTCGATTATATGAGGtt                                                                                                                                                                                                                        >  2:43121/1‑250 (MQ=255)
                                           tttATCCCTCATATTTATAGGTGAAAGTACCTGTAAAGTAGTTTCTAATTACCCGCCTCGCAAACTCCAGTGTTCAAATATCTATTTCAGGTCAAAAGATTGAATTATTCACTGGAGACGATTCGTAAATTCGCTGGATTTTTGACTAAAGTTTCTTTAGGATTTGTCTCATCGTATGATTAATTACACCAATCGTTCAACGTCTTGCGGATAAGTGATGGCGATTATATGACGTTAGCAATGCCTTCTg                                                                                                                                                                                                          <  2:104276/250‑1 (MQ=255)
                                                      tatTTATAGGGTAAATTACCTGTAAAGTCGTGTCTAATTACCCGCCTCGCAAACTCAAGTTTTCAAATATCTCTTTCAGGTCAAAAGATTGAATTATTCACTGGAGACGATTCGTAAATTCGCTGGATTTTTGACTAAAGTTTCTTTAGGATTTGTCTCATCGTATGATTAGTTACACCAATCGATCAACGTCTTGCTGTTAAGTGATGGCGTTTATATGAGGTTAGCAATGCCTTGTGGCAATCAGGatc                                                                                                                                                                                              >  2:558773/1‑249 (MQ=255)
                                                                                                                     cAAATATCTATTTCAGGTCAAAAGATTGAATTATTCACTGGAGACGATTCGTAAATTCGCTGGATTTTTGACTAAAGTTTCTTTAGGATTTGTCTCATCGTATGATTAATTACACCAATCGTTCAACGTTTTGCTGATAAGTGATGGCGATTATAGGAGGTTAGCAATGCCTTCTGGAAATCAGGAACCCCGCCGCGATCCTGAATTAAAACGTAAAGCCTGGCTGGCGGTTTTTCTTGGTTCTGCACtt                                                                                                                                <  1:305129/250‑1 (MQ=255)
                                                                                                                            cTATTTCAGGTCAAAAGATTGAATTATTCACTGGAGACGATTCGTTAATTCGCTGGATTTTTGACTAAAGTTTCTTTAGGATTTGTCTCATCGTATGATTAATTACACCAATCGTTCAACGTTTTGCTGATAAGTGATGGCGATTATATGAGGTTAGCAATGCCTTCTGGAAATCAGGAACCCCGCCGCGATCCTGAATTAAAACGTAAAGCCTGGCTGGCGGTTTTTCTTTGTTCTGCACTTTTCttggg                                                                                                                        >  1:90444/1‑247 (MQ=255)
                                                                                                                                                                                          gACTAAAGTTTCTTTAGGATTTTTCTCATCGTATGATTAATTACACCAATCGTTCAACGTTTGGCTGATAAGTGATGGCGATTATAGGAGGTTAGCAATGCCTTCTGGAAATCAGGAACCCCGCCGCGATCCTGAATTAAAACGTAAAGCCTGGCTGGCGGTTTTTCTTGGTTCTGCACTTTTCTGGGTGGTTGTCGCACTGCTGATTTGGAAAGTGTGGGGATAACTATGGTCGGTCAGGAGCAACTgg                                                           <  1:233286/250‑1 (MQ=255)
                                                                                                                                                                                                      tttAGGATTTGTCTCATCGTATGATTAAATACACCAATCGTTCAACGTTTTGCTGATAAGTGATTGCGATTATATGAGGTAAGAAATGCCTTCTGGAACTCAGGAACACCTCCGCGATCCTGAATTAAAACGTAAAGCCTGGCTGGCGGTGTTTCTTGGTTTTGCAATTATATGGGAGGATGTCGCACTGGTGGTTTAGAAAGAGTGGGGAAAACTGTGGACGGTTAGGAGCAACAGGAGTAGTGACCAtt                                              >  2:14824/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  aCGTTTTGCTGATAAGTGAGGGCGATTATTTGGGTTTAGCAAGGTCTTCGGGAAATCAGAAATGCCGCCGCGATCCTGAATTAAAACGTAAAGCCTGGCTGGCGGGTTTTCTTGGTTCTGCACTTTTCTGGGTGGTTGTCGCACTGCTGATTTGGAAAGTGTGGGGATAACTATGGTCGGTGAGGAGCAACTGGAGTCGTCACCATTATGCCAGCATAGTGACAATGAAACAGAAACGAAACGGGAATGtt  <  1:369711/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      
TATATTCTCTTTTGTTACTTATTTCTATACCGTGCACAATATTTTTATCCCTCATATTTATAGGGTAAATTACCTGTAAAGTAGTGTCTAATTACCCGCCTCGCAAACTCCAGTTTTCAAATATCTATTTCAGGTCAAAAGATTGAATTATTCACTGGAGACGATTCGTAAATTCGCTGGATTTTTGACTAAAGTTTCTTTAGGATTTGTCTCATCGTATGATTAATTACACCAATCGTTCAACGTCTTGCTGATAAGTGATGGCGATTATATGAGGTTAGCAATGCCTTCTGGAAATCAGGAACCCCGCCGCGATCCTGAATTAAAACGTAAAGCCTGGCTGGCGGTTTTTCTTGGTTCTGCACTTTTCTGGGTGGTTGTCGCACTGCTGATTTGGAAAGTGTGGGGATAACTATGGTCGGTCAGGAGCAACTGGAGTCGTCACCATTATGCCAGCATAGTGACAATGAAACAGAAACGAAACGGGAATGTT  >  NC_000913/1334877‑1335369

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