Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,642,370 C→G 11.9% intergenic (+159/+119) ynfS → / ← ydfT Qin prophage; protein YnfS/Qin prophage; putative antitermination protein Q

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,642,3700CG11.9% 111.0 / 9.8 42intergenic (+159/+119)ynfS/ydfTQin prophage; protein YnfS/Qin prophage; putative antitermination protein Q
Reads supporting (aligned to +/- strand):  ref base C (21/16);  new base G (3/2);  total (24/18)
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.24e-01

TTTAGGGACATTGTACTGGAAGAAAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTATCTCGTCACCTTCCTCAATT  >  NC_000913/1642229‑1642511
                                                                                                                                             |                                                                                                                                             
tttAGGGACATTGTACTGGAAGAAAACATTTTACACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCtca                                                                                                                                              <  2:712543/143‑1 (MQ=255)
    gggACATTGTACTGGAAGAAAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAg                                                                                                                                          >  1:101375/1‑143 (MQ=255)
       aCATTGTACTGGAAGAAAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCgg                                                                                                                                       >  1:544617/1‑143 (MQ=255)
       aCATTGTACTGGAAGAAAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCgg                                                                                                                                       <  2:44374/143‑1 (MQ=255)
       aCATAGTACTGGAAGAAAACATTTTAAACATCAGGCAAATAACCTGCGCGTAACCGTAGTTGTAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCgg                                                                                                                                       <  2:166080/143‑1 (MQ=255)
           tGTAATGGAAGGTAACATTTGAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACACACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTGATGAGCGGGttt                                                                                                                                   <  2:583014/143‑1 (MQ=255)
              aCTGGAAGAAAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTGATGAGCGGGttttgt                                                                                                                                <  1:280979/143‑3 (MQ=255)
               cTGGAAGAAAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTa                                                                                                                               <  2:705901/143‑1 (MQ=255)
                      aaaaCATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTc                                                                                                                        <  1:435688/143‑1 (MQ=255)
                      aaaaCATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTc                                                                                                                        >  1:734167/1‑143 (MQ=255)
                              ttAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAACCGGTTTTCTACTTTTTCTTAACGTc                                                                                                                >  2:368122/1‑143 (MQ=255)
                               tAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAACACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCg                                                                                                               <  2:133210/143‑1 (MQ=255)
                                       aGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACa                                                                                                       >  1:750138/1‑143 (MQ=255)
                                       aGGCAAATAACCAAGTCACCAGCTAAATAAATAGTTAACAGACAGCAGTCACGGGATGAGATTCAACATTACCATTGCCCCAGTTAAAGCACAAAACCCGCTCATCAGCGGGTTGTCTACTTTTTCTTAACGTCGGGTATACa                                                                                                       >  2:730406/1‑143 (MQ=255)
                                               aaCCAAGTCACCAGCTAAATAATAAGTTAACAGACATTAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCGCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATTCAGAGCCCAt                                                                                               >  2:201645/1‑143 (MQ=255)
                                                  cAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGt                                                                                            >  2:363166/1‑143 (MQ=255)
                                                   aaGTCACCAGCTAAAGAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGtt                                                                                           <  2:464448/143‑1 (MQ=255)
                                                          cAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGaaaaaa                                                                                    <  2:620176/143‑1 (MQ=255)
                                                                   aataaGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTACCGTCGGGGATACAAAGCCCATAGTCGCAAACACTTTCTCCa                                                                           >  2:520195/1‑143 (MQ=255)
                                                                    ataaGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGCGTTTTCTACTTTTTCTTAACGTCGCGCATACACAGCCCATCGTCGAAACAAtttcaactac                                                                          >  2:770662/1‑136 (MQ=255)
                                                                           aaCAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATCtttttt                                                                   >  1:761198/1‑143 (MQ=255)
                                                                                  aTGAGTCCCGGGAGGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAAt                                                            <  1:86767/143‑1 (MQ=255)
                                                                                       tCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCCTCGTTGAAAAAAATTTATCCATATTTTTTGAACAACGCaaa                                                       >  2:701329/1‑143 (MQ=255)
                                                                                       gcccGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCaaa                                                       <  1:625601/142‑1 (MQ=255)
                                                                                          cGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAAcat                                                    >  1:234852/1‑143 (MQ=255)
                                                                                           gggATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAAcatc                                                   <  1:156212/143‑1 (MQ=255)
                                                                                                       aaCATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCt                                       <  1:244237/143‑1 (MQ=255)
                                                                                                        aCATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCtt                                      <  1:542152/143‑1 (MQ=255)
                                                                                                          aTTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATAAAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGCCGCAATATTCa                                    >  2:647/1‑143 (MQ=255)
                                                                                                            tACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGc                                  <  2:251519/143‑1 (MQ=255)
                                                                                                                aTTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTCTAAAAAAATGCAAACATCATGTCGCCATCTTCAGCaaca                              >  1:742849/1‑141 (MQ=255)
                                                                                                                 ttGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATACTTTTTTAAAAATGCAAACATCATGTCGCCATCTTCAGCaaaaa                             <  2:169293/143‑1 (MQ=255)
                                                                                                                    ccccATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAACATCa                          >  2:770604/1‑143 (MQ=255)
                                                                                                                     cccATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCAt                         >  1:637369/1‑143 (MQ=255)
                                                                                                                      ccATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCAtt                        >  1:467449/1‑143 (MQ=255)
                                                                                                                         tttAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTGAAAGCTGGTATACAAAGCCCATCGTTGAAAAAATTTTCTCCATATTTTTTGAAAAATGCAAACATCATGTCGCCACTCTCAGCAACAATCATTTAt                     >  2:168527/1‑143 (MQ=255)
                                                                                                                          ttAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTAtc                    <  2:362039/143‑1 (MQ=255)
                                                                                                                                    aaaCCCGCTGATGAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTATCTCGTCACCtt          >  1:580889/1‑143 (MQ=255)
                                                                                                                                    aaaCCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATTGTTGTATAAATTTTATCCATATTTTTTGAAACATGCAAACATCATGTTGCCATCCTCAGCAACAATCATTTATCTCGTCACCtt          >  2:727260/1‑143 (MQ=255)
                                                                                                                                     aaCCCGCTGATGAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTATCTCGTCACCTTc         >  1:632326/1‑143 (MQ=255)
                                                                                                                                     aaCCCGCTGATGAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTATCTCGTCACCTTc         >  1:639400/1‑143 (MQ=255)
                                                                                                                                      aCCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTATCTCGTCACCTTcc        <  1:730406/143‑1 (MQ=255)
                                                                                                                                          gCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGCAAAAATTTTATCAATATTTTTTGAAAAATGCAAACATCATGTCGCCGTCTTCAGGCGAAATCATTTATCACGTCACCTTTCTCaa    >  1:295400/1‑143 (MQ=255)
                                                                                                                                            tcagcaGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTATCTCGTCACCTTCCTCAAtt  <  2:751254/143‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                             
TTTAGGGACATTGTACTGGAAGAAAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTATCTCGTCACCTTCCTCAATT  >  NC_000913/1642229‑1642511

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