Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,272,271 Δ1 bp coding (481/885 nt) garR ← tartronate semialdehyde reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,272,2680C.97.1% 159.8 / ‑1.6 35coding (484/885 nt)garRtartronate semialdehyde reductase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base . (16/18);  total (16/19)
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.16e-01

CAGTCCACCGCGAATTGCCTGATAAACCAGGTCCGGGTTAACGCCCGCTTTAGTTGCCAGCGTTAACGCTTCTGACATCGCGGCAATATTCAGCGCCACAATGACCTGATTTGCCAGTTTGGTGACGTTACCTGCACCGATTTCCCCGGTATGCACCACGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGTCGAAAATAGCCTTGTCGCCGCCCACCATCACTGACAGCGTACCGTCGATGGCTTTCGGTTCACCGCCGCTCACCGGAGCATCCAGCAT  >  NC_000913/3272125‑3272409
                                                                                                                                               |                                                                                                                                             
cAGTCCACCGCGAATTGCCTGATAAACCAGGTCCGGGTTAACGCCCGCTTTAGTTGCCAGCGTTAACGCTTCTGACATCGCGGCAATATTCAGCGCCACAATGACCTGATTTGCCAGTTTGGTGACGTTACCTGCACCGATTTCCCGGTAt                                                                                                                                       <  2:360336/151‑1 (MQ=255)
   tCCACCGCGAATTGCCTGATAAACCAGGTCCGGGTTAACGCCCGCTTTAGTTGCCAGCGTTAACGCTTCTGACATCGCGGCAATATTCAGCGCCACAATGACCTGATTTGCCAGTTTGGTGACGTTACCTGCACCGATTTCCCGGTATGc                                                                                                                                     >  1:118737/1‑150 (MQ=255)
      aCCGCGAATTGCCTGATAAACCAGGTCCGGGTTAACGCCCGCTTTAGTTGCCAGCGTTAACGCTTCTGACATCGCGGCAATATTCAGCGCCACAATGACCTGATTTGCCAGTTTGGTGACGTTACCTGCACCGATTTCCCGGTATGcacca                                                                                                                                 <  2:460653/151‑1 (MQ=255)
           gAATTGCCTGATAAACCAGGTCCGGGTTAACGCCCGCTTTAGTTGCCAGCGTTAACGCTTCTGACATCGCGGCAATATTCAGCGCCACAATGACCTGATTTGCCAGTTTGGTGACGTTACCTGCACCGATTTCCCGGTATGCACCACGGaa                                                                                                                            <  1:457515/151‑1 (MQ=255)
                    gATAAACCAGGTCCGGGTTAACGCCCGCTTTAGTTGCCAGCGTTAACGCTTCTGACATCGCGGCAATATTCAGCGCCACAATGACCTGATTTGCCAGTTTGGTGACGTTACCTGCACCGATTTCCCGGTATGCACCACGGAACCCGCCa                                                                                                                     <  2:260948/149‑1 (MQ=255)
                    gATAAACCAGGTCCGGGTTAACGCCCGCTTTAGTTGCCAGCGTTAACGCTTCTGACATCGCGGCAATATTCAGCGCCACAATGACCTGATTTGCCAGTTTGGTGACGTTACCTGCACCGATTTCCCGGTATGCACCACGGAACCCGCCa                                                                                                                     <  2:370045/149‑1 (MQ=255)
                           cAGGTCCGGGTTAACGCCCGCTTTAGTTGCCAGCGTTAACGCTTCTGACATCGCGGCAATATTCAGCGCCACAATGACCTGATTTGCCAGTTTGGTGACGTTACCTGCACCGATTTCCCGGTATGCACCACGGAACCCGCCATCGCTTtca                                                                                                            <  2:391408/151‑1 (MQ=255)
                                ccGGGTTAACGCCCGCTTTAGTTGCCAGCGTTAACGCTTCTGACATCGCGGCAATATTCAGCGCCACAATGACCTGATTTGCCAGTTTGGTGACGTTACCTGCACCGATTTCCCGGTATGCACCACGGAACCCGCCATCGCTTTCATCaa                                                                                                        <  2:554670/150‑1 (MQ=255)
                                    gTTAACGCCCGCTTTAGTTGCCAGCGTTAACGCTTCTGACATCGCGGCAATATTCAGCGCCACAATGACCTGATTTGCCAGTTTGGTGACGTTACCTGCACCGATTTCCCGGTATGCACCACGGAACCCGCCATCGCTTTCATCAAATCAt                                                                                                   >  2:546610/1‑151 (MQ=255)
                                      tAACGCCCGCTTTAGTTGCCAGCGTTAACGCTTCTGACATCGCGGCAATATTCAGCGCCACAATGACCTGATTTGCCAGTTTGGTGACGTTACCTGCACCGATTTCCCGGTATGCACCACGGAACCCGCCATCGCTTTCATCAAATCAt                                                                                                   >  2:390209/1‑149 (MQ=255)
                                       aaCGCCCGCTTTAGTTGCCAGCGTTAACGCTTCTGACATCGCGGCAATATTCAGCGCCACAATGACCTGATTTGCCAGTTTGGTGACGTTACCTGCACCGATTTCCCGGTATGCACCACGGAACCCGCCATCGCTTTCATCAAATCATAGt                                                                                                <  1:146676/151‑1 (MQ=255)
                                          gCCCGCTTTAGTTGCCAGCGTTAACGCTTCTGACATCGCGGCAATATTCAGCGCCACAATGACCTGATTTGCCAGTTTGGTGACGTTACCTGCACCGATTTCCCGGTATGCACCACGGAACCCGCCATCGCTTTCATCAAATCATAGTAtt                                                                                             <  2:505691/151‑1 (MQ=255)
                                                  tAGTTGCCAGCGTTAACGCTTCTGACATCGCGGCAATATTCAGCGCCACAATGACCTGATTTGCCAGTTTGGTGACGTTACCTGCACCGATTTCCCGGTATGCACCACGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGTCGaa                                                                                      <  1:73219/150‑1 (MQ=255)
                                                   aGTTGCCAGCGTTAACGCTTCTGACATCGCGGCAATATTCAGCGCCACAATGACCTGATTTGCCAGTTTGGTGACGTTACCTGCACCGATTTCCCGGTATGCACCACGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGTCGaaa                                                                                     >  2:443101/1‑150 (MQ=255)
                                                      tGCCAGCGTTAACGCTTCTGACATCGCGGCAATATTCAGCGCCACAATGACCTGATTTGCCAGTTTGGTGACGTTACCTGCACCGATTTCCCCGGTATGCACCACGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGTCGAAAAt                                                                                   <  2:473425/150‑1 (MQ=255)
                                                      tGCCAGCGTTAACGCTTCTGACATCGCGGCAATATTCAGCGCCACAATGACCTGATTTGCCAGTTTGGTGACGTTACCTGCACCGATTTCCCGGTATGCACCACGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGTCGAAAATa                                                                                  <  2:242911/150‑1 (MQ=255)
                                                       gCCAGCGTTAACGCTTCTGACATCGCGGCAATATTCAGCGCCACAATGACCTGATTTGCCAGTTTGGTGACGTTACCTGCACCGATTTCCCGGTATGCACCACGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGTCGAAAATAGc                                                                                <  2:599489/151‑1 (MQ=255)
                                                           gCGTTAACGCTTCTGACATCGCGGCAATATTCAGCGCCACAATGACCTGATTTGCCAGTTTGGTGACGTTACCTGCACCGATTTCCCGGTATGCACCACGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGTCGAAAATAGCCTTg                                                                            >  2:362003/1‑151 (MQ=255)
                                                            cGTTAACGCTTCTGACATCGCGGCAATATTCAGCGCCACAATGACCTGATTTGCCAGTTTGGTGACGTTACCTGCACCGATTTCCCGGTATGCACCACGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGTCGAAAATAGCCTTGt                                                                           <  1:221193/151‑1 (MQ=255)
                                                                      tCTGACATCGCGGCAATATTCAGCGCCACAATGACCTGATTTGCCAGTTTGGTGACGTTACCTGCACCGATTTCCCGGTATGCACCACGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGTCGAAAATAGCCTTGTCGCCGCcca                                                                  <  2:136251/150‑1 (MQ=255)
                                                                        tGACATCGCGGCAATATTCAGCGCCACAATGACCTGATTTGCCAGTTTGGTGACGTTACCTGCACCGATTTCCCGGTATGCACCACGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGTCGAAAATAGCCTTGTCGCCGCccacc                                                                >  2:309067/1‑150 (MQ=255)
                                                                               gcgGCAATATTCAGCGCCACAATGACCTGATTTGCCAGTTTGGTGACGTTACCTGCACCGATTTCCCGGTATGCACCACGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGTCGAAAATAGCCTTGTCGCCGCCCACCATCACTg                                                         >  2:299968/1‑150 (MQ=255)
                                                                                        ttCAGCGCCACAATGACCTGATTTGCCAGTTTGGTGACGTTACCTGCACCGATTTCCCGGTATGCACCACGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGTCGAAAATAGCCTTGTCGCCGCCCACCATCACTGACAGCGTAcc                                               >  1:473686/1‑151 (MQ=255)
                                                                                            gcgcCACAATGACCTGATTTGCCAGTTTGGTGACGTTACCTGCACCGATTTCCCGGTATGCACCACGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGTCGAAAATAGCCTTGTCGCCGCCCACCATCACTGACAGCGTACCGTCg                                           <  2:108680/151‑1 (MQ=255)
                                                                                                cacaATGACCTGATTTGCCAGTTTGGTGACGTTACCTGCACCGATTTCCCGGTATGCACCACGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGTCGAAAATAGCCTTGTCGCCGCCCACCATCACTGACAGCGTACCGTCGATgg                                       >  1:112123/1‑151 (MQ=255)
                                                                                                cacaATGACCTGATTTGCCAGTTTGGTGACGTTACCTGCACCGATTTCCCGGTATGCACCACGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGTCGAAAATAGCCTTGTCGCCGCCCACCATCACTGACAGCGTACCGTCGATgg                                       >  2:59007/1‑151 (MQ=255)
                                                                                                 acaATGACCTGATTTGCCAGTTTGGTGACGTTACCTGCACCGATTTCCCGGTATGCACCACGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGTCGAAAATAGCCTTGTCGCCGCCCACCATCACTGACAGCGTACCGTCGATGGc                                      >  1:270012/1‑151 (MQ=255)
                                                                                                  caATGACCTGATTTGCCAGTTGGGTGACGTTACCTGCACCGATTTCCCGGTATGCACCACGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGTCGAAAATAGCCTTGTCGCCGCCCACCATCACTGACAGCGTACCGTCGATGGCt                                     <  2:59140/151‑1 (MQ=255)
                                                                                                      gACCTGATTTGCCAGTTTGGTGACGTTACCTGCACCGATTTCCCGGTATGCACCACGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGTCGAAAATAGCCTTGTCGCCGCCCACCATCACTGACAGCGTACCGTCGATGGCTTTCg                                 >  1:456547/1‑151 (MQ=255)
                                                                                                               tGCCAGTTTGGTGACGTTACCTGCACCGATTTCCCGGTATGCACCACGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGTCGAAAATAGCCTTGTCGCCGCCCACCATCACTGACAGCGTACCGTCGATGGCTTTCGGTTCAccgc                        >  2:528868/1‑151 (MQ=255)
                                                                                                                      ttGGTGACGTTACCTGCACCGATTTCCCGGTATGCACCACGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGTCGAAAATAGCCTTGTCGCCGCCCACCATCACTGACAGCGTACCGTCGATGGCTTTCGGTTCACCGCCGCTCAc                 >  1:91471/1‑151 (MQ=255)
                                                                                                                                 aCCTGCACCGATTTCCCGGTATGCACCACGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGTCGAAAATAGCCTTGTCGCCGCCCACCATCACTGACAGCGTACCGTCGATGGCTTTCGGTTCACCGCCGCTCACCGGAGCATcc       <  2:112123/150‑1 (MQ=255)
                                                                                                                                  ccTGCACCGATTTCCCGGTATGCACCACGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGTCGAAAATAGCCTTGTCGCCGCCCACCATCACTGACAGCGTACCGTCGATGGCTTTCGGTTCACCGCCGCTCACCGGAGCATcc       >  1:264098/1‑149 (MQ=255)
                                                                                                                                  ccTGCACCGATTTCCCGGTATGCACCACGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGTCGAAAATAGCCTTGTCGCCGCCCACCATCACTGACAGCGTACCGTCGATGGCTTTCGGTTCACCGCCGCTCACCGGAGCATCCa      <  1:270259/150‑1 (MQ=255)
                                                                                                                                     gCACCGATTTCCCGGTATGCACCACGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGTCGAAAATAGCCTTGTCGCCGCCCACCATCACTGACAGCGTACCGTCGATGGCTTTCGGTTCACCGCCGCTCACCGGAGCATCCAGCat  >  1:206972/1‑151 (MQ=255)
                                                                                                                                               |                                                                                                                                             
CAGTCCACCGCGAATTGCCTGATAAACCAGGTCCGGGTTAACGCCCGCTTTAGTTGCCAGCGTTAACGCTTCTGACATCGCGGCAATATTCAGCGCCACAATGACCTGATTTGCCAGTTTGGTGACGTTACCTGCACCGATTTCCCCGGTATGCACCACGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGTCGAAAATAGCCTTGTCGCCGCCCACCATCACTGACAGCGTACCGTCGATGGCTTTCGGTTCACCGCCGCTCACCGGAGCATCCAGCAT  >  NC_000913/3272125‑3272409

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