Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,178,139 G→A 20.0% intergenic (‑53/+137) PP_5410 ← / ← glmU DeoR family transcriptional regulator/bifunctional N‑acetyl glucosamine‑1‑phosphate uridyltransferase/glucosamine‑1‑phosphate acetyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,178,1390GA20.0% 37.0 / 4.6 20intergenic (‑53/+137)PP_5410/glmUDeoR family transcriptional regulator/bifunctional N‑acetyl glucosamine‑1‑phosphate uridyltransferase/glucosamine‑1‑phosphate acetyl transferase
Reads supporting (aligned to +/- strand):  ref base G (9/7);  new base A (2/2);  total (11/9)
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.00e+00

CTTACCTCACCCTGCTCGCTGAGCAGAGCCAGAATATTGTGTCGCCGTTGTGGTGTGTTTCGTTTCGACATGTGCGTTTAAGTTTCGATTCGAAAGTTAATGGTTGCAATCAAAACCTAAGACGCGCGATTCGTCAAGCTGTCGCTGACCTGTGGATTGTTTGCTTCTGCTTTGAGGCGGTCCTATAGGCGCAGCCGCGCCAATAAAAAAGCC  >  NC_002947/6178016‑6178228
                                                                                                                           |                                                                                         
cTTACCTCACCCTGCTCGCTGAGCAGAGCCAGAATATTGTGTCGCCGTTGTGGTGTGTTTCGTTTCGACATGTGCGTTTAAGTTTCGATTCGAAAGTTAATGGTTGCAATCAAAACCTAAGACGCGCGATTCGTCa                                                                               <  1:452554/136‑1 (MQ=255)
      tCACCCTGCTCGCTGAGCAGAGCCAGAATATTGTGTCGCCGTTGTGGTGTGTTTCGTTTCGACATGTGCGTTTAAGTTTCGATTCGAAAGTTAATGGTTGCAATCAAAACCTAAGACACGCGATTCGTCAAGCTGt                                                                         >  2:177061/1‑136 (MQ=255)
      tCACCCTGCTCGCTGAGCAGAGCCAGAATATTGTGTCGCCGTTGTGGTGTGTTTCGTTTCGACATGTGCGTTTAAGTTTCGATTCGAAAGTTAATGGTTGCAATCAAAACCTAAGACACGCGATTCGTCAAGCTGt                                                                         >  2:174763/1‑136 (MQ=255)
           cTGCTCGCTGAGCAGAGCCAGAATATTGTGTCGCCGTTGTGGTGTGTTTCGTTTCGACATGTGCGTTTAAGTTTCGATTCGAAAGTTAATGGTTGCAATCAAATCCTAAGACGCGCGATTCGTCAAGCTGTCGCTg                                                                    >  1:209150/1‑136 (MQ=255)
            tGCTCGCTGAGCAGAGCCAGAATATTGTGTCGCCGTTGTGGTGTGTTTCGTTTCGACATGTGCGTTTAAGTTTCGATTCGAAAGTTAATGGTTGCAATCAAAACCTAAGACGCGCGATTCGTCAAGCTGTCGCTGa                                                                   >  1:291538/1‑136 (MQ=255)
            tGCTCGCTGAGCAGAGCCAGAATATTGTGTCGCCGTTGTGGTGTGTTTCGTTTCGACATGTGCGTTTAAGTTTCGATTCGAAAGTTAATGGTTGCAATCAAAACCTAAGACGCGCGATTCGTCAAGCTGTCGCTGa                                                                   >  1:492192/1‑136 (MQ=255)
                cGCTGAGCAGAGCCAGAATATTGTGTCGCCGTTGTGGTGTGTTTCGTTTCGACATGTGCGTTTAAGTTTCGATTCGAAAGTTAATGGTTGCAATCAAAACCTAAGACGCGCGATTCGTCAAGCTGTCGCTGACCtg                                                               >  1:223167/1‑136 (MQ=255)
                cGCTGAGAAGAGCCAGAATATTGTGTCGCCGTTGTGGTGTGTTTCGTTTCGACATGTGCGTTTAAGTTTCGATTCGAAAGTTAATGGTTGAAATCAAAACCTAAGACGCTCGATTCGTCAAGATGTCGATGACCtg                                                               >  1:297786/1‑136 (MQ=255)
                 gCTGAGCAGAGCCAGAATATTGTGTCGCCGTTGTTTTGTGTTTCGTTTCGACATTTGCGTTTAAGTTTCGATTCGAAAGTTAATGGTTGCAATCAAAACCTAAGACGCGCTATTCGTCAAGCTGTCGCTGACCtgt                                                              >  1:58991/1‑136 (MQ=255)
                           gCCAGAATATTGTGTCGCCGTTGTGGTGTGTTTCGTTTCGACATGTGCGTTTAAGTTTCGATTCGAAAGTTAATGGTTGCAATCAAAACCTAAGACGCGCGATTCGTCAAGCTGTCGCTGACCTGTGGATTGTTTg                                                    <  1:106245/136‑1 (MQ=255)
                                       tgtCGCCGTTGTGGTGTGTTTCGTTTCGACATGTGCGTTTAAGTTTCGATTCGAAAGTTAATGGTTGCAATCAAAACCTAAGACGCGCGATTCGTCAAGCTGTCGCTGACCTGTGGATTGTTTGCTTCTGCTTTGa                                        <  1:113150/136‑1 (MQ=255)
                                       tgtCGCCGTTGTGGTGTGTTTCGTTTCGACATGTGCGTTTAAGTTTCGATTCGAAAGTTAATGGTTGCAATCAAAACCTAAGACGCGCGATTCGTCAAGCTGTCGCTGACCTGTGGATTGTTTGCTTCTGCTTTGa                                        <  1:303213/136‑1 (MQ=255)
                                       tgtCGCCGTTGTGGTGTGTTTCGTTTCGACATGTGCGTTTAAGTTTCGATTCGAAAGTTAATGGTTGCAATCAAAACCTAAGACGCGCGATTCGTCAAGCTGTCGCTGACCTGTGGATTGTTTGCTTCTGCTTTGa                                        <  1:344624/136‑1 (MQ=255)
                                            ccGTTGTGGTGTGTTTCGTTTCGACATGTGCGTTTAAGTTTCGATTCGAAAGTTAATGGTTGCAATCAAAACCTAAGACGCGCGATTCGTCAAGCTGTCGCTGACCTGTGGATTGTTTGCTTCTGCTTTGAGGCgg                                   >  1:440852/1‑136 (MQ=255)
                                             cGTTGTGGTGTGTTTCGTTTCGACATGTGCGTTTAAGTTTCGATTCGAAAGTTAATGGTTGCAATCAAAACCTAAGACGCGCGATTCGTCAAGCTGTCGCTGACCTGTGGATTGTTTGCTTCTGCTTTGAGGCGGt                                  >  1:355086/1‑136 (MQ=255)
                                             cGTTGTGGTGTGTTTCGTTTCGACATGTGCGTTTAAGTTTCGATTCGAAAGTTAATGGTTGCAATCAAAACCTAAGACGCGCGATTCGTCAAGCTGTCGCTGACCTGTGGATTGTTTGCTTCTGCTTTGAGGCGGt                                  >  1:357690/1‑136 (MQ=255)
                                                           tcgtttcgACATGTGCGTTTAAGTTTCGATTCGAAAGTTAATGGTTGCAATCAAAACCTAAGACACGCGATTCGTCAAGCTGTCGCTGACCTGTGGATTGTTTGCTTCTGCTTTGAGGCGGTCCTATAGGCGCAGc                    <  1:177061/136‑1 (MQ=255)
                                                           tcgtttcgACATGTGCGTTTAAGTTTCGATTCGAAAGTTAATGGTTGCAATCAAAACCTAAGACACGCGATTCGTCAAGCTGTCGCTGACCTGTGGATTGTTTGCTTCTGCTTTGAGGCGGTCCTATAGGCGCAGc                    <  1:174763/136‑1 (MQ=255)
                                                                             ttATGTTTCGATTCGAAAGTTAATGGTTGCAATCAAAACCTAAGACGCGCGATTCGTCAAGCTGTCGCTGACCTGGGGATTGTTTGCTTCTGGTTTGTGGCGGTCCTATAGGCGCAGCCGCGCCAATAAAAAAGcc  <  2:291538/136‑1 (MQ=255)
                                                                             ttAAGTTTCGATTCGAAAGTTAATGGTTGCAATCAAAACCTAAGACGCGCGATTCGTCAAGCTGTCGCTGACCTGTGGATTGTTTGCTTCTGCTTTGAGGCGGTCCTATAGGCGCAACCGCGCCAATAAAAAAGcc  <  2:492192/136‑1 (MQ=255)
                                                                                                                           |                                                                                         
CTTACCTCACCCTGCTCGCTGAGCAGAGCCAGAATATTGTGTCGCCGTTGTGGTGTGTTTCGTTTCGACATGTGCGTTTAAGTTTCGATTCGAAAGTTAATGGTTGCAATCAAAACCTAAGACGCGCGATTCGTCAAGCTGTCGCTGACCTGTGGATTGTTTGCTTCTGCTTTGAGGCGGTCCTATAGGCGCAGCCGCGCCAATAAAAAAGCC  >  NC_002947/6178016‑6178228

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: