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

GCAATTTCAGTTTGCCCTTTTCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGATAGCGCCCATTAAGTACGCCGGGATTTACTGCCGTCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTGAGCTTTGTTTTGCGCGGGTAACGTCGCTAATGAATC  >  NC_000913/2328846‑2329329
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                      
gCAATTTCAGTTTGCCCTTTTCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGATAGCGCCCATTAAGTACGCCGGGATTTACTGCCGTCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGCCGCATAACggg                                                                                                                                                                                                                                             <  1:250855/249‑1 (MQ=255)
 cAATTTCAGTTTGCCCTTTTCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCAGTAGTAAACGTCACACGATAGCGCCCATTAAGTGCGCCGGGATTTACTGCCGTCGTACTCTTTTCCGCTGTGATTTTCTTTAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCGGCGTCGGCACGGGCAAAACCCGATTAAGCACCGCCGCATAACGGGTTa                                                                                                                                                                                                                                          <  1:184522/251‑1 (MQ=255)
       cAGTTTGCCCTTTTCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGATAGCGCCCATTAAGTACGCCGGGATTTACTGCCGTCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGCCGCATAACGGGTTAGTACCg                                                                                                                                                                                                                                    <  1:40802/251‑1 (MQ=255)
                          aGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGATAGCGCCCATTAAGTACGCCGGGATTTACTGCCGTCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCCCCTCTACACATTGCCCCCTGGCAACCGGCGTCGGCACTGGCAAACCCCGATTAAGCACCGCCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGcc                                                                                                                                                                                                                 <  2:68145/251‑1 (MQ=255)
                          aGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGATAGCGCCCATTAAGTACGCCGGGATTTACTGCCGTCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGCCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGcc                                                                                                                                                                                                                 >  1:274948/1‑251 (MQ=255)
                                  cgtcTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGATAGCGCCCATTAAGTACGCCGGGATTTACTGCCGTCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGCCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAAccc                                                                                                                                                                                                         <  1:313193/251‑1 (MQ=255)
                                       cACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGATAGCGCCCATTAAGTACGCCGGGATTTACTGCCGTCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGCCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAg                                                                                                                                                                                                     <  1:131118/250‑1 (MQ=255)
                                        aCCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGATAGCGCCCATTAAGTACGCCGGGATTTACTGCCGTCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGCCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCCCAGCg                                                                                                                                                                                                   >  1:106232/1‑251 (MQ=255)
                                                 gACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGATAGCGCCCATTAAGTACGCCGGGATTTACTGCCGTCGTTCTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCTCCTCTACCCATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGCCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCAt                                                                                                                                                                                          <  1:231707/251‑1 (MQ=255)
                                                                   ccctttAGTAACCGTCACACGATATCGCCCATTAAGTACGCCGGGATTTACTGCCGTCGTACTCTTTCCCGCTGTGATTTTCTTCCGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGCCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCaa                                                                                                                                                                        <  2:221800/247‑1 (MQ=255)
                                                                          gTAAACGTCACCCGATAGCGCCCATTAAGTACGCCGGGATTTACTGCCGCCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAACAAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGCCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAaccac                                                                                                                                                                  <  2:226836/250‑1 (MQ=255)
                                                                            aaaCGGCACACGATAGCGCCCATTAAGTACGCCGGGATTTACTGCCGTCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCCCTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGCCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCCCAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCAccc                                                                                                                                                                <  1:306252/250‑1 (MQ=255)
                                                                               cGTCACACGATAGCGCCCATTAAGTACGCCGGGATTTACTGCCGTCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGAATAAGCACCGCCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCAAAGCGGTGTACCATCGGGTAACCAGCCACCACAAACACCCTGAc                                                                                                                                                            >  2:169439/1‑251 (MQ=255)
                                                                                        aTAGCGCTCATTAAGTACGCCGCGATTTACTGCCGTCGTACCCTTTCCCGCTGTGATTTTCTTCGGCGGATATCGGGCAAACAGCTCCACCTCCCCACTATTCCCCCTGGCAACCTGCGTCTGCACTGGAAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTTCCATTGCCAGTAACCCACAGCGGTGTACCATCGGTTAAGCAGCCAGCAAAACCACCCTGACGGAttttct                                                                                                                                                   <  2:52791/251‑4 (MQ=255)
                                                                                                 aTTAAGTACGCCGGGATTTACTGCCGTCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTCCACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCAACGCCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACAGCCAGTAACCCAAAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACAACCCTGAAGGATTTCCTGTTTatcat                                                                                                                                          >  1:301648/1‑251 (MQ=255)
                                                                                                              ggATTTACTGCCGTCTGAATCTTTTCCGCGTTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAGCCCGATTAAGCACCGCCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCGCCCTGACGGAGTTCCTGTTTATCATCGGCAAACCAGCt                                                                                                                             <  2:131624/251‑1 (MQ=255)
                                                                                                                                                           gaaTAGCGGGCAAACAGCTCCCCCTCTACCCAGTTCCCACTGGCACCCTGCGTCGGCATTGGCAAAACCCGATTAAGCACCGCCGCATAACGGGTTAGGACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCTGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGAACCCGACTACCGAGCATGGATGCGACAcc                                                                                <  2:18130/249‑1 (MQ=255)
                                                                                                                                                                                        ccaTTGCCCACTGGCAACCGGCGGCGGCACTGGCAAAACCCGATTAAGCACCGCCTCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGGAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGGACCCGACTACCTAGCATGGATGCGACACCGATTTTCGCCGCGTCCAGCACCACATCCa                                                   <  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)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                      
GCAATTTCAGTTTGCCCTTTTCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGATAGCGCCCATTAAGTACGCCGGGATTTACTGCCGTCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTGAGCTTTGTTTTGCGCGGGTAACGTCGCTAATGAATC  >  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: