Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 650,848 A→T 14.3% intergenic (‑137/+43) aroQ‑I ← / ← dsbD‑I type II 3‑dehydroquinate dehydratase/thiol/disulfide interchange protein 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947650,8480AT14.3% 85.0 / 6.9 35intergenic (‑137/+43)aroQ‑I/dsbD‑Itype II 3‑dehydroquinate dehydratase/thiol/disulfide interchange protein 1
Reads supporting (aligned to +/- strand):  ref base A (19/11);  new base T (3/2);  total (22/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

GATTCCTTGGATTTGTAGGGCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATCAGCGGTAATTTCGCTGCCGTTCGCAGC  >  NC_002947/650717‑650995
                                                                                                                                   |                                                                                                                                                   
gATTCCTTGGATTTGTAGGGCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCAttt                                                                                                                                     >  7:173528/1‑148 (MQ=255)
      ttGGATTTGTAGGGCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATTTAAATATCCGGTCAca                                                                                                                                          >  8:72965/1‑136 (MQ=255)
      ttGGATTTGTAGGGCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATTTAAATATCCGGTCAca                                                                                                                                          <  7:72965/137‑2 (MQ=255)
      ttGGATTTGTAGGGCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGAtt                                                                                                                              >  3:5607/1‑149 (MQ=255)
          aTTTGTAGGGCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATTTAAATATCCGGTc                                                                                                                                             >  8:8671/1‑130 (MQ=255)
          aTTTGTAGGGCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATTTAAATATCCGGTc                                                                                                                                             <  7:8671/130‑1 (MQ=255)
              gTAGGGCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGa                                                                                                                       <  8:173528/148‑1 (MQ=255)
                  ggCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTcg                                                                                                                  >  1:126471/1‑149 (MQ=255)
                  ggCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACACATTTCAGGATTTTTCCGAGCTcg                                                                                                                  >  3:130466/1‑149 (MQ=255)
                             acaGCGCGGGGAGTGTGCCGTAAAGCAACGGAAAGTGTCCAGTTCCCGGCAATAGGCAGCACGATGTCCGAGGTTCGCGCTAATTATGTGACCGGATATTTAAATCGGGTCACTCTTTTCAGGATTTTTCCGAGCTCGCGGGAAGtttt                                                                                                       <  6:10669/149‑3 (MQ=255)
                                            tgCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCcg                                                                                        <  2:33036/149‑1 (MQ=255)
                                                        aCGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGt                                                                             >  2:32872/1‑148 (MQ=255)
                                                                  tCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTc                                                                                                                   >  5:30691/1‑100 (MQ=255)
                                                                  tCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTc                                                                                                                   <  6:30691/100‑1 (MQ=255)
                                                                         cccGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAAtt                                                            >  7:9909/1‑148 (MQ=255)
                                                                            ggCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCAcc                                                        <  4:5607/149‑1 (MQ=255)
                                                                             gCAATAGCCAGGACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCTCACGCGCGTCAGAATCTGCGCGAATTCACCg                                                       <  2:126471/149‑1 (MQ=255)
                                                                               aaTAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACACATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGc                                                     <  4:130466/149‑1 (MQ=255)
                                                                                    ccAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTa                                                >  1:517/1‑149 (MQ=255)
                                                                                            aTGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGcc                                         >  8:146161/1‑148 (MQ=255)
                                                                                                        cgcgcgAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCCAGGGAAGTTATAGACAGACTTTGCCGCGCACGCGCGTCAGAATCTCCGCGAAGTCACCGGCATTTAACTCGCCAATCACCCCAt                              >  4:48940/1‑147 (MQ=255)
                                                                                                         gcgcgAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCAc                                                         >  5:53714/1‑119 (MQ=255)
                                                                                                         gcgcgAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCAc                                                         <  6:53714/119‑1 (MQ=255)
                                                                                                                  aTGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATCAGCGGTAAttt                  >  3:80218/1‑149 (MQ=255)
                                                                                                                    gtgACCGGATTTAAATATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATCAGCGGTAATTTCg                >  2:137869/1‑149 (MQ=255)
                                                                                                                     tgACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATCAGCGGTAATTTCGc               >  5:214382/1‑149 (MQ=255)
                                                                                                                      gACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATCAGCGGTAATTTCGCt              >  7:169368/1‑149 (MQ=255)
                                                                                                                      gACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATCAGCGGTAATTTCGCt              >  5:131197/1‑149 (MQ=255)
                                                                                                                        ccGGATATTTAAATCCGGTTACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATCAGCGGTAATTTCGCTGc            >  5:23307/1‑149 (MQ=255)
                                                                                                                         cGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGTTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCgg                                                      >  4:10945/1‑106 (MQ=255)
                                                                                                                         cGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGTTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCgg                                                      <  3:10945/106‑1 (MQ=255)
                                                                                                                          ggATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATCAGCGGTAATTTCGCTGCCg          >  6:209370/1‑149 (MQ=255)
                                                                                                                           gATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATCAGCGGTaa                     <  4:35605/137‑1 (MQ=255)
                                                                                                                           gATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATCAGCGGTaa                     >  3:35605/1‑137 (MQ=255)
                                                                                                                               tttAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATCAGCGGTAATTTCGCTGCCGTTCg      <  6:131197/148‑1 (MQ=255)
                                                                                                                                  aaaTCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATCAGCGGTAATTTCGCTGCCGTTCgcag   >  3:246340/1‑148 (MQ=255)
                                                                                                                                   aaTCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATCAGCGGTAATTTCGCTGCCGTTCgcagc  <  6:214382/148‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                                   
GATTCCTTGGATTTGTAGGGCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATCAGCGGTAATTTCGCTGCCGTTCGCAGC  >  NC_002947/650717‑650995

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