Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,329,418 T→C Q127R (CAA→CGA)  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,4180TC96.2% 85.8 / NA 26Q127R (CAA→CGA) yfaQtandem DUF2300 domain‑containing protein YfaQ
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (12/13);  minor base G (0/1);  total (12/14)
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

TGTTTTGCGCGGGTAACGTCGCTAATGAATCAAGCAGCGATTTAACCGTTACCGAGGTTTCTGGCTTCAGCGTAGTCAGAGACGCTAACCACGCAGGTGCCTGACGCTGTTGCCAGTATTGCCCCCACACATCTGCGCCAATATGCAACCGCTGCGGCGCAAAATAGAGCCCACAGGAACGTACCAGCGCGGTATCGCGCGTTATGCTCTCCCCAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGAC  >  NC_000913/2329299‑2329553
                                                                                                                       |                                                                                                                                       
tGTTTTGCGCGGGTAACGTCGCTAATGAATCAAGCAGCGATTTAACCGTTACCGAGGTTTCTGGCTTCAGCGTAGTCAGAGACGCTAACCACGCAGGTGCCTGACGCTGTTGCCAGTATCGCCCCcac                                                                                                                                 <  2:784864/128‑1 (MQ=255)
tGTTTTGCGCGGGTAACGTCGCTAATGAATCAAGCAGCGATTTAACCGTTACCGAGGTTTCTGGCTTCAGCGTAGTCAGAGACGCTAACCACGCAGGTGCCTGACGCTGTTGCCAGTATCGCCCCcac                                                                                                                                 >  1:784864/1‑128 (MQ=255)
tGTTTTGCGCGGGTAACGTCGCTAATGAATCAAGCAGCGATTTAACCGTTACCGAGGTTTCTGGCTTCAGCGTAGTCAGAGACGCTAACCACGCAGGTGCCTGACGCTGCTGCCAGTATCGCCCCCACACATCTGCGCCa                                                                                                                     >  1:508052/1‑140 (MQ=255)
                      tAATGAATCAAGCAGCGATTTAACCGTTACCGAGGTTTCTGGCTTCAGCGTAGTCAGAGACGCTAACCACGCAGGTGCCTGACGCTGTTGCAAGTATCGCCCCCCCACGTCTGCGCCCATATACAACCGCTGCGGCGCaa                                                                                               >  2:349815/1‑140 (MQ=255)
                       aaTGAATCAAGCAGCGATTTAACCGTTACCGAGGTTTCTGGCTTCAGCGTAGTCAGAGACGCTAACCACGCAGGTGCCTGACGCTGTTGCCAGTATCGCCCCCACACATCTGCGCCAATATGCAACCGCTGCGGCGCaaa                                                                                              <  2:39854/140‑1 (MQ=255)
                        aTGAATCAAGCAGCGATTTAACCGTTACCGAGGTTTCTGGCTTCAGCGTAGTCAGAGACGCTAACCACGCAGGTGCCTGACGCTGTTGCCAGTATCGCCCCCACACATCTGCGCCAATATGCAACCGCTGCGGCGCaaaa                                                                                             <  1:635458/140‑1 (MQ=255)
                            aTCAAGCAGCGATTTCACCGTTTCCCAGGCTTCTGGCCTCAGCGTAGTCAGAGACGCTAACCACGCAGGTGCCTGACGCTGTTGCCAGTATCGCCCCCACACATCTGCGCCAATATGCAACCGCTGCGGCGCAAAATaga                                                                                         <  1:218428/140‑1 (MQ=255)
                             tCAAGCAGCGATTTAACCGTTACCGAGGTTTCTGGCTTCAGCGTAGTCAGAGACGCTAACCACGCAGGTGCCTGACGCTGCTGCCAGTATCGCCCCCACACATCTGCGCCAATATGCAACCGCTGCGGCGCAAAATagag                                                                                        <  2:508052/140‑1 (MQ=255)
                               aaGCAGCGATTTAACCGTTACCGAGGTTTGTGGCTTCAGCGTAGTCAGAGACGCTAACCACGCAGGTGCCTGACGCTGTTGCCAGTATCGCCCCCACACATCTGCGCCAATATGCAACCGCTGCGGCGCAAAATAGAGcc                                                                                      <  1:349815/140‑1 (MQ=255)
                                agcagcGATTTAACCGTTACCGAGGTTTCTGGCTTCAGCGTAGTCAGAGACGCTAACCACGCAGGTGCCTGACGCTGTTGCCAGTATCGCCCCCACACATCTGCGCCAATATGCAACCGCTGCGGCGCAAAATaga                                                                                         >  1:34851/1‑136 (MQ=255)
                                agcagcGATTTAACCGTTACCGAGGTTTCTGGCTTCAGCGTAGTCAGAGACGCTAACCACGCAGGTGCCTGACGCTGTTGCCAGTATCGCCCCCACACATCTGCGCCAATATGCAACCGCTGCGGCGCAAAATaga                                                                                         <  2:34851/136‑1 (MQ=255)
                                 gcagcGATTTAACCGTTACCGAGGTTTCTGGCTTCAGCGTAGTCAGAGACGCTAACCACGCAGGTGCCTGACGCTGTTGCCAGTATCGCCCCCACACATCTGCGCCAATATGCAACCGCTGCGGCGCAAAATAGAGCCca                                                                                    >  2:735783/1‑140 (MQ=255)
                                            aCCGTTACCGAGGTTTCTGGCTCCACCGTAGTCAGAGACGCTAACCACGCAGGTGCCTGACGCTGTTGCCAGTATCGCCCCCACACATCTGCGCCAATATGCAACCGCTGCGGCGCAAAATAGAGCCCACAGGAACGTAc                                                                         <  1:305481/140‑1 (MQ=255)
                                                 tACCGAGGTTTCTGGCTTCAGCGTAGTCAGAGACGCTAACCACGCAGGTGCCTGACGCTGTTGCCAGTATCGCCCCCACACATCTGCGCCAATATGCAACCGCTGCGGCGCAAAATAGAGCCCACAGGAACGTACCAgcg                                                                    >  1:767318/1‑140 (MQ=255)
                                                 tACCGAGGTTTCTGGCTTCAGCGTAGTCAGAGACGCTAACCACGCAGGTGCCTGACGCTGTTGCCAGTATCGCCCCCACACATCTGCGCCAATATGCAACCGCTGCGGCGCAAAATAGAGCCCACAGGAACGTACCAgcg                                                                    <  2:576998/140‑1 (MQ=255)
                                                  aCCGAGGTTTCTGGCTTCAGCGTAGTCAGAGACGCTAACCACGCAGGTGCCTGACGCTGTTGCCAGTATCGCCCCCACACATCTGCGCCAAGATGCAACCACTGGGGCCCAAAATAGCGCCCACAGGAACGTACCAgcgc                                                                   >  2:305481/1‑140 (MQ=255)
                                                   ccGAGGTTTCTGGCTTCAGCGTAGTCAGAGACGCTAACCACGCAGGTGCCTGACGCTGTTGCCAGTATCGCCCCCACACATCTGCGCCAATATGCAACCGCTGCGGCGCAAAATAGAGCCCACAGGAACGTACCAGcgcg                                                                  <  1:22793/140‑1 (MQ=255)
                                                        gTTTCTGGCTTCAGCGTAGTCAGAGACGCTAACCACGCAGGTGCCTGACGCTGTTGCCAGTATCGCCCCCACACATCTGCGCCAATATGCAACCGCTGCGGCGCAAAATAGAGCCCACAGGAACGTACCAGCGCGGTATc                                                             <  2:467116/140‑1 (MQ=255)
                                                                   cAGCGTAGTCAGAGACGCTAACCACGCAGGTGCCTGACGCTGTTGCCAGTATCGCCCCCACACATCTGCGCCAATATGCAACCGCTGCGGCGCAAAATAGAGGCCACAGGGACGGACGAGCGCGGTATCGCGCGTTATGc                                                  >  1:230228/1‑140 (MQ=255)
                                                                    aGCGTAGTCAGAGACGCTAACCACGCAGGTGCCTGACGCTGTTGCCAGTATCGCCCCCACACATCTGCGCCAATATGCAACCGCTGCGGCGCAAAATAGAGCCCACAGGAACGTACCAGCGCGGTATCGCGCGTTATGct                                                 >  2:764422/1‑140 (MQ=255)
                                                                    aGCGTAGTCAGAGACGCTAACCACGCAGGTGCCTGACGCTGTTGCCAGTATCGCCCCCACACATCTGCGCCAATATGCAACCGCTGCGGCGCAAAATAGAGCCCACAGGAACGTACCAGCGCGGTATCGCGCGTTATGc                                                  >  2:2381/1‑139 (MQ=255)
                                                                    aGCGTAGTCAGAGACGCTAACCACGCAGGTGCCTGACGCTGTTGCCAGTATCGCCCCCACACATCTGCGCCAATATGCAACCGCTGCGGCGCAAAATAGAGCCCACAGGAACGTACCAGCGCGGTATCGCGCGTTATGc                                                  <  1:2381/139‑1 (MQ=255)
                                                                               agaCGCTAACCACGCAGGTGCCTGACGCTGTTGCCAGTATCGCCCCCACACATCTGCGCCAATATGCAACCGCTGCGGCGCAAAATAGAGCCCACAGGAACGGACCAGCGTGGTATCGCGCGTTATGCTCTCCCCAGCCt                                      >  1:309006/1‑140 (MQ=255)
                                                                                  cGCTAACCACGCAGGGGCCTGACGCGGTTGCCAGTATGGCCCCCACACATCTGCGCCAATATGCAACCGCTGCGGCGCAAAATAGAGCCCACAGGAACGTACCAGCGCGGTATCGCGCGTTATGCTCTCCCCAGCCTGac                                   <  1:285044/140‑1 (MQ=255)
                                                                                            gCAGGTGCCTGACGCTGTTGCCAGTATCGCCCCCACACATCTGCGCCAATATGCAACCGCTGCGGCGCAAAATAGAGCCCACAGGAACGTACCAGCGCGGTATCGCGCGTTATGCTCTCCCCAGCCTGACAACAATAAAc                         >  1:482554/1‑140 (MQ=255)
                                                                                                                   gTATCGCCCCCACACATCTGCGCCAATATGCAACCGCTGCGGCGCAAAATAGAGCCCACAGGAACGTACCAGCGCGGTATCGCGCGTTATGCTCTCCCCAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGac  <  1:735783/140‑1 (MQ=255)
                                                                                                                       |                                                                                                                                       
TGTTTTGCGCGGGTAACGTCGCTAATGAATCAAGCAGCGATTTAACCGTTACCGAGGTTTCTGGCTTCAGCGTAGTCAGAGACGCTAACCACGCAGGTGCCTGACGCTGTTGCCAGTATTGCCCCCACACATCTGCGCCAATATGCAACCGCTGCGGCGCAAAATAGAGCCCACAGGAACGTACCAGCGCGGTATCGCGCGTTATGCTCTCCCCAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGAC  >  NC_000913/2329299‑2329553

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