Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,637,431 A→G 14.7% K129E (AAA→GAA)  ydfO → Qin prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,637,4310AG14.7% 73.3 / 3.4 27K129E (AAA→GAA) ydfOQin prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (13/10);  new base G (2/2);  total (15/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TACATTTACTATTTAGCCACAGAGAATGTTCATATTGTATTAAAAAATGATAATACAGTGTTATTAAAGGGCCTAAAAAACATTGTGTCTGTCAAATTTTCAAAGGATAGGCATCTTATAGAAACGACCTCTAATAAGCTGAAATCCAGAGAGATCACATTTCAGGAATACAGAAGAAACCTTGCTAAAGCAGGAGTTTTTCGGTGGGTTACAAATATCCACGAACAAAAAAGATATTACTATACCTTTGATAATTCATTACTATTTACTGAAAGCATCCAGAAAACTACACAGATCTTACCACGCTAAACCATAACGTCCGGCTTCTCTCACTCCTGAGCCGGACTGCATTGGTTTAATAAAAACCATCAACAATTGTGATTTAGATATTCGGAACCATTCAAATATAACAAAACCCCGTAAAAACGAGGTTTATGGATAAATTTTATTATTGAATACATCAGATTAAATTAATCTTGACATCATAGCTTTCAAGACCCGTCATTTTTTCCCGTGCGGTAAACTGAATACTGGTAACTTCTTTCCCGGTCTTTTTCTTAAGTTCAA  >  NC_000913/1637149‑1637715
                                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                            
tacattTACTATTTAGCCACAGAGAATGTTCATATTGTATTAAAAAATGATAATACAGTGTTATTAAAGGGCCTAAAAAACATTGTGTCTGTCAAATTTTCAAAGGATAGGCATCTTATAGAAACGACCTCTAATAAGCTGAAATCCAGAGAGATCACATTTCAGGAATACAGAAGAAACCTTGCTAAAGCAGGAGTTTTTCGGTGGGTTACAAATATCCACGAACAAAAAAGATATTACTATACCTTTGATAATTCATTACTATTTACTGAAAGCATCCAGAAAACt                                                                                                                                                                                                                                                                                         >  1:197369/1‑288 (MQ=255)
            ttAGCCACAGAGAATGTTCATATTGTATTAAAAAATGATAATACAGTGTTATTAAAGGGCCTAAAAAACATTGTGTCTGTCAAATTTTCAAAGGATAGGCATCTTATAGAAACGACCTCTAATAAGCTGAAATCCAGAGAGATCACATTTCAGGAATACAGAAGAAACCTTGCTAAAGCAGGAGTTTTTCGGTGGGTTACAAATATCCACGAACAAAAAAGATATTACTATACCTTTGATAATTCATTACTATTTACTGAAAGCATCCAGAAAACTACACAGATCTTa                                                                                                                                                                                                                                                                             >  2:726489/1‑288 (MQ=255)
             tAGCCACAGAGAATGTTCATATTGTATTAAAAAATGATAATACAGTGTTATTAAAGGGCCTAAAAAACATTGTGTCTGTCAAATTTTCAAAGGATAGGCATCTTATAGAAACGACCTCTAATAAGCTGCAATCCAGAGAGAGCACATTTAAGGAATACAGAAGAAACCTTGCTAAAGCAGGAGTTTGTGGGTGGGTTACAAATATCCACGACCAACAAAGATATAACTATCCCTTTGATAATTAATTACTATTTACTGAAAGCATCAAGGAAAATACACAGATATCAcc                                                                                                                                                                                                                                                                           >  2:732746/1‑289 (MQ=255)
             tAGCCACAGAGAATGTTCATATTGTATTAAAAAATGATAATACAGTGTTATTAAAGGGCCTAAAAAACATTGTGTCTGTCAAATTTTCAAAGGATAGGCATCTTATAGAAACGACCTCTAATAAGCTGAAATCCAGAGAGATCACATTTCAGGAATACAGAAGAAACCTTGCTAAAGCAGGAGTTTTTCGGTGGGTTACAAATATCCACGAACAAAAAAGATATTACTATACCTTTGATAATTCATTACTATTTACTGAAAGCATCCAGaaaa                                                                                                                                                                                                                                                                                           >  1:179829/1‑273 (MQ=255)
             tAGCCACAGAGAATGTTCATATTGTATTAAAAAATGATAATACAGTGTTATTAAAGGGCCTAAAAAACATTGTGTCTGTCAAATTTTCAAAGGATAGGCATCTTATAGAAACGACCTCTAATAAGCTGAAATCCAGAGAGATCACATTTCAGGAATACAGAAGAAACCTTGCTAAAGCAGGAGTTTTTCGGTGGGTTACAAATATCCACGAACAAAAAAGATATTACTATACCTTTGATAATTCATTACTATTTACTGAAAGCATCCAGaaaa                                                                                                                                                                                                                                                                                           <  2:179829/273‑1 (MQ=255)
             tAGCCACAGAGAATGTTCATATTGTATTAAAAAATGATAATACAGTGTTATTAAAGGGCCTAAAAAACATTGTGTCTGTCAAATTTTCAAAGGATAGGCATCTTATAGAAACGACCTCTAATAAGCTGAAATCCAGAGAGATCACATTTCAGGAATACAGAAGAAACCTTGCTAAAGCAGGAGTTTTTCGGTGGGTTACAAATATCCACGAACAAAAAAGATATTACTATACCTTTGATAATTCATTACTATTTACTGAAAGCATCCAGGAAACTACACAGATCTTAca                                                                                                                                                                                                                                                                           >  2:732767/1‑288 (MQ=255)
                 cacaGAGAATGTTCATATTGTATTAAAAAATGATAATACAGTGTTATTAAAGGGCCTAAAAAACATTGTGTCTGTCAAATTTTCAAAGGATAGGCATCTTATAGAAACGACCTCTAATAAGCTGAAATCCAGAGAGATCACATTTCAGGAATACAGAAGAAACCTTGCTAAAGCAGGAGTTTTTCGGTGGGTTACAAATATCCACGAACAAAAAAGATATTACTATACCTTTGATAATTCATTACTATTTACTGAAAGCATCCAGAAAACTACACAGATCTTACCACg                                                                                                                                                                                                                                                                        <  1:726489/288‑1 (MQ=255)
                                         aaaaaaTGATAATACAGTGTTATTAAAGGGCCTAAAAAACATTGTGTCTGTCAAATTTTCAAAGGATAGGCATCTTATAGAAACGACCTCTAATAAGCTGAAATCCAGAGAGATCACATTTCAGGAATACAGAAGAAACCTTGCTAAAGCAGGAGTTTTTCGGTGGGTTACAAATATCCACGAACAAAAAAGATATTACTATACCTTTGATAATTCATTACTATTTACTGAAAGCATCCAGAAAACTACACAGATCTTACCACGCTAAACCATAACGTCCGGCTtctct                                                                                                                                                                                                                                               <  2:20922/289‑1 (MQ=255)
                                                          tgtTATTAAAGGGCCTAAAAAACATTGTGTCTGTCAAATTTTCAAAGGATAGGCATCTTATAGAAACGACCTCTAATAAGCTGAAATCCAGAGAGATCACATTTCAGGAATACAGAAGAAACCTTGCTAAAGCAGGAGTTTTTCGGTGGGTTACAAATATCCACGAACAAAAAAGATATTACTATACCTTTGATAATTCATTACTATTTACTGAAAGCATCCAGAAAACTACACAGATCTTACCACGCTAAACCATAACGTCCGGCTTCTCTCACTCCTGAGCCgg                                                                                                                                                                                                                                 >  2:700291/1‑286 (MQ=255)
                                                              attaAAGGGCCTAAAAAACATTGTGTCTGTCAAATTTTCAAAGGATAGGCATCTTATAGAAACGACCTCTAATAAGCTGAAATCCAGAGAGATCACATTTCAGGAATACAGAAGAAACCTTGCTAAAGCAGGAGTTTTTCGGTGGGTTACAAATATCCACGAACAAAAAAGATATTACTATACCTTTGATAATTCATTACTATTTACTGAAAGCATCCAGAAAACTACACAGATCTTACCACGCTAAACCATAACGTCCGGCTTCTCTCACTCCTGAGCCGGACTGCAt                                                                                                                                                                                                                          <  1:700291/289‑1 (MQ=255)
                                                                    gggCCTAAAAAACATTGTGTCTGTCAAATTTTCAAAGGATAGGCATCTTATAGAAACGACCTCTAATAAGCTGAAATCCAGAGAGATCACATTTCAGGAATACAGAAGAAACCTTGCTAAAGCAGGAGTTTTTCGGTGGGTTACAAATATCCACGAACAAAAAAGATATTACTATACCTTTGATAATTCATTACTATTTACTGAAAGCATCCAGAAAACTACACAGATCTTACCACGCTAAACCATAACGTCCGGCTTCTCTCACTCCTGAGCCGGACTGCAt                                                                                                                                                                                                                          <  2:611447/283‑1 (MQ=255)
                                                                    gggCCTAAAAAACATTGTGTCTGTCAAATTTTCAAAGGATAGGCATCTTATAGAAACGACCTCTAATAAGCTGAAATCCAGAGAGATCACATTTCAGGAATACAGAAGAAACCTTGCTAAAGCAGGAGTTTTTCGGTGGGTTACAAATATCCACGAACAAAAAAGATATTACTATACCTTTGATAATTCATTACTATTTACTGAAAGCATCCAGAAAACTACACAGATCTTACCACGCTAAACCATAACGTCCGGCTTCTCTCACTCCTGAGCCGGACTGCAt                                                                                                                                                                                                                          >  1:611447/1‑283 (MQ=255)
                                                                                        ctgtcAAATTTTCAAAGGATAGGCATCTTATAGAAACGACCTCTAATAAGCTGAAATCCAGAGAGATCACATTTCAGGAATACAGAAGAAACCTTGCTAAAGCAGGAGTTTTTCGGTGGGTTACAAATATCCACGAACAAAAAAGATATTACTATACCTTTGATAATTCATTACTATTTACTGAAAGCATCCAGAAAACTACACAGATCTTACCACGCTAAACCATAACGTCCGGCTTCTCTCACTCCTGAGCCGGACTGCATTGGTTTAATAAAAACCATCAACAAtt                                                                                                                                                                                                <  2:197369/289‑1 (MQ=255)
                                                                                                                                          cTGAAATCCAGAGAGATCACATTTCAGGAATACAGAAGAAACCTTGCTAAAGCAGGAGTTTTTCGGTGGGTTACAAATATCCACGAACAAAAAAGATATTACTATACCTTTGATAATTCATTACTATTTACTGAAAGCATCCAGAAAACTACACAGATCTTACCACGCTAAACCATAACGTCCGGCTTCTCTCACTCCTGAGCCGGACTGCATTGGTTTAATAAAAACCATCAACAATTGTGATTTAg                                                                                                                                                                                       >  2:535523/1‑248 (MQ=255)
                                                                                                                                          cTGAAATCCAGAGAGATCACATTTCAGGAATACAGAAGAAACCTTGCTAAAGCAGGAGTTTTTCGGTGGGTTACAAATATCCACGAACAAAAAAGATATTACTATACCTTTGATAATTCATTACTATTTACTGAAAGCATCCAGAAAACTACACAGATCTTACCACGCTAAACCATAACGTCCGGCTTCTCTCACTCCTGAGCCGGACTGCATTGGTTTAATAAAAACCATCAACAATTGTGATTTAg                                                                                                                                                                                       <  1:535523/248‑1 (MQ=255)
                                                                                                                                                                        tACAGAAGAAACCTTGCTAAAGCAGGAGTTTTTCGGTGGGTTACAAATATCCACGAACAAAAAAGATATTACTATACCTTTGATAATTCATTACTATTTACTGAAAGCATCCAGaa                                                                                                                                                                                                                                                                                             <  1:191306/116‑1 (MQ=255)
                                                                                                                                                                        tACAGAAGAAACCTTGCTAAAGCAGGAGTTTTTCGGTGGGTTACAAATATCCACGAACAAAAAAGATATTACTATACCTTTGATAATTCATTACTATTTACTGAAAGCATCCAGaa                                                                                                                                                                                                                                                                                             >  2:191306/1‑116 (MQ=255)
                                                                                                                                                                                 aaCCTTGCTAAAGCAGTAGTTTTTCGGTGGGTTACAAATATCCACGAACAACAAAGATATTACTATACCTTTGATAATTCATTACTATTTACTGAAAGCATCCAGGAAACTACACAGATCTTACCACGCTAAACCATAACGGCCGGCTTCTCTCACTCCTGAGCCGGACTGCATTGGTTTAATAAAAACCATCAACAATTGTGATTTAGATATTCGGAACCATTCAAATATAACAAAACCCCGTAAAAACGAGGTTTATGGATAAATTTTATTATTGAATACATCAGat                                                                                                       <  1:732746/289‑1 (MQ=255)
                                                                                                                                                                                 aaCCTTGCTAAAGCAGGAGTTTTTCGGTGGGTTACAAATATCCACGAACAAAAAAGATATTACTATACCTTTGATAATTCATTACTATTTACTGAAAGCATCCAGGAAACTACACAGATCTTACCACGCTAAACCATAACGTCCGGCTTCTCTCACTCCTGAGCCGGACTGCATTGGTTTAATAAAAACCATCAACAATTGTGATTTAGATATTCGGAACCATTCAAATATAACAAAACCCCGTAAAAACGAGGTTTATGGATAAATTTTATTATTGAATACATCAGat                                                                                                       <  1:732767/289‑1 (MQ=255)
                                                                                                                                                                                 aaCCTTGCTAAAGCAGGAGTTTTTCGGTGGGTTACAAATATCCACGAACAAAAAAGATATTACTATACCTTTGATAATTCATTACTATTTACTGAAAGCATCCAGAAAACTACACAGATCTTACCACGCTAAACCATAACGTCCGGCTTCTCTCACTCCTGAGCCGGACTGCATTGGTTTAATAAAAACCATCAACAATTGTGATTTAGATATTCGGAACCATTCAAATATAACAAAACCCCGTAAAAACGAGGTTTATGGATAAATTTTATTATTGAATACATCAGa                                                                                                        <  1:793092/288‑1 (MQ=255)
                                                                                                                                                                                                 gAGTTTTTCGGTGGGTTACAAATATCCACGAACAAAAAAGATATTACTATACCTTTGATAATTCATTACTATTTACTGAAAGCATCCAGAAAACTACACAGATCTTACCACGCTAAACCATAACGTCCGGCTTCTCTCACTCCTGAGCCGGACTGCATTGGTTTAATAAAAACCATCAACAATTGTGATTTAGATATTCGGAACCATTCAAATATAACAAAACCCCGTAAAAACGAGGTTTATGGATAAATTTTATTATTGAATACATCAGATTAAATTAATCTTGAca                                                                                       >  1:339871/1‑289 (MQ=255)
                                                                                                                                                                                                                          ccACGAACAAAAAAGATATTACTATACCTTTGATAATTCATTACTATTTACTGAAAGCATCCAGAAAACTACACAGATCTTACCACGCTAAACCATAACGTCCGGCTTCTCTCACTCCTGAGCCGGACTGCATTGGTTTAATAAAAACCATCAACAATTGTGATTTAGATATTTGGAACCATTCAAATATAACAAAACCCCGTAAAAACg                                                                                                                                             <  2:458540/210‑1 (MQ=255)
                                                                                                                                                                                                                          ccACGAACAAAAAAGATATTACTATACCTTTGATAATTCATTACTATTTACTGAAAGCATCCAGAAAACTACACAGATCTTACCACGCTAAACCATAACGTCCGGCTTCTCTCACTCCTGAGCCGGACTGCATTGGTTTAATAAAAACCATCAACAATTGTGATTTAGATATTCGGAACCATTCAAATATAACAAAACCCCGTAAAAACg                                                                                                                                             >  1:458540/1‑210 (MQ=255)
                                                                                                                                                                                                                                                       ttGATAATTCATTACTATTTACTGAAAGCATCCAGAAAACTACACAGATCTTACCACGCTAAACCATAACGTCCGGCTTCTCTCACTCCTGAGCCGGACTGCATTGGTTTAATAAAAACCATCAACAATTGTGATTTAGATATTCGGAACCATTCAAATATAACAAAACCCCGTAAAAACGAGGTTTATGGATAAATTTTATTATTGAATACATCAGATTAAATTAATCTTGACATCATAGCTTTCAAGACCCGTCATTTTTTCCCGTGCGGTAAACTGAATACTGGTa                                 >  1:270459/1‑289 (MQ=255)
                                                                                                                                                                                                                                                       ttGATAATTCATTACTATTTACTGAAAGCATCCAGAAAACTACACAGATCTTACCACGCTAAACCATAACGTCCGGCTTCTCTCACTCCTGAGCCGGACTGCATTGGTTTAATAAAAACCATCAACAATTGTGATTTAGATATTCGGAACCATTCAAATATAACAAAACCCCGTAAAAACGAGGTTTATGGATAAATTTTATTATTGAATACATCAGATTAAATTAATCTTGACATCATAGCTTTCAAGACCCGTCATTTTTTCCCGTGCGGTAAACTGAATACTGGTa                                 >  2:76863/1‑289 (MQ=255)
                                                                                                                                                                                                                                                         gATAATTCATTACTATTTACTGAAAGCATCCAGAAAACTACACAGATCTTACCACGCTAAACCATAACGTCCGGCTTCTCTCACTCCTGAGCCGGACTGCATTGGTTTAATAAAAACCATCAACAATTGTGATTTAGATATTCGGAACCATTCAAATATAACAAAACCCCGTAAAAACGAGGTTTATGGATAAATTTTATTATTGAATACATCAGATTAAATTAATCTTGACATCATAGCTTTCAAGACCCGTCATTTTTTCCCGTGCGGTAAACTGAATACTGGTa                                 >  1:578230/1‑287 (MQ=255)
                                                                                                                                                                                                                                                                   tACTATTTACTGAAAGCATCCAGAAAACTACACAGATCTTACCACGCTAAACCATAACGTCCGGCTTCTCTCACTCCTGAGCCGGACTGCATTGGTTTAATAAAAACCATCAACAATTGTGATTTAGATATTCGGAACCATTCAAATATAACAAAACCCCGTAAAAACGAGGTTTATGGATAAATTTTATTATTGAATACATCAGATTAAATTAATCTTGACATCATAGCTTTCAAGACCCGTCATTTTTTCCCGTGCGGTAAACt                                            >  1:830533/1‑266 (MQ=255)
                                                                                                                                                                                                                                                                   tACTATTTACTGAAAGCATCCAGAAAACTACACAGATCTTACCACGCTAAACCATAACGTCCGGCTTCTCTCACTCCTGAGCCGGACTGCATTGGTTTAATAAAAACCATCAACAATTGTGATTTAGATATTCGGAACCATTCAAATATAACAAAACCCCGTAAAAACGAGGTTTATGGATAAATTTTATTATTGAATACATCAGATTAAATTAATCTTGACATCATAGCTTTCAAGACCCGTCATTTTTTCCCGTGCGGTAAACt                                            <  2:830533/266‑1 (MQ=255)
                                                                                                                                                                                                                                                                                    aTCCAGAAAACTACACAGATCTTACCACGCTAAACCATAACGTCCGGCTTCTCTCACTCCTGAGCCGGACTGCATTGGTTTAATAAAAACCATCAACAATTGTGATTTAGATATTCGGAACCATTCAAata                                                                                                                                                                  <  1:72231/131‑1 (MQ=255)
                                                                                                                                                                                                                                                                                    aTCCAGAAAACTACACAGATCTTACCACGCTAAACCATAACGTCCGGCTTCTCTCACTCCTGAGCCGGACTGCATTGGTTTAATAAAAACCATCAACAATTGTGATTTAGATATTCGGAACCATTCAAata                                                                                                                                                                  >  2:72231/1‑131 (MQ=255)
                                                                                                                                                                                                                                                                                      ccAGAAAACTACACAGATCTTACCACGCTAAACCATAACGTCCGGCTTCTCTCACTCCTGAGCCGGACTGCATTGGTTTAATAAAAACCATCAACAATTGTGATTTAGATATTCGGAACCATTCAAATATAACAAAACCCCGTAAAAACGAGGTTTATGGATAAATTTTATTATTGAATACATCAGATTAAATTAATCTTGACATCATAGCTTTCAAGACCCGTCATTTTTTCCCGTGCGGTAAACTGAATACTGGTAACTTCTTTCCCGGTCTTTTTCTTAAGTTCaa  >  1:173902/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                            
TACATTTACTATTTAGCCACAGAGAATGTTCATATTGTATTAAAAAATGATAATACAGTGTTATTAAAGGGCCTAAAAAACATTGTGTCTGTCAAATTTTCAAAGGATAGGCATCTTATAGAAACGACCTCTAATAAGCTGAAATCCAGAGAGATCACATTTCAGGAATACAGAAGAAACCTTGCTAAAGCAGGAGTTTTTCGGTGGGTTACAAATATCCACGAACAAAAAAGATATTACTATACCTTTGATAATTCATTACTATTTACTGAAAGCATCCAGAAAACTACACAGATCTTACCACGCTAAACCATAACGTCCGGCTTCTCTCACTCCTGAGCCGGACTGCATTGGTTTAATAAAAACCATCAACAATTGTGATTTAGATATTCGGAACCATTCAAATATAACAAAACCCCGTAAAAACGAGGTTTATGGATAAATTTTATTATTGAATACATCAGATTAAATTAATCTTGACATCATAGCTTTCAAGACCCGTCATTTTTTCCCGTGCGGTAAACTGAATACTGGTAACTTCTTTCCCGGTCTTTTTCTTAAGTTCAA  >  NC_000913/1637149‑1637715

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