Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 634,722 A→G A379A (GCA→GCG ybdL → methionine‑oxo‑acid transaminase, PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913634,7220AG92.9% 78.4 / ‑1.7 28A379A (GCA→GCGybdLmethionine‑oxo‑acid transaminase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (13/13);  total (14/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 = 9.78e-01

TGTTGAGTTTTGCCAGTGGCTGACGCAGGAGCACGGCGTAGCGGCGATTCCGCTGTCGGTGTTTTGCGCCGATCCCTTCCCACATAAACTGATTCGTCTCTGTTTTGCCAAGAAGGAATCGACGTTGCTGGCAGCAGCTGAACGCCTGCGCCAGCTTTAGCTATTTAACCGTCCAGGCGCGGGAATATTGACGGTCGGCAAACAACTCTTGCAGGCCGTTAATTTGCTTCAGGCGCAGCACTTCGTCGCTGTCCATTCCCAGCTCTTTGC  >  NC_000913/634587‑634856
                                                                                                                                       |                                                                                                                                      
tgTTGAGTTTTGCCAGTGGCTGACGCAGGAGCACGGCGTAGCGGCGATTCCGCTGTCGGTGTTTTGCGCCGATCCCTTCCCACATAAACTGATTCGTCTCTGTTTTGCCAAGAAGGAATCGACGTTGCTGGCAGCGGCtt                                                                                                                                    >  2:480128/1‑139 (MQ=255)
tgTTGAGTTTTGCCAGTGGCTGACGCAGGAGCACGGCGTAGCGGCGATTCCGCTGTCGGTGTTTTGCGCCGATCCCTTCCCACATAAACTGATTCGTCTCTGTTTTGCCAAGAAGGAATCGACGTTGCTGGCAGCAGCTg                                                                                                                                    >  1:533825/1‑140 (MQ=255)
     aGTTTTGCCAGTGGCTGACGCAGGAGCACGGCGTAGCGGCGATTCCGCTGTCGGTGTTTTGCGCCGATCCCTTCCCACATAAACTGATTCGTCTCTGTTTTGCCAAGAAGGAATCGACGTTGCTGGCAGCGGCTGAACGc                                                                                                                               >  1:310996/1‑140 (MQ=255)
     aGTTTTGCCAGTGGCTGACGCAGGAGCACGGCGTAGCGGCGATTCCGCTGTCGGTGTTTTGCGCCGATCCCTTCCCACATAAACTGATTCGTCTCTGTTTTGCCAAGAAGGAATCGACGTTGCTGGCAGCGGCTGAACGc                                                                                                                               >  2:282957/1‑140 (MQ=255)
      gTTTTGCCAGTGGCTGACGCAGGAGCACGGCGTAGCGGCGATTCCGCTGTCGGTGTTTTGCGCCGATCCCTTCCCACATAAACTGATTCGTCTCTGTTTTGCCAAGAAGGAATCGACGTTGCTGGCAGCGGCTGAACGcc                                                                                                                              >  2:248057/1‑140 (MQ=255)
           gCCAGTGGCTGACGCAGGAGCACGGCGTAGCGGCGATTCCGCTGTCGGTGTTTTGCGCCGATCCCTTCCCACATAAACTGATTCGTCTCTGTTTTGCCAAGAAGGAATCGACGTTGCTGGCAGCGGCTGAACGCCTgcgc                                                                                                                         <  2:76644/140‑1 (MQ=255)
           gCCAGTGGCTGACGCAGGAGCACGGCGTAGCGGCGATTCCGCTGTCGGTGTTTTGCGCCGATCCCTTCCCACATAAACTGATTCGTCTCTGTTTTGCCAAGAAGGAATCGACGTTGCTGGCAGCGGCTGAACGCCTgcgc                                                                                                                         <  1:454325/140‑1 (MQ=255)
                tGGCTGACCCAGGAGCACGGCGTAGCGGCGATTCCGCTGTCGGTGTTTTGCGCCGATCCCTTCCCACATAAACTGATTCGTCTCTGTTTTGCCAAGAAGGAATCGACGTTGCTGGCAGCGGCTGAACGCCTGCGCCAGCt                                                                                                                    <  2:210614/140‑1 (MQ=255)
                  gCTGACGCAGGAGCACGGCGTAGCGGCGATTCCGCTGTCGGTGTTTTGCGCCGATCCCTTCCCACATAAACTGATTCGTCTCTGTTTTGCCAAGAAGGAATCGACGTTGCTGGCAGCGGCTGAACGCCTGCGCCAGCttt                                                                                                                  <  1:282957/140‑1 (MQ=255)
                  gCTGACGCAGGAGCACGGCGTAGCGGCGATTCCGCTGTCGGTGTTTTGCGCCGATCCCTTCCCACATAAACTGATTCGTCTCTGTTTTGCCAAGAAGGAATCGACGTTGCTGGCAGCGGCTGAACGCCTGCGCCAGCttt                                                                                                                  <  2:451874/140‑1 (MQ=255)
                  gCTGACGCAGGAGCACGGCGTAGCGGCGATTCCGCTGTCGGTGTTTTGCGCCGATCCCTTCCCACATAAACTGATTCGTCTCTGTTTTGCCAAGAAGGAATCGACGTTGCTGGCAGCGGCTGAACGCCTGCGCCAGCttt                                                                                                                  <  2:746532/140‑1 (MQ=255)
                    tGACGCAGGAGCACGGCGTAGCGGCGATTCCGCTGTCGGTGTTTTGCGCCGATCCCTTCCCACATAAACTGATTCGTCTCTGTTTTGCCAAGAAGGAATCGACGTTGCTGGCAGCGGCTGAACGCCTGCGCCAGCTTTAg                                                                                                                >  2:56810/1‑140 (MQ=255)
                                aCGGCGTAGCGGCGATTCCGCTGTCGGTGTTTTGCGCCGATCCCTTCCCACATAAACTGATTCGTCTCTGTTTTGCCAAGAAGGAATCGACGTTGCTGGCAGCGGCTGAACGCCTGCGCCAGCTTTAGCTATTTAACCGt                                                                                                    >  1:716381/1‑140 (MQ=255)
                                aCGGCGTAGCGGCGATTCCGCTGTCGGTGGTTTGCGCCGATCCCTTCCCACATAAACTGATTCGTCTCTGTTTTGCCAAGAAGGAATCGACGTTGCTGGCAGCGGCTGAACGCCTGCGCCAGCTTTAGCTATTTAACCGt                                                                                                    >  2:199386/1‑140 (MQ=255)
                                 cGGCGTAGCGGCGATTCCGCTGTCGGTGTTTTGCGCCGATCCCTTCCCACATAAACTGATTCGTCTCTGTTTTGCCAAGAAGGAATCGACGTTGCTGGCAGCGGCTGAACGCCTGCGCCAGCTTTAGCTATTTAACCGTc                                                                                                   >  1:105095/1‑140 (MQ=255)
                                      tAGCGGCGATTCCGCTGTCGGTGTTTTGCGCCGATCCCTTCCCACATAAACTGATTCGTCTCTGTTTTGCCAAGAAGGAATCGACGTTGCTGGCAGCGGCTGAACGCCTGCGCCAGCTTTAGCTATTTAACCGTCCAGgc                                                                                              >  2:17285/1‑140 (MQ=255)
                                       aGCGGCGATTCCGCTGTCGGTGTTTTGCGCCGATCCCTTCCCACATAAACTGATTCGTCTCTGTTTTGCCAAGAAGGAATCGACGTTGCTGGCAGCGGCTGAACGCCTGCGCCAGCTTTAGCTATTTAACCGTCCAGgcg                                                                                             <  2:698353/140‑1 (MQ=255)
                                         cggcgATTCCGCTGTCGGTGTTTTGCGCCGATCCCTTCCCACATAAACTGATTCGTCTCTGTTTTGCCAAGAAGGAATCGACGTTGCTGGCAGCGGCTGAACGCCTGCGCCAGCTTTAGCTATTTAACCGTCCAGGcgcg                                                                                           <  1:56810/140‑1 (MQ=255)
                                          ggcgATTCCGCTGTCGGTGTTTTGCGCCGATCCCTTCCCACATAAACTGATTCGTCTCTGTTTTGCCAAGAAGGAATCGACGTTGCTGGCAGCGGCTGAACGCCTGCGCCAGCTTTAGCTATTTAACCGTCCAGGCGCgg                                                                                          <  2:310996/140‑1 (MQ=255)
                                                    cTGTCGGTGTTTTGCGCCGATCCCTTCCCACATAAACTGATTCGTCTCTGTTTTGCCAAGAAGGAATCGACGTTGCTGGCAGCGGCTGAACGCCTGCGCCAGCTTTAGCTATTTAACCGTCCAGGCGCGGGAATATTGAc                                                                                >  1:93368/1‑140 (MQ=255)
                                                        cGGTGTTTTGCGCCGATCCCTTCCCACATAAACTGATTCGTCTCTGTTTTGCCAAGAAGGAATCGACGTTGCTGGCAGCGGCTGAACGCCTGCGCCAGCTTTAGCTATTTAACCGTCCAGGCGCGGGAATATTGACGGTc                                                                            <  1:541633/140‑1 (MQ=255)
                                                         ggTGTTTTGCGCCGATCCCTTCCCACATAAACTGATTCGTCTCTGTTTTGCCAAGAAGGAATCGACGTTGCTGGCAGCGGCTGAACGCCTACGCCAGCTTTAGCTATTTAACCGTCCAGGCGCGGGAATATTGACGGTCg                                                                           <  1:289728/140‑1 (MQ=255)
                                                           tgtTTTGCGCCGATCCCTTCCCACATAACTGATGCGGCTCTGTTTTGCCAAGAAGGAATCGACGTTGCTGGCAGCGGCTGAACGCCTGCGCCAGCTTTAGCTAGTTAACCGTCCAGGCGCGGGAATATTGACGGTCGGCa                                                                        <  2:304039/140‑1 (MQ=255)
                                                                 gcgcCGATCCCTTCCCACATAAACTGATTCGTCTCTGTTTTGCCAAGAAGGAATCGACGTTGCTGGCAGCAGCTGAACGCCTGCGCCAGCTTTAGCTATTTAACCGTCCAGGCGCGGGAATATTGACGGTCGGCAaacaa                                                                   <  2:533825/140‑1 (MQ=255)
                                                                            ttCCCACATAAACTGATTCGTCTCTGTTTTGCCAAGAAGGAATCGACGTTGCTGGCAGCGGCTGAACGCCTGCGCCAGCTTTAGCTATTTAACCGTCCAGGCGCGGGAATATTGACGGTCGGCAAACAACTCTTGCAGGc                                                        >  1:727696/1‑140 (MQ=255)
                                                                                   aTAAACTGATTCGTCTCTGTTTTGCCAAGAAGGAATCGACGTTGCTGGCAGCGGCTGAACGCCTGCGCCAGCTTTAGCTATTTAACCGTCCAGGCGCGGGAATATTGACGGTCGGCAAACAACTCTTGCAGGCCGTTAAt                                                 >  2:98564/1‑140 (MQ=255)
                                                                                                           ccAAGAAGGAATCGACGTTGCTGGCAGCGGCTGAACGCCTGCGCCAGCTTTAGCTATTTAACCGTCCAGGCGCGGGAATATTGACGGTCGGCAAACAACTCTTGCAGGCCGTTAATTTGCTTCAGGCGCAGCACTtcgtc                         >  2:471774/1‑140 (MQ=255)
                                                                                                             aagaagGAATCGACGTTGCTGGCAGCGGCTGAACGCCTGCGCCAGCTTTAGCTATTTAACCGTCCAGGCGCGGGAATATTGACGGTCGGCAAACAACTCTTGCAGGCCGTTAATTTGCTTCAGGCGCAGCACTTCGTCGc                       <  2:716381/140‑1 (MQ=255)
                                                                                                              agaagGAATCGACGTTGCTGGCAGCGGCTGAACGCCTGCGCCAGCTTTAGCTATTTAACCGTCCAGGCGCGGGAATATTGACGGTCGGCAAACAACTCTTGCAGGCCGTTAATTTGCTTCAGGCGCAGCACTTCGTCGCt                      >  2:108883/1‑140 (MQ=255)
                                                                                                                  ggAATCGACGTTGCTGGCAGCGGCTGAACGCCTGCGCCAGCTTTAGCTATTTAACCGTCCAGGCGCGGGAATATTGACGGTCGGCAAACAACTCTTGCAGGCCGTTAATTTGCTTCAGGCGCAGCACTTCGTCGCTGTcc                  <  1:17285/140‑1 (MQ=255)
                                                                                                                             tGCTGGCAGCGGCTGAACGCCTGCGCCAGCTTTAGCTATTTAACCGTCCAGGCGCGGGAATATTGACGGTCGGCAAACAACTCTTGCAGGCCGTTAATTTGCTTCAGGCGCAGCACTTCGTCGCTGTCCATTCCCAGctc       >  2:664278/1‑140 (MQ=255)
                                                                                                                                 ggCAGCGGCTGAACGCCTGCGCAAGCTTTAGCTATTTAACCGTCCAGGCGCGGGAATATTGACGGTCTGCAAACAACTCTTGCAGGCCGTTAATTTGCTTCAGGCGCAGCACTTCGTCGCTGTCCATTCCCAGCGCTTTg   <  2:727696/140‑1 (MQ=255)
                                                                                                                                  gcagcgGCTGAACGCCTGCGCCAGCTTTAGCTATTTAACCGTCCAGGCGCGGGAATATTGACGGTCGGCAAACAACTCTTGCAGGCCGTTAATTTGCTTCAGGCGCAGCACTTCGTCGCTGTCCATTCCCAGCTCTTTGc  <  2:301079/140‑1 (MQ=255)
                                                                                                                                  gcagcgGCTGAACGCCTGCGCCAGCTTTAGCTATTTAACCGTCCAGGCGCGGGAATATTGACGGTCGGCAAACAACTCTTGCAGGCCGTTAATTTGCTTCAGGCGCAGCACTTCGTCGCTGTCCATTCCCAGCTCTTTGc  <  1:772882/140‑1 (MQ=255)
                                                                                                                                  gcagcgGCTGAACGCCTGCGCCAGCTTTAGCTATTTAACCGTCCAGGCGCGGGAATATTGACGGTCGGCAAACAACTCTTGCAGGCCGTTAATTTGCTTCAGGCGCAGCACTTCGTCGCTGTCCATTCCCAGCTCTTTGc  <  1:108883/140‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                      
TGTTGAGTTTTGCCAGTGGCTGACGCAGGAGCACGGCGTAGCGGCGATTCCGCTGTCGGTGTTTTGCGCCGATCCCTTCCCACATAAACTGATTCGTCTCTGTTTTGCCAAGAAGGAATCGACGTTGCTGGCAGCAGCTGAACGCCTGCGCCAGCTTTAGCTATTTAACCGTCCAGGCGCGGGAATATTGACGGTCGGCAAACAACTCTTGCAGGCCGTTAATTTGCTTCAGGCGCAGCACTTCGTCGCTGTCCATTCCCAGCTCTTTGC  >  NC_000913/634587‑634856

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