Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,668,063 G→A 13.4% intergenic (+81/‑86) prpF → / → prpD aconitate isomerase/2‑methylcitrate dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,668,0630GA13.4% 91.2 / 5.8 39intergenic (+81/‑86)prpF/prpDaconitate isomerase/2‑methylcitrate dehydratase
Reads supporting (aligned to +/- strand):  ref base G (19/13);  new base A (3/2);  total (22/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.35e-01

TCATGAGCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCAC‑TGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTA‑ATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTA  >  NC_002947/2667918‑2668199
                                                                                                                                                   |                                                                                                                                        
tCATGAGCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCAC‑TGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTA‑ATCACAGGGttt                                                                                                                                       >  3:69783/1‑149 (MQ=255)
  aTGAGCCGTAGTGCTTGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCAC‑TGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTA‑ATCACAGGGTTTGa                                                                                                                                     <  5:383413/149‑1 (MQ=255)
      gCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCAC‑TGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTA‑ATCACAGGGttt                                                                                                                                       >  1:232604/1‑143 (MQ=255)
      gCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCAC‑TGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTA‑ATCACAGGGttt                                                                                                                                       <  2:232604/143‑1 (MQ=255)
        cGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCAC‑TGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTA‑ATCACAGGGTTTGAATGCtg                                                                                                                               <  6:381921/149‑1 (MQ=255)
          tAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCAC‑TGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTA‑ATCACAGGGTTTGAATGCtgtg                                                                                                                             <  4:69783/149‑1 (MQ=255)
                tCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCAC‑TGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTA‑ATCACAGGGTTTGAATGCTGTGATATcc                                                                                                                       >  5:341566/1‑149 (MQ=255)
                       cTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCAC‑TGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTA‑ATCACAGGGTTTGAATGCTGTGatat                                                                                                                         >  7:181381/1‑140 (MQ=255)
                       cTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCAC‑TGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTA‑ATCACAGGGTTTGAATGCTGTGatat                                                                                                                         <  8:181381/140‑1 (MQ=255)
                        tGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCAC‑TGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGAT‑TACCTGTgg                                                                                                                 >  7:355327/1‑147 (MQ=18)
                         gATGGCGGGCTGGGTTCGGGTGCCTGGCTATAGCTTCTAAACCACCAC‑TGGCGCTCGGCTTCAAGTGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTA‑ATCACAGGGTTTGAATGCtgt                                                                                                                              <  2:190065/133‑1 (MQ=38)
                         gATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCAC‑TGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTA‑ATCACAGGGTTTGAATGCtgt                                                                                                                              >  7:173407/1‑133 (MQ=255)
                         gATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCAC‑TGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTA‑ATCACAGGGTTTGAATGCtgt                                                                                                                              <  8:173407/133‑1 (MQ=255)
                         gATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCAC‑TGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTA‑ATCACAGGGTTTGAATGCtgt                                                                                                                              >  1:190065/1‑133 (MQ=255)
                             gAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCAC‑TGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGAT‑TACCTGTGGGAGCgg                                                                                                           <  8:355327/148‑1 (MQ=18)
                                          gcgTGCCTGGCGATAGCTTCTAAAGCACCACGTGGCGCGCGGCTTCAAGAGGGGCCGCTTTGCGCCCCTTCCCGGGCTTGCCCCCTCCCACAGGTA‑ATAACAGGGTTTGAATGGTGTTATATCCCGGTGGGCGCGGGCAAGCCcgcg                                                                                                >  1:39018/1‑147 (MQ=255)
                                              gCCTGGCGATAGCTTCTAAAGCACCAC‑TGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACCCGAA‑ATAACAGGGGTTGTATGGGGTGGTATGCCGTGGGGAGCGGGCAAGGCCGCGCAGGGCt                                                                                         >  8:291950/1‑149 (MQ=37)
                                                 tGGCGATAGCTTCTAAAGCACCAC‑TGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGAT‑TACCTGTgg                                                                                                                 >  3:64917/1‑122 (MQ=14)
                                                 tGGCGATAGCTTCTAAAGCACCAC‑TGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGAT‑TACCTGTgg                                                                                                                 <  4:64917/122‑1 (MQ=14)
                                                   gCGATAGCTTCTAAAGCACCAC‑TGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTA‑ATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGca                                                                                      <  6:183962/147‑1 (MQ=255)
                                                           ttCTAAAGCACCAC‑TGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTA‑ATGACAGGGTTTTTATGCTGTGATATGCCTGTGGGTGCGGGCAAGCCCGCGCAGGGCTGGACAGCCGcccc                                                                            >  3:117795/1‑149 (MQ=37)
                                                                aaGCACCAC‑TGGCGCTCGGCTTCAAGAGGGGCCGCTTTGAGGCCCTTCGCGGGCATGCCCCCTCCCACAGGTA‑ATCACACGGTTTGAATGGTGTGATAACCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGGACAGCAGCCCCGGCaa                                                                       >  6:164039/1‑149 (MQ=38)
                                                                           ggCGGTCGGCTTCAAGTGGGGCCGCTTTGCGCCCCTTGGCGGGGTTCCCCGCTCCCACAGTTA‑ATCCCCGGTTTTTAAGGCTGTTATATCCTTGTGGGAGGGGGAAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAggg                                                             <  6:341566/149‑1 (MQ=255)
                                                                            gcgcgcGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGCAA‑ATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAgg                                                              <  7:291950/147‑1 (MQ=255)
                                                                                tCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTTCCCGCTCCCACAGGTA‑ATAACAGGGGTTGAATGCTGTGATATACCTGTGGGAGCGGGCAACCCCGCGAAGGGCTGCACAGCAGCCCCCGCAATTGTGAAGGGAGcag                                                        >  8:8035/1‑147 (MQ=37)
                                                                                    cTTCAAGAGGGGACGCTTTGCGGCCCTTCGCGGGCTTTCCCGCTCCCACACGTA‑ATCAAACGGTTTTAAAGCTGTGATATACCTGAGGGAGCGGGCAACCCCCCCAAGGGCTGGCCAGCAGCCCCCGCAATTAAGAAGTGAGCAGGat                                                     >  6:88947/1‑148 (MQ=255)
                                                                                       cAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTA‑ATCACAGGGGTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGCAGGGCTGGACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGatat                                                   >  3:61682/1‑147 (MQ=255)
                                                                                           agGGGCCGCTTGGGGGCCGTGCGCGGGGTGGCCGGCTCCCACTGGTA‑ATCACTGGGTTTGAATGCTGTTATATCCGTGGGGGAGCGGGAAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAgcgc                                             <  7:111457/149‑1 (MQ=255)
                                                                                                     tggCGGCCCTTCGCGGGCTTCCCCGCTCCCCCCGTTA‑ATCACCGGTTTTTAATGCTGTCATATCCCGGGGGGTGGGGGAAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGa                                    <  4:61682/146‑1 (MQ=255)
                                                                                                         ggCCCTTCGCGGGCTTGCCCCCTACCACAGGTA‑ATAACAAGGTTTTAATTGTGTGATATACCTGTGGGAGCGGGCAAGCCCGCGAAGGGCAGCCCAGCAGCCCCGGCCAATCTGAAGGGAGCTGGATATTAGCGCCAACGTAGACCtg                                >  2:315859/1‑147 (MQ=25)
                                                                                                             cTTCGCGGGCTTGCCCCCTCCCAAAGGTA‑ATCACAGGGTTTGAATGCTGTGATATACCAGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTAAGAAAGGAGCTGGAAATGAGCGCCAACGTAGACCTAAACg                            >  2:111337/1‑148 (MQ=37)
                                                                                                                         ccccGCTCCCACCGTTA‑ATCCCCGGGTTTCAATGCTGTTATAACCCTGGGGGTGCGGGAAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGAtt                <  2:39018/147‑1 (MQ=255)
                                                                                                                          cccGCTCCCACAGGTA‑ATCACAGGGTTTGATTGCTGTCATATCCCTGTGGGAGCGGGAAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTAc              <  7:8035/149‑1 (MQ=255)
                                                                                                                           ccGCTCCCACAGGTA‑ATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTAc              >  4:5565/1‑148 (MQ=255)
                                                                                                                           ccGCTCCCACAGGTA‑ATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTAc              >  8:339574/1‑148 (MQ=255)
                                                                                                                               tCCCACAGGGATATCACAGCATTCAAACCCTGTGAT‑TACCTGTGGGAGCGGGCAAGCCCGCGAAGGGCCGCAAAGCGGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCg         >  5:150592/1‑149 (MQ=17)
                                                                                                                                       gTA‑ACCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTa  >  1:212624/1‑148 (MQ=255)
                                                                                                                                              acaGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGa                                                      >  7:286749/1‑90 (MQ=255)
                                                                                                                                              acaGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGa                                                      <  8:286749/90‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                        
TCATGAGCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCAC‑TGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTA‑ATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTA  >  NC_002947/2667918‑2668199

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: ‑