Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,776,100 G→A E172K (GAA→AAA)  ydiF → putative acetate‑CoA transferase YdiF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,776,1000GA96.8% 90.3 / NA 31E172K (GAA→AAA) ydiFputative acetate‑CoA transferase YdiF
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (13/17);  minor base C (1/0);  total (14/17)
Fisher's exact test for biased strand distribution p-value = 4.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.26e-01

ATTATGTGGTCACTGGGGACAATCGCCGCGTATTTCTGAACTCGCAGAACAAAATAAAATTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATATTGTGGTGGT  >  NC_000913/1775859‑1776341
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                  
catATGTGGTCACTGGGGACAATCGCCGCGTATTTCTGAACTCGCAGAACAAAATAAAAGTATTGCGTATAACTACCCAAAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAAAATATctct                                                                                                                                                                                                                                           <  2:176228/249‑1 (MQ=255)
   aTGTGGTCACGTGGGACAATCTCCGCGTATTTCTTAACTCGCAGAACAAAATAAAATTATTGCTTATAACTACCCACAAGGTTTACTTACACCAACCCTACGCTCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGTTCGAGTTTGATAACAAAAAATATCTCTAtt                                                                                                                                                                                                                                        <  2:454963/251‑1 (MQ=255)
    tgtgGTCACTGGGGACAATCGCCGCGTATTTCTGAACTCGCAGAACAAAATAAAATTATTGATTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGCATTATTAGTGATATTGGCATCGGGACATTTGTAGATCCACGCCAGCAATGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGCTTAAACTGGGCGAGTTTGATAACAAAAAATATCTCTAtt                                                                                                                                                                                                                                        >  2:383338/1‑250 (MQ=255)
              tgggACAATCGCCGCGTCTTTCTGACGTCGAAGAACAAATTAAAATTATTTATTATAACTACCCACAAGGTGTACTTACACAAACCTTACGAGCCGCCGCAGCCCACCAGCCTGGAATTATTAGAAATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAAATGGTCGAGATTGATAACAAACAAGATATATATTAAAAAGAGAtt                                                                                                                                                                                                                             >  2:78204/2‑251 (MQ=255)
              ggggACAATCGCCGCGTATTTCTGAACTCGCAGAACAAAATAAAATTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAAAATATCTCTATTACAAAGCGAtt                                                                                                                                                                                                                             <  1:46079/251‑1 (MQ=255)
                   cAATCGCCGCGTATTTCTGAACTCGCAGAACAAAATAAAATTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAAAATATCTCTATTACAAAGCGATTGCGcc                                                                                                                                                                                                                        >  1:192222/1‑251 (MQ=255)
                       cgccgtggaTTTCTGAACGCGCAGAACAAAATAAAATTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAAAATATCTCTATTACAAAGCGATTGCGCCAGat                                                                                                                                                                                                                    <  2:511761/243‑1 (MQ=255)
                                            cAGAACAAAATAAAATTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAAAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTacc                                                                                                                                                                                               <  1:130125/251‑1 (MQ=255)
                                            cAGAACAAAATAAAATTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAAAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTacc                                                                                                                                                                                               <  1:113673/251‑1 (MQ=255)
                                                            tattGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAAAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGa                                                                                                                                                                                 >  1:506415/1‑249 (MQ=255)
                                                                 cTTATAACTACCCACACGGGGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGCTCCACGCCAGCAAGGCGGCAAACTGCATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAAAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCTCGCTACCACCTTCGACAGTGAAGGcaac                                                                                                                                                                          >  2:276889/1‑248 (MQ=255)
                                                                      aaCTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAAAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCAc                                                                                                                                                                     >  1:83361/1‑251 (MQ=255)
                                                                        cTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAAAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTc                                                                                                                                                                                                           <  1:478012/211‑1 (MQ=255)
                                                                        cTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAAAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTc                                                                                                                                                                                                           >  2:478012/1‑211 (MQ=255)
                                                                                       aCTTACACAATCCTGACGCGCCGCCGCAACCCAACAGCCTGGTAGTATTAGTGATATTGTCATCGGGTCATTTGTCGATCCACGCCAGCAAGGCGGCTAACTGAATGAAGTCACTTAAGAATACCTGATTAAACTGGTCGATTTTGATAAAAAAAAATATCCCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCCCCTGCGACAGTGAAGGCTACTCCAcctattaatctgaggtg                                                                                                                                                     <  2:446628/250‑17 (MQ=255)
                                                                                          tACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGCGTTTGATAACAAAAAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGATACCACCTGCGACAGTGAGGGCTACGCCACTTTTGAAGATGAGGTGCTGt                                                                                                                                                 >  2:198605/1‑251 (MQ=255)
                                                                                             aacaacctttccggccgccgcAGACCCACAGCCTGGTATTATTAGTGATATTGGCATTGGTACATTTGTCGATCCAAGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAATACATGATTAAACTGGTCGAGTTTGATAACAAAAAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCCCTTTTGAAGATGAGGTGATGTAtc                                                                                                                                              <  2:127976/238‑1 (MQ=255)
                                                                                                       cgcgCAGCCGCAGCCCACCAGCCTGGTATTATTAGAGATATTGGCATCGGGACCTTTGTCGATCCACGCCAGAAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGGCGAGTTTGATAACAAAAAATATCTCTATTACAAAGAGATTGTGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGCCAGTGAAGGCTACGCCACTTTTGACGATGAGGTAATGTATCTAGAAGCAtt                                                                                                                                    >  2:363737/1‑251 (MQ=255)
                                                                                                              ccgcAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAAAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTAtt                                                                                                                             >  1:467978/1‑251 (MQ=255)
                                                                                                                 cAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAAAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGcc                                                                                                                          <  1:48470/251‑1 (MQ=255)
                                                                                                                      ccaccaGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAAAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAgg                                                                                                                      <  1:383590/251‑1 (MQ=255)
                                                                                                                                ggTATTATTAGTGATATTGGCATCGGGACATTTTGCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATGAAACTGGTCGAGTTTGATAACAAAAAATATCCCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGTTGATGTATCTCGACGCATTTGTTATTGCCCAGGCGGTGcac                                                                                                              <  2:192222/249‑1 (MQ=255)
                                                                                                                                           tGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAAAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGt                                                                                                  >  1:25200/1‑250 (MQ=255)
                                                                                                                                             atatTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAAAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGGGATGTATCTCGACGCATTGGTTATTGCCCCGGCGGAGCACAATAAcggcggggt                                                                                                >  2:162580/1‑247 (MQ=255)
                                                                                                                                                  tgCATCGGGACAGTTTTCTAGCCACTCCAGCAAGGCGGCAAACTTAATGAAGACACTAAAGAAGACCTGATTAAACTGGACGAGTTTGATAACAAAAAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAATGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCAATGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGtga                                                                                           <  2:85504/249‑1 (MQ=255)
                                                                                                                                                     aTCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAAAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGc                                                                                                                     >  1:414692/1‑221 (MQ=255)
                                                                                                                                                     aTCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAAATGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAAAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGc                                                                                                                     <  2:414692/221‑1 (MQ=255)
                                                                                                                                                                                                                 gaCCTGATTAAACTGGTCGAGTTTTATAACAAAAAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCgtg                          <  2:87695/251‑3 (MQ=255)
                                                                                                                                                                                                                    cTGATTAAACTGGTCGAGTTTGATAACAAAAAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTTGGATGATGCAGGGGGCGCAAAAGGGTTAGAAAGCCACCCCGCATCCTCAATATGTTCGCATTTCggggt                        >  1:279761/1‑248 (MQ=255)
                                                                                                                                                                                                                         tAAACTGGTCGACATGGAGAAAAAAAAATATCGCTACTAAAAAGCGAATGCTCCAGATAGTGCCTTCATACGCGATACCACATGCGATAGTGAAGGCTACGCCACTGATTAAGATGAGGTGATGTATCTCGACGGATTGGTTATTGGCCAGGCGGTGCACAAAAACGGCGGTATTGGGATAATGCAGTTGTAGAAAATGTGTAAGAAAGCCACGCTGCATCCTAAATCTGTGCGTATGCCGGGTTATCtg                   <  2:237663/250‑1 (MQ=255)
                                                                                                                                                                                                                                         gATAACAAAAAATATCTCTAGTACAATGCGATTGCGCCAGATAATGCCTTCATTCGCGCTACCGACTGCGACTGTGAAGGCTACGCAACTTTTGAAGATGAGGTGATGTACTTCGACGCATTGGTGATTGCCCAGGCTGTTCCCAAGATCGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGTTGCATCCTAATTCTGTCCGTATTCCGGGTTACCTGGTGGATCATgtggttga  <  2:200952/251‑4 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                  
ATTATGTGGTCACTGGGGACAATCGCCGCGTATTTCTGAACTCGCAGAACAAAATAAAATTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATATTGTGGTGGT  >  NC_000913/1775859‑1776341

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