Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,344,480 A→G 14.4% intergenic (+55/‑123) pgi‑II → / → acsA‑II glucose 6‑phosphate isomerase/acetyl‑CoA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,344,4800AG14.4% 60.9 / 4.0 34intergenic (+55/‑123)pgi‑II/acsA‑IIglucose 6‑phosphate isomerase/acetyl‑CoA synthetase
Reads supporting (aligned to +/- strand):  ref base A (17/9);  new base G (2/2);  total (21/13)
Fisher's exact test for biased strand distribution p-value = 6.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.77e-01

CCAGCGCTTGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTC  >  NC_002947/5344335‑5344608
                                                                                                                                                   |                                                                                                                                 
ccAGCGCTTGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCTGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTTCCAGGCAACGCGTTACCt                                                                                                                                >  1:173284/1‑149 (MQ=255)
 cAGCGCTTGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCTCGGTTGATCCGAGCTTGTGCCTGCCTCTTCGAGGGTAATCCCC‑‑CCGACAAGGCAACGCGTTACATc                                                                                                                               >  4:37910/1‑149 (MQ=25)
 cAGCGCTTGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCCGGTAAACCCC‑‑CCGACCAGGCAACCCCTTACCTc                                                                                                                               >  3:144700/1‑149 (MQ=37)
                  ccTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGTCCTCTTCGCGGGTAAACCCGCTCCTACGAGGTAACGCGTTGCCTGGTAG‑‑GCGGGTTTAc                                                                                                               <  2:11064/148‑1 (MQ=21)
                    tGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCACTTCGCGGGTTAACCCC‑‑CCAACCAGGGCACCGGTTTACTCGGAGGAGCGGGTTTTCCcg                                                                                                            >  3:61933/1‑149 (MQ=25)
                      gAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCcgc                                                                                                           <  7:150346/148‑1 (MQ=255)
                        aGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCCGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCATCGCCGGTAAAACCG‑‑CCAGCCACGCAACTCGTTACCTCGTAGGAGCTGGTTAACCcccg                                                                                                          >  1:62215/1‑147 (MQ=37)
                          aCAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAaga                                                                                                      <  8:71229/149‑1 (MQ=255)
                            aGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGTAACGCGTTGCCTGGTAG‑‑GCGGGTTTAc                                                                                                               >  1:197775/1‑138 (MQ=21)
                            aGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGTAACGCGTTGCCTGGTAG‑‑GCGGGTTTAc                                                                                                               <  2:197775/138‑1 (MQ=21)
                            aGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAgg                                                                                                    <  3:269532/149‑1 (MQ=255)
                               gcCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAggccg                                                                                                <  6:282612/149‑1 (MQ=255)
                                      aCGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCCgg                                                                                         >  5:486/1‑149 (MQ=255)
                                                           acaacaaGTTCCGCGGCCGGCAGCCTGGTTTACCCGAGTTTGTGCCGCGCTCTTCGTTGTTAAACCCG‑‑CCTACCATGCATCCGGTATCCTCTTAGGTCCGTGTTTACCCGCTTTTTGGCCCGGCCCAGGTTGAACACCGCTCCCAGCTa                                                                     <  2:62215/148‑1 (MQ=255)
                                                              actactTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCcgcg                                                                                                          >  8:50527/1‑109 (MQ=255)
                                                              actactTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCcgcg                                                                                                          <  7:50527/109‑1 (MQ=255)
                                                                  ctTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAAAGCGTTTCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCCGGTTGCACACCGCTCCCCGCTACCCAGTCt                                                             >  7:178157/1‑149 (MQ=38)
                                                                  ctTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCGACCACGCAACGCGGTACCTCGTAGGAGCCGGTTTAACCGCGGAGACGCCGGCCCAGGTTGAACACCCCGCCCAGCTACACTGTcg                                                             >  8:248564/1‑148 (MQ=38)
                                                                    tCCGCGGCCGTCACCGGGTTTTATCCTTGCTTGTGCCGGCCTTTTCGCGGGTAAACCCG‑‑CCTCCCCGCCAACGCGTTACCTCGTTGGTGCGGTTTTACCCGCTAAGTGGCCGGCCCAGGTTGAac                                                                                   <  6:78835/125‑1 (MQ=25)
                                                                    tCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCGCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAac                                                                                   >  5:78835/1‑125 (MQ=255)
                                                                                  cggggTTGATCGGAGCTTGGGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTTGGAGGGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTg                                             <  7:248564/146‑1 (MQ=38)
                                                                                           tCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGGCTTATCGAATAGCCGTTGCCGTCCCTc                                    >  3:222381/1‑149 (MQ=255)
                                                                                           tCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACAACCCTCCCAGCTACACTGTCTTATCGAATAGCCGTTTCAGTCCa                                      >  6:86749/1‑147 (MQ=255)
                                                                                             cGAGCTTGTGCCGGCCTCTTAGCGGGTAAACCCG‑‑CCGACCCGGCAACGCGGTACCGCGTAGGCGCGGGGTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCGGCCAGCTACACTGGCTTATCGAATAGCCGTTGCACTGCATCGc                                  >  2:22238/1‑149 (MQ=37)
                                                                                                 ctggggccgccCTCTTCGCGGCAAACCCCG‑‑CCAACCCGCCACCCCGTTCCCTCTCCGGCGGGGTTTTACCCGCTACGTGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCt                              <  4:222381/144‑1 (MQ=255)
                                                                                                               ttCGCGGGTAAACCCGCTCCTACGAGGTAACGCGTTGCCTGGTAG‑‑GCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGac                >  2:3324/1‑149 (MQ=21)
                                                                                                                           cccG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGt    >  6:25200/1‑149 (MQ=255)
                                                                                                                            ccG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTg                                             <  3:244149/107‑1 (MQ=255)
                                                                                                                            ccG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGc                     >  4:209200/1‑131 (MQ=255)
                                                                                                                            ccG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGc                     <  3:209200/131‑1 (MQ=255)
                                                                                                                            ccG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGtt   >  1:155310/1‑149 (MQ=255)
                                                                                                                             cG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGTCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTc  >  6:51883/1‑149 (MQ=255)
                                                                                                                              cccccTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTg                                             >  4:244149/4‑107 (MQ=255)
                                                                                                                                        ggCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAg                   >  8:19989/1‑123 (MQ=255)
                                                                                                                                                   |                                                                                                                                 
CCAGCGCTTGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTC  >  NC_002947/5344335‑5344608

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: