Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,329,083 T→C T239A (ACG→GCG)  yfaQ ← tandem DUF2300 domain‑containing protein YfaQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,329,0830TC95.7% 70.2 / ‑6.9 23T239A (ACG→GCG) yfaQtandem DUF2300 domain‑containing protein YfaQ
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (7/15);  total (7/16)
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

GCAATTTCAGTTTGCCCTTTTCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGATAGCGCCCATTAAGTACGCCGGGATTTACTGCCGTCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACAC‑ATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTGAGCTTTGTTTTGCGCGGGTAACGTCGCTAATGAATC  >  NC_000913/2328846‑2329329
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                      
gCAATTTCAGTTTGCCCTTTTCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGATAGCGCCCATTAAGTACGCCGGGATTTACTGCCGTCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACAC‑ATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGCCGCATAACggg                                                                                                                                                                                                                                             <  1:250855/249‑1 (MQ=255)
 cAATTTCAGTTTGCCCTTTTCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCAGTAGTAAACGTCACACGATAGCGCCCATTAAGTGCGCCGGGATTTACTGCCGTCGTACTCTTTTCCGCTGTGATTTTCTTTAGCGGATAGCGGGCAAACAGCTCCACCTCTACAC‑ATTGCCCACTGGCAACCGGCGTCGGCACGGGCAAAACCCGATTAAGCACCGCCGCATAACGGGTTa                                                                                                                                                                                                                                          <  1:184522/251‑1 (MQ=255)
       cAGTTTGCCCTTTTCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGATAGCGCCCATTAAGTACGCCGGGATTTACTGCCGTCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACAC‑ATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGCCGCATAACGGGTTAGTACCg                                                                                                                                                                                                                                    <  1:40802/251‑1 (MQ=255)
                          aGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGATAGCGCCCATTAAGTACGCCGGGATTTACTGCCGTCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCCCCTCTACAC‑ATTGCCCCCTGGCAACCGGCGTCGGCACTGGCAAACCCCGATTAAGCACCGCCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGcc                                                                                                                                                                                                                 <  2:68145/251‑1 (MQ=255)
                          aGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGATAGCGCCCATTAAGTACGCCGGGATTTACTGCCGTCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACAC‑ATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGCCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGcc                                                                                                                                                                                                                 >  1:274948/1‑251 (MQ=255)
                                  cgtcTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGATAGCGCCCATTAAGTACGCCGGGATTTACTGCCGTCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACAC‑ATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGCCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAAccc                                                                                                                                                                                                         <  1:313193/251‑1 (MQ=255)
                                       cACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGATAGCGCCCATTAAGTACGCCGGGATTTACTGCCGTCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACAC‑ATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGCCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAg                                                                                                                                                                                                     <  1:131118/250‑1 (MQ=255)
                                        aCCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGATAGCGCCCATTAAGTACGCCGGGATTTACTGCCGTCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACAC‑ATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGCCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCCCAGCg                                                                                                                                                                                                   >  1:106232/1‑251 (MQ=255)
                                                 gACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGATAGCGCCCATTAAGTACGCCGGGATTTACTGCCGTCGTTCTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCTCCTCTACCC‑ATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGCCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCAt                                                                                                                                                                                          <  1:231707/251‑1 (MQ=255)
                                                                   ccctttAGTAACCGTCACACGATATCGCCCATTAAGTACGCCGGGATTTACTGCCGTCGTACTCTTTCCCGCTGTGATTTTCTTCCGCGGATAGCGGGCAAACAGCTCCACCTCTACAC‑ATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGCCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCaa                                                                                                                                                                        <  2:221800/247‑1 (MQ=255)
                                                                          gTAAACGTCACCCGATAGCGCCCATTAAGTACGCCGGGATTTACTGCCGCCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAACAAGCTCCACCTCTACAC‑ATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGCCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAaccac                                                                                                                                                                  <  2:226836/250‑1 (MQ=255)
                                                                            aaaCGGCACACGATAGCGCCCATTAAGTACGCCGGGATTTACTGCCGTCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACAC‑ATTGCCCCCTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGCCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCCCAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCAccc                                                                                                                                                                <  1:306252/250‑1 (MQ=255)
                                                                               cGTCACACGATAGCGCCCATTAAGTACGCCGGGATTTACTGCCGTCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACAC‑ATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGAATAAGCACCGCCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCAAAGCGGTGTACCATCGGGTAACCAGCCACCACAAACACCCTGAc                                                                                                                                                            >  2:169439/1‑251 (MQ=255)
                                                                                        aTAGCGCTCATTAAGTACGCCGCGATTTACTGCCGTCGTACCCTTTCCCGCTGTGATTTTCTTCGGCGGATATCGGGCAAACAGCTCCACCTCCCCACTATT‑CCCCCTGGCAACCTGCGTCTGCACTGGAAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTTCCATTGCCAGTAACCCACAGCGGTGTACCATCGGTTAAGCAGCCAGCAAAACCACCCTGACGGAttttct                                                                                                                                                   <  2:52791/251‑4 (MQ=255)
                                                                                                 aTTAAGTACGCCGGGATTTACTGCCGTCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTCCAC‑ATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCAACGCCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACAGCCAGTAACCCAAAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACAACCCTGAAGGATTTCCTGTTTatcat                                                                                                                                          >  1:301648/1‑251 (MQ=255)
                                                                                                              ggATTTACTGCCGTCTGAATCTTTTCCGCGTTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACAC‑ATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAGCCCGATTAAGCACCGCCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCGCCCTGACGGAGTTCCTGTTTATCATCGGCAAACCAGCt                                                                                                                             <  2:131624/251‑1 (MQ=255)
                                                                                                                                                           gaaTAGCGGGCAAACAGCTCCCCCTCTACCC‑AGTTCCCACTGGCACCCTGCGTCGGCATTGGCAAAACCCGATTAAGCACCGCCGCATAACGGGTTAGGACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCTGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGAACCCGACTACCGAGCATGGATGCGACAcc                                                                                <  2:18130/249‑1 (MQ=255)
                                                                                                                                                                                        cc‑aTTGCCCACTGGCAACCGGCGGCGGCACTGGCAAAACCCGATTAAGCACCGCCTCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGGAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGGACCCGACTACCTAGCATGGATGCGACACCGATTTTCGCCGCGTCCAGCACCACATCCa                                                   <  2:130283/250‑1 (MQ=255)
                                                                                                                                                                                            ttGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGCCGCATAACGGGTTAGAACCGTTTTGCTGGTGACACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACAACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGATTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGt                                                 >  1:302636/1‑250 (MQ=255)
                                                                                                                                                                                                                        ggCAAAACCCGATTAAGCACCGCCGCATAACGGGTTAGTACCGTTTGGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCTTCGGTTAACCAGCCAGCAAAACCTCCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCTACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTGATCTTTGTTTTGCGCgg                     <  2:274948/250‑1 (MQ=255)
                                                                                                                                                                                                                                      aaGCCCCGCCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGATGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTGTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTGAGCTTTGTTTTGCGCGGGTAACGTCGCTAATg      <  1:315007/251‑1 (MQ=255)
                                                                                                                                                                                                                                       aGCACCGCCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTGAGCTTTGTTTTGCGCGGGTAACGTCGCTAATg      >  1:119057/1‑250 (MQ=255)
                                                                                                                                                                                                                                          aCCGCCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTGAGCTTTGTTTTGCGCGGGTAACGTCGCTAATGAATc  >  1:154606/1‑251 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                      
GCAATTTCAGTTTGCCCTTTTCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGATAGCGCCCATTAAGTACGCCGGGATTTACTGCCGTCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACAC‑ATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTGAGCTTTGTTTTGCGCGGGTAACGTCGCTAATGAATC  >  NC_000913/2328846‑2329329

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