Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_pae_tpiA 4,111,596 G→A L179L (CTG→TTG)  yiiQ ← DUF1454 family putative periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_tpiA4,111,5960GA90.0% 17.8 / ‑6.1 10L179L (CTG→TTG) yiiQDUF1454 family putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (6/3);  minor base T (0/1);  total (6/4)
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

GCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTG  >  NC_000913_3_pae_tpiA/4111475‑4111714
                                                                                                                         |                                                                                                                      
cgatctggTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCATCTTAATCGTTTCAACATCTAAGGTCAACCCCTTTACGCCGt                                                                                                          <  1:84294/130‑1 (MQ=255)
                     ataGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTAGACGTTGTCCGCGACAACATagatcg                                                                                     >  2:84294/1‑131 (MQ=255)
                          gttgcATCGGTAAAAAGTCTTTGTTGTTTGACATTTTTTCAAACCTTCAAGCTATCCAGATTGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCATCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCa                                                                                <  2:135353/132‑1 (MQ=255)
                                               tgctttTGGTCGTTTAGTCAAACGTTTAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGGCAACCCCTTTTCGCCGTTGTCCGAGACGACATAACGCGGTGCACCTTCTGTCTCGCTCtcataa                                                           >  2:121873/1‑136 (MQ=255)
                                                     ttGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACgtttg                                                     >  1:149136/1‑136 (MQ=255)
                                                              tttCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTTTGTCTCGGTGTTATATCGTTTGTTTTTcccc                                            <  1:149593/135‑1 (MQ=255)
                                                              aTTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGacaac                                                                                             <  1:176323/87‑1 (MQ=255)
                                                              aTTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGacaac                                                                                             >  2:176323/1‑87 (MQ=255)
                                                                                   aGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGATCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGt                       >  1:108443/1‑136 (MQ=255)
                                                                                                        ttCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTg  >  2:41738/1‑136 (MQ=255)
                                                                                                                         |                                                                                                                      
GCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTG  >  NC_000913_3_pae_tpiA/4111475‑4111714

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