Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,473,265 G→A intergenic (‑56/+75) yjgH ← / ← bdcA RutC family protein YjgH/c‑di‑GMP‑binding biofilm dispersal mediator protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,473,2650GA90.9% 22.6 / ‑5.9 12intergenic (‑56/+75)yjgH/bdcARutC family protein YjgH/c‑di‑GMP‑binding biofilm dispersal mediator protein
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (6/4);  total (6/6)
Fisher's exact test for biased strand distribution p-value = 4.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGCAAATTATCAAATGCCAGTCTGACTTGTTGCTGAAAATCGGGTTCTGGTGTTCCGTCCTCTCGACTTCCTACTTGCCCGGAAACAAACAGCAAATCGCCGGAACGAATAGCCGCAGAATAACGATGCTCAGCATATAGTGAATGTCGGCCAGCAGGGAAAACAGCGGTTCTTTCTACCATTTGGTTATCCTCAAGATTTACGACATGAACAGAAGATTTCTCTTTACCGGGAGCCGCTTTTAGCGGACGACGTGAGTAAACAAAACCCAGACATCATGGATAATGGCTGGGCTTAATTGAGCGTAGTCGGTTATGCGCCAAACGCGCCATCAATGGTATGCATCGCGCCGGTAACAAAACTGGCTTCTGGCCCTGCTAACCATGCGACCATACCAGCGACCTCTTCCGGTTGCCCATGTCTTTTGATAGCCATCAAACTATGCAACATATCGCGCATTGGCC  >  NC_000913/4473027‑4473491
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                  
aTGCAAATTATCAAATGCCAGTCTGACTTGTTGCTGAAAATCGGGTTCTGGTGTTCCGTCCTCTCGACTTCCTACTTGCCCGGAAACAAACAGCAAATCGCCGGAACGAATAGCCGCAGAATAACGATGCTCAGCATATAGTGAATGTCGGCCAGCAGGGAAAACAGCGGTTCTTTCTACCATTTGGTTATCCACAAGATTTACGACATGAACAGAAGATTTCTCTTTACCGGGAGCCACTTTTAGCGgac                                                                                                                                                                                                                        <  1:234046/251‑1 (MQ=255)
            aaaTGCCAGTCTGACTTGTTGCTGAAAATCGGGTTCTGGTGTTCCGTCCTCTCGACTTCCTACTTGCCCGGAAACAAACAGCAAATCGCCGGAACGAATAGCCGCAGAATAACGATGCTCAGCATATAGTGAATGTCGGCCAGCAGGGAAAACAGCGGTTCTTTCTACCATTTGGTTATCCTCAAGATTTACGACATGAACAGAAGATTTCTCTTTACCGGGAGCCACTTTTAGCgn                                                                                                                                                                                                                          <  1:158380/237‑2 (MQ=255)
                                                     ttCCGNNTTNNNGANNACCTACTTGCCCGGAAACAAACANCAAATCGCCGGGACGAATAGCAGCAGAATNANGNTGCCNNGAANATAGTGAATGNCGGCCCACAGGGAAAACAGCGGGTCTTTCTACCATTTGGTTATCCTCANGNTTTACGACCTGAGCAGAAGATTTCTCTTTACNGGGAGCCACTTTTAGCGGCCGAAGNGAGTnnaanannnnnnnnn                                                                                                                                                                                                >  2:183484/1‑210 (MQ=255)
                                                                       ntACTTGCCCGGAAACAAACAGCAAATCGCCGGAACGAATAGCCGCAGAATAACGATGCTCAGCATATAGTGAATGTCGGCCAGCCGGGAAAACAGCGGTTCTTTCTACCATTTGGTTATCCTCAAGATTTACGACATGAACAGAAGATTTCTCTTTACCGGGAGCCACTTTTAGCGGACGACGTGAGTAAACAAAACCCAGACATCATGGATAATGGCTGGGCTTAATTGAGCGTAGTCGGTTATGCGcc                                                                                                                                                 >  1:66614/2‑251 (MQ=255)
                                                                         acttacccgtcaacaaacaGCAAATCGCCGGAACGAATAGCCGCAGAATAACGATGCTCATCATATAGTGAATGTCGGCCAGCAGGGAAAACAGCGGTTCTTTCTATCATTTTGTTATCCTCAAGATTTACGACATGAACAGAAGATGTCTCTTTACCGGGAGCCGCTTTTAGCGCACCACGTCAGTAATCAAATCCTAAACCTCATGGATAATGGCTGGGCTTAATTGAGCGTAGTCGGTTATGCGcc                                                                                                                                                 <  1:183484/238‑1 (MQ=255)
                                                                               ccGGAAACAAACAGCAAATCGCCGGAACGAATAGCCGCAGAATAACGATGCTCAGCATATAGTGAATGTCGGCCAGCAGGGAAAACAGCGGTTCTTTCTACCATTTGGTTATCCTCAAGATTTACGACATGAACAGAAGATTTCTCTTTACCGTGAGCCTCTTTTAGCGGACGACGTGAGTAAACAAAACCCAGACATCATGGATAATGGCTGGGCTTAATTGAGCGTAGTCGGTTATGCGCCAAACgcgc                                                                                                                                         <  1:147653/251‑1 (MQ=255)
                                                                               ccGGAAACAAACAGCAAATCGCCGGAACGAATAGCCGCAGAATAACGATGCTCAGCATATAGTGAATGTCGGCCAGCAGGGAAAACAGCGGTTCTTTCTACCATTTGGTTATCCTCAAGATTTACGACATGAACAGAAGATTTCTCTTTACCGGGAGCCACTTTTAGCGGACGACGTGAGTAAACAAAACCCAGACATCATGGATAATGGCTGGGCTTAATTGAGCGTAGTCGGTTATGCGCCAAACgcgc                                                                                                                                         >  1:124455/1‑251 (MQ=255)
                                                                                cGGAAACAAACAGCAAATCGCCGGAACGAATAGCCGCAGAATAACGATGCTCAGCATATCGTGAATGTCGGCCAGCAGGGAAAACAGCGGTTCTTTCTACCATTTGGTTATCCTCAAGATTTACGACATGAACAGAAGATTTCTCTTTACCGGGAGCCACTTTTAGCGGACGACGTGAGTAAACAAAACCCAGACATCATGGATAATGGCTGGGCTTAATTGAGCGTAGTCGGTTATGCGCCAAACGCGcc                                                                                                                                        >  1:163874/1‑251 (MQ=255)
                                                                                          caGCAAATCGCCGGAACGAATAGCCGCAGAATAACGATGCTCAGCATATAGTGAATGTCGGCCAGCAGGGAAAACAGCGGTTCTTTCTACCATTTGGTTATCCTCAAGATTTACGCCATGAACAGAAGATTTCTCTTTACCGGGAGCCACTTTTAGCGGACGACGTGAGTAAACAAAACCCAGACATCATGGATAATGGCTGGGCTTAATTGAGCGTAGTCGGTTATGCGCCAAACgcgc                                                                                                                                         >  1:173065/1‑240 (MQ=255)
                                                                                                                                                                                              tCCTCNNGANNNACNNCATGAACAGAAGATTTCTCTTTANCGGGAGCCACTTTTAGCGGACGACGTGAGNANANAAAANNCAGNCATCATGGATNATGGCTGGGCTTAATTGAg                                                                                                                                                                   >  2:189903/1‑114 (MQ=255)
                                                                                                                                                                                              tCCTCAAGATTTACGACATGAACAGAAGATTTCTCTTTACCGGGAGCCACTTTTAGCGGACGACGTGAGTAAACAAAACCCAGACATCATGGATAATGGCTGGGCTTAATTGAg                                                                                                                                                                   <  1:189903/114‑1 (MQ=255)
                                                                                                                                                                                                                      gaagaTTTCTCTTTACCGGGAGCCACTTTTAGCGGACGACGTGAGTAAACAAAACCCAGACATCATGGATAATGGCTGGGCTTAATTGAGCGTAGTCGGTTATGCGCCAAACGCGCCATCAATGGTATGCATCGCGCCGGTAACAAAACTGGCTTCTGGCCCTGCTAACCATGCGACCATACCAGCGACCTCTTCCGGTTGCCCATGTCTTTTGATAGCCATCAAACTATGCAACATATCGCGCATTGGcc  <  1:125775/251‑1 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                  
ATGCAAATTATCAAATGCCAGTCTGACTTGTTGCTGAAAATCGGGTTCTGGTGTTCCGTCCTCTCGACTTCCTACTTGCCCGGAAACAAACAGCAAATCGCCGGAACGAATAGCCGCAGAATAACGATGCTCAGCATATAGTGAATGTCGGCCAGCAGGGAAAACAGCGGTTCTTTCTACCATTTGGTTATCCTCAAGATTTACGACATGAACAGAAGATTTCTCTTTACCGGGAGCCGCTTTTAGCGGACGACGTGAGTAAACAAAACCCAGACATCATGGATAATGGCTGGGCTTAATTGAGCGTAGTCGGTTATGCGCCAAACGCGCCATCAATGGTATGCATCGCGCCGGTAACAAAACTGGCTTCTGGCCCTGCTAACCATGCGACCATACCAGCGACCTCTTCCGGTTGCCCATGTCTTTTGATAGCCATCAAACTATGCAACATATCGCGCATTGGCC  >  NC_000913/4473027‑4473491

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