Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 3,823,901 +T coding (1248/1293 nt) kguT ← 2‑ketogluconate transporter, putative

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,823,9011.T90.5% 52.2 / 1.5 21coding (1248/1293 nt)kguT2‑ketogluconate transporter, putative
Reads supporting (aligned to +/- strand):  ref base . (2/0);  new base T (9/10);  total (11/10)
Fisher's exact test for biased strand distribution p-value = 4.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.97e-01

CAGGTCACCTCGACGCGATCTTGCAGGCGTGCCATCAAGTCGTCGGACAGGCGTTTATACAGGACGATTCGCTTCTTCATTGCAGGTTTCCTAATGGTTCATGGCCAGACGCGAGCTATCGCGCTTTGCCCGGCCAGAGGTCTGGGAAGTATTGAGCGCAGCGGTAAGCGCCACTGCAGTCAGCAGCGCACCACACATGAACAGGTAAGAGGCGCCGGGCCCGCCCGTAGCCCCATTGAGGTAGCCCACCAGCCACGAGCCGCCAAA  >  NC_002947/3823766‑3824032
                                                                                                                                        |                                                                                                                                   
cAGGTCACCTCGACGCGATCTTGCAGGCGTGCCATCAAGTCGTCGGACAGGCGTTTATACAGGACGATTCGCTTCTTCATTGCAGGTTTCCTAATGGTTCATGGCCAGACGCGAGCTATCGCGCTTTGCCCGGCCa                                                                                                                                      <  2:5672/136‑1 (MQ=255)
cAGGTCACCTCGACGCGATCTTGCAGGCGTGCCATCAAGTCGTCGGACAGGCGTTTATACAGGACGATTCGCTTCTTCATTGCAGGTTTCCTAATGGTTCATGGCCAGACGCGAGCTATCGCGCTTTGCCCGGCCa                                                                                                                                      <  1:362780/136‑1 (MQ=255)
     cACCTCGACGCGATCTTGCAGGCGTGCCATCAAGTCGTCGGACAGGCGTTTATACAGGACGATTCGCTTCTTCATTGCAGGTTTCCTAATGGTTCATGGCCAGACGCGAGCTATCGCGCTTTGCCCGGCCatgagg                                                                                                                                 >  2:170481/1‑131 (MQ=255)
        cTCGACGCGATCTTGCAGGCGTGCCATCAAGTCGTCGGACAGGCGTTTATACAGGACGATTCGCTTCTTCATTGCAGGTTTCCTAATGGTTCATGGCCAGACGCGAGCTATCGCGCTTTGCCCGGCCATGAGGTCt                                                                                                                              <  1:141140/136‑1 (MQ=255)
        cTCGACGCGATCTTGCAGGCGTGCCATCAAGTCGTCGGACAGGCGTTTATACAGGACGATTCGCTTCTTCATTGCAGGTTTCCTAATGGTTCATGGCCAGACGCGAGCTATCGCGCTTTGCCCGGCCATGAGGTCt                                                                                                                              <  2:414768/136‑1 (MQ=255)
        cTCGACGCGATCTTGCAGGCGTGCCATCAAGTCGTCGGACAGGCGTTTATACAGGACGATTCGCTTCTTCATTGCAGGTTTCCTAATGGTTCATGGCCAGACGCGAGCTATCGCGCTTTGCCCGGCCATGAGGTCt                                                                                                                              <  1:318760/136‑1 (MQ=255)
           gACGCGATCTTGCAGGCGTGCCATCAAGTCGTCGGACAGGCGTTTATACAGGACGATTCGCTTCTTCATTGCAGGTTTCCTAATGGTTCATGGCCAGACGCGAGCTATCGCGCTTTGCCCGGCCATGAGGTCTggg                                                                                                                           >  2:376203/1‑136 (MQ=255)
           gACGCGATCTTGCAGGCGTGCCATCAAGTCGTCGGACAGGCGTTTATACAGGACGATTCGCTTCTTCATTGCAGGTTTACTAATGGTTCATGGCCAGACGCGAGCTATCGCGCTTTGCCCGGCCATGAGGTCTggg                                                                                                                           >  2:77844/1‑136 (MQ=255)
                    ttGCAGGCGTGCCATCAAGTCGTCGGACAGGCGTTTATACAGGACGATTCGCTTCTTCATTGCAGGTTTCCTAATGGTTCATGGCCAGACGCGAGCTATCGCGCTTTGCCAGGCCATGAGGTCTGGGAAGTATTGa                                                                                                                  <  1:437024/136‑1 (MQ=255)
                                   cAAGTCGTCGGACAGGCGTTTATCCCGGACGATGCGCTTCTTCATTGCAGGTTTCCTAATGGTTCATGGCCAGACGCGAGCTATCGCGCTTTGCCCGGCCATGAGGTCTGGGAAGTATTGAGCGCAGCGGTAAgcg                                                                                                   <  1:440383/136‑1 (MQ=255)
                                     aGTCGTCGGACAGGCGTTTATACAGGACGATTCGCTTCTTCATTGCAGGTTTCCTAATGGTTCATGGCCAGACGCGAGCTATCGCGCTTTGCCCGGCCATGAGGTCTGGGAAGTATTGAGCGCAGCGGTAAGCGcc                                                                                                 <  1:136372/136‑1 (MQ=255)
                                     aGTCGTCGGACAGGCGTTTATACAGGACGATTCGCTTCTTCATTGCAGGTTTCCTAATGGTTCATGGCCAGACGCGAGCTATCGCGCTTTGCCCGGCCATGAGGTCTGGGAAGTATTGAGCGCAGCGGTAAGCGcc                                                                                                 <  1:189592/136‑1 (MQ=255)
                                                                 gATTCGCTTCTTCATTGCAGGTTTCCTAATGGTTCATGGCCAGACGCGAGCTATCGCGCTTTGCCCGGCCATGAGGTCTGGGAAGTATTGAGCGCAGCGGTAAGCGCCACTGCAGTCAGCAGCGCACCACACATGa                                                                     >  1:396938/1‑136 (MQ=255)
                                                                        ttcttcATTGCAGGTTTCCTAATGGTTCATGGCCAGACGCGAGCTATCGCGCTTTGCCCGGCCATGAGGTCTGGGAAGTATTGAGCGCAGCGGTAAGCGCCACTGCAGTCAGCAGCGCACCACACATGAACAGGTa                                                              <  1:170481/136‑1 (MQ=255)
                                                                              aTTGCAGGTTTCCTAATGGTTCATGGCCAGACGCGAGCTATCGCGCTTTGCCCGGCCATGAGGTCTGGGAAGTATTGAGCGCAGCGGTAAGCGCCACTGCAGTCAGCAGCGCACCACACATGAACAGGTAAGAGgc                                                        >  1:376633/1‑136 (MQ=255)
                                                                                       ttCCTAATGGTTCATGGCCAGACGCGAGCTATCGCGCTTTGCCCGGCCATGAGGTCTGGGAAGTATTGAGCGCAGCGGTAAGCGCCACTGCAGTCAGCAgcg                                                                                 <  1:367163/102‑1 (MQ=255)
                                                                                       ttCCTAATGGTTCATGGCCAGACGCGAGCTATCGCGCTTTGCCCGGCCATGAGGTCTGGGAAGTATTGAGCGCAGCGGTAAGCGCCACTGCAGTCAGCAgcg                                                                                 >  2:367163/1‑102 (MQ=255)
                                                                                             aTGGTTCATGGCCAGACGCGAGCTATCGCGCTTTGCCCGGCCATGAGGTCTGGGAAGTATTGAGCGCAGCGGTAAGCGCCACTGCAGTCAGCAGCGCACCACACATGAACAGGTAAGAGGCGCCGGGCCCGCCCGt                                         <  1:97016/136‑1 (MQ=255)
                                                                                                      ggCCAGACGCGAGCTATCGCGCTTTGCCCGGCCATGAGGTCTGGGAAGTATTGAGCGCAGCGGTAAGCGCCACTGCAGTCAGCAGCGCACCACACATGAACAGGTAAGAGGCGCCGGGCCCGCCCGTAGCCCCAtt                                >  2:73487/1‑136 (MQ=255)
                                                                                                      ggCCAGACGCGAGCTATCGCGCTTTGCCCGGCCATGAGGTCTGGGAAGTATTGAGCGCAGCGGTAAGCGCCACTGCAGTCAGCAGCGCACCACACATGAACAGGTAAGAGGCGCCGGGCCCGCCCGTAGCCCCAtt                                >  1:129047/1‑136 (MQ=255)
                                                                                                              gcgAGCTATCGCGCTTTGCCAGGCCATGAGGTCTGGGAAGTATTGAGCGCAGCGGTAAGCGCCACTGCAGTCAGCAGCGCACCACACATGAACAGGTAAGAGGCGCCGGGCCCGCCCGTAGCCCCATTGAGGTAGc                        >  1:357187/1‑136 (MQ=255)
                                                                                                                         cgcTTTGCCCGGCCATGAGGTCTGGGAAGTATTGAGCGCAGCGGTAAGCGCCACTGCAGTCAGCAGCGCACCACACATGAACAGGTAAGAGGCGCCGGGCCCGCCCGTAGCCCCATTGAGGTAGCCCACCAGCCAc             >  1:172542/1‑136 (MQ=255)
                                                                                                                                   ggCCAGAGGTCTGGGAAGTATTGAGCGCAGCGGTAAGCGCCACTGCAGTCAGCAGCGCACCACACATGAACAGGTAAGAGGCGCCGGGCCCGCCCGTAGCCCCATTGAGGTAGCCCACCAGCCACGAGCCGCCaaa  >  1:42775/1‑136 (MQ=255)
                                                                                                                                   ggCCAGAGGTCTGGGAAGTATTGAGCGCAGCGGTAAGCGCCACTGCAGTCAGCAGCGCACCACACATGAACAGGTAAGAGGCGCCGGGCCCGCCCGTAGCCCCATTGAGGTAGCCCACCAGCCACGAGCCGCCaaa  >  1:367462/1‑136 (MQ=255)
                                                                                                                                        |                                                                                                                                   
CAGGTCACCTCGACGCGATCTTGCAGGCGTGCCATCAAGTCGTCGGACAGGCGTTTATACAGGACGATTCGCTTCTTCATTGCAGGTTTCCTAATGGTTCATGGCCAGACGCGAGCTATCGCGCTTTGCCCGGCCAGAGGTCTGGGAAGTATTGAGCGCAGCGGTAAGCGCCACTGCAGTCAGCAGCGCACCACACATGAACAGGTAAGAGGCGCCGGGCCCGCCCGTAGCCCCATTGAGGTAGCCCACCAGCCACGAGCCGCCAAA  >  NC_002947/3823766‑3824032

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