Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,769,648 (C)6→5 intergenic (+162/‑55) yfjQ → / → yfjR CP4‑57 prophage, uncharacterized protein/CP4‑57 prophage, putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,769,6430C.90.0% 31.6 / ‑1.7 10intergenic (+157/‑60)yfjQ/yfjRCP4‑57 prophage, uncharacterized protein/CP4‑57 prophage, putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base . (5/4);  minor base G (0/1);  total (5/5)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

GTTAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCAGCCGTCTGGTTGCGGGT  >  NC_000913/2769507‑2769777
                                                                                                                                        |                                                                                                                                      
gttAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGccc                                                                                                                                      <  1:17512/139‑1 (MQ=255)
gttAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGccc                                                                                                                                      <  1:139278/139‑1 (MQ=255)
gttAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTGTCCATACCCTGTCCGcccc                                                                                                                                     <  2:3980/139‑1 (MQ=255)
               caccTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGcccc                                                                                                                                     <  1:338459/125‑1 (MQ=255)
               caccTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGcccc                                                                                                                                     >  2:338459/1‑125 (MQ=255)
                                           tgCCTCAATCCCCCGGTTATAGCTTTTAACCCCCATTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCTCTTTAAAGTAAGCACATCATTTTCAGTCAGTTAACTTTcc                                                                                          <  2:648868/139‑1 (MQ=255)
                                             ccttagtCCCCCGGTTATAGCTTCTAACGCGCGATACATCTGGCTTTTGCAGATATAAAAAATAGTTCCTGCGTTGTCCATGCCCTGTCGGCCCCCTCTTTAGCGTAATCACATCATTTTCAGTCAGTTAACTTTCCTg                                                                                        <  1:74835/133‑1 (MQ=18)
                                             ccTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCTCTTTAGAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTg                                                                                        >  2:634936/1‑139 (MQ=255)
                                             ccTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTg                                                                                        >  2:565066/1‑139 (MQ=255)
                                                cAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCTCTTTAAAGTAATCACATCATTTtc                                                                                                          <  2:724743/118‑1 (MQ=255)
                                                cAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCTCTTTAAAGTAATCACATCATTTtc                                                                                                          >  1:724743/1‑118 (MQ=255)
                                                       ccGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGccc                                                                                                                                      <  2:163231/84‑1 (MQ=255)
                                                       ccGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGccc                                                                                                                                      <  2:352359/84‑1 (MQ=255)
                                                       ccGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGccc                                                                                                                                      >  1:352359/1‑84 (MQ=255)
                                                       ccGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGccc                                                                                                                                      >  1:163231/1‑84 (MQ=255)
                                                                    ttAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCTCTTTAAAGTAATCACATCATTTtca                                                                                                         <  2:281671/99‑1 (MQ=255)
                                                                    ttAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCTCTTTAAAGTAATCACATCATTTtca                                                                                                         >  1:281671/1‑99 (MQ=255)
                                                                           ccGTTACATCTGGCTTTTGCAGATATAAAAAATAGTTTCAGCGTTGTCCATACCCTGTCCGccccctct                                                                                                                                 <  1:583111/69‑4 (MQ=21)
                                                                           ccGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGccccctct                                                                                                                                 >  2:583111/1‑66 (MQ=255)
                                                                                                  aaTAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTggctgtcagg                                   <  1:565066/139‑1 (MQ=255)
                                                                                                                                   gTCCGCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCAGCCGTCTGGTTGCGGGt  >  2:341846/1‑139 (MQ=255)
                                                                                                                                        |                                                                                                                                      
GTTAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCAGCCGTCTGGTTGCGGGT  >  NC_000913/2769507‑2769777

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