Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,668,054:1 +T 11.9% intergenic (+72/‑95) 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,0541.T11.9% 152.4 / 17.0 42intergenic (+72/‑95)prpF/prpDaconitate isomerase/2‑methylcitrate dehydratase
Reads supporting (aligned to +/- strand):  ref base . (20/17);  new base T (3/2);  total (23/19)
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.44e-01

AGGCAATCATGAGCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCAC‑TGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTA‑ATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTA  >  NC_002947/2667912‑2668199
                                                                                                                                                |                                                                                                                                                 
aGGCAATCATGAGCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCAC‑TGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTA‑ATcaca                                                                                                                                             <  5:30223/149‑1 (MQ=255)
      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‑ATCACAgg                                                                                                                                           >  8:79827/1‑137 (MQ=255)
              cGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCAC‑TGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTA‑ATCACAgg                                                                                                                                           <  7:79827/137‑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)
                                              tCGCGTGCCTGGGGTTAGCTTCTAATGAACCAC‑TGGCGCTGGGCTTCAAGCGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTA‑ATCACAgg                                                                                                                                           <  2:106615/105‑1 (MQ=37)
                                              tCGCGTGCCTGGCGATAGCTTCTAAAGCACCAC‑TGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTA‑ATCACAgg                                                                                                                                           >  1:106615/1‑105 (MQ=255)
                                                gcgTGCCTGGCGATAGCTTCTAAAGCACCACGTGGCGCGCGGCTTCAAGAGGGGCCGCTTTGCGCCCCTTCCCGGGCTTGCCCCCTCCCACAGGTA‑ATAACAGGGTTTGAATGGTGTTATATCCCGGTGGGCGCGGGCAAGCCcgcg                                                                                                >  1:39018/1‑147 (MQ=255)
                                                    gCCTGGCGATAGCTTCTAAAGCACCAC‑TGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACCCGAA‑ATAACAGGGGTTGTATGGGGTGGTATGCCGTGGGGAGCGGGCAAGGCCGCGCAGGGCt                                                                                         >  8:291950/1‑149 (MQ=37)
                                                       tGGCGATAGCTTCTAAAGCACCAC‑TGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGAT‑TACCTGTgg                                                                                                                 <  4:64917/122‑1 (MQ=14)
                                                       tGGCGATAGCTTCTAAAGCACCAC‑TGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGAT‑TACCTGTgg                                                                                                                 >  3:64917/1‑122 (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              >  8:339574/1‑148 (MQ=255)
                                                                                                                                 ccGCTCCCACAGGTA‑ATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTAc              >  4:5565/1‑148 (MQ=255)
                                                                                                                                     tCCCACAGGGATATCACAGCATTCAAACCCTGTGAT‑TACCTGTGGGAGCGGGCAAGCCCGCGAAGGGCCGCAAAGCGGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCg         >  5:150592/1‑149 (MQ=17)
                                                                                                                                             gTA‑ACCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTa  >  1:212624/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
AGGCAATCATGAGCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCAC‑TGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTA‑ATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTA  >  NC_002947/2667912‑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: ‑