Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,642,373 C→G 11.9% intergenic (+162/+116) 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,3730CG11.9% 119.2 / 9.8 42intergenic (+162/+116)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 = 3.37e-01

GGGACATTGTACTGGAAGAAAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTATCTCGTCACCTTCCTCAATTGCG  >  NC_000913/1642233‑1642514
                                                                                                                                            |                                                                                                                                             
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:734167/1‑143 (MQ=255)
                  aaaaCATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTc                                                                                                                           <  1:435688/143‑1 (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:639400/1‑143 (MQ=255)
                                                                                                                                 aaCCCGCTGATGAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTATCTCGTCACCTTc            >  1:632326/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)
                                                                                                                                           tcaGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTATCTCGTCACCTTCCTCAATTgcg  >  1:472159/1‑143 (MQ=255)
                                                                                                                                            |                                                                                                                                             
GGGACATTGTACTGGAAGAAAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTATCTCGTCACCTTCCTCAATTGCG  >  NC_000913/1642233‑1642514

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: