Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 476,337 C→T G12D (GGC→GAC)  ybaZ ← DNA base‑flipping protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913476,3370CT94.7% 55.6 / NA 19G12D (GGC→GAC) ybaZDNA base‑flipping protein
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (11/7);  total (12/7)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGATTAACCACCCGGTGCCAGGGTAAGGTGCTGCCTTCAGGGAGACGCTTTAACACACCGCCCACCTGGCGCGCGGCGCGGGGCGATCCCGCCAGTTTCGCCACATCACCGTAAGTGGTGACATAGCCTTCGGGAATAGCGGCGACGATTTGCCAGACGCGTTGGGGAAATGAATCTTCTTTTTCCATCTTTTCTTCCTGAGGTAATTTTTCAGCATAATCTGGAAAAACGCCCGAGTGAAGTCGCATTGCGCAAGAAACCAGCATCTGGCACGCGATGGGTTGCAATTAGCCGGGGCAGCAGTGATAATGCGCCTGCGCGTTGGTTCTCAACGCTCTCAATGGGGGCTCTGTTGGTTCTCCCGCAACGCTACTCTGTTTACCAGGTCAGGTCCGGAAGGAAGCAGCCAAGGCAGATGACGCGTGTGCCGGGATGTAGCTGGCAGGGCCCCCACCCATTTCTGCCTCCCACCGT  >  NC_000913/476106‑476579
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  
cGATTAACCACCCGGTGCCAGGGTAAGGTGCTGCCTTCAGGGAGACGCTTTAACACACCGCCCACCTGGCGCGCGGCGCGGGGCGATCCCGCCAGTTTCGCCACATCACCGTAAGTGGTGACATAGCCTTCGGGAATAGCGGCGACGATTTGCCAGACGCGTTGGGGAAATGAATCTTCTTTTTCCATCTTTTCTTCCTGAGGTAATTTTTCAGCATAATCTGGAAAAACGTCCGAGTGAAGTCGCATTgc                                                                                                                                                                                                                                 <  1:49874/251‑1 (MQ=255)
                gCCAGGGTAAGGTGCTGCCTGCAGGGAGACGCTTTAACACACCGCCCACCTGGCGCGCGGCGCGGGGCGATCCCGCCAGTTTCGCCACATCACCGTAAGTGGTGACATAGCCTTCGGGAATAGCGGCGACGATTTGCCAGACGCGTTGGGGAAATGAATCTTCTTTTTCCATCTTTTCTTCCTGAGGTAATTTTTCAGCATAATCTGGAAAAACGTCCGAGTGAAGTCGCATTGCGCAAGAAACCAGCATc                                                                                                                                                                                                                 <  2:84622/251‑1 (MQ=255)
                                     cAGGGAGACGCTGTAACACACCGCCCACCTGGCGCGCGGCGCGGGGCGATCCCGCCAGTTTCGCCACATCCCCGTAAGTGGTGACATAGCCTTCGGGAATAGCGGCGACGATTTGCCAGACGCGTTGGGGAAATGAATCTTCTTTTTCCATCCTTTCTTCCTGAGGTAATTTTTCAGCATAATCTGGAAAAACGTCCGAGTGAAGTCGCATTGCGCAAGAAACCAGCATCTGGCACGCGATGGGTTGCAAt                                                                                                                                                                                            >  1:114909/1‑251 (MQ=255)
                                                          cGCCCACCTAGCGCGCGGCGCGGGGCGATCCCGCCAGTTTCGCCACATCACCGTAAGCGGTGACATAGCCTTCGGGAATAGCGGCGACGATTTGCCAGACGCGTAGGGGAAATGAATCTTCTTTTTCAATCTTTTCTTCATGAGGTAATTTTTCAGCATAATCTGGAAAAACGCCCGAGTGAAGTCGCATTGCGCAAGAAACAAGCATCTGGCACGCGATGGGTTGCAATTACACGGGGCAGCAGTGATa                                                                                                                                                                        >  2:20088/1‑250 (MQ=255)
                                                           cccctccgggcgcgtttgctggggCGTTCCCGCCAGTTTCGCCTCATCACCGTAAGTGGTGACATAGCCTTCGGGAATAGCGGCGACGATTTGCCAGACGCGTTGGGGAAATGAATCTTCTTGTTCCATCTTTTCTTCCTGAGGTAATTTTTCATCATAATCTGGAAAAACGTCCGAGTGAAGTCGCATTGCGCAAGAAACCAGCATCTGGCACGCGATTTGTTGCAATTATGCGGTGCAGCAGTGATAAt                                                                                                                                                                      <  2:162417/231‑1 (MQ=255)
                                                                          ggcgcggGGCGATCCCGCCAGTTTCGCCACATCACCGTAAGTGGTGACATAGCCTTCGGGAATAGCGGCGACGATTTGCCAGACGCGTTGGGGAAATGAATCTTCTTTTTCCATCTTTTCTTCCTGAGGTAATTTTTCAGCATAATCTGGAAAAACGTCCGAGTGAAGTCGCATTGCGCAAGAAACCAGCATCTGGCACGCGATGGGTTGCAATTAGCCGGGGCAGCAGTGATAATGCGCCTGCGCGTTgg                                                                                                                                                       >  1:227130/1‑251 (MQ=255)
                                                                                gggCGATCCCGCCAGTTTCGCCACATCACCGTAAGTGGTGACATAGCCTTCGGGAATAGCGGCGACGATTTGCCAGACGCGTTGGGGAAATGAATCTTCTTTTTCCATCTTTTCTTCCTGAGGTAATTTTTCAGCATAATCTGGAAAAACGTCCGAGTGAAGTCGCATTGCGCAAGAAACCAGCATCTGGCACGCGATGGGTTGCAATTAGCCGGGGCAGCAGTGATAATGCGCCTGCGCGTTGGTtctc                                                                                                                                                  >  2:86953/1‑250 (MQ=255)
                                                                                    gATCCCGCCAGTTTCGCCACATCACCGTAAGTGGTGACATAGCCTTCGGGAATAGCGGCGACGATTTGCCAGACGCGTTGGGGAAATGAATCTTCTTTTTCCATCTTTTCTTCCTGAGGTAATTTTTCAGCATAATCTGGAAAAACGTCCGAGTGAAGTCGCATTGCGCAAGAAACCAGCATCTGGCACGCGATGGGTTGCAATTAGCCGGGGCAGCAGTGATAATGCGCCTGCGCGTTGGTTCTCAACGc                                                                                                                                             >  1:181130/1‑251 (MQ=255)
                                                                                                                  gtggagACATAGCCTTCGGGAATAGCGGCGACGATTTGCCAGACGCGTTGGGGAAATGAATCTTCTTTTTCCATCTTTTCTTCCTGAGGTAATTTTTCAGCATAATCTGGAAAAAAGTCCGAGTGAAGTCGCATTGCGCAAGAAACCAGCATCTGGCACGCGATGGGTTGCAATTAGCCGGGGCAGCAGTGATAATGCGCCTGCGCGTTGGTTCTCAACGCTCTCAATGGGGGCTCTGTTGGTTCTCCCGc                                                                                                               >  1:220168/1‑251 (MQ=255)
                                                                                                                                  cGGGAATAGCGGCGACGATTTGCCAGACGCGTTGGGGAAATGAATCTTCTTTTTCCATCTTTTCTTCCTGAGGTAATTTTTCAGCATAATCTGGAAAAACGTCCGAGTGAAGTCGCATTGCGCAAGAAACCAGCATCTGGCACGCGATGGGTTGCAATTAGCCGGGGCAGCAGTGATAATGCGCCTGCGCGTTGGTTCTCAACGCTCTCAATGGGGGCTCTGTTGGTTCTCCCGCACCGCTACTCTGttt                                                                                                >  2:38635/1‑250 (MQ=255)
                                                                                                                                                       gCCAGACGCGTTGGGGAAATGAATCTTCTTTTTCCATCTTTTCTTCCTGAGGTAATTTTTCAGCATAATCTGGAAAAACGTCCGAGAGAAGTCGCATTGCGCAAGAAACCAGCATCTGGCACGCGATGGGTTGCAATTAGCCGGGGCAGCAGTGATAATGCGCCTGCGCGTTGGTTCTCAACGCTCTCAATGGGGGCTCTGTTGGTTCTCCCGCAACGCTACTCTGTTTACCAGGTCAGGTCCggaagga                                                                           >  1:248909/1‑250 (MQ=255)
                                                                                                                                                                      gAAATGAATCTTCTTTTTCCATCTTTTCTTCCTGAGGTAATTTTTCAGCATAATCTGGAAAAACGTCCGAGTGAAGTCGCATTGCGCAAGAAACCAGCATCTGGCACGCGATGGGTTGCAATTAGCCGGGGCAGCAGTGATAATGCGCCTGCGCGTTGGTTCTCAACGCTCTCAATGGGGGCTCTGTTGGTTCTCCCGCAACGCTACTCTGTTTACCAGGTCAGGTCCGGAAGGAAGCAGCCAAGGc                                                               >  2:78043/1‑247 (MQ=255)
                                                                                                                                                                            aaTCTTCTTTTTCCATCTTTTCTTCCTGAGGTAATTTTTCAGCATAATCTGGAAAAACGTCCGAGTGAAGTCGCATTGCGCAAGAAACCAGCATCTGGCACGCGATGGGTTGCAATTAGCCGGGGCAGCAGTGATAATGCGCCTGCGCGTTGGTTCTCAACGCTCTCAATGGGGGCTCTGTTGGTTCTCCCGCAACGCTACTCTGTTTACCAGGTCAGGTCCGGAAGGAAGCAGCCAAGGCAGATGAcgcg                                                     >  1:249705/1‑251 (MQ=255)
                                                                                                                                                                                               ttcttcCTGAGGTAATTTTTCAGCATAATCTGGAAAAACGTCCGCGTGAAGTCGCATTGCGCAAGAAACCAGCATCTGGCACGCGATGGGTTGCAATTAGCCGGGGCAGCAGTGATAATGCGCCTGCGCGTTGGTTCTAACCGCTCTCAATGGGGGCTCTGTTGGTTCTCCCGCAACGCTACTCTGTTTACCAGGTCAGGTCCGGAAGGAAGCAGCCAAGGCAGATGACGCGTGTGCCGGGATGTAGCTgg                                  >  1:32894/1‑251 (MQ=255)
                                                                                                                                                                                                tcttcCTGAGGTAAGTTTTCAGCATAATTAGGCAAAACGTCCTTGTGAAGGCGCAGTGGGCAAGAAACCAGCATCTGGCACGCGAGGGGTTTCAATTAGCCGGGGCAGCAGTGATAATGCGCCTGCGCGTTGGTTCTCCACGCTCTCAATGGGGGCTCTGTTGGTTCTCGCGCAACGCTACTCTGTTTACCAGGTCAGGTCCGGAATGAAGCAGCCAAGGCAGATGACGCGTGTGCCGGGATGTAGCttt                                  <  2:140202/250‑3 (MQ=255)
                                                                                                                                                                                                                 ttCAGCATAAGCTGGAAAAACGTCCGAGTGAAGTCGCATTGCGCAAGAAACCAGCATCTGGCACGCGATGGGTTGCAATGAGCCGGGGCAGCAGTGATAATGCGCCTGCGCGTTGGTTCTCAACGCTCTGAATGGGGGCTCTGTTGGTTCTCCCGCAACGCTACTCTGTTTACCAGGTCAGGTCCGGAAGGAAGCAGCCAAGGCAGATGACGCGTGTGCCGGGATGTAGCTGGCAGGGCCCCCACCCAtt                 <  2:249705/250‑1 (MQ=255)
                                                                                                                                                                                                                    aGCATAATCTGGAAAAACGTCCGAGTTAAGTCGCATTGCGCAAGAAACCAGCATCTGGCACTCGATGGGTTGCAATTAGCCGGGGCAGCAGTGATAATGCGCCTGCGCGTTGGTTCTCAACGCTCTCAATGGGGGCTCTGTTGGTTCTCCCGCAACGCTACTCTGTTTACCAGGTCAGGTCCGGAAGGAAGCAGCCAAGGCAGATGACGCGTGTGCCGGGATGTAGCTGGCAGGGCCCCCACCCATTTCTg             <  1:214592/251‑1 (MQ=255)
                                                                                                                                                                                                                             tGGAAAAACGTCCGAGTGAAGTCGCATTGCGCAAGAAACCAGCATCTGGCACGCGATGGGTTGCAATTAGCCGGGGCAGCAGTGATAATGCGCCTGCGCGTTGGTTCTCAACGCTCTCAATGGGGGCTCTGTTGGTTCTCCCGCAACGCTACTCTGTTTACCAGGTCAGGTCCGGAAGGAAGCAGCCAAGGCAGATGACGCGTGTGCCGGGATGTAGCTGGCAGGGCCCCCACCCATTTCTGCCTCCCAcc    >  1:58151/1‑251 (MQ=255)
                                                                                                                                                                                                                               gAAAAACGTCCGAGTGAAGTCGCATTGCGCAAGAAACCAGCATCTGGCACGCGATGGGTTGCAATTAGCCGGGGCAGCAGTGATAATGCGCCTGCGCGTTGGTTCTCAACGCTCTCAATGGGGGCTCTGTTGGTTCTCCCGCAACGCTACTCTGTTTACCAGGTCAGGTCCGGAAGGAAGCAGCCAAGGCAGATGACGCGTGTGCCGGGATGTAGCTGGCAGGGCCCCCACCCATTTCTGCCTCCCACCGt  <  2:135979/251‑1 (MQ=255)
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  
CGATTAACCACCCGGTGCCAGGGTAAGGTGCTGCCTTCAGGGAGACGCTTTAACACACCGCCCACCTGGCGCGCGGCGCGGGGCGATCCCGCCAGTTTCGCCACATCACCGTAAGTGGTGACATAGCCTTCGGGAATAGCGGCGACGATTTGCCAGACGCGTTGGGGAAATGAATCTTCTTTTTCCATCTTTTCTTCCTGAGGTAATTTTTCAGCATAATCTGGAAAAACGCCCGAGTGAAGTCGCATTGCGCAAGAAACCAGCATCTGGCACGCGATGGGTTGCAATTAGCCGGGGCAGCAGTGATAATGCGCCTGCGCGTTGGTTCTCAACGCTCTCAATGGGGGCTCTGTTGGTTCTCCCGCAACGCTACTCTGTTTACCAGGTCAGGTCCGGAAGGAAGCAGCCAAGGCAGATGACGCGTGTGCCGGGATGTAGCTGGCAGGGCCCCCACCCATTTCTGCCTCCCACCGT  >  NC_000913/476106‑476579

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