Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,463,057 G→T 10.5% E16* (GAA→TAA)  yfdF → uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,463,0570GT10.5% 126.7 / 5.9 38E16* (GAA→TAA) yfdFuncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (20/14);  new base T (2/2);  total (22/16)
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.73e-01

CTCAAATTAAAGATAGGTTCTAAATAAATGAGCGTTTTTTGATAGTCTATTTCATTAGGTAATATATATTTGTAACAAATCAATCAAAATGGAATAAAATCATGCTACCATCTATTTCAATCAACAATACCAGCGCAGCTTACCCAGAATCCATCAATGAAAATAACAATGATGAAGTTAATGGATTAGTACAAGAGTTCAAAAACCTTTTTAATGGTAAGGAAGGAATAAGCACCTGTATTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCC  >  NC_000913/2462911‑2463202
                                                                                                                                                  |                                                                                                                                                 
cTCAAATTAAAGATAGGTTCTAAATAAATGAGCGTTTTTTGATAGTCTATTTCATTAGGTAATATATATTTGTAACAAATCAATCAAAATGGAATAAAATCATGCTACCATCTATTTCAATCAACAATACCAGCGCAGCTTACCCAGaa                                                                                                                                                 <  2:471398/149‑1 (MQ=255)
 tCAAATTAAAGATAGGTTCTAAATAAATGAGCGTTTTTTGATAGTCTATTTCATTAGGTAATATATATTTGTAACAAATCAATCAAAATGGAATAAAATCATGCTACCATCTATTTCAATCAACAATACCAGCGCAGCTTACCCAGAAt                                                                                                                                                >  3:91611/1‑149 (MQ=255)
                 ttCTAAATAAATGAGCGTTTTTTGATAGTCTATTTCATTAGGTAATATATATTTGTAACAAATCAATCAAAATGGAATAAAATCATGCTACCATCTATTTCAATCAACAATACCAGCGCAGCTTACCCAGAATCCATCAATGAAAATaa                                                                                                                                >  3:173809/1‑149 (MQ=255)
                        taaatGAGCGTTTTTTGATAGTCTATTTCATTAGGTAATATATATTTGTAACAAATCAATCAAAATGGAATAAAATCATGCTACCATCTATTTCAATCAACAATACCAGCGCAGCTTACCCAGAATCCATCAATGAAAATAACAatgat                                                                                                                         >  4:116350/1‑149 (MQ=255)
                           atGAGCGTTTTTTGATAGTCTATTTCATTAGGTAATATATATTTGTAACAAATCAATCAAAATGGAATAAAATCATGCTACCATCTATTTCAATCAACAATACCAGCGCAGCTTACCCAGAATCCATCAATGAAAATAACAATGATGaa                                                                                                                      <  4:173809/149‑1 (MQ=255)
                                 gTTTTTTGATAGTCTATTTCATTAGGTAATATATATTTGTAACAAATCAATCAAAATGGAATAAAATCATGCTACCATCTATTTCAATCAACAATACCAGCGCAGCTTACCCAGAATCCATCAATGAAAATAACAATGATGAAGTTAAt                                                                                                                <  3:75645/149‑1 (MQ=255)
                                           aGTCTATTTCATTAGGTAATATATATTTGTAACAAATCAATCAAAATGGAATAAAATCATGCTACCATCTATTTCAATCAACAATACCAGCGCAGCTTACCCAGAATCCATCAATGAAAATAACAATGATGAAGTTAATGGATTAGTAc                                                                                                      <  2:190928/149‑1 (MQ=255)
                                             tCTATTTCATTAGGTAATATATATTTGTAACAAATCAATCAAAATGGAATAAAATCATGCTACCATCTATTTCAATCAACAATACCAGCGCAGCTTACCCAGAATCCATCAATGAAAATAACAATGATGaa                                                                                                                      <  2:25739/131‑1 (MQ=255)
                                             tCTATTTCATTAGGTAATATATATTTGTAACAAATCAATCAAAATGGAATAAAATCATGCTACCATCTATTTCAATCAACAATACCAGCGCAGCTTACCCAGAATCCATCAATGAAAATAACAATGATGaa                                                                                                                      >  1:25739/1‑131 (MQ=255)
                                             tCTATTTCATTAGGTAATATATATTTGTAACAAATCAATCAAAATGGAATAAAATCATGCTACCATCTATTTCAATCAACAATACCAGCGCAGCTTACCCAGAATCCATCAATGAAAATAACAATGATGAAGTTAATGGATTAGTACaa                                                                                                    >  3:510386/1‑149 (MQ=255)
                                               tATTTCATTAGGTAATATATATTTGTAACAAATCAATCAAAATGGAATAAAATCATGCTACCATCTATTTCAATCAACAATACCAGCGCAGCTTACCCAGAATCCATCAATGAAAATAACAATGATGAAGTTAATGGATTAGTACAaga                                                                                                  >  2:110061/1‑149 (MQ=255)
                                                                tatatTTGTAACAAATCAATCAAAATGGAATAAAATCATGCTACCATCTATTTCAATCAACAATACCAGCGCAGCTTACCCAGAATCCATCAATGAAAATAACAATGATGAAGTTAATGGATTAGTACAAGAGTTCAAAAACCTTTTTa                                                                                 >  3:512400/1‑149 (MQ=255)
                                                                       gTAACAAATCAATCAAAATGGAATAAAATCATGCTACCATCTATTTCAATCAACAATACCAGCGCAGCTTACCCAGAATCCATCAATGAAAATAACAATGATGAAGTTAATGGATTAGTACAAGAGTTCAAAAACCTTTTTAATGGTaa                                                                          <  4:512400/149‑1 (MQ=255)
                                                                             aatcaatcAAAATGGAATAAAATCATGCTACCATCTATTTCAATCAACAATACCAGCGCAGCTTACCCAGAATCCATCAATGAAAATAACAATGATGAAGTTAATGGATTAGTACAAGAGTTCAAAAACCTTTTTAATGGTaaggaagg                                                                    >  4:111609/1‑149 (MQ=255)
                                                                                      aaaTGGAATAAAATCATGCTACCATCTATTTCAATCAACAATACCAGCGCAGCTTACCCAGAATCCATCAATGAAAATAACAATGATGAAGTTAATGGATTAGTACAAGAGTTCAAAAACCTTTTTAATGGTAAGGAAGGAATAAGCAc                                                           >  4:83057/1‑149 (MQ=255)
                                                                                          ggAATAAAATCATGCTACCATCTATTTCAATCAACAATACCAGCGCAGCTTACCCAGAATCCATCAATGAAAATAACAATGATGAAGTTAATGGATTAGTACAAGAGTTCAAAAACCTTTTTAATGGTAAGGAAGGAATAAGCACCTGt                                                       <  3:83057/149‑1 (MQ=255)
                                                                                             aTAAAATCATGCTACCATCTATTTCAATCAACAATACCAGCGCAGCTTACCCAGAATCCATCAATGAAAATAACAATGATGAAGTTAATGGATTAGTACAAGAGTTCAAAAACCTTTTTAATGGTAAGGAAGGAATAAGCACCTGTAtt                                                    >  4:251706/1‑149 (MQ=255)
                                                                                                  aTCATGCTACCATCTATTTCAATCAACAATACCAGCGCAGCTTACCCAGAATCCATCAATGAAAATAACAATGATGAAGTTAATGGATTAGTACAAGAGTTCAAAAACCTTTTTAATGGTAAGGAAGGAATAAGCACCTGTATTAAACa                                               >  4:103473/1‑149 (MQ=255)
                                                                                                          aCCATCTATTTCAATCAACAATACCAGCGCAGCTTACCCAGAATCCATCAATGAAAATAACAATGATGAAGTTAATGGATTAGTACAAGAGTTCAAAAACCTTTTTAATGGTAaggaagga                                                                   <  2:53038/121‑1 (MQ=255)
                                                                                                          aCCATCTATTTCAATCAACAATACCAGCGCAGCTTACCCAGAATCCATCAATGAAAATAACAATGATGAAGTTAATGGATTAGTACAAGAGTTCAAAAACCTTTTTAATGGTAaggaagga                                                                   >  1:53038/1‑121 (MQ=255)
                                                                                                             aTCTATTTCAATCAACAATACCAGCGCAGCTTACCCAGAATCCATCAATGAAAATAACAATGATGAAGTTAATGGATTAGTACAAGAGTTCAAAAACCTTTTTAATGGTAAGGAAGGAATAAGCACCTGTATTAAACATCTACTTGAGc                                    >  2:591086/1‑149 (MQ=255)
                                                                                                             aTCTATTTCAATCAACAATACCAGCGCAGCTTACCCAGAATCCATCAATGAAAATAACAATGATGAAGTTAATGGATTAGTACAAGAGTTCAAAAACCTTTTTAATGGTAAGGAAGGAATAAGCACCTGTATTAAACATCTACTTGAGc                                    <  3:116350/149‑1 (MQ=255)
                                                                                                                  tttCAATCAACAATACCAGCGCAGCTTACCCAGAATCCATCAATGAAAATAACAATGATGAAGTTAATGGATTAGTACAAGAGTTCAAAAACCttt                                                                                    <  4:568781/96‑1 (MQ=255)
                                                                                                                  tttCAATCAACAATACCAGCGCAGCTTACCCAGAATCCATCAATGAAAATAACAATGATGAAGTTAATGGATTAGTACAAGAGTTCAAAAACCttt                                                                                    >  3:568781/1‑96 (MQ=255)
                                                                                                                  tttCAATCAACAATACCAGCGCAGCTTACCCAGAATCCATCAATGAAAATAACAATGATGAAGTTAATGGATTAGTACAAGAGTTCAAAAACCTTTTTAATGGTAAGGAAGGAATAAGCACCTGTATTAAACATCTACTTGAGCTtata                               >  4:69346/1‑149 (MQ=255)
                                                                                                                      aatcaaCAATACCAGCGCAGCTTACCCAGAATCCATCAATGAAAATAACAATGATGAAGTTAATGGATTAGTACAAGAGTTCAAAAACCttt                                                                                    >  1:303604/1‑92 (MQ=255)
                                                                                                                      aatcaaCAATACCAGCGCAGCTTACCCAGAATCCATCAATGAAAATAACAATGATGAAGTTAATGGATTAGTACAAGAGTTCAAAAACCttt                                                                                    <  2:303604/92‑1 (MQ=255)
                                                                                                                      aatcaaCAATACCAGCGCAGCTTACCCAGAATCCATCAATGAAAATAACAATGATGAAGTTAATGGATTAGTACAAGAGTTCAAAAACCTTTTTAATGGTAAGGAAGGAATAAGCACCTGTATTaaa                                                 >  2:94032/1‑127 (MQ=255)
                                                                                                                      aatcaaCAATACCAGCGCAGCTTACCCAGAATCCATCAATGAAAATAACAATGATGAAGTTAATGGATTAGTACAAGAGTTCAAAAACCTTTTTAATGGTAAGGAAGGAATAAGCACCTGTATTaaa                                                 <  1:94032/127‑1 (MQ=255)
                                                                                                                      aatcaaCAATACCAGCGCAGCTTACCCAGAATCCATCAATGAAAATAACAATGATGAAGTTAATGGATTAGTACAAGAGTTCAAAAACCTTTTTAATGGTAAGGAAGGAATAAGCACCTGTATTAAACATCTACTTGAGCTTATaaaaa                           <  4:91611/149‑1 (MQ=255)
                                                                                                                        tcaaCAATACCAGCGCAGCTTACCCATAATCCATCAATGAAAATAACAATGATGAAGTTAATGGATTAGTACAAGAGTTCAAAAACCTTTTTAATGGTAAGGAAGGAATAAGCACCTGTATTAAACATCTACTTGAGc                                    <  2:4221/138‑1 (MQ=255)
                                                                                                                        tcaaCAATACCAGCGCAGCTTACCCATAATCCATCAATGAAAATAACAATGATGAAGTTAATGGATTAGTACAAGAGTTCAAAAACCTTTTTAATGGTAAGGAAGGAATAAGCACCTGTATTAAACATCTACTTGAGc                                    <  2:330278/138‑1 (MQ=255)
                                                                                                                        tcaaCAATACCAGCGCAGCTTACCCATAATCCATCAATGAAAATAACAATGATGAAGTTAATGGATTAGTACAAGAGTTCAAAAACCTTTTTAATGGTAAGGAAGGAATAAGCACCTGTATTAAACATCTACTTGAGc                                    >  1:4221/1‑138 (MQ=255)
                                                                                                                        tcaaCAATACCAGCGCAGCTTACCCATAATCCATCAATGAAAATAACAATGATGAAGTTAATGGATTAGTACAAGAGTTCAAAAACCTTTTTAATGGTAAGGAAGGAATAAGCACCTGTATTAAACATCTACTTGAGc                                    >  1:330278/1‑138 (MQ=255)
                                                                                                                                   aGCGCAGCTTACCCAGAATCCATCAATGAAAATAACAATGATGAAGTTAATGGATTAGTACAAGAGTTCAAAAACCTTTTTAATGGTAAGGAAGGAATAAGCACCTGTATTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGt              <  4:510386/149‑1 (MQ=255)
                                                                                                                                               ccAGAATCCATCAATGAAAATAACAATGATGAAGTTAATGGATTAGTACAAGAGTTCAAAAACCTTTTTAATGGTAAGGAAGGAATAAGCACCTGTATTAAACATCTACTTGAGCTTATAAACAACGCCATACGAGTAAACGACGATcc  >  1:379904/1‑149 (MQ=255)
                                                                                                                                               ccAGAATCCATCAATGAAAATAACAATGATGAAGTTAATGGATTAGTACAAGAGTTCAAAAACCTTTTTAATGGTAAGGAAGGAATAAGCACCTGTATTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATcc  >  2:60706/1‑149 (MQ=255)
                                                                                                                                               ccAGAATCCATCAATGAAAATAACAATGATGAAGTTAATGGATTAGTACAAGAGTTCAAAAACCTTTTTAATGGTAAGGAAGGAATAAGCACCTGTATTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATcc  >  4:513657/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
CTCAAATTAAAGATAGGTTCTAAATAAATGAGCGTTTTTTGATAGTCTATTTCATTAGGTAATATATATTTGTAACAAATCAATCAAAATGGAATAAAATCATGCTACCATCTATTTCAATCAACAATACCAGCGCAGCTTACCCAGAATCCATCAATGAAAATAACAATGATGAAGTTAATGGATTAGTACAAGAGTTCAAAAACCTTTTTAATGGTAAGGAAGGAATAAGCACCTGTATTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCC  >  NC_000913/2462911‑2463202

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: