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

AGGCGGGTGGGTATCCACAGAAGGTCGATTGTCAGCCTCAACCTTAACCCGCGTGCGCTGAGCCGCATTCACCACAACGCTAAGGATTCTCTCTGGTTGAAAATACTTAGCTGTTATGTGCCTGTCTTTTCACCACTTCAGGCTCGGTGGTATGCTGGAGTTCTCACACAGCCAGCAAGCAAGGAAACTTAATGAACCAGTTTTATGTTCACGTTCGTCTATTTGAAGACACAGCCGAACAGACCAAAAAATTTGAAGAATTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACATTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAACAAAGCGCATGCTTACTTCCGTCTGCTTCCTTAGTTGGTAACGATT  >  NC_000913/1201306‑1201791
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    
aGGCGGGTGGGTATCCACAGAAGGTCGATTGTCAGCCTCAACCTTAACCCGCGTGCGCTGAGCCGCATTCACCACAACGCTAAGGATTCTCTCTGGTTGAAAATACTTAGCTGTTATGTGCCTGTCTTTTCACCACTTCAGGCTCGGTGGTATGCTGGAGTTCTCACACAGCCAGCAAGCAAGGAAACTTAATGAACCAGTTTTATGTTCACGTTCGTCTATTTGAAGACACAGCCGAACAGCCaaaaaa                                                                                                                                                                                                                                             >  2:419674/1‑250 (MQ=255)
       ggggtATCCACAGAAGGTCGATTGTCAGCCTCAACCTTAACCCGCGTGCGCTGAGCCGCATTCACCACAACGCTAAGGATTCTCTCTGGTTGAAAATACTTAGCTGTTATGTGCCTGTCTTTTCACCACTTCAGGCTCGGTGGTATGCTGGAGTTCTCACACAGCCAGCAAGCAAGGAAACTTAATGAACCAGTTTTATGTTCACGTTCGTCTATTTGAAGACACAGCCGAACAGCCAAAAAATTTgaaga                                                                                                                                                                                                                                     <  2:218808/250‑1 (MQ=255)
         ggtATCCACAGAAGGTCGATTGTCAGCCTCAACCTTAACCCGCGTGCGCTGAGCCGCATTCACCACAACGCTAAGGATTCTCTCTGGTTGAAAATACTTAGCTGTTATGTGCCTGTCTTTTCACCACTTCAGGCTCGGTGGTATGCTGGAGTTCTCACACAGCCAGCAAGCAAGGAAACTTAATGAACCAGTTTTATGTTCACGTTCGTCTATTTGAAGACACAGCCGAACAGCCAAAAAAATTTgaagaa                                                                                                                                                                                                                                    <  2:160465/251‑1 (MQ=255)
                 caGAAGGTCGATTGTCAGCCTCCACCTTAACCCGAGTGCGCTGAGCCGCATTCACCACACCGCTAAGGATTCTCTCTGTTTGAAAATACTTAGCTGTTATGTGCCTGTCTTTTCACCACTTCAGGCTCGGTGGTATGCTGGAGTTCTCACACAGCCAGCAAGCAAGGAAACTTAATGAACCCGTTTTATGTTCACGTTCGTCTATTTGAAGCCACCGCCGAACAGCCAAAACAATTTGAAGAATTACTGCt                                                                                                                                                                                                                            >  2:564582/1‑251 (MQ=255)
                              gTCAGCCTCAACCTTAACCCGCGTGCGCTGAGCCGCATTCACCACAACGCTAAGGATTCTCTCTGGTTGAAAATACTTAGCTGTTATGTGCCTGTCTTTTCACCACTTCAGGCTCGGTGGTATGCTGGAGTTCTCACACAGCCAGCAAGCAAGGAAACTTAATGAACCAGTTTTATGTTCACGTTCGTCTATTTGAAGACACAGCCGAACAGCCAAAAAAATTTGAAGAATTAATGCTTAACTTTCTGTAc                                                                                                                                                                                                               <  2:155394/251‑1 (MQ=255)
                                    cTCAACCTTAACCCGCGTGCGCTGAGCCGCATTCACCACAACGCTAAGGATTCTCTCTGGTTGAAAATACTTAGCTGTTATGTGCCTGTCTTTTCACCACTTCAGGCTCGGTGGTATGCTGGAGTTCTCACACAGCCAGCAAGCAAGGAAACTTAATGAACCAGTTTTATGTTCACGTTCGTCTATTTGAAGACACAGCCGAACAGCCAAAAAAATTTGAAGAATTAATGCTTAACTTTCTGTACCAGaa                                                                                                                                                                                                          >  2:183597/1‑250 (MQ=255)
                                              aCCCGCGTGCGCTGAGCCGCATTCACCACAACGCTAAGGATTCTCTCTGGTTGAAAATACTTAGCTGTTATGTGCCTGTCTTTTCACCACTTCAGGCTCGGTGGTATGCTGGAGTTCTCACACAGCCAGCAAGCAAGGAAACTTAATGAACCAGTTTTATGTTCACGTTCGTCTATTTGAAGACACAGCCGAACAGCCAAAAAAATTTGAAGAATTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAag                                                                                                                                                                                               <  1:574413/251‑1 (MQ=255)
                                                ccGCGTGCGCTGAGCCGCATTCACCACAACGCTAAGGATTCTCTCTGGTTGAAAATACTTAGCTGTTATGTGCCTGTCTTTTCACCACTTCAGGCTCGGTGGTATGCTGGAGTTCTCACACAGCCAGCAAGCAAGGAAACTTAATGAACCAGTTTTATGTTCACGTTCGTCTATTTGAAGACACAGCCGAACAGCCAAAAAAATTTGAAGAATTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAagag                                                                                                                                                                                             <  1:352867/251‑1 (MQ=255)
                                                      gcgcTGAGCCGCATTCATCACAACGCTAAGGATTCTCTCTGGTTGAAAATACTTAGCTGTTATGTGCCTGTCTTTTCACCACTTCAGGCTCGGTGGTATGCTGGAGTTCTCACACAGCCAGCAAGCAAGGAAACTTAATGAACCAGTTTTATGTTCACGTTCGTCTATTTGAAGACACAGCCGAACAGCCAAAAAAATTTGAAGAATTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGAc                                                                                                                                                                                       <  1:515861/251‑1 (MQ=255)
                                                       cgcTGAGCCGCATTCACCACAACGCTAAGGATTCTCTCTGGTTGAAAATACTTAGCTGTTATGTGCCTGTCTTTTCACCACTTCAGGCTCGGTGGTATGCTGGAGTTCTCACACAGCCAGCAAGCAAGGAAACTTAATGAACCAGTTTTATGTTCACGTTCGTCTATTTGAAGACACAGCCGAACAGCCAAAAAAATTTGAAGAATTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACg                                                                                                                                                                                      <  1:261649/251‑1 (MQ=255)
                                                                             cgctacgggttaactctGGTTCAACATACATAGCTGTTATGTGCCTGTCTTTTCACCACTTCAGGCTCGGTGGTATGCTGGAGTTCTAACACAGCCAGCAAGCAAGGACACTTAATGAACCAGTTTTATGTTCACGTTCGTCTATTTGAAGACACAGCCGAACAGCCAAAAAAATTTGAAGAATTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTcc                                                                                                                                                                >  2:36145/14‑251 (MQ=255)
                                                                                  aGGATTCTCTCTGGTTGAAAATACTTAGCTGTTATGTGCCTGTCTTTTCACCACTTCAGGCTCGGTGGTATGCTGGAGTTCTCACACAGCCAGCAAGCAAGGAAACTTAATGAACCAGTTTTATGTTCACGTTCGTCTATTTGAAGACACAGCCGAACAGCCAAAAAAATTTGAAGAATTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAgg                                                                                                                                                           <  1:572829/251‑1 (MQ=255)
                                                                                             tGGTTGAAAATACTTAGCTGTTATGTGCCTGTCTTTTCACCACTTCAGGCTCGGTGGTATGCTGGAGTTCTCACACAGCCAGCAAGCAAGGAAACTTAATGAACCAGTTTTATGTTCACGTTCGTCTATTTGAAGACACAGCCGAACAGCCAAAAAAATTTGAAGAATTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGCCTGATTCCAGAGGGATATATCCTc                                                                                                                                                >  2:572933/1‑251 (MQ=255)
                                                                                              ggTTGAAAATACTTAGCTGTTATGTGCCTGTCTTTTCACCACTTCAGGCTCGGTGGTATGCTGGAGTGCTCACACAGCCAGCAAGCAAGGAAACTTAATGAACCAGTTTTATGTTCACGTTCGTCTATTTGAAGACACAGCCGAACAGCCAAAAAAATTTGAAGAATTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGGCGGATTCCAGAGGGAGATATCCt                                                                                                                                                 <  2:611719/249‑1 (MQ=255)
                                                                                                        aCTTAGCTGGTATGTGCCTGTCTTTTCACCTCTTCAGGCTCGGTGGTATGCTGGAGTTCTCACACAGCCAGCAAGCAAGGAAACTTAATGAACCAGTTCTATGGTCACGTGCGTCTATTTGAAGACACAGCCGAACAGCCAAAAAAATTTGAAGAATGAATTCGTAACTTGCTGTGCCGGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAAt                                                                                                                                     <  1:371816/251‑1 (MQ=255)
                                                                                                         cTTAGCTGTTATGTGCCTGTCTTTTCACCACTTCAGGCTCGGTGGAATGCTGGAGTTCTCACACAGCCAGCAAGCAAGGAAACTTAATGAACCAGTTTTATGTTCACGTTCGTCTATTTGAAGACACAGCCGAACAGCAAAAAAAATTTGAAGAATTAATGCTTACCTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATg                                                                                                                                    >  1:26817/1‑251 (MQ=255)
                                                                                                                gTTATGTGCCTGTCTTTTCACCACTTCAGGCTCGGTGGTATGCTGGAGTTCTCACACAGCCAGCAAGCAAGGAAACTTAATGAACCAGTTTTATGTTCACGTTCGTCTATTTGAAGACACAGCCGAACAGCCAAAAAAATTTGAAGAATTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGcc                                                                                                                             >  1:495743/1‑251 (MQ=255)
                                                                                                                       gCCTGTCTTTTCACCACTTAAGGCTCGGTGGTATGCTGGAGTTCTCACACAGCCAGCAAGCAAGGAAACTTAATGAACCAGTTTTATGTTCACGTTCGTCTATTTGAAGACACAGCCGAACAGCCAAAAAAATTTGAAGAATTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAAt                                                                                                                      >  2:576043/1‑251 (MQ=255)
                                                                                                                                                cGGTGGTATGCTGGAGTTCTCACACAGCCAGCAAGCAAGTAAACTTTATGAACCAGTTTTATGTTCACGTTCGTCTATTTGAAGACACAGCCGAACAGCCAAAAAAATTTGAAGAATTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACATTTTCAATTg                                                                                               <  2:220221/249‑1 (MQ=255)
                                                                                                                                                cGGTGGTATGCTGGAGTTCTCACACAGCCAGCAAGCAAGGAAACTTAATGAACCAGTTTTATGTTCACGTTCGTCTATTTGAAGACACAGCCGAACAGCCAAAAAAATTTGAAGAATTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAAGGAACTGCCAACAAATCCTTGATCAAACATTTTCAATTg                                                                                               <  2:305404/249‑1 (MQ=255)
                                                                                                                                                           tGGAGTTCTCACACAGCCAGCAAGCAAGGAAACTTAATGAACCAGTTTTATGTTCACGTTCGTCTATTTGAAGACACAGCCGAACAGCCAAAAAAATTTGAAGAATTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACATTTTCAATTGCTAACAGTGCCgg                                                                                  >  2:374718/1‑251 (MQ=255)
                                                                                                                                                            ggAGTTCTCACACAGCCAGCAAGCAAGGAAACTTAATGAACCAGTTTTATGTTCACGTTCGTCTATTTGAAGACACAGCCGAACAGCCAAAAAAATTTGAAGAATTTATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACATTTTCAATTGCTAACAGTGCCGgt                                                                                 <  1:216125/251‑1 (MQ=255)
                                                                                                                                                                           ccagcAAGCAAGGAAACTTAATGAACCAGTTTTATGTTCACGTTCGTCTATTTGAAGACACAGCCGAACAGCCAAAAAAATTTGAAGAATTAAGGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACATTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAatata                                                                  <  2:432618/251‑1 (MQ=255)
                                                                                                                                                                                        aaaCTTAATGAACCAGTTTTATGTTCACGTTCGTCTATTTGAAGACACAGCCGAACAGCCAAAAAAATTTGAAGAATTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAGCAAATCCTTGATCAATCATTTTGAATTGCTAACAGTGCCGGTGGTGACGCAAATATATGTGTCTGTAAAt                                                     <  2:405044/251‑1 (MQ=255)
                                                                                                                                                                                                 gAACCAGTTTTATGTTCACGTTAGTCTATTTGACGACACAGCCGAACAGCCAAAAAAATTTGAAGAATTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACATTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAACaa                                             >  1:178561/1‑250 (MQ=255)
                                                                                                                                                                                                            aTGTGCACGTTCGTCTATTTGAAGACACAGCCGAACAGCCAAAAAAATGTGAAGAATTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGGACAATGAACTGCCAACAAATCCTTGATCAAACATTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAACAAAGCGCATGCtt                                  <  2:466395/250‑1 (MQ=255)
                                                                                                                                                                                                                 cACGTTCGTCTATTTGAAGACACAGCCGAACAGCCAAAAAAATTTGAAGAATTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACATTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAACAAAGCGCATGCTTACTTcc                            <  1:442203/251‑1 (MQ=255)
                                                                                                                                                                                                                                          ccGAACAGCCAAAAAAATTTGAAGAATTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACATTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAACAAAGCGCATGCTTACTTCCGTCTGCTTCCTTAGTTGGTAACGAt   >  1:411013/1‑251 (MQ=255)
                                                                                                                                                                                                                                           cGAACAGCCAAAAAAATTTGAAGAATTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACATTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAACAAAGCGCATGCTTACTTCCGTCTGCTTCCTTAGTTGGTAACGAtt  <  1:419674/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    
AGGCGGGTGGGTATCCACAGAAGGTCGATTGTCAGCCTCAACCTTAACCCGCGTGCGCTGAGCCGCATTCACCACAACGCTAAGGATTCTCTCTGGTTGAAAATACTTAGCTGTTATGTGCCTGTCTTTTCACCACTTCAGGCTCGGTGGTATGCTGGAGTTCTCACACAGCCAGCAAGCAAGGAAACTTAATGAACCAGTTTTATGTTCACGTTCGTCTATTTGAAGACACAGCCGAACAGACCAAAAAATTTGAAGAATTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACATTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAACAAAGCGCATGCTTACTTCCGTCTGCTTCCTTAGTTGGTAACGATT  >  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: