Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,894,331 G→A V67I (GTT→ATT)  yoaC → DUF1889 domain‑containing protein YoaC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,894,3310GA81.2% 20.6 / 3.5 16V67I (GTT→ATT) yoaCDUF1889 domain‑containing protein YoaC
Reads supporting (aligned to +/- strand):  ref base G (3/0);  new base A (9/4);  total (12/4)
Fisher's exact test for biased strand distribution p-value = 5.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.77e-01

GGATTTCATTGGTGCCATGGATGTATCAGCGCCAACACCAAGTTCGATGAATGAAAGCACGGCGAAGGGAATCTTTAAATATTTAAAAGAACTGGGAGTACCCGCGAGTGCCGCTGATATTACGGCGCGAGCTGACCAGGAAGGCTGGAACCCTGGGTTCACGGAAAAAATGGTTGGATGGGCAAAAAAAATGGAGACAGGTGAACGTTCTGTGATTAAAAATCCTGAATACTTTTCAACATATATGCAGGAAGAACTAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGTCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTACTCGTCTTCATCCTC  >  NC_000913/1894159‑1894569
                                                                                                                                                                            |                                                                                                                                                                                                                                              
ggATTTCATTGGTGCCATGGATGTATCAGCGCCAACACCAAGTTCGATGAATGAAAGCACGGCGAAGGGAATCTTTAAATATTTAAAAGAACTGGGAGTTCCCGCGATTGCCGCTGATATTACGGCGCGAGCTGACCAGGAAGGCTGGAACCCTGGGTTCACGGAAAAAATTATTGGATGGGCAAAAAAAATGGAGACAGGTGAACGGTCTGTGATTAAAAATCCTGAATACTTTTCAACatata                                                                                                                                                                        <  1:149021/245‑1 (MQ=255)
                             cgcCAACACCAAGTTCGATGAATGAAAGCACGGCGAAGGGCATCTTTAAATATTTAAAAGAACTGGGAGTCCCCGCGAGTGCCGCTGATATTACGGCGCGAGCTGACCAGGAAGGCTGGAACCCTGGGTTCACGGAAAAAATGATTGGGTGGGCAAAAAAAATg                                                                                                                                                                                                                            >  1:165939/1‑164 (MQ=255)
                                              atgaatgaAAGCACGGCGAAGGGAATCTTTAAATATTTAAAAGAACTGGGAGTACCCGCGAGTGCCGCTGATATGACGGCGCGAGCTGACCAGGAAGGCTGGAACCCTGGGTTCACGGAAAAAATGATTGGATGGGCAAAAAAAATGGAGACAGGTGAACGTTCTGTGATTAAAAATCCTGAATACTTTTCAACATATATGCAgg                                                                                                                                                                  <  1:163577/205‑1 (MQ=255)
                                                                           tAAATATTTAAAAGAACTGGGAGTACCCGCGAGTGCCGCTGATATTACGGCGCGAGCTGACAAGGAAGGCTGGAACCCTGGGTTAACGGAAAAAATGATTGGATGGGCAAAAAAAATGGAGACAGGTGAACGTTCTGTGATTAAAAATCCTGAATACTTTTCAACATATACGCAGGAAGAACTAAAAGCACTGGTCTGAGTTAAATTTATATCAg                                                                                                                           >  1:187807/1‑215 (MQ=255)
                                                                                    naaaGAACTGGGAGTACCCGCGAGTGCCGCTGATATTACGGCGCGAGCTGACCAGGAAGGCTGGAACCCTGGGTTCAAGGAAAAAATGATTGGATGGGCAAACAAAATGGAGACAGGTGAACGTTCTGTGATTAAAAATCCTGAATACTTTTCAACATATATGCAGGAAGAACTAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGAc                                                                                                          >  1:59503/2‑223 (MQ=255)
                                                                                                    cccGCGAGTGCCGCTGATATTACGGCGCGAGCTGACCAGGAAGGCTGGAACCCTGGGTTCACGGAAAAAATGATTGGATGGGCAAAAAAAATGGAGACAGGTGAACGTTCTGTGATTAAAAATCCTGAATACTTTTCAACATATATGCAGGAAGAACTAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGTCCTTTTTCTCACAATCaa                                                               >  1:251545/1‑250 (MQ=255)
                                                                                                      cgcgAGTGCCGCTGATATTACGGCGCGAGCTGACCAGGAAGGCTGGAACCCTGGGTTCACGGAACAAATGGTTGGATGGGCAAAAAAAATGGAGACAGGTGAACGTTCTGTGATTAAAAATCCTGAATACTTTTCAACATATATGCAGGAAGAACTAAAAGCACTGGGCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGTCCTTTTTCTCACAATCAAATg                                                            >  1:247796/1‑251 (MQ=255)
                                                                                                         gAGTGCCGCTGATATTACGGCGCGAGCTGACCAGGAAGGCTGGAACCCTGGGTTCACGGAAAAAATGATTGGATGGGCAAAAAAAATGGAGACAGGTGAACGTTCTGTGATTAAAAATCCTGAATACTTTTCAACATATATGCAGGAAGAACTAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGTCCTTTTTCTCACAATCAAATGatc                                                         >  1:76297/1‑251 (MQ=255)
                                                                                                                                   cTGACCAGGAAGGCTGGAACCCTGGGTTCACGGAAAAACTGGTTGGATGGGCAAACAAACTGGAGACAGGCGAACGTTCTGTGATTAAAAATCCTGAATACTTTTCAACATATATGCAGGAAGAACTAAAAGCACTGGTCTg                                                                                                                                            >  1:131969/1‑142 (MQ=255)
                                                                                                                                   cTGACCAGGAAGGCTGGAACCCTGGGTTCACGGAAAAAATGATTGGATGGGCAAAAAAAATGGAGACAGGTGAACGTTCTGTTATTAAAAATCCTGAAGACTGTTCAACATATATGCAGGAAGAACTAAAAGCACTGGTCTGAGTTAAAtt                                                                                                                                   <  1:171151/151‑1 (MQ=255)
                                                                                                                                       ncAGGAAGGCTGGAACCCTGGGTTCACGGAAAAAATGGTTGGATGGGCAAAAAAAATGGAGACAGGTGAACGTTCTGTGATTAAAAATCCTGAATACTTTTCAACATATATGCAGGAAGAACTAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACAtata                                                                                         >  1:83678/2‑189 (MQ=255)
                                                                                                                                       ncAGGAAGGCTGGAACCCTGGGTTCACGGAAAAAATGATTGGATGGGCAAAAAAAATGGAGACAGGTGAACGTTCTGTGATTAAAAATCCTGAATACTTTTCAACATATATGCAGGAAGAACTAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGTCCTTTTTCTcac                                                                     >  1:201883/2‑209 (MQ=255)
                                                                                                                                                cTGGAACCCTGGGTTCACGGAAAAAATGATTGGATGGGCAAAAAAAATGGAGACAGGTGAACGTTCTGTGATTAAAAATCCTGAATACTTTTCAACATATATGCAGGAAGAACTAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAAtn                                                                                                                   <  1:195275/154‑2 (MQ=255)
                                                                                                                                                cTGGAACCCTGGGTTCACGGAAAAAATGAGTGGATGGGCAAAAACAATGGAGACAGGTGAACGTTCTGTGATTAAATATCCTGAATACTTTTCAACATATATGCAGGAAGAACTAATTGCACTGGTCTGAGTTNATTNNNNNTcnnnnnnnnnn                                                                                                                   >  2:195275/1‑144 (MQ=255)
                                                                                                                                                            gTTCACGGAAAAAATGATTGGATGGGCAAAAAAAATGGAGAAAGGTGAACGTTCTGTGATTAAAAATCCTGAATACTTTTCCACATATATGCAGGAAGAACTAACAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGTCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTACTCGTCTTCAt      >  1:213223/1‑251 (MQ=255)
                                                                                                                                                                 cGGAAAAAATGATTGGATGGGCAAAAAAAATGGAGACAGGTGAACGTTCTGTGATTAAAAATCCTGAATACTTTTCAACATATATGCAGGAAGAACTAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATACAAAGCGTCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGACATCGATACGGATTAGTGTGTTACTCGTCTTCATCCTc  >  1:219848/1‑250 (MQ=255)
                                                                                                                                                                            |                                                                                                                                                                                                                                              
GGATTTCATTGGTGCCATGGATGTATCAGCGCCAACACCAAGTTCGATGAATGAAAGCACGGCGAAGGGAATCTTTAAATATTTAAAAGAACTGGGAGTACCCGCGAGTGCCGCTGATATTACGGCGCGAGCTGACCAGGAAGGCTGGAACCCTGGGTTCACGGAAAAAATGGTTGGATGGGCAAAAAAAATGGAGACAGGTGAACGTTCTGTGATTAAAAATCCTGAATACTTTTCAACATATATGCAGGAAGAACTAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGTCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTACTCGTCTTCATCCTC  >  NC_000913/1894159‑1894569

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