Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,429,767 C→T A240V (GCG→GTG)  stfR → putative prophage side tail fiber protein StfR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,429,7670CT95.8% 69.4 / NA 25A240V (GCG→GTG) stfRputative prophage side tail fiber protein StfR
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base T (10/13);  minor base G (1/0);  total (12/13)
Fisher's exact test for biased strand distribution p-value = 4.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACAGGCCGCGTCGTCGGCTCAGTCAGCGTCTTCCAGCGCAGGAACGGCATCAACAAAGGCCACTGAAGCATCAAAAAGTGCTGCCGCTGCAGAGTCCTCAAAAAGCGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCAGCGTCACTACAATCAGCAGCCACATCTGCAT‑CCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTC  >  NC_000913/1429527‑1429959
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                
gACAGGCTGCATCGTCAGCTCAGGAAGCGTCCTCCGGCGCAGAAGCGGCATCAGCAAAGGCCACTGAAGCGGAAAAAAGGGCCGCAGCCGCAGAGTCCTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCT‑CCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCGCaaaa                                                                                                                                                                                        <  1:290812/251‑1 (MQ=255)
               cagctaaggaagCGTCCTCCGGCGCAGAAGCGGCATCAGCAAAGGCCACTGAAGCGGAAAAAAGTGCCGCAGCCGCAGAGTCCTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCT‑CCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCCGCAAAatc                                                                                                                                                                          >  2:292519/11‑250 (MQ=255)
                                      gcagaagcGGCATCAGCAAAGGCCACTGAAGCGGAAAAAAGTGCCGCAGCCGCAGAGTCCTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCT‑CCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGc                                                                                                                                                           <  1:140396/235‑1 (MQ=255)
                                      gcagaagcGGCATCAGCAAAGGCCACTGAAGCGGAAAAAAGTGCCGCAGCCGCAGAGTCCTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCT‑CCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGc                                                                                                                                                           >  2:140396/8‑242 (MQ=255)
                                      gcAGGAGCGGCATCAGCAAAGGCCACTGAGACGAAAAAAAGTGCCCCAGCCTCAGAGCCCTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCGGCAGCGTCACAACAATCACCCGCACCGTCTGCCC‑CCACCGCTGCCACGAAAGCGGCAGAGGCCGCCACGTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAATCGAACGCATCATCAAg                                                                                                                                                  <  2:228919/251‑1 (MQ=255)
                                        agaagcGGCATCAGCAAAGGCCACTGAAGCGGAAAAAAGTGCCGCAGCCGCAGAGTCCTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCT‑CCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAACTCATCAGAAACGAACGCATCATCAAgtg                                                                                                                                                >  2:282569/6‑249 (MQ=255)
                                                   cagcaaAGGCCACTGAAGCGGACAAAAGCGCCGCAGCCGCCGAGTCCTCAAAAAACGCGGCGGCCACCAGTGCCGGCGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGCCACAACAATCAGCCGCCACGTCTGCCTACCAC‑GCGGCCACGCAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGGGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGccggtcgggc                                                                                                                                      >  2:187912/4‑242 (MQ=255)
                                                        aaGGCCACTGAAGCGGAAAAAAGTGCCGCAGCCGCAGAGTCCTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCT‑CCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTc                                                                                                                                >  2:124821/1‑251 (MQ=255)
                                                           gCCACTGAAGCGGAAAAAAGTGCCGCAGCCGCAGAGTCCTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCT‑CCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTc                                                                                                                             <  1:282569/251‑1 (MQ=255)
                                                                     cggaaaaaaGTGCCGCAGCCGCAGAGTCCTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCT‑CCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAAcggcg                                                                                                                   >  2:25873/5‑251 (MQ=255)
                                                                      ggaaaaaagtcccgcagccgCAGAGTCCTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCT‑CCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAAcggcgg                                                                                                                  <  1:90255/232‑1 (MQ=255)
                                                                            aaGTGCCGCAGCCGCAGAGTCCTCAAAAAACGCGGCGGCCACCAGGGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCT‑CCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAgaaa                                                                                                            <  1:387432/251‑4 (MQ=255)
                                                                             agtgccgcagccgCAGAGTCCTCAAAAAACGCCGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAACGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCT‑CCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGCGGCCTCAAAAGAGGCAGCAAAATCCTCCGAAACGAACGCATCATCAAGCGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAgaaa                                                                                                            >  1:148141/13‑247 (MQ=255)
                                                                                    cagccgCAGAGTCCTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCT‑CCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTctgcc                                                                                                    <  2:49892/246‑5 (MQ=255)
                                                                                      gccgCAGAGTCCTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCT‑CCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTctgccag                                                                                                  <  1:25873/248‑7 (MQ=255)
                                                                                                ccTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGACAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCT‑CCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGAAGCGGTGGCCTCAAAAGAGGCAGCACAATCATCAGAAACGACCGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTATGCCAGGGCGGCaaa                                                                                         >  1:193819/1‑250 (MQ=255)
                                                                                                        aacgcgGCGGCCCCCAGGGCCGGGGGGGCGAAAACGTCGGAAACGAATGCTGCAGCGTCACAACAAGCAGCCGCCACGTCTGCCT‑CCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCga                                                                                <  2:109587/249‑1 (MQ=255)
                                                                                                            cggcggCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCT‑CCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTATGCCAGGGCGGCAAAAACGTCCGAGACg                                                                            >  1:378681/1‑251 (MQ=255)
                                                                                                                      caGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCT‑CCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCCGCTTCCTCGGCAACGGCGGCAGACAATTCTGCAAGGGCGGCAAAAACGACCGAGACGAATGCCAggg                                                                  >  2:332566/1‑250 (MQ=255)
                                                                                                                               ggtggagAAAACGTAAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCT‑CCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGaaa                                                         <  1:346502/245‑1 (MQ=255)
                                                                                                                                gcgccgAAAAGGCCAGAGACGAAGGCTGCAGCGTACCAACAACCAGCCGCCACGTCTGCCT‑CCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCGCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAGAAACGTCCGATACGAATGCCAGTGCATCTGAAAc                                                        <  2:267531/247‑1 (MQ=255)
                                                                                                                                  ggcgAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCT‑CCCCCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAAcag                                                      <  1:218591/251‑3 (MQ=255)
                                                                                                                                                           gcAGCGCCACAACAATCAGCCGCCACGTCTGCCC‑CCCCCGCGGCCACGAAAGCTTCAGAGGACGCCACTTCAGCACGAGATGTGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGTTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAGCGCCTctgccg                              <  1:359887/249‑6 (MQ=255)
                                                                                                                                                                   acaaCAATCAGCCGCCACGTCTGCCT‑CCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACCGCAGCGGCACGGAGCGACTCTGCCGCGGCAgacg                     >  2:362131/4‑248 (MQ=255)
                                                                                                                                                                                        tgcct‑CCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAGCGCCTCTGCCGCGGCAGACGCAAAAACAgcggcggcggg  >  2:288181/5‑247 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                
GACAGGCCGCGTCGTCGGCTCAGTCAGCGTCTTCCAGCGCAGGAACGGCATCAACAAAGGCCACTGAAGCATCAAAAAGTGCTGCCGCTGCAGAGTCCTCAAAAAGCGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCAGCGTCACTACAATCAGCAGCCACATCTGCAT‑CCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTC  >  NC_000913/1429527‑1429959

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