Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,889,791 A→T 17.1% K14* (AAG→TAG) ‡ glxR → 2‑hydroxy‑3‑oxopropionate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,889,7910AT17.1% 49.9 / 5.5 24K14* (AAG→TAG) ‡glxR2‑hydroxy‑3‑oxopropionate reductase
Reads supporting (aligned to +/- strand):  ref base A (10/9);  new base T (2/2);  total (13/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.48e-01

GGCAAAGCCGGTGCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGCCCATGGCTCAGAACCTGCAAAAGGCAGGTCACAGCCTGTTCATTTCCACCCACCACGACGCCGCGCCAGCCGACCTGATCGCTGCCGGTGCCGTGGCCCTGGCCAACCCGAAAGAGGTTGCCCAGGAAGCCGAATTCATCATCGTC  >  NC_002947/4889645‑4889940
                                                                                                                                                  |                                                                                                                                                     
ggCAAAGCCGGTGCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGc                                                                                                                                                    <  2:359577/150‑1 (MQ=255)
           tGCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGCCCATGGCTCAg                                                                                                                                         <  1:363614/150‑1 (MQ=255)
                                          gAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGCCCATGGCTCAGAACCTGCAAAAGGCAGGTCACAGCCTGTTCa                                                                                                          >  2:458902/1‑150 (MQ=255)
                                              ccGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCGCATGGCTAAAATCGGTTTCATCGGCACCCGACTAATACGGCTTCCCATTGCTCAGAACCTGCAAAAGGCAGGTCACAGCCTGTTCATTTc                                                                                                      >  2:173635/1‑150 (MQ=255)
                                                gCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGCCCATGGCTCAGAACCTGCAAAAGGCAg                                                                                                                         >  1:173786/1‑129 (MQ=255)
                                                gCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGCCCATGGCTCAGAACCTGCAAAAGGCAg                                                                                                                         <  2:173786/129‑1 (MQ=255)
                                                gCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGCCCATGGCTCAGAACCTGCAAAAGGAAGGTCACAGCCTGTTCATTTcca                                                                                                    >  1:62652/1‑150 (MQ=255)
                                                     cacaGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCTTGCCCATgatg                                                                                                                                            >  2:324230/1‑102 (MQ=255)
                                                     cacaGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCTTGCCCATgatg                                                                                                                                            <  1:324230/105‑4 (MQ=255)
                                                     cacaGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGCCCATGGCTCAGAACCTGCAAAAGGCAGGTCACAGCCTGTTCATTTccacccac                                                                                               <  2:103628/150‑1 (MQ=255)
                                                              ccAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGCCCATGGCTCAGAACCTGCAAAAGGCAGGTCACAGCCTGTTCATTTCCACCCACCACGAcgcc                                                                                      >  2:600493/1‑150 (MQ=255)
                                                                                 tCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCACCATGGGCTTGCCCATga                                                                                                                                              <  2:250242/75‑2 (MQ=255)
                                                                                 tCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCACCATGGGCTTGCCCATga                                                                                                                                              >  1:250242/1‑74 (MQ=255)
                                                                                      cttACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGCCCATGGCTCAGAACCTGCAAAAGGCAGGTCACAGCCTGTTCATTTCCACCCACCCCGACGCCGCGCCAGCCGACCTGATCGCTGcc                                                              <  1:173635/149‑1 (MQ=255)
                                                                                                           aTGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGCCCATGGCTCAGAACCTGCAAAAGTCAGGTCACAGCCTGTTCATTTCCACCCACCTCGACGCCGCGCCAGCCGACCTGATCGCTGCCGGTGCCGTGGCCCTGGCCAAc                                         <  2:62652/150‑1 (MQ=255)
                                                                                                                tAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGCCCATGGCTCCGCAAATGCAACAAGCAGGTCACAGCCTGTTCATTTCCACCCACCACGACGCCGCGCCAGCCGACCTGATCGCTGCCGCTGCCGTGGCCCTGGCCAACCCGaa                                    >  2:152354/1‑150 (MQ=255)
                                                                                                                tAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGCCCATGGCTCAGAACCTGCAAAAGGCAGGTCACAGCCTGTTCATTTCCACCCACCACGACGCCGCGCCAGCCGACCTGATCGCTGCCGGTGCCGTGGCCCTGGCCAACCCGaa                                    >  2:152678/1‑150 (MQ=255)
                                                                                                                     tCGGTTTCATCGGCACCGGCATCATGGGCAAGCCCATGGCTCAGAACCTg                                                                                                                                   <  1:277999/50‑1 (MQ=255)
                                                                                                                     tCGGTTTCATCGGCACCGGCATCATGGGCAAGCCCATGGCTCAGAACCTg                                                                                                                                   >  2:277999/1‑50 (MQ=255)
                                                                                                                      cGGTTTCATCGGCACCGGCATCATGGGCAAGCCCATGGCTCAGAACCTGCAAAAGGCAGGTCACAGCCTGTTCATTTCCACCCACCACGACGCCGCGCCAGCCGACCTGATCGCTGCCGGTGCCGTGGCCCTGGCCAACCCGAAAGAGGt                              <  1:152678/150‑1 (MQ=255)
                                                                                                                             aTCGGCACCGGCATCATGGGCAAGCCCATGGCTCAGAACCTGCAAAAGGCAGGTCACAGCCTGTTCATTTCCACCCACCACGACGCCGCGCCAGCCGACCTGATCGCTGCCGGTGCCGTGGCCCTGGCCAACCCGAAAGAGGTTGCCCAg                       >  1:123144/1‑150 (MQ=255)
                                                                                                                                    ccGGCATCATGGGCAAGCCCATGGCTCAGAACCTGCAAAAGGCAGGTCACCGCCTGTTCATTTCCACCCACCACGACGCCGCGCCAGCCGACCTGa                                                                      <  2:87821/96‑1 (MQ=255)
                                                                                                                                    ccGGCATCATGGGCAAGCCCATGGCTCAGAACCTGCAAAAGGCAGGTCACAGCCTGTTCATTTCCACCCACCACGACGCCGCGCCAGCCGACCTGa                                                                      >  1:87821/1‑96 (MQ=255)
                                                                                                                                     cGGCATCATGGGCAAGCCCATGGCTCAGAACCTGCAAAAGGCAGGTCACAGCCTGTTCATTTCCACCCACCACAACGCCGCGCCAGCCGACCTGATCGCTGCCGGTGCCGTGGCCCTGGCCAACCCGAAAGAGGTTGCCCAGGAAGCCGa               >  1:41824/1‑150 (MQ=255)
                                                                                                                                                  aaGCCCATGGCTCAGAACCTGCAAAAGGCAGGTCACAGCCTGTTCATTTCCACCCACCACGACGCCGCGCCAGCCGACCTGATCGCTGCCGGTGCCGTGGCCCTGGCCAACCCGAAAGAGGTTGCCCAGGAAGCCGAATTCATCATCGTc  >  2:123816/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GGCAAAGCCGGTGCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGCCCATGGCTCAGAACCTGCAAAAGGCAGGTCACAGCCTGTTCATTTCCACCCACCACGACGCCGCGCCAGCCGACCTGATCGCTGCCGGTGCCGTGGCCCTGGCCAACCCGAAAGAGGTTGCCCAGGAAGCCGAATTCATCATCGTC  >  NC_002947/4889645‑4889940

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