Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,026,760 C→T A420T (GCT→ACT)  dgcQ ← diguanylate cyclase DgcQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,026,7600CT90.0% 21.2 / NA 10A420T (GCT→ACT) dgcQdiguanylate cyclase DgcQ
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (4/5);  minor base G (0/1);  total (4/6)
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

CTCCGCAGCCTCCGTCAGACTCGCGCCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTAAGGTGTCGTGCCACGCCTGCCACTGCAACGAGCTTTGCAGAACATACATGTTGCTAACCATCCGGCGAATCACATACCAGGAGATG  >  NC_000913/2026533‑2026896
                                                                                                                                                                                                                                   |                                                                                                                                        
cTCCGCAGCCTCCGTCAGACTCGCGCCTGGCAGACTCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGATAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGTGAGCGGACGGGCTTTTTCGAAc                                                                                                                    >  1:114890/1‑250 (MQ=255)
           ccGTCAGACTCGCGCCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAATGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAAGGGTCAAGAGCGACCTGAATGACAGAAAAAGGATGTTGTTGCGTCTGACACAATTTAGTGAGCGGACGGGCTTTTTCGAACAGTGCGCCACgg                                                                                                        <  1:41645/251‑1 (MQ=255)
                                                    ccaccGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACGCGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATGGCTTTAAAAGGGTCAAGAGCGACCTGAATGACAGAAAAAGGATGTTGGTGCGGCTCACACAGTTTAGGGGGCGGAGGGGCGTTGGCGAACAGGGCGCCACGGTTATATAAAGGCGTGAAGGTGTCGTGCCACGCCTGCCACTg                                                               <  1:127539/251‑1 (MQ=255)
                                                          ncccgcccgGCAACGTCCTGCGCACGCAAGGAACTACTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGTGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTAAGGTGTCGTGCCACGCCTGCCACTGCAACg                                                          >  1:66056/2‑250 (MQ=255)
                                                                           cNGCNNACGCNATGAACTGCTAATTAATCCGGCAGNATGAGAAAGAACACGGTCGCCCGCCTGATGCNCNAAGCGGTCATTAATCGCTTTAAAAGGGTCAAGATCGACCTGAATGACANAAAAAGGATGNTGGNNCGTCNGNCNCAATTTAGTGAGCGGACGGGCTTTTTCGANCAGTGCGCCACGGTTATATAAACNNGTNNNGGNNTCGTg                                                                              <  2:240144/213‑1 (MQ=255)
                                                                                          aCTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGTGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTAAGGTGTCGTGCCACGCCTGCCACTGCAACGAGCTTTGCAGAACATACATGTTGCTAACCATcc                         <  1:11967/251‑1 (MQ=255)
                                                                                                       nnnnnngcnNNNNAAANAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGTGAGCGGACGGGCTTTTTCGAACAGt                                                                                                                 <  2:204907/144‑1 (MQ=255)
                                                                                                                                                            aaTCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGTGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTAAGGTGTCGTGCCACGCCTGCCACTGCAACGAGCTTTGCAGAACATACATGTTGCTAACCATCCGGCGAATCACATACCAGGAGATg  <  1:186049/208‑1 (MQ=255)
                                                                                                                                                                     aaaaTGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTCGTGCGTCTGACACCATTTAGTGAGCGGACGGGCTTTTTCGAACAGTGCGACACGGTACTATAAACGCGTTAAGGTGTAGTGCCACGc                                                                        >  1:283431/1‑129 (MQ=255)
                                                                                                                                                                                gATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACCCAATTTAGTGAGCGGACGGGCTTTTTc                                                                                                                        >  1:69162/1‑70 (MQ=255)
                                                                                                                                                                                                                                   |                                                                                                                                        
CTCCGCAGCCTCCGTCAGACTCGCGCCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTAAGGTGTCGTGCCACGCCTGCCACTGCAACGAGCTTTGCAGAACATACATGTTGCTAACCATCCGGCGAATCACATACCAGGAGATG  >  NC_000913/2026533‑2026896

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