Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,403,881 T→C T229A (ACC→GCC)  abgA ← p‑aminobenzoyl‑glutamate hydrolase subunit A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,403,8810TC91.7% 28.9 / ‑3.6 12T229A (ACC→GCC) abgAp‑aminobenzoyl‑glutamate hydrolase subunit A
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (6/5);  total (7/5)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.76e-01

GGCGATTGCATGCAGTGCAAGAGTGGCTTGTGCTGCCGCCAACAAGGCATTGTGACCGTCTTCTGGTTTTGCGCCTGCGTGAGCGGCGGTACCGGTGAAGTGCGCGTCAAATTTGGTGGTTGCCATAAAATTATCACTGCCGCACACCACGGTGCCCGCAGGTACGCCAGTGCCAATGTGCACGGCAGTAAAATAATCAACATCATCTACGACACCTGCATCGACCATCGCCCGCGCGCCACGC  >  NC_000913/1403765‑1404008
                                                                                                                     |                                                                                                                               
ggCGATTGCATGCAGTGCAAGAGTGGCTTGTGCTGCCGCCCACAAGGCATTGTGACCGTCTTCTGGGTTTTCGCCTGCGTGAGCGGCGGGGCCGGGGGAGGGGGCGGCAAAATTTGGGGGTGGCCTAAAAATATTACtcc                                                                                                          >  2:225929/1‑138 (MQ=255)
                                   ccgccAACAAGGCATTGTGACCGTCTTCTGGTTTTGCGCCTGCGTGAGCGGCGGTACCGGTGAAGTGCGCGTCAAATTTGGCGGTTGCCATAAAATTATCACTGCCGCACACCACGGTGCCCGCAGGTACGCCAGTGCCa                                                                       >  1:335964/1‑140 (MQ=255)
                                   ccgccAACAAGGCATTGTGACCGTCTTCTGGTTTTGCGCCTGCGTGAGCGGCGGTACCGGTGAAGTGCGCGTCAAATTTGGCGGTTGCCATAAAATTATCACTGCCGCACACCACGGTGCCCGCAGGTACGCCAGTGCCa                                                                       <  2:335964/140‑1 (MQ=255)
                                             ggCATTGTGACCGTCTTCTGGTTTTGCGCCTGCGTGAGCGGCGGTACCGGTGAAGTGCGCGTCAAATTTGGCGGTTGCCATAAAATTATCACTGCCGCACACCACGGTGCCCGCAGGTACGCCAGTGCCAATGTGCACgg                                                             <  2:197152/140‑1 (MQ=255)
                                                    tgACCGTCTTCTGGTTTTGCGCCTGCGTGAGCGGCGGTACCGGTGAAGTGCGCGTCAAATTTGGCGGTTGCCATAAAATTATCACTGCTGCACACCACGGTGCCCGCAGGTACGCCAGTGCCAATGTGCACGGCAGTaaa                                                      <  1:9298/140‑1 (MQ=255)
                                                         gTCTTCTGGTTTTGCGCCTGCGTGAGCGGCGGTACCGGTGAAGTGCGCGTCAAATTTGGCGGTTGCCATAAAATTATCACTGCCGCACACCACGGTGCCCGCAGGTACGCCAGTGCCAATGTGCACGGCAGTAAaataat                                                 <  1:324050/140‑1 (MQ=255)
                                                              ctGGTTTTGCGCCTGCGTGAGCGGCGGTACCGGTGAAGTGCGCGTCAAATTTGGCGGTTGCCAAAAAAATATCACTGCCGCACACCCCGGGGCCCGCAGGGACCCCAGTGCCAATGTGCACGGGAGGAAAAAAATCAAca                                            >  2:37899/1‑140 (MQ=255)
                                                                     tGCGCCTGCGTGAGCGGCGGTACCGGTGAAGTGCGCGTCAAATTTGGCGGTTGCCATAAAATTATCACTGCCGCACACCACGGGGCCCGCAGGGACGCCAGTGCCAATGTGCACGGCAGTAAAATAAACAACATCATCtc                                     >  2:324050/1‑139 (MQ=255)
                                                                     tGCGCCTGCGTGAGCGGCGGTACCGGTGAAGTGCGCGTCAAATTTGGCGGGTGCCATAAAATTATCACTGCCGCACACCACGGTGCCCGCAGGGACGCCAGTGCCAATGTGCACGGCAGTAAAATAATCAACATCATCTa                                     >  2:269704/1‑140 (MQ=255)
                                                                                   cggcggTACCGGTGAAGTGCGCGTCAAATTTGGCGGTTGCCATAAAATTATCACTGCCGCACACCACGGTGCCCGCAGGGACGCCAGTGCCAATGTGCACGGCAGTAAAATAATCAACATCATCTACGACACCTGCATCg                       >  2:334227/1‑140 (MQ=255)
                                                                                     gcggTACCGGTGAAGTGCGCGTCAAATTTGGCGGTTGCCATAAAATTATCACTGCCGCACACCACGGTGCCCGCAGGTACGCCAGTGCCAATGTGCACGGCAGTAAAATAATCAACATCATCTACGACACCTGCATCGAc                     <  1:183931/140‑1 (MQ=255)
                                                                                        gTACCGGGGAAGTGCGCGTCAAATTTGGCGGGTGCCATAAAATTATCACCGCCGCACACCACGGGGCCCGCAGGGACGCCCGTGCCAATGTGCACGGCAGTAAAATAATCAACATCATCTACGACACCTGCATCGACCAt                  >  2:117135/1‑140 (MQ=255)
                                                                                                   gTGCGCGTCAAATTTGGCGGTTGCCATAAAATTATCACTGCCGCACACCACGGTGCCCGCAGGTACGCCAGTGCCAATGTGCACGGCAGTAAAATAATCAACATCATCTACGACACCTGCATCGACCATCGCCCgcgcgc       >  1:48971/1‑140 (MQ=255)
                                                                                                        cgTCAAATTTGGTGGTTGCCATAAAATTATCACTGCCGCACACCACGGTGCCCGCAGGTACGCCAGTGCCAATGTGCACGGCAGTAAAATAATCAACATCATCTACGACACCTGCATCGACCATCGCCCGCGCGCCAcgc  >  2:242101/1‑140 (MQ=255)
                                                                                                                     |                                                                                                                               
GGCGATTGCATGCAGTGCAAGAGTGGCTTGTGCTGCCGCCAACAAGGCATTGTGACCGTCTTCTGGTTTTGCGCCTGCGTGAGCGGCGGTACCGGTGAAGTGCGCGTCAAATTTGGTGGTTGCCATAAAATTATCACTGCCGCACACCACGGTGCCCGCAGGTACGCCAGTGCCAATGTGCACGGCAGTAAAATAATCAACATCATCTACGACACCTGCATCGACCATCGCCCGCGCGCCACGC  >  NC_000913/1403765‑1404008

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