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.93.3% 114.3 / 2.1 30intergenic (+157/‑60)yfjQ/yfjRCP4‑57 prophage, uncharacterized protein/CP4‑57 prophage, putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base . (16/12);  total (17/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.47e-01

TGTTAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCAGCCGTCTGGTTGCGGGT  >  NC_000913/2769506‑2769777
                                                                                                                                         |                                                                                                                                      
tgttAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGcc                                                                                                                                       <  2:956233/139‑1 (MQ=255)
 gttAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGccc                                                                                                                                      <  2:84004/139‑1 (MQ=255)
  ttAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGcccc                                                                                                                                     <  2:616730/139‑1 (MQ=255)
  ttAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCGTTTAACCCCCGTTACAGCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGcccc                                                                                                                                     <  2:1085298/139‑1 (MQ=255)
            acacCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGcccc                                                                                                                                     >  2:733301/1‑129 (MQ=255)
            acacCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGcccc                                                                                                                                     <  1:733301/129‑1 (MQ=255)
            acacCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCTCTTTAAAg                                                                                                                          >  2:799542/1‑139 (MQ=255)
                caccTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCTCTTTAAAGTAAt                                                                                                                      <  2:433291/139‑1 (MQ=255)
                  ccTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCTCTTTAAAGTAATca                                                                                                                    >  1:412389/1‑139 (MQ=255)
                     gTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCTCTTTAAAGTAATCAcat                                                                                                                 <  1:633060/139‑1 (MQ=255)
                        cGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCTCTTTAAAGTAAt                                                                                                                      <  2:622650/131‑1 (MQ=255)
                        cGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCTCTTTAAAGTAAt                                                                                                                      >  1:622650/1‑131 (MQ=255)
                          cATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCTCTTTAAAGTAATCACATCAttt                                                                                                            >  1:68263/1‑139 (MQ=255)
                                  tgtgCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCTCTTTAAAGTAATCACATCATTTtcagtcag                                                                                                    >  1:417077/1‑139 (MQ=255)
                                       tCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAAc                                                                                               >  1:1247/1‑139 (MQ=255)
                                       tCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAAc                                                                                               >  2:409460/1‑139 (MQ=255)
                                         gtgtGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACtt                                                                                             >  1:730869/1‑139 (MQ=255)
                                                 cAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGgag                                                                                     <  2:322958/139‑1 (MQ=255)
                                                   aTCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGaa                                                                                   <  1:587716/139‑1 (MQ=255)
                                                     cccccGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAAtc                                                                                 <  2:412389/139‑1 (MQ=255)
                                                        ccGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGccccctct                                                                                                                                 >  1:131642/1‑86 (MQ=255)
                                                        ccGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGccccctct                                                                                                                                 <  2:131642/89‑4 (MQ=255)
                                                              ataGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGac                                                                        >  1:542375/1‑139 (MQ=255)
                                                                     ttAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCa                                                                 <  2:542375/139‑1 (MQ=255)
                                                                                      tGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGac                                                                        >  1:631045/1‑115 (MQ=255)
                                                                                      tGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGac                                                                        <  2:631045/115‑1 (MQ=255)
                                                                                       ggCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGccc                                                                                                                                      >  1:641795/1‑53 (MQ=255)
                                                                                       ggCTTTTCCAGAAATAAAAAATATTTTCTGCGTTGTCCATACCCTGTCCGccc                                                                                                                                      <  2:641795/53‑1 (MQ=18)
                                                                                                  aaaTAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCTCTTTAAAGTAATc                                                                                                                     >  1:203879/1‑58 (MQ=255)
                                                                                                  aaaTAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCTCTTTAAAGTAATc                                                                                                                     <  2:203879/58‑1 (MQ=255)
                                                                                                   aaTAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCC‑‑CTTTAAAGTAATCACATCATTTTCAGTCAGTTaa                                                                                                <  1:113840/77‑1 (MQ=17)
                                                                                                   aaTAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCC‑‑CTTTAAAGTAATCACATCATTTTCAGTCAGTTaa                                                                                                >  2:113840/1‑77 (MQ=17)
                                                                                                   aaTAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTggctgtcagg                                   <  1:799542/139‑1 (MQ=255)
                                                                                                   aaTAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTggctgtcagg                                   <  2:1247/139‑1 (MQ=255)
                                                                                                                                    gTCCGCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCggc                                              >  2:861982/1‑95 (MQ=255)
                                                                                                                                    gTCCGCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCggc                                              <  1:861982/95‑1 (MQ=255)
                                                                                                                                    gTCCGCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCAGCCGTCTGGTTGCGGGt  >  2:461892/1‑139 (MQ=255)
                                                                                                                                    gTCCGCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCAGCCGTCTGGTTGCGGGt  >  1:464459/1‑139 (MQ=255)
                                                                                                                                    gTCCGCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCAGCCGTCTGGTTGCGGGt  >  2:860063/1‑139 (MQ=255)
                                                                                                                                    gTCCGCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCAACCGTCTGGTTGCGGGt  >  1:661952/1‑139 (MQ=255)
                                                                                                                                         |                                                                                                                                      
TGTTAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCAGCCGTCTGGTTGCGGGT  >  NC_000913/2769506‑2769777

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