Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,366,277 G→A G243D (GGC→GAC)  puuE → 4‑aminobutyrate aminotransferase PuuE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,366,2770GA90.0% 14.2 / ‑5.9 10G243D (GGC→GAC) puuE4‑aminobutyrate aminotransferase PuuE
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (5/4);  minor base C (0/1);  total (5/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 = 1.00e+00

GGTTCGGTGTATCACGTACCTTATCCGTCAGATTTACACGGCATTTCAACACAGGACTCCCTCGACGCCATCGAACGCTTGTTTAAATCAGACATCGAAGCGAAGCAGGTGGCGGCGATTATTTTCGAACCGGTGCAGGGCGAGGGCGGTTTCAACGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTGCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCTGTTTGCCATGGATCATTACGCCGATAAGCCGGATTTAATGACGATGGCGAAAAGCCTCGCGGGCGGGATGCCGCTTTCGGGCGTGGTCGGTAACGCGAATATTATGGACGCACCCGCGCCGGGCGGGCTTGGCGGCACCTACGCCGGTAACCCGCTGGCGGTGGCTGCCGCGCACGCGGTGCTCAACATTATCGACAAAGAATCACTCTG  >  NC_000913/1366036‑1366508
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                       
ggTTCGGTGTATCACGTACCTTATCCGTCAGATTTACACGGCATTTCAACACAGGACTCCCTCGACGCCATCGAACGCTTGTTTAAATCAGACATCGAAGCGAAGCAGGTGGCGGCGATTATTTTCGAACCGGTGCAGGGCGAGGGCGGTTTCAACGGTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTGCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGACTTTGcgcg                                                                                                                                                                                                                                <  1:248712/251‑1 (MQ=255)
                              gATTTACACGGCATTTCAACACAGGACTCCCTCGACGCCATCGAACGCTTGTTTAAATCAGACATTGAAGCGAAGCAGGTGGCGGCGATTATTTTCGACCAGGTGCAGGGCGCGGGCGGTTTCAACGTTGCGCCAAAAGGGCTGGTTGCCGCTATTCGCCGCCTGTGCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGACTTTGCGCGTACCGGTAAGCTGTTTGCCATGGATCAtt                                                                                                                                                                                                   >  1:467796/1‑250 (MQ=255)
                                                                                            cATCGAAGCGAAGCAGGTGGCGGCGATTATTTTCGAACCGGTGCAGGGCGAGGGCGGTTTCAACGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTGCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGACTTTGCGCGTACCGGTAAGCTGTTTGCCATGGATCATTACGCCGATAAGCCGGATTTAATGACGATGGCGa                                                                                                                                                                  >  1:474277/1‑221 (MQ=255)
                                                                                                                   cgATTATTTTCGAACCGGTGCAGGGCGAGGGCGGTTTCAACGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTGCGAAGCGCACGGTATTGCGATGATTGCTGATGAAGTGCAAAGCGACTTTGCGCGTACCGGTAAGCTGTTTGCCATGGATCATTACGCCGATAAGCCGGATTTAATGACGATGGCGAAAAGCCTCGCGGGCGGGATGCCGCTTTCGGGCGTGGTCGGTAACGCGAAtatt                                                                                                             >  1:289396/1‑251 (MQ=255)
                                                                                                                                        agggcgagggcgGTTTCAACGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTGCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAGGCGACTTGGCTCTTTCCGGTAAGCTGTTTGCCATGGATCATTACGCCGATAAGCCGGAGTTAATGACGATGGCGAAAAGCCTCGCGGGCGGGATGc                                                                                                                                             <  1:252286/198‑1 (MQ=255)
                                                                                                                                                                                    gCTATTCGCCGCCTGTGCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAGGCGACTTTGCGCGTACCGGTAAGCTGTTTGCCATGGATCATTACGCCGATACGCCGGATTTAATGACGATGGCGAAAAGCCTCGCGGGCGGGCTGCCGCTTTCGGGCGTGTTCGGTAACGCGAATATTATGGACGCACCCGCGCCGGGCGGGCTTgg                                                                                >  1:464554/1‑215 (MQ=255)
                                                                                                                                                                                        ttCGCCGCCTGTGCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGACTTTGCGCGTACCGGTAAGCTGTTTGCCATGGATCATTACGCCGATAAGCCGGATTTAATGACGATGGCGAAAAGCCTCGCGGGCGGGATGCCGCTTTCgg                                                                                                                                    <  1:260902/159‑1 (MQ=255)
                                                                                                                                                                                                                              gCTGATGAAGTGCAAAGCGACTTTGCGCGTACCTGTAAGCTGTTTGCCATGGATCATTACGCCGATAAGCCGGATTTAATGACGATGGCGAAAAGCCTCGCGGGCGGGATGCCGCTTTCGGGCGTGGTCGGTAACGCGAATATTATGGACGCACCCGCGCCGGGCGGGCTTGGCGACACCTACGCCGGTAACCCGCTGGCGGGGGCTGCCGCGCACGCGGTGCTCAACATTATCGACAAAGAATCACTCTg  <  1:330391/251‑1 (MQ=255)
                                                                                                                                                                                                                                        tGCAAAGCGACTTTGCGCGTACCGGTAAGCTGTTTGCCATGGATCATTACGCCGATAAGCAGGATTTAATGAAGATGGCGAAAAGCCTCGCGGGCGGGATGCCGCTTTCGGGCGTGGTCGGTaa                                                                                                                       >  1:64582/1‑124 (MQ=255)
                                                                                                                                                                                                                                          cAAAGCGCCTTTGCGCGTACCGGTAAGCTGTTTGCCATGGATCATTACGCCGATAAGCCGGATTTAATGACGAGGGCGAAAAGCCTCGCGGGCGGGATGCCGCTTTCGGGCGTGGTCGGTAACGCGAATATTATGGACGCACCCGCGCCGGGCGGGCTTGGCGACACCTACGCCGGTAACCCGCTGGCGGTGGCTGCcgc                                         <  1:363576/200‑1 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                       
GGTTCGGTGTATCACGTACCTTATCCGTCAGATTTACACGGCATTTCAACACAGGACTCCCTCGACGCCATCGAACGCTTGTTTAAATCAGACATCGAAGCGAAGCAGGTGGCGGCGATTATTTTCGAACCGGTGCAGGGCGAGGGCGGTTTCAACGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTGCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCTGTTTGCCATGGATCATTACGCCGATAAGCCGGATTTAATGACGATGGCGAAAAGCCTCGCGGGCGGGATGCCGCTTTCGGGCGTGGTCGGTAACGCGAATATTATGGACGCACCCGCGCCGGGCGGGCTTGGCGGCACCTACGCCGGTAACCCGCTGGCGGTGGCTGCCGCGCACGCGGTGCTCAACATTATCGACAAAGAATCACTCTG  >  NC_000913/1366036‑1366508

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