Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,023,436 G→A R234H (CGC→CAC)  tatC → twin arginine protein translocation system ‑ TatC protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,023,4360GA96.0% 74.3 / NA 25R234H (CGC→CAC) tatCtwin arginine protein translocation system ‑ TatC protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (9/15);  minor base T (1/0);  total (10/15)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTTGGTGCATTCGTTGTCGGGATGTTGCTGACGCCGCCGGATGTCTTCTCGCAAACGCTGTTGGCGATCCCGATGTACTGTCTGTTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATG  >  NC_000913/4023326‑4023556
                                                                                                              |                                                                                                                        
ggttggtGCATTCGTTGTCGGGATGTTGCTGACGCCGCCGGATGTCTTCTCGCAAACGCTGTTGGCGATCCCGATGTACTGTCTGTTTGAAATCGGTGTCTTCTTCTCACACTTTTACGTTGGTAAAGGGCGAAATCggg                                                                                             <  2:444309/140‑1 (MQ=255)
 gttggtGCATTCGTTGTCGGGATGTTGCTGACGCCGCCGGATGTCTTCTCGCAAACGCTGTTGGCGATCCCGATGTACTGTCTGTTTGAAAACGGTGTCTACTTCTGACTCTTTTACGTTGGACAAGGGCGCAATCGGGa                                                                                            >  2:359263/1‑140 (MQ=255)
   tggtGCATTCGTTGTCGGGATGTTGCTGACGCCGCCGGATGTCTTCTCGCAAACGCTGTTGGCGATCCCGATGTACTGTCTGTTTGACATCGGTGTCTTCTTCTCACACTTTTACGTTGGTAAAGGGCGAAATCGGGAag                                                                                          <  2:430278/140‑1 (MQ=255)
      agCATTCGGTGTCGGGATGTTGCTGACGCCGCCGGATGTCTTGTCGCAAACGCTGTTGGCGGTCCCGATGCACTGTCTGTTTGAAATCGGTGTCTTCTTCTCACACTTTTACGTTGGTAAAGGGCGAAATCGGGAAgggg                                                                                       >  2:578411/2‑140 (MQ=255)
             gTTGTCGGGATGTTGCTGACGCCGCAGGATGTCTTCTAGCAAACGCTGTTGGCGATCCCGATGTACTGTCTGTTTGAAATCGGTGTCTTCTTCTCACACTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAacga                                                                                <  1:327761/140‑1 (MQ=255)
               tGTCGGGATGTTGCTGACGCCGCCGGATGTCTTCTCGCAAACGCTGTTGGCGATCCCGATGTACTGTCTGTTTGAAATCGGTGTCTTCTTCTCACACTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAacgacg                                                                              <  1:65271/140‑1 (MQ=255)
                                       ggATGTCTTCTCGCAGACGCTGGTGCCGATCCCGATGCACTGTCTGTTTGAAATCGGTGTCTTCTTCTCACACTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTg                                                      <  1:359263/140‑1 (MQ=255)
                                                        cGCTGTTGGCGATCCCGATGTACTGTCTGTTTGAAATCGGTGTCTTCTTCTCACACTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCg                                     >  2:306974/1‑140 (MQ=255)
                                                                     cccGATGTACTGTCTGTTTGAAATCGGTGTCTTCTTCTCACACTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCgg                        <  1:306974/140‑1 (MQ=255)
                                                                     cccGATGTACTGTCTGTTTGAAATCGGTGTCTTCTTCTCACACTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCgg                        <  1:65623/140‑1 (MQ=255)
                                                                     cccGATGTACTGTCTCTTTGAACTCCGTGACTTCTTCCCACACTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCgg                        <  1:237061/140‑1 (MQ=255)
                                                                        agacagctgtctgtTAGAAATCCGTATCTCCTTCTCATACTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTg                     <  1:240347/134‑1 (MQ=255)
                                                                          tGTACTGTCTGTTTGAAATCGGTGTCTTCTTCTCACACTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGc                          <  2:108513/133‑1 (MQ=255)
                                                                          tGTACTGTCTGTTTGAAATCGGTGTCTTCTTCTCACACTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGc                          >  1:108513/1‑133 (MQ=255)
                                                                             aCTGTCTGTTTGAAATCGGTGTCTTCTTCTCACACTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCata                <  1:733650/140‑1 (MQ=255)
                                                                               tgtctgtTTGAAATCGGTGTCTTCTTCTCACACTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATg              <  2:28238/140‑1 (MQ=255)
                                                                                 tctgtTTGAAATCGGTGTCTTCTTCTCACACTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGa            >  2:385399/1‑140 (MQ=255)
                                                                                      ttGAAATCGGTGTCTTCTTCTCACACTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTc                   <  1:590622/128‑1 (MQ=255)
                                                                                      ttGAAATCGGTGTCTTCTTCTCACACTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTc                   >  2:590622/1‑128 (MQ=255)
                                                                                      ttGAAATCGGTGTCTTCTTCTCACACTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCata                >  1:696317/1‑131 (MQ=255)
                                                                                      ttGAAATCGGTGTCTTCTTCTCACACTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCata                <  2:696317/131‑1 (MQ=255)
                                                                                      ttGAAATCGGTGTCTTCTTCTCACACTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACa       >  2:502921/1‑140 (MQ=255)
                                                                                         aaaTCGGTGTCTTCTTCTCACACTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAAGCGCGCGTCAGGGGGGTTGTCATATGGAGtctgtct    >  2:240347/1‑134 (MQ=255)
                                                                                           aTCGGTGTCTTCTTCTCACACTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCata                <  1:77339/126‑1 (MQ=255)
                                                                                           aTCGGTGTCTTCTTCTCACACTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCata                >  2:77339/1‑126 (MQ=255)
                                                                                           aTCGGTGTCTTCTTCTCACACTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATg  >  1:366520/1‑140 (MQ=255)
                                                                                                              |                                                                                                                        
GGTTGGTGCATTCGTTGTCGGGATGTTGCTGACGCCGCCGGATGTCTTCTCGCAAACGCTGTTGGCGATCCCGATGTACTGTCTGTTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATG  >  NC_000913/4023326‑4023556

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