Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,366,433 G→A G295D (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,4330GA96.8% 89.5 / NA 31G295D (GGC→GAC) puuE4‑aminobutyrate aminotransferase PuuE
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (15/15);  total (15/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

CAACGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTGCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCTGTTTGCCATGGATCATTACGCCGATAAGCCGGATTTAATGACGATGGCGAAAAGCCTCGCGGGCGGGATGCCGCTTTCGGGCGTGGTCGGTAACGCGAATATTATGGACGCACCCGCGCCGGGCGGGCTTGGCGGCACCTACGCCGGTAACCCGCTGGCGGTGGCTGCCGCGCACGCGGTGCTCAACATTATCGACAAAGAATCACTCTGCGAACGCGCGAATCAACTGGGCCAGCGTCTCAAAAACACGTTGATTGATGCCAAAGAAAGCGTTCCGGCCATTGCTGCGGTACGCGGCCTGGGGTCGATGATTGCGGTAGAGTTTAACGATCCGCAAACGGGCGAGCCGTCA  >  NC_000913/1366188‑1366650
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                         
cAACGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTGCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCTGTTTGCCATGGATCATTACGCCGATAAGCCGGATTTAATGACGATGGCGAAAAGCCTCGCGGGCGGGATGCCGCTTTCGGGCGTGTTCGGTAACGCGAATATTATGGACGCACCCGCGCCGGGCGGGCTTGGCGACAcc                                                                                                                                                                                                                       <  1:300248/250‑1 (MQ=255)
       gcgcCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTGCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCTGTTTGCCATGGATCATTACGCCGATAAGCCGGATTTAATGACGATGGCGAAAAGCCTCGCGGGCGGGATGCCGCTTTCGGGCGTGGTCGGTAACGCGAATCTTATGGACGCACCCGCGCCGGGCGGGCTTGGCGACACCTACGCCgg                                                                                                                                                                                                               >  1:56014/1‑251 (MQ=255)
                   cTGGTTGCCGCTATTCGCCGCCTGTGCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCTGTTTGCCATGGATCATTACGCCGATAAGCCGGATTTAATGACGATGGCGAAAAGCCTCGCGGGCGGGATGCCGCTTTCGGGCGTGGTCGGTAACGCGAATATTATGGACGCACCCGCGCCGGGCGGGCTTGGCGACACCTACGCCGGTAACCCGCTGGc                                                                                                                                                                                                   >  1:283976/1‑251 (MQ=255)
                                       ccTGTGCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCTGTTTGCCATGGATCATTACGCCGATAAGCCGGATTTAATGACGATGGCGAAAAGCCTCGCGGGCGGGATGCCGCTTTCGGGCGTGGTCGGTAACGCGAATATTATGGACGCACCCGCGCCGGGCGGGCTTGGCGACACCTACGCCGGTAAccc                                                                                                                                                                                                         <  2:532816/225‑1 (MQ=255)
                                       ccTGTGCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCTGTTTGCCATGGATCATTACGCCGATAAGCCGGATTTAATGACGATGGCGAAAAGCCTCGCGGGCGGGATGCCGCTTTCGGGCGTGGTCGGTAACGCGAATATTATGGACGCACCCGCGCCGGGCGGGCTTGGCGACACCTACGCCGGTAAccc                                                                                                                                                                                                         >  1:532816/1‑225 (MQ=255)
                                                     acggtatgtagatgatTGTTGATAAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGATGGGTGCAATGGATCATGACGCCGATAAGCCGGATTGAATGACGATGGCGAAAAGCCTCGCGGGCTGTATGCCGCTTTCGTGCGGGGTCGGTAACGCGAATAATATGGACGCACCCGCGCCGGGCGGGCTGGGCGGTACCTACGCCGGTATACCGCTGGCGGTGGCTGCTGTGCACGAGGTGCTTAACATTAt                                                                                                                                                                  <  2:354924/240‑1 (MQ=255)
                                                       ggtctagtgagaatTGATGATGCAGTGCAAATCGGCTTTGCCCGGACAGGTAAGCTGATTGCCATGGATCATTAACCCGCTAAACCGGATTTAATAACGATGACGAAAAGCCTCACGGGAGGGATGCAGCTTTCGGGCGTGGTCGGTAACGCGAATATTCTGGACGCACCCGCGCCGGGCGGGCTTGGCGACAACTACGCAGGTAACCCGCTGCCGGTGCGTGACGCGGACGCGGTGCTCAACaataacg                                                                                                                                                                >  2:213226/13‑244 (MQ=255)
                                                                atgatTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCTGTTTGCCATGGATCATTACGCCGATAAGCCGGATTTAATGACGATGGCGAAAAGCCTCGCGGGCGGGATGCCGCTTTCGGGCGTGGTCGGTAACGCGAATATTATGGACGCACCCGCGCCGGGCGGGCTTGGCGACACCTACGCCGGTAACCCGCTGGCGGTGGCTGCCGCGCACGCGGTGCTCAACATTATCGACAAAGAATc                                                                                                                                                      <  1:28636/251‑1 (MQ=255)
                                                                                 gCAAAGCGGCTTTGCGCGTACCGGTAAGCTGTTTGCCATGGATCATTACGCCGAGAAGCCGGATTTAATGACGATGGCGAAAAGCCTCGCGGGCGGGATGCCGCTTTCGGGCGTGGTCGGTAACGCGAATATTATGGACGCACCCGCGCCGGGCGGGCTTGGCGACACCTACGCCGGTAACCCGCTGGCGGTGGCTGCCGCGCACGCGGTGCTCAACATTATCGACAAAGAATCACTCTGCGAAcgcgcg                                                                                                                                      <  1:353664/250‑1 (MQ=255)
                                                                                   aaaGCGGCTTTGCGCGTACCGGTAAGCTGTTTGCCATGGATCATTACGCCGATAAGCCGGATTTAATGACGATGGCGAAAAGCCTCGCGGCCGGGATGCCGCTTTCGGGCGTGGACGGTACCGCGAATATTATGGACGCACCCGCGCCGGGCGGGCTTGGCGACACCTACGCCGGTAACCCGCTGGCGGTGGCTGCCGCGCACGCGGTGCTCAACATTATCGACAAAGAATCCATCTGCGAACGCGGGaac                                                                                                                                   >  1:505043/1‑250 (MQ=255)
                                                                                                      cGGTAAGCTGTTTGCCATGGCTCATGACGCCGATAAGCCGGATTTAATGACGATGGCGAAAAGCCTCGCGGGCGTGATGCCGCTTTCGGGCGTGGTCGGTAACGCGAATATTATGGACGCACCCGCGCCGGGCGGGCTTGGCGACTCCTACGCCGGTAACCCGCTGGCGGTGGCTGCCGCGCACGCGGTGCTCAACATTATCGACAAAGAAGCACTCTGCGAACGCGCGAATCAACTGGGCCAGCGTCTCa                                                                                                                <  2:297151/251‑1 (MQ=255)
                                                                                                       ggTAAGCTGTTGGCCATGGATCATTACGCAGATAAGCCGGATTTAATGACGATGGCGAAAAGCCTTGCGGGCGGGTCGCCGCTTTCGGGCGTGGTCGGTAACGCGAATATTATGGACGCACCCGCGCCGGGCGGGCTTGGCGACACCTACGCAGGTAACCCGATGGCGGTGGCTGCCGCGCACGCGGTGCTCAACATTATCGACAAAGCATCACTCTGCGAACGCGCGAAAAAAATGGGCCAGCGtctcc                                                                                                                >  2:7921/1‑249 (MQ=255)
                                                                                                              ggTTTGCCATGGATCATTACGCCGATAAGCCGGATTTAATGACGATGGCGAAAAGCCTCGCGGGCGGGATGCCGCTTTCGGGGGTGGTCGGTAACGCGAATATTATGGACGCACCCGCGCCGGGCGGGCTTGGCGACACCTACGCCGGTAACCCGCTGGCGGTGGCTGCCGCGCACGCGGTGCTCAACATTATCGACAAAGAATCACTCTGCGAACGCGCGAATCAACTGGGCCAGCGTCTCAAAAACACg                                                                                                        <  2:115962/250‑1 (MQ=255)
                                                                                                               ggttGCCATGGATCATTACGCCGATAAGCCGTATTGAATGGCGATGGCGAAAAGCCTCGCGGGCGGGGTGCCGCTTTCGGGCGGGGTCGGTAACGCGAATATTAGGGACGCACCCGCGCCGGGCGGGCGTGGCGACACCTACGCCGGTAACCCGCTGGCGGTGTCTGCCGCGCACGCGGTTCTCAACATTATCGACAAAGAATCACTCTGCGAGCGTGCGAATCAACTGGGCCTGCGGCTCAAAAACACg                                                                                                        <  2:505043/248‑1 (MQ=255)
                                                                                                               gTTTGCCATGGATCATTACGCCGATAAGCCGGATTTAATGACGATGGCGAAAAGCCTCGCGGGCGGGATGCCGCTTTCGGGCGTGGTCGGTAACGCGAATATTATGGACGCACCCGCGCCGGGCGGGCTTGGCGACACCTACGCCGGTAACCCGCTGGCGGTGGCTGCCGCGCACGCGGTGCTCAACATTATCGACAAAGAATCACTCTGCGAACGCGCGAATCAACTGGGCCAGCGTCTCAAAAACACGt                                                                                                       <  1:386363/251‑1 (MQ=255)
                                                                                                                         gATCATTACGCCGATAAGCCGGATTTAATGACGATGGCGAAAAGCCTCGCGGGCGGGATGCCGCTTTCGGGCGTGGTCGGTAACGCGAATATTATGGACGCACCCGCGCCGGGCGGGCTTGGCGACACCTACGCCGGTAACCCGCTGGCGGTGGCTGCCGCGCACGCGGCGCTCAACATTATCGACAAAGAATCACTCTGCGCACGCGCGAATCAACTGGGCCAGCGGCTCAAAAACACGTTGATTTATGc                                                                                             >  2:117340/1‑251 (MQ=255)
                                                                                                                            cATTACGCCGATAAGCCGGATTTAATGACGATGGCGAAAAGCCTCGCGGGCGGGATGCCGCTTTCGGGCGTGGTCGGTAACGCGAATATTATGGACGCACCCGCGCCGGGCGGGCTTGGCGACACCTACGCCGGTAACCCGCTGGCGGTGGCTGCCGCGCACGCGGTGCTCAACATTATCGACAAAGAATCACTCTGCGAACGCGCGAATCAACTGGGCCAGCGTCTCAAAAACACGTTGATTAATGCCa                                                                                           <  1:153009/250‑1 (MQ=255)
                                                                                                                             aTTACGCCGATAAGCCGGATTTAATGACGATGGCGAAAAGCCTCGCGGGCGGGATGCCGCTTTCGGGCGTGTTCGGTAACGCGAATATTATGGACGCACCCGCGCCGGGCGGGCTTGGCGACACCTACGCCGGTAACCCGCTGGCGGTCGATGCCTCGCACGCCGTGCTCAACATTATCGCCCAAGAATCACTCTGCGAACGCGCGCATCAACTAGGCCCGCCTCTCCAAACCACCTTTATTGATGccca                                                                                          >  2:247591/1‑248 (MQ=255)
                                                                                                                                aCGCCGATAAGCCGGATTTAATGACGATGGCGAAAAGCCTCGCGGGCGGGATGCCGCTTTCGGGCGTGGTCGGTAACGCGAATATTATGGACGCATCCGCGCCGGGCGGGCTTGGCGACACCTACGCCGGTAACCCGCTGGCGGTGGCTGCCGCGCACGCGGGGCTCAACATTATCGACAAAGAATCACTCTGCGAACGCGCGAATAAACTGGGCCAGCGTCCCACAAACACGTTGATTGATGCCaaagc                                                                                       >  2:453395/1‑249 (MQ=255)
                                                                                                                                 cGCCGATAAGCCGGAGTGAATGACGATGGCGAAAAGCCACGCGGGCGGGATGCCGCTTGCGGGCGTGGTCGGTAACGCTAATATTATGGACGCACCCGCGCCGGGCGGTCTTGGCGACACCTACGCCGGTAACCCGCTGGCGGGGGCTGCCGCGCACGCGGTGCTCAACATTATCGACAAAGAATCACTCGGCGAACGCGCGAATCAACTGTTCCTGCTTCTCAAATACACGTTGATTGATGCgaaagaaa                                                                                     <  2:56014/251‑1 (MQ=255)
                                                                                                                                     gATAAGCCGGATTTAATGACGATGGCGAAAAGCCTCGCGGGCGGGATGCCGCTTTCGGGCGTGGTCGGTAACGCGAATATTATGGACGCACCCGCGCCGGGCGGGCTTGGCGACACCTACGCCGGTAACCCGCTGGCGGTGGCTGCCGCGCACGCGGGGCTCAACATTATCGACAAAGAATCACTCTGCGAACGCGCGAATCAACTGGGCCAGCGTCTCAAAAACACGTTGATTGATGCCAAAGAAAGCGt                                                                                 >  1:444931/1‑251 (MQ=255)
                                                                                                                                        aaGCCGGATTTAATGACGACGGCGAAAAGCCTCGCGGGCGGGCTGCCGCATTCGGGCGTGGTAGGTAACGCGACTATTATGGACGCACCCGCGCCGGGCGGGCTTGGCGACACCTACGCCGGTAACCCGCTGGCGGTGGCTGCCGCGCACGCGGTGCTCAACCTTCTCGACCAAGAATCACTCTGCGCACGCGCGAATCAACTGGGCCAGCGTCTCACAAACACGCTGATTGATGCCAAAGAAAGCGTTcc                                                                              >  2:225157/1‑251 (MQ=255)
                                                                                                                                          gCCGGAGTTAATGACGATGGCGCAAAGACTCGCGGGCGGGCTGCCGCCTTCGGGCGTGGTCGGTAACGCGAATATTATTGACGCACCCGCGCCGGGCGGGCTTGGCGACACCTACGCCGGTAACCCGCTGGCGGTGGCCGCCGGGCCCGCGGGGCTCAACACTAACGACAAAGACTCACCCTGCGAACGCGCGGATCAACTGGGACAACGTCTCAAAAACACGGTGATTGAAGCCaaagaaagaggtcagg                                                                            >  2:310781/1‑243 (MQ=255)
                                                                                                                                                         gATGGCGAAAAGCCTCGCGGGCGGGATGCCGCTTTCGGGCGTGGTCGGTAACGCGAATATTATGGACGCACCCGCGCCGGGCGGGCTTGGCGACACCTACGCCGGTAACCCGCTGGCGGTGGCTGCCGCGCACGCGGTGCTCAACATTATCGACAAAGAATCACTCTGCGAACGCGCGAATCAACTGGGCCAGCGTCTCAAAAACACGTTGATTGATGCCAAAGAAAGCGTTCCGGCCATTGCTGCGGTAc                                                             >  1:242615/1‑251 (MQ=255)
                                                                                                                                                           tGGCGAAAAGCCTCGCGGGCGGGATGCCGCTTTCGGGCGTGGTCGGTAACGCGAATATTATGGACGCACCCGCGCCGGGCGGGCTTGGCGACACCTACGCCGGTAACCCGCTGGCGGTGGCTGCCGCGCACGCGGTGCTCAACATTATCGACAAAGAATCACTCTGCGAACGCGCGAATCAACTGGGCCAGCGTCTCAAAAACACGTTGATTGATGCCAAAGAAAGCGTTCCGGCCATTGCTGCGGTAcgc                                                           >  1:195679/1‑251 (MQ=255)
                                                                                                                                                                                  gtGCCGCTTTCGGGGGTGGTCGGTAACGCGAATAGTATGGACGCCCCGGCGCCGGGCTGGCTTGGCGACACCTACGCGGGTAACCCGCTGGCGGTGGCTGCCGCGCACGCGGTGCTAAACATTATCGACAAAGAATCACTCTGCGAACGCGCGAATCGACTGGGCCAGCGTCTCAAAAACACGTTGATTGATGCCAAAGAAAGCGTTCCGGCCATTGCTGCGGTACGCGGCCTGGGGTCGATGATTGCGGt                                    <  2:283976/250‑1 (MQ=255)
                                                                                                                                                                                  aTGCCGCTTTCGGGCGTGGTCGGTAACGCGAATATTATGGACGCACCCGCGCCGGGCGGGCTTGGCGACACCTACGCCGGTAACCCGCTGGCGGTGGCTGCCGCGCACGCGGTGCTCAACATTATCGACAAAGAATCACTCTGCGAACGCGCGAATCAACTGGGCCAGCGTCTCAAAAACACGTTGATTGATGCCAAAGAAAGCGTTCCGGCCATTGCTGCGGTACGCGGCCTGGGGTCGATGATTGCGGt                                    <  1:131716/251‑1 (MQ=255)
                                                                                                                                                                                                 gTGGTCGGTAACGCGAATATTATGGACGCACCCGCGCCGGGCGGGCTTGGCGACACCTACGCCGGTAACCCGCTGGCGGTGGCTGCCGCGCACGCGGTGCTCAACATTATCGACAAAGAATCACTCTGCGAACGCGCGAATCAACTGGGCCAGCGTCTCAAAAACACGTTGATTGATGCCAAAGAAAGCGTTCCGGCCATTGCTGCGGTACGCGGCCTGGGGTCGATGATTGCGGTAGAGTTTAACGATcc                     <  2:14171/251‑1 (MQ=255)
                                                                                                                                                                                                 gTGGTCGGTAACGCGAATATTATGGACGCACCCGCGCCGGGCGGGCTTGGCGACACCTACGCCGGTAACCCGCTGGCGGTGGCTGCCGCGCACGCGGTGCTCAACATTATCGACAAAGAATCACTCTGCGAACGCGCGAATCAACTGGGCCAGCGTCTCAAAAACACGTTGATTGATGCCAAAGAAAGCGTTCCGGCCATTGCTGCGGTACGCGGCCTGGGGTCGATGATTGCGGTAGAGTTTAACGATcc                     <  2:242615/251‑1 (MQ=255)
                                                                                                                                                                                                              cgAATATTATGGACGCACCCGCGCCGGGCGGGCTTGGCGACACCTACGCCGGTAACCCGCTGGCGGTGGCTGCCGCGCACGCGGTGCTCAAAATTATCGACCAAGAATCACTCTGCGAACGCGCGAATCAACTGGGCCAGCGTCTCAAAAACACGTTGATTGATGCCAAAGAAAGCGTTCCGGCCATTGCTGCGGCACGCGGCCTGGGGTCGATGCTTGCGGTAGAGTTTAACGATCCGCAAACGGGCg          >  1:21456/1‑249 (MQ=255)
                                                                                                                                                                                                                     tagggACGCACCCGCGCCGGGCGGGCTTGGCGACACCTACGCCGGTAACCCGCTGGCGGTGGCTGCCGCGCACGCGGTGCTCAACATTATCGACAAAGAATCACTCTGCGAACGCGCGAATCAACTGGGCCAGCGTCTCAAAAACACGTTGATTGATGCCAAAGAAAGCGTGCCGGCCATTGCGGCGGTACGCGGCCTGGGGTCGATGATTGCGGTAGAGTTTAACGATCCGCAAACGGGCGAGCCGTCa  <  2:159989/247‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                         
CAACGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTGCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCTGTTTGCCATGGATCATTACGCCGATAAGCCGGATTTAATGACGATGGCGAAAAGCCTCGCGGGCGGGATGCCGCTTTCGGGCGTGGTCGGTAACGCGAATATTATGGACGCACCCGCGCCGGGCGGGCTTGGCGGCACCTACGCCGGTAACCCGCTGGCGGTGGCTGCCGCGCACGCGGTGCTCAACATTATCGACAAAGAATCACTCTGCGAACGCGCGAATCAACTGGGCCAGCGTCTCAAAAACACGTTGATTGATGCCAAAGAAAGCGTTCCGGCCATTGCTGCGGTACGCGGCCTGGGGTCGATGATTGCGGTAGAGTTTAACGATCCGCAAACGGGCGAGCCGTCA  >  NC_000913/1366188‑1366650

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