Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,271,825 T→A 8.7% intergenic (‑55/+42) garK ← / ← garR glycerate 2‑kinase 1/tartronate semialdehyde reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,271,8250TA 8.7% 114.8 / 4.1 47intergenic (‑55/+42)garK/garRglycerate 2‑kinase 1/tartronate semialdehyde reductase
Reads supporting (aligned to +/- strand):  ref base T (21/21);  new base A (2/2);  total (24/23)
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.91e-01

TCCTTTTTCTATCGCCTGCGCAACCTCGCTGGCAGATAAACTTTCTTTATAAGAGTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAGGGCGCTATGATCCGCCGTTCCTAAACCATCTGCTCGCAGTGCCTGCATCATC  >  NC_000913/3271693‑3271966
                                                                                                                                    |                                                                                                                                             
tctctttgCTAGCGCCTGCGCAACTTCGCTGGCAGATAAACTTTCTTTATAAGAGTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGAATCCAGGccga                                                                                                                                     <  2:523344/139‑3 (MQ=255)
tCCTTTTTCTATCGCCTGCGCAACCTCGCTGGCAGATAAACTTTCTTTATAAGAGTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCg                                                                                                                                     <  2:180025/143‑1 (MQ=255)
tCCTTTTTCTATCGCCGGCGCAACCTCGCTGGCAGATAATCTTTCTTGATAAGAGTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACATGCCGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCg                                                                                                                                     <  1:147341/143‑1 (MQ=255)
   tttttCTATCGCCTGCGCAACCTCGCTGGCAGATAAACTTTCTTTATAAGAGTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGAt                                                                                                                                  >  1:150953/1‑143 (MQ=255)
    ttttCTATCGCCTGCGCAACCTCGCTGGCAGATAAACTTTCTTTATAAGAGTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCGCGGAtt                                                                                                                                 >  1:535741/1‑143 (MQ=255)
     tttCTATCGCCTGCGCAACCTCGCTGGCAGATAACCATTCACTCTGTCATTGTGTGGCGCTTACGTTTTTCATACCTATGCCTGTTACATCGTGTCTTCTTTTGGCGTTTCTTTTATTCGGCCTGGATTCCAGGCCCGGATTg                                                                                                                                <  1:840741/143‑1 (MQ=255)
     tttCTACTGCCTGCGCAACCTCGCTGGCAGATAAACTTTCTTTATAAGAGTCTGGGGCGATTACGATTTTCATACCTATGCTTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGAATCCAGGccgaataag                                                                                                                                <  1:261754/143‑8 (MQ=255)
     gcggtATCGCCTGCGCAACCTCGCTGGCAGATAAACTTTCTTTATAAGAGTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTg                                                                                                                                <  1:206607/139‑1 (MQ=255)
         tATCGCCTGCGCAACCTCGCTGGCAGATAAACTTTCTTTATAAGAGTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACTACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGAGTCCAGGCCCGGATTGCaag                                                                                                                            >  1:765226/1‑142 (MQ=255)
           tCGCCTGCGCAACCTCGCTGGCAGATAAACTTTCTTTATAAGAGTCTGGGGCGGTTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTCTTCGGCCTGGATTCCAGGCCCGGATTGCAATAc                                                                                                                          >  2:804149/1‑143 (MQ=255)
              ccTGCGCAACCTCGCTGGCAGATAAACTTTCTCTATAAGAGTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATCCCAGGCCCGGATTGCAATACGcc                                                                                                                       >  2:494686/1‑143 (MQ=255)
              ccTGCGCAACCTCGCTGGCAAATCGACTGTCTTTATAAGAGTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAACACAcc                                                                                                                       >  2:573494/1‑143 (MQ=255)
                tGCGCAACCTCGCTGGCAGATAAACTTTCTTTATAAGAGTCTGGGGCGATTACGATTTTCATACCTATGCCTGGTACCACATAACGCCGGAGGGAGTTTCTCTTATTAGGCCTGGATTCCAGGCCCGGCTTGCAATCCGcctt                                                                                                                     >  2:55625/1‑141 (MQ=255)
                    cAACCTCGCTGGCAGATAAACTTTCTTTATAAGAGTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCgg                                                                                                                 <  1:55706/143‑1 (MQ=255)
                    cAACCTCGCTGCTAGATAAACTTTCTTTATAAGAGTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACAGGACGCCGGAGGGCGTTTCTCTTATGCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCgg                                                                                                                 <  1:376404/143‑1 (MQ=255)
                       ccTCGCTGGCAGATAAACTTTCTTTATAAGAGTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCgcgcc                                                                                                              >  2:486033/1‑142 (MQ=255)
                         tCGCTGGCAGATAAACTTTCTTTATAAGAGTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATGGCAATACGCCATGCGGGCagg                                                                                                            >  1:808052/1‑141 (MQ=255)
                         tCGCTGGCAGATAAACTTTCTTTATAAGAGTCTGGGGCGATTACGATTTTAATACCTATGCGTGTTACCGCGGGACGCCGGAGGGCGTTACACTTATGCGGCCTGGATTACAGGCCCGGATTGGAATACGCGATCAGGGCacg                                                                                                            >  2:903385/1‑143 (MQ=255)
                             tGGCAGATAAACTTTCTTTATAAGAGTCTGGGGCGATTACGATTGTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTGATTCGGCCTGGATTCCAGGCCCGGATTGCACTACGCCATCCGGGCacgacgg                                                                                                        >  1:363235/1‑142 (MQ=255)
                                cAGATAAACTTTCTTTATAAGAGTCTGGGGCGATTACGCTAGACATACCTGTGCCTGTTACCACATGACACCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCAt                                                                                                     <  1:682544/143‑1 (MQ=255)
                                 aGATAAACTTTCTTTATAAGAGTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCACCCAGGCCCGAAGTCagt                                                                                                    >  2:604013/1‑141 (MQ=255)
                                    tAAACTTTCTTTATAAGAGTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAAc                                                                                                 <  2:365518/143‑1 (MQ=255)
                                        ctttctttATAAGAGTCTGGGGCGCTTACGATTTTCATACCTATTCCTGTTACCGCATGACGCCGGAGGGCGATCATCTTAATCGGCCTGGATACGAGGCCCGGATCACAATACGCCACCCGGGCACGACGTCATTAACgtac                                                                                             >  2:707178/1‑140 (MQ=255)
                                                 taAGAGTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGAc                                                                                    >  1:336252/1‑143 (MQ=255)
                                                 taAGAGTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATGGCAATACGCCACCCGGGCACGACGTCATTAACGAGTGACTTCGAc                                                                                    >  2:452647/1‑143 (MQ=255)
                                                               gATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCACTAACGAGTAACTTCGACTTTAGCCAGtcttt                                                                      >  2:809020/1‑143 (MQ=255)
                                                                aTTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTc                                                                     <  1:522624/143‑1 (MQ=255)
                                                                    cgtttttCACATCTCTGTCTGTTACCACATGACGCCTGAGGGCGTTTCTTTTATTCGGTCTGGATTCCAGGCCCGGATTGAAATCCGCTATTCGGGCTCGACGTCATTAACGAGTTACTTAGATTTTCGCCAGTTTTTCgtag                                                                 <  2:117051/140‑1 (MQ=255)
                                                                     gATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGCAACTCCGACTTTCGCCAGTTCTTCgcagt                                                                >  2:621748/1‑143 (MQ=255)
                                                                            aTACCTCTGCCTGTTACCTCATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGc                                                         <  2:313422/143‑1 (MQ=255)
                                                                                         ttACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAGGGCGCTATgc                                            >  2:267005/1‑142 (MQ=255)
                                                                                         tcaCCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAGGGCGCTATGa                                            <  1:459304/141‑1 (MQ=255)
                                                                                              acaTGACGCCGGAGGGCGTTTCTCTTATGCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAGGGCGCTATGATccgc                                       <  1:881210/143‑1 (MQ=255)
                                                                                                  ggcGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCCTCCTGGCTCGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAGGGCGCTATGATCCGCCGtt                                   <  2:818569/141‑1 (MQ=255)
                                                                                                      ccGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAGGGCGCTATGATCCGCCGTTCCTa                               <  1:602246/143‑1 (MQ=255)
                                                                                                            ggCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAGGGCGCTATGATCCGCCGTTCCTAAACCAt                         <  1:252707/143‑1 (MQ=255)
                                                                                                              cccttCTCTTATTCGGCCTGGATTCCAGGCCCCGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAGGGCGCTATGATCCGCCGTTCCTAAACCATCt                       <  1:17744/140‑1 (MQ=255)
                                                                                                              cGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAGGGCGCTATGATCCGCCGTTCCTAAACCATCt                       >  1:511654/1‑143 (MQ=255)
                                                                                                                tttCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTCCGTAGTAGCACGCCAGGGCGCTCTGCTCCGCCGTTCCTAACCCTTTTGc                     >  1:728173/1‑143 (MQ=255)
                                                                                                                    tctTATTCGGCCTGGATTCCAGGCCCGGATCGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAGGGCGCTATGATCCGCCGTTCCTAAACCATCTGCTCGc                 <  2:150953/143‑1 (MQ=255)
                                                                                                                     ctTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAGGGCGCTATGATCCGCCGTTCCTAAACCATCTGCCCGCa                >  1:520722/1‑143 (MQ=255)
                                                                                                                      aatccgggCCTGGAATCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTGGCACGCCAGGGCGCTATGATCGGCCGTTCCTCAACCATCTGCTCGCAg               >  1:909224/7‑143 (MQ=255)
                                                                                                                           gggCCTGGAATCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAGGGCGCTATGATCCGCCGGTCCTAAACCATCTGCTCGCAGTGCCt          >  1:50546/2‑143 (MQ=255)
                                                                                                                            ggCCTGGATGCCAGCCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAGGGCGCTATGATCCGCCGTTCCTAAACCATCTGCTCGCAGTGCCTg         <  2:30195/143‑1 (MQ=255)
                                                                                                                               cTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTCTTCATAGTACCACGCCAGGGCGCGATGATCCGCCGTTCCCAAACCATCGGCGCGCAGTGCCCGcct      >  2:810146/1‑141 (MQ=255)
                                                                                                                                  gATCCCAGGCCCGGATTGCAATGCGCCATCCGGGCACGACGTCATTACCGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGCAGCACGCCAGGGCGCTATGATCCGCCGTTCCTAAACCATCTGCTCGCAGTGCCTGcatcat   <  2:583949/143‑1 (MQ=255)
                                                                                                                                   aTTTCATGCCCGGATTGCAATACGCCGTCCGGGCACGACGTGATTAACGAGGAACTGCGACGTTCGCCAGTTTTTGGTAGTAGCACGCCAGGGCGCTATGATCCGCCGTTCCTAAACCATCTGCTCGCAGTGCCTGCatcatc  <  2:30798/143‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                             
TCCTTTTTCTATCGCCTGCGCAACCTCGCTGGCAGATAAACTTTCTTTATAAGAGTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAGGGCGCTATGATCCGCCGTTCCTAAACCATCTGCTCGCAGTGCCTGCATCATC  >  NC_000913/3271693‑3271966

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