Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000964 3,770,066 (A)8→9 intergenic (‑613/+38) ywzE ← / ← csbD conserved protein of unknown function (putative phage interference protein)/stress response protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009643,770,0661.A95.2% 60.7 / ‑3.2 21intergenic (‑613/+38)ywzE/csbDconserved protein of unknown function (putative phage interference protein)/stress response protein
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base A (9/11);  total (10/11)
Fisher's exact test for biased strand distribution p-value = 4.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.92e-01

TTGTGTCGATATCAACTGTTTTCTCGCACTCCCGTGGCTTCCGACACTTTCTTCGGTGCAGTACATTCCGTCATCAAAACCCCTTTGTTTGATTTCTTCAAATGATGTTATAAAAGTTTACAACTTATGTCATCATTATATAGAAGAGGCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGATTAATCTTTATCTGAAAATTTATCCTTGGCTTTCCCGATATCCTTTTGGATTTCGCCTTTTGCCTTATCTTTTTTGCCTTCAGCCTGCATGTCCGTTCTGTCAGCCATATCTCCGACTTTATCCTTCACTTCTCCTTTGGCTTTATTGAAGCCGCCCTTCATTTTGTCTTTTACACTATCGTTACCCATGTTAAT  >  NC_000964/3769820‑3770298
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                        
ttGTGTAGATATCAACTGTTTTCTCGCACTCCCGTGGCTTCCGAAACTTTCTTCGGTGCAGTACATTCCGTCATCAAAACCCCTTTGTTTGATTTCTTCAAATGATGTTATAAAAGTTTACAACTTATGTCATCATTATATAGAAGAGGCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGaaaaaaaaagac                                                                                                                                                                                                                                       >  2:342832/1‑247 (MQ=255)
                   tttCTCGCACTCCCGTGGCTTCCGACACTTTCTTCGGTGCAGTACATTCCGTCATCAAAACCCCTTTGTTTGATTTCTTCAAATGATGTTATAAAAGTTTACAACTTATGTCATCATTATATAGAAGAGGCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGCAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCt                                                                                                                                                                                                                    >  2:275858/1‑251 (MQ=255)
                                cGTGGCTTCCGACACTTTCTTCGGTGCAGTACATTCCGTCATCAAAACCCCTTTGTTTGATTTCTTCAAATGATGTTATAAAAGTTTACAACTTATGTCATCATTATATAGAAGAGGCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTAc                                                                                                                                                                                                       <  2:142946/251‑1 (MQ=255)
                                       tCCGACACTTTCTTCGGTGCAGTACATTCCGTCATCAAAACCCCTTTGTTTGATTTCTTCAAATGATGTTATAAAAGTTTACAACTTATGTCATCATTATATAGAAGAGGCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGATTa                                                                                                                                                                                                  <  1:278116/249‑1 (MQ=255)
                                                         tcAGTACATTCCGTCATCCAAAACCCTTTGTTTGATTTATTCAAATGATGTTATAAAAGTTTACACCTTATGTCATCATTATATAGAAGAGGCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCATACGATTAATCTTTATCTGAAAAtttct                                                                                                                                                                              >  2:52221/2‑249 (MQ=255)
                                                                                cccTTTGTTTGATTTCTTCAAATGATGTTATAAAAGTTTACAACTTATGTCATCATTATATAGAAGAGGCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAATAACCCTTTCGTTTCCTGCCTCTTCCTTACGATTAATCTTTATCTGAAAATTTATCCTTGGCTTTCCCGATATCCttt                                                                                                                                                       <  1:278720/251‑1 (MQ=255)
                                                                                                   aaaTGATGTTATAAAAGTTTACAACTTATGTCATCATTATATAGAAGAGGCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGATTAATCTTTATCTGAAAATTTATCCTTGGCTTTCCCGATATCCTTTTGGATTTCGCCTTTTGCCtt                                                                                                                                   >  2:177616/1‑251 (MQ=255)
                                                                                                             ataAAAGTTTACAACTTATGTCATCATTATATAGAAGAGGCCTATTAATTAAAAGAATATTTAAAAAATTTAGTTAGCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATGAGAACAAGCGAAAAAAAAAGACAAGAAACATTTCGTTCCCGGCATCTGCCGTACGATTAATCTTTATCTGAAAATTTATCCTTGGCTTACCCGCTAGCCTGTTGGATTGCGCCTTTTGCCTTATCttttt                                                                                                                           <  2:111679/250‑1 (MQ=255)
                                                                                                                aaaGTTTACAACTTATGTCATCATTATATAGAAGAGGCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGATTAATCTTTATCTGAAAATTTATCCTTGGCTTTCCCGATATCCTTTTGGATTTCGCCTTTTGCCTTATCTTTTTTGcc                                                                                                                       >  1:304750/1‑251 (MQ=255)
                                                                                                                   gTTTACAACTTATGTCATCATTATATAGAAGAGGCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGATTAATCTTTATCTGAAAATTTATCCTTGGCTTTCCCGATATCCTTTTGGATTTCGCCTTTTGCCTTATCTTTTTTGCCTTc                                                                                                                    <  1:192990/251‑1 (MQ=255)
                                                                                                                         aaCTTATGTCATCATTATATAGAAGAGGCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTGTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGATTAATCTTTATCTGAAAATTTATCCTTGGCTTTCCCGATATCCTTTTGGATTTCGCCTTTTGCCTTATCTTTTTTGCCTTCAGCCTg                                                                                                              <  2:225185/251‑1 (MQ=255)
                                                                                                                                                              tCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGATTAATCTTTATCTGAAAATTTATCCTTGGCTTTCCCGATATCCTTTTGGATTTCGCCTTTTGCCTTATCTTTTTTGCCTTCAGCCTGCATGTCCGTTCTGTCAGCCATATCTCCGACTTTATcc                                                                         >  1:153275/1‑251 (MQ=255)
                                                                                                                                                                  aGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGATTAATCTTTATCTGAAAATTTATCCTTGGCTTTCCCGATATCCTTTTGGATTTCGCCTTTTGCCTTATCTTTTTTGCCTTCAGCCTGCATGTCCGTTCTGTCAGCCATATCTCCGACTTTATCCTTcc                                                                     >  1:327895/1‑250 (MQ=255)
                                                                                                                                                                               aaTTTAGTTACCTTGACCAAATCAGCTGCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGATTAATCTTTATCTGAAAATTTATCCTTGGCTTTCCCGATATCCTTTTGGATTTCGCCTTTTGCCTTATCTTTTTTGCCTTCAGCCTGCATGTCCGTTCTGTCAGCCATATCTCCGACTTTATCCTTCACTGCGCCTTTGGc                                                        <  2:221172/251‑1 (MQ=255)
                                                                                                                                                                                         ccTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGATTAATCTTTATCTGAAAATTTATCCTTGGCTTTCCCGATATCCTTTTGGATTTCGCCTTTTGCCTTATCTTTTTTGCCTTCAGCCTGCATGTCCGTTCTGTCAGCCATATCTCCGACTTTATCCTTCACTTCTCCTTTGGCTTTATTGAAg                                              <  2:159781/251‑1 (MQ=255)
                                                                                                                                                                                                 aaaTCAGCGTCTTTTTTTGTGATATTTATACCCATCAAAAAAAGCAAAAAAAAAAGACAAGAAACAGTTCTTTTCCTGCCTCTTCCGTACGATTAATCTTTATCTGAAAATTTATCCTTGGCTTTCCCGATATCCTTTTGTATTTCGCCTTTTGCCTTATCTTTTTTGCCTTCAGCCTGCATGTCCGTTCAGTCAGCCATATCTCCGACTGTATCCTTCACTTCTCCTTTGGCTTTATTCAAgccgccgtt                                      <  2:208283/251‑4 (MQ=255)
                                                                                                                                                                                                  aaTCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGATTAATCTTTATCTGAAAATTTATCCTTGGCTTTCCCGATATCCTTTTGGATTTCGCCTTTTGCCTTATCTTTTTTGCCTTCAGCCTGCATGTCCGTTCTGTCAGCCATATCTCCGACTTTATCCTTCACTTCTCCTTTGGCTTTATTGAAGCCGCCCTTc                                     >  2:44056/1‑251 (MQ=255)
                                                                                                                                                                                                                 ttGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGATTAATCTTTATCTGAAAATTTATCCTTGGCTTTCCCGATATCCTTTTGGATTTCGCCTTTTGCCTTATCTTTTTTGCCTTCAGCCTGCATGTCCGTTCTGTCAGCCATATCTCCGACTTTATCCTTCACTTCTCCTTTGGCTTTATTGAAGCCGCCCTTCATTTTGTCTTTTac                       >  2:120013/1‑250 (MQ=255)
                                                                                                                                                                                                                        atttattCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGATTAATCTTTATCTGAAAATTTATCCTTGGCTTTCCCGATATCCTTTTGGATTTCGCCTTTTGCCTTATCTTTTTTGCCTTCAGCCTGCATGTCCGTTCTGTCAGCCATATCTCCGACTTTATCCTTCACTTCTCCTTTGGCTTTATTGAAGCCGCCCTTCATTTTGTCTTTTACACTATCGt               >  1:330825/1‑251 (MQ=255)
                                                                                                                                                                                                                                    ctgAAAAAGCGAAAAAAAAATACAAGAAACATTTCTTTTCCTGCCTCTTCCGTTCGTTGAATCTTTATCTGAAAATTTTTCCTTGTCTTTCCCGATATCCTTTTGGTTTTCGCCTTGTGCCTTATCTTTTTTGCCTGCAGCCTGCATGTCCGTTCTGTCAACCATATCTCCTACTTTATCCTCCTCTTCTCCTTTGGCTTTATTTAAGCCGCCCTTCATTTTGTCTTTTACACTATCGTTACCCATGtt     <  2:129058/247‑1 (MQ=255)
                                                                                                                                                                                                                                    cAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGATTAATCTTTATCTGAAAATTTATCCTTGGCTTTCCCGATATCCTTTTGGATTTCGCCTTTTGCCTTATCTTTTTTGCCTTCAGCCTGCATGTCCGTTCTGTCAGCCATATCTCCGACTTTATCCTTCACTTCTCCTTTGGCTTTATTGAAGCCGCCCTTCATTTTGTCTTTTACACTATCGTTACCCATGTTa    >  2:114377/1‑250 (MQ=255)
                                                                                                                                                                                                                                     aGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGATTAATCTTTATCTGAAAATTTATCCTTGGCTTTCCCGATATCCTTTTGGATTTCGCCTTTTGCCTTATCTTTTTTGCCTTCAGCCTGCATGTCCGTTCTGTCAGCCATATCTCCGACTTTATCCTTCACTTCTCCTTTGGCTTTATTGAAGCCGCCCTTCATTTTGTCTTTTACACTATCGTTACCCATGTTAAt  <  1:259703/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                        
TTGTGTCGATATCAACTGTTTTCTCGCACTCCCGTGGCTTCCGACACTTTCTTCGGTGCAGTACATTCCGTCATCAAAACCCCTTTGTTTGATTTCTTCAAATGATGTTATAAAAGTTTACAACTTATGTCATCATTATATAGAAGAGGCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGATTAATCTTTATCTGAAAATTTATCCTTGGCTTTCCCGATATCCTTTTGGATTTCGCCTTTTGCCTTATCTTTTTTGCCTTCAGCCTGCATGTCCGTTCTGTCAGCCATATCTCCGACTTTATCCTTCACTTCTCCTTTGGCTTTATTGAAGCCGCCCTTCATTTTGTCTTTTACACTATCGTTACCCATGTTAAT  >  NC_000964/3769820‑3770298

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