Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,201,548 A→C T18P (ACC→CCA)  ymfI → uncharacterized protein YmfI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,201,5480AC95.0% 53.7 / ‑2.4 20coding (52/342 nt)ymfIuncharacterized protein YmfI
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base C (9/10);  minor base . (1/0);  total (10/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.53e-01

TTGTCAGCCTCAACCTTAACCCGCGTGCGCTGAGCCGCATTCACCACAACGCTAAGGATTCTCTCTGGTTGAAAATACTTAGCTGTTATGTGCCTGTCTTTTCACCACTTCAGGCTCGGTGGTATGCTGGAGTTCTCACACAGCCAGCAAGCAAGGAAACTTAATGAACCAGTTTTATGTTCACGTTCGTCTATTTGAAGACACAGCCGAACAGACCAAAAAATTTGAAGAATTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACATTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAACAAAGCGCATGCTTAC  >  NC_000913/1201334‑1201761
                                                                                                                                                                                                                      |                                                                                                                                                                                                                      
ttGTCAGCCTCAACCTTAACCCGCGTGCGCTTAGCCGCATTCACCACAACGCTAAGGATTCTCTCTGGTTGAAAATACTTAGCTGTTATGTGCCTGTCTTTTCACCACTTCAGGCTCGGTGGTATGCTGGAGTTCTCACACAGCCAGCAAGCAAGGAAACTTAATGAACCAGTTTTATGTTCACGTTCGTCTATTTGAAGACACAGCCGAACAGCCAAAAAAATTTGAAGAATTAATGCTTAACTTTCTg                                                                                                                                                                                    <  1:349419/250‑1 (MQ=255)
            aCCTTAAACCGCGTGCGCTGCGCCGCATTCACCACACCGCTAAGGATTCTCTCTGGTTGAAAATACTTAGCTGTTATGTGCCTGTCTTTTCACCACTTCAGGCTCGGTGGTATGCTTGAGTTCTCACACAGCCAGCAAGCAAGGAAACTTAATGAACCAGTTTTATGTTCACGTTCGTCTATTTGAAGACACAGCCGAACAGCCAAAAAAATTTGAAGAATTAATGCTTAACTTTCTGTACCAGAAAACAg                                                                                                                                                                       <  2:396651/251‑1 (MQ=255)
                              ggAGCCGCACTCACCACAACGCTAAGAATTCCCACTGGGTGAAAAGACTTAGCTGTTATGGGCCTGTCTTTTCAACACGTCAGGCTCGGTGGTATGCTGGAGTTCGCACACAGCCAGCAAGCAAGGAAACTTAATGAACCAGTTTTATGTTCACGTTCGTCTATTTGAAGACACAGCCGAACAGCCAAAAAAATTTGAAGAATTAATGCTTAACGGTCTGTACCAGAAATCAGTTAAAGAGTCTGCCGATa                                                                                                                                                     <  2:173953/250‑1 (MQ=255)
                                                                  ggTTGAAAATACGTAGCTGTTATGTGCCTGTCTTTTCACCACTTCAGGCTCGGTGGTATGCTGGAGTTCTCACACAGCCAGCAAGCAAGGAAACTTAATGAACCAGTTTTATGTTCACGTTCGTCTATTTGAAGACACAGCCGAACAGCCAAAAAAATTTGAAGAATTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTc                                                                                                                  <  2:182790/250‑1 (MQ=255)
                                                                       aaaaaaCTTAGCTGTTATGTGCCTGTCTTTTCACCACTTCAGGCTCGATGGTATGCTGGAGTTCTCACCCAGCCAGCAAGCCAGGAAACTTAATGAACCAGTTTTCTGTTCACGCTAGTCTATTTGCAGACACAGCCGAAACGCCAAAAAAATTTGAAGAATTAATGCTTAACTTTCTGTACCAGAACACAGTTAAAGCGTCTGACCATAGCAGCTGCAGACTGATTCCAGAGGGTTATATCCTCaaa                                                                                                               >  2:384857/1‑248 (MQ=255)
                                                                        aaaTACTTAGCTGTTATGTGCCTGTCTTTTCACCACTTCAGGCTCGGTGGTATGCTGGAGTTCTCACACAGCCAGCAAGCAAGGAAACTTAATGAACCAGTTTTATGTTCACGTTCGTCTATTTGAAGACACAGCCGAACAGCCAAAAAAATTTGAAGAATTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGCCGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGt                                                                                                            >  1:374299/1‑250 (MQ=255)
                                                                             cTTAGCTGTTATGTGCCTGTCTTTTCACCACTTCAGGCTCGGTGGTATGATGGAGTTCTCACACAGCCAGCAAGCAAGGAAACTTAATGAACCAGTTTTATGTTCACGTTCGTCTATTTGAAGACACAGCCGAACAGCCAAAAAAATTTGAAGAATTACTGCTTAACTTTCTGTACCAGAACACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCacaacaacaatg                                                                                                      >  2:198927/1‑240 (MQ=255)
                                                                                       aTGTGCCTGTCTTTTCACCACTTCAGGCTCGGTGGTATGCTGGAGTTCTCACACAGCCAGCAAGCAAGGAAACTTAATGAACCAGTTTTATGTTCACGTTCGTCTATTTGAAGACACAGCCGAACAGCCAAAAAAATTTGAAGAATTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCcaac                                                                                            <  1:222396/251‑1 (MQ=255)
                                                                                         gtgCCAGTCTTTTCACCACTTCAGGCTCGGTGGTATGATGGAGTTCTCACACAGCCAGCAAGCAAGGAAACTTAATGAACCAGTTTTATGTTCACGTTAGTCTATTTGAAGACACAGCCGAACAGCCAAAAAATTTTGAAGAATTAATGATTAACTTTCTGTAACAGAAAAAAGTAAAAGAGTCAGACGATATTTGCTTCAGACTGATTCCAGAGGAATATATCCTCAAAAGTACAATAAACTGCcacca                                                                                           >  2:381615/1‑247 (MQ=255)
                                                                                                    ttCACCACTTCAGGCTCGGTGGTATGCTGGAGTTCTCACACAGCCAGCAAGCAAGGAAACTTAATGAACCAGTTTTATGTTCACGTTCGTCTATTTGAAGACACAGCCGAACAGCCAAAAAAATTTGAAGAATTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGCTTCCAGAGGGATATATCCTCAAAAGAACAATGAACTGCCAACAAATCCTTGATc                                                                                >  2:102456/1‑250 (MQ=255)
                                                                                                     tCACCACTTCAGGCTCGGTGGTATGCTGGAGTTCTCACACAGCCAGCAAGCAAGGAAACTTAATGAACCAGTTTTATGTTCACGTTCGTCTATTTGAAGACACAGCCGAACAGCCAAAAAAATTTGAAGAATTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCa                                                                               >  1:311606/1‑250 (MQ=255)
                                                                                                               aGGCTCGGTTGTATGCTGGAGTTCTCACACAGCCAGCAAGCAAGGAAACTTAATGAACCAGTTTTATGTTCACGTTCGTCTATTTGAAGACACAGCCGAACAGCCAAAAAAATTTGAAGAATTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACATTTTCaa                                                                    <  2:30099/251‑1 (MQ=255)
                                                                                                                    cGGTGGTATGCTGGAGTTCTCACACAGCCAGCAAGCAAGGAAACTTAATGAACCAGTTTTATGTTCACGTTCGTCTATTTGAAGACACAGCCGAACAGCCAAAAAAATTTGAAGAATTAATGCTTAACTTGCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACATTTTCAATTGc                                                                <  2:109571/250‑1 (MQ=255)
                                                                                                                                        cacacaGCCAGCAAGCAAGGAAACTTAATGAACCAGTTTTATGTTCACGTTCGTCTATTTGAAGACACAGCCGAACAGCCAAAAAAATTTGAAGAATTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACATTTTCAATTGCTAACAGTGCCGGTGTTGAc                                             <  1:391137/249‑1 (MQ=255)
                                                                                                                                             agccagcAAGCAAGGAAACTTAATGAACCAGTTTTATGTTCACGTTCGTCTATTTGAAGACACAGCCGAACAGCCAAAAAAATTTGAAGAATTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACATTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAat                                       >  1:158601/1‑250 (MQ=255)
                                                                                                                                               ccagcAAGCAAGGAAACTTAATGAACCAGTTTTATGTTCACGTTCGTCTATTTGAAGACACAGCCGAACAGCCAAAAAAATTTGAAGAATTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACATTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAatata                                    <  1:134308/251‑1 (MQ=255)
                                                                                                                                                        aagGAAACTTAATGAACCAGTTTTATGTTCACGTTCGTCTATTTGAAGACACAGCCGAACAGCCAAAAAAATTTGAAGAATTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACATTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTATGt                           >  1:148786/1‑251 (MQ=255)
                                                                                                                                                                             tttATGTTCACGTTCGTCTATTTGAAGACACAGCCGAACAGCCAAAAAAATTTGAAGAATTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACATTTTCAATTGCTAACAGTGCCGGTGTTGAc                                             <  1:445237/212‑1 (MQ=255)
                                                                                                                                                                             tttATGTTCACGTTCGTCTATTTGAAGACACAGCCGAACAGCCAAAAAAATTTGAAGAATTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACATTTTCAATTGCTAACAGTGCCGGTGTTGAc                                             >  2:445237/1‑212 (MQ=255)
                                                                                                                                                                                 tGTTCACGTTCGTCTATTTGAAGACACAGCCGAACAGCCAAAAAAATTTGAAGAATTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACATTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAACAAAGCGCATGCTTAc  >  2:396121/1‑251 (MQ=255)
                                                                                                                                                                                                                      |                                                                                                                                                                                                                      
TTGTCAGCCTCAACCTTAACCCGCGTGCGCTGAGCCGCATTCACCACAACGCTAAGGATTCTCTCTGGTTGAAAATACTTAGCTGTTATGTGCCTGTCTTTTCACCACTTCAGGCTCGGTGGTATGCTGGAGTTCTCACACAGCCAGCAAGCAAGGAAACTTAATGAACCAGTTTTATGTTCACGTTCGTCTATTTGAAGACACAGCCGAACAGACCAAAAAATTTGAAGAATTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACATTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAACAAAGCGCATGCTTAC  >  NC_000913/1201334‑1201761

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