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,5480AC89.7% 76.2 / 5.5 29coding (52/342 nt)ymfIuncharacterized protein YmfI
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base C (9/17);  minor base . (2/1);  total (11/18)
Fisher's exact test for biased strand distribution p-value = 5.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.82e-01

AGGCGGGTGGGTATCCACAGAAGGTCGATTGTCAGCCTCAACCTTAACCCGCGTGCGCTGAGCCGCATTCACCACAACGCTAAGGATTCTCTCTGGTTGAAAATACTTAGCTGTTATGTGCCTGTCTTTTCACCACTTCAGGCTCGGTGGTATGCTGGAGTTCTCACACAGCCAGCAAGCAAGGAAACTTAATGAACCAGTTTTATGTTCACGTTCGTCTATTTGAAGACACAGCCGAACAGACCAAAA‑AATTTGAAGAATTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACATTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAACAAAGCGCATGCTTACTTCCGTCTGCTTCCTTAGTTGGTAACGATT  >  NC_000913/1201306‑1201791
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    
aGGCGGGTGGGTATCCACAGAAGGTCGATTGTCAGCCTCAACCTTAACCCGCGTGCGCTGAGCCGCATTCACCACAACGCTAAGGATTCTCTCTGGTTGAAAATACTTAGCTGTTATGTGCCTGTCTTTTCACCACTTCAGGCTCGGTGGTATGCTGGAGTTCTCACACAGCCAGCAAGCAAGGAAACTTAATGAACCAGTTTTATGTTCACGTTCGTCTATTTGAAGACACAGCCGAACAG‑CCaaaa‑aa                                                                                                                                                                                                                                             >  2:419674/1‑250 (MQ=255)
       ggggtATCCACAGAAGGTCGATTGTCAGCCTCAACCTTAACCCGCGTGCGCTGAGCCGCATTCACCACAACGCTAAGGATTCTCTCTGGTTGAAAATACTTAGCTGTTATGTGCCTGTCTTTTCACCACTTCAGGCTCGGTGGTATGCTGGAGTTCTCACACAGCCAGCAAGCAAGGAAACTTAATGAACCAGTTTTATGTTCACGTTCGTCTATTTGAAGACACAGCCGAACAG‑CCAAAA‑AATTTgaaga                                                                                                                                                                                                                                     <  2:218808/250‑1 (MQ=255)
         ggtATCCACAGAAGGTCGATTGTCAGCCTCAACCTTAACCCGCGTGCGCTGAGCCGCATTCACCACAACGCTAAGGATTCTCTCTGGTTGAAAATACTTAGCTGTTATGTGCCTGTCTTTTCACCACTTCAGGCTCGGTGGTATGCTGGAGTTCTCACACAGCCAGCAAGCAAGGAAACTTAATGAACCAGTTTTATGTTCACGTTCGTCTATTTGAAGACACAGCCGAACAGCCAAAAA‑AATTTgaagaa                                                                                                                                                                                                                                    <  2:160465/251‑1 (MQ=255)
                 caGAAGGTCGATTGTCAGCCTCCACCTTAACCCGAGTGCGCTGAGCCGCATTCACCACACCGCTAAGGATTCTCTCTGTTTGAAAATACTTAGCTGTTATGTGCCTGTCTTTTCACCACTTCAGGCTCGGTGGTATGCTGGAGTTCTCACACAGCCAGCAAGCAAGGAAACTTAATGAACCCGTTTTATGTTCACGTTCGTCTATTTGAAGCCACCGCCGAACAG‑CCAAAACAATTTGAAGAATTACTGCt                                                                                                                                                                                                                            >  2:564582/1‑251 (MQ=255)
                              gTCAGCCTCAACCTTAACCCGCGTGCGCTGAGCCGCATTCACCACAACGCTAAGGATTCTCTCTGGTTGAAAATACTTAGCTGTTATGTGCCTGTCTTTTCACCACTTCAGGCTCGGTGGTATGCTGGAGTTCTCACACAGCCAGCAAGCAAGGAAACTTAATGAACCAGTTTTATGTTCACGTTCGTCTATTTGAAGACACAGCCGAACAGCCAAAAA‑AATTTGAAGAATTAATGCTTAACTTTCTGTAc                                                                                                                                                                                                               <  2:155394/251‑1 (MQ=255)
                                    cTCAACCTTAACCCGCGTGCGCTGAGCCGCATTCACCACAACGCTAAGGATTCTCTCTGGTTGAAAATACTTAGCTGTTATGTGCCTGTCTTTTCACCACTTCAGGCTCGGTGGTATGCTGGAGTTCTCACACAGCCAGCAAGCAAGGAAACTTAATGAACCAGTTTTATGTTCACGTTCGTCTATTTGAAGACACAGCCGAACAGCCAAAAA‑AATTTGAAGAATTAATGCTTAACTTTCTGTACCAGaa                                                                                                                                                                                                          >  2:183597/1‑250 (MQ=255)
                                              aCCCGCGTGCGCTGAGCCGCATTCACCACAACGCTAAGGATTCTCTCTGGTTGAAAATACTTAGCTGTTATGTGCCTGTCTTTTCACCACTTCAGGCTCGGTGGTATGCTGGAGTTCTCACACAGCCAGCAAGCAAGGAAACTTAATGAACCAGTTTTATGTTCACGTTCGTCTATTTGAAGACACAGCCGAACAGCCAAAAA‑AATTTGAAGAATTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAag                                                                                                                                                                                               <  1:574413/251‑1 (MQ=255)
                                                ccGCGTGCGCTGAGCCGCATTCACCACAACGCTAAGGATTCTCTCTGGTTGAAAATACTTAGCTGTTATGTGCCTGTCTTTTCACCACTTCAGGCTCGGTGGTATGCTGGAGTTCTCACACAGCCAGCAAGCAAGGAAACTTAATGAACCAGTTTTATGTTCACGTTCGTCTATTTGAAGACACAGCCGAACAGCCAAAAA‑AATTTGAAGAATTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAagag                                                                                                                                                                                             <  1:352867/251‑1 (MQ=255)
                                                      gcgcTGAGCCGCATTCATCACAACGCTAAGGATTCTCTCTGGTTGAAAATACTTAGCTGTTATGTGCCTGTCTTTTCACCACTTCAGGCTCGGTGGTATGCTGGAGTTCTCACACAGCCAGCAAGCAAGGAAACTTAATGAACCAGTTTTATGTTCACGTTCGTCTATTTGAAGACACAGCCGAACAGCCAAAAA‑AATTTGAAGAATTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGAc                                                                                                                                                                                       <  1:515861/251‑1 (MQ=255)
                                                       cgcTGAGCCGCATTCACCACAACGCTAAGGATTCTCTCTGGTTGAAAATACTTAGCTGTTATGTGCCTGTCTTTTCACCACTTCAGGCTCGGTGGTATGCTGGAGTTCTCACACAGCCAGCAAGCAAGGAAACTTAATGAACCAGTTTTATGTTCACGTTCGTCTATTTGAAGACACAGCCGAACAGCCAAAAA‑AATTTGAAGAATTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACg                                                                                                                                                                                      <  1:261649/251‑1 (MQ=255)
                                                                             cgctacgggttaactctGGTTCAACATACATAGCTGTTATGTGCCTGTCTTTTCACCACTTCAGGCTCGGTGGTATGCTGGAGTTCTAACACAGCCAGCAAGCAAGGACACTTAATGAACCAGTTTTATGTTCACGTTCGTCTATTTGAAGACACAGCCGAACAGCCAAAAA‑AATTTGAAGAATTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTcc                                                                                                                                                                >  2:36145/14‑251 (MQ=255)
                                                                                  aGGATTCTCTCTGGTTGAAAATACTTAGCTGTTATGTGCCTGTCTTTTCACCACTTCAGGCTCGGTGGTATGCTGGAGTTCTCACACAGCCAGCAAGCAAGGAAACTTAATGAACCAGTTTTATGTTCACGTTCGTCTATTTGAAGACACAGCCGAACAGCCAAAAA‑AATTTGAAGAATTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAgg                                                                                                                                                           <  1:572829/251‑1 (MQ=255)
                                                                                             tGGTTGAAAATACTTAGCTGTTATGTGCCTGTCTTTTCACCACTTCAGGCTCGGTGGTATGCTGGAGTTCTCACACAGCCAGCAAGCAAGGAAACTTAATGAACCAGTTTTATGTTCACGTTCGTCTATTTGAAGACACAGCCGAACAGCCAAAAA‑AATTTGAAGAATTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGCCTGATTCCAGAGGGATATATCCTc                                                                                                                                                >  2:572933/1‑251 (MQ=255)
                                                                                              ggTTGAAAATACTTAGCTGTTATGTGCCTGTCTTTTCACCACTTCAGGCTCGGTGGTATGCTGGAGTGCTCACACAGCCAGCAAGCAAGGAAACTTAATGAACCAGTTTTATGTTCACGTTCGTCTATTTGAAGACACAGCCGAACAGCCAAAAA‑AATTTGAAGAATTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGGCGGATTCCAGAGGGAGATATCCt                                                                                                                                                 <  2:611719/249‑1 (MQ=255)
                                                                                                        aCTTAGCTGGTATGTGCCTGTCTTTTCACCTCTTCAGGCTCGGTGGTATGCTGGAGTTCTCACACAGCCAGCAAGCAAGGAAACTTAATGAACCAGTTCTATGGTCACGTGCGTCTATTTGAAGACACAGCCGAACAGCCAAAAA‑AATTTGAAGAATGAATTCGTAACTTGCTGTGCCGGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAAt                                                                                                                                     <  1:371816/251‑1 (MQ=255)
                                                                                                         cTTAGCTGTTATGTGCCTGTCTTTTCACCACTTCAGGCTCGGTGGAATGCTGGAGTTCTCACACAGCCAGCAAGCAAGGAAACTTAATGAACCAGTTTTATGTTCACGTTCGTCTATTTGAAGACACAGCCGAACAGCAAAAAA‑AATTTGAAGAATTAATGCTTACCTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATg                                                                                                                                    >  1:26817/1‑251 (MQ=255)
                                                                                                                gTTATGTGCCTGTCTTTTCACCACTTCAGGCTCGGTGGTATGCTGGAGTTCTCACACAGCCAGCAAGCAAGGAAACTTAATGAACCAGTTTTATGTTCACGTTCGTCTATTTGAAGACACAGCCGAACAGCCAAAAA‑AATTTGAAGAATTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGcc                                                                                                                             >  1:495743/1‑251 (MQ=255)
                                                                                                                       gCCTGTCTTTTCACCACTTAAGGCTCGGTGGTATGCTGGAGTTCTCACACAGCCAGCAAGCAAGGAAACTTAATGAACCAGTTTTATGTTCACGTTCGTCTATTTGAAGACACAGCCGAACAGCCAAAAA‑AATTTGAAGAATTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAAt                                                                                                                      >  2:576043/1‑251 (MQ=255)
                                                                                                                                                cGGTGGTATGCTGGAGTTCTCACACAGCCAGCAAGCAAGTAAACTTTATGAACCAGTTTTATGTTCACGTTCGTCTATTTGAAGACACAGCCGAACAGCCAAAAA‑AATTTGAAGAATTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACATTTTCAATTg                                                                                               <  2:220221/249‑1 (MQ=255)
                                                                                                                                                cGGTGGTATGCTGGAGTTCTCACACAGCCAGCAAGCAAGGAAACTTAATGAACCAGTTTTATGTTCACGTTCGTCTATTTGAAGACACAGCCGAACAGCCAAAAA‑AATTTGAAGAATTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAAGGAACTGCCAACAAATCCTTGATCAAACATTTTCAATTg                                                                                               <  2:305404/249‑1 (MQ=255)
                                                                                                                                                           tGGAGTTCTCACACAGCCAGCAAGCAAGGAAACTTAATGAACCAGTTTTATGTTCACGTTCGTCTATTTGAAGACACAGCCGAACAGCCAAAAA‑AATTTGAAGAATTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACATTTTCAATTGCTAACAGTGCCgg                                                                                  >  2:374718/1‑251 (MQ=255)
                                                                                                                                                            ggAGTTCTCACACAGCCAGCAAGCAAGGAAACTTAATGAACCAGTTTTATGTTCACGTTCGTCTATTTGAAGACACAGCCGAACAGCCAAAAA‑AATTTGAAGAATTTATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACATTTTCAATTGCTAACAGTGCCGgt                                                                                 <  1:216125/251‑1 (MQ=255)
                                                                                                                                                                           ccagcAAGCAAGGAAACTTAATGAACCAGTTTTATGTTCACGTTCGTCTATTTGAAGACACAGCCGAACAGCCAAAAA‑AATTTGAAGAATTAAGGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACATTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAatata                                                                  <  2:432618/251‑1 (MQ=255)
                                                                                                                                                                                        aaaCTTAATGAACCAGTTTTATGTTCACGTTCGTCTATTTGAAGACACAGCCGAACAGCCAAAAA‑AATTTGAAGAATTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAGCAAATCCTTGATCAATCATTTTGAATTGCTAACAGTGCCGGTGGTGACGCAAATATATGTGTCTGTAAAt                                                     <  2:405044/251‑1 (MQ=255)
                                                                                                                                                                                                 gAACCAGTTTTATGTTCACGTTAGTCTATTTGACGACACAGCCGAACAGCCAAAAA‑AATTTGAAGAATTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACATTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAACaa                                             >  1:178561/1‑250 (MQ=255)
                                                                                                                                                                                                            aTGTGCACGTTCGTCTATTTGAAGACACAGCCGAACAGCCAAAAA‑AATGTGAAGAATTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGGACAATGAACTGCCAACAAATCCTTGATCAAACATTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAACAAAGCGCATGCtt                                  <  2:466395/250‑1 (MQ=255)
                                                                                                                                                                                                                 cACGTTCGTCTATTTGAAGACACAGCCGAACAGCCAAAAA‑AATTTGAAGAATTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACATTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAACAAAGCGCATGCTTACTTcc                            <  1:442203/251‑1 (MQ=255)
                                                                                                                                                                                                                                          ccGAACAGCCAAAAA‑AATTTGAAGAATTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACATTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAACAAAGCGCATGCTTACTTCCGTCTGCTTCCTTAGTTGGTAACGAt   >  1:411013/1‑251 (MQ=255)
                                                                                                                                                                                                                                           cGAACAGCCAAAAA‑AATTTGAAGAATTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACATTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAACAAAGCGCATGCTTACTTCCGTCTGCTTCCTTAGTTGGTAACGAtt  <  1:419674/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    
AGGCGGGTGGGTATCCACAGAAGGTCGATTGTCAGCCTCAACCTTAACCCGCGTGCGCTGAGCCGCATTCACCACAACGCTAAGGATTCTCTCTGGTTGAAAATACTTAGCTGTTATGTGCCTGTCTTTTCACCACTTCAGGCTCGGTGGTATGCTGGAGTTCTCACACAGCCAGCAAGCAAGGAAACTTAATGAACCAGTTTTATGTTCACGTTCGTCTATTTGAAGACACAGCCGAACAGACCAAAA‑AATTTGAAGAATTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACATTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAACAAAGCGCATGCTTACTTCCGTCTGCTTCCTTAGTTGGTAACGATT  >  NC_000913/1201306‑1201791

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