Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,671,109 T→A 8.1% P403P (CCA→CCT)  acnB ← aconitate hydratase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,671,1090TA 8.1% 140.8 / 3.0 49P403P (CCA→CCT) acnBaconitate hydratase B
Reads supporting (aligned to +/- strand):  ref base T (22/23);  new base A (2/2);  total (24/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.96e-01

CGGATAGGCCGCGGTGTGGCAGAACGACTGCATGACCAGGTCAGCAGAGAAGCCCAGGCATGCCA‑GGTCTTTCAGTTCGTCGCGCGTCATCGGGCCAGTGGTGTCCTGGGAGCCGACGGTGGTCATTTTTGGCTCGCAGTAGGCGCCTGGGCGCACGCCCTGGCCTTCTGGCAGGCCACAGGCACGGCCGACCATTTTCTGCGCCAGGGTGAAGCCCTTGCCGGAATCGGCAGGCTGCTCTGGTTTCTTGAACAGGTCGGAAGCACCCAGGCCCAGCTCGGCACGGGCTTTTTC  >  NC_002947/2670962‑2671255
                                                                                                                                                    |                                                                                                                                                  
cGGATAGGCCGCGGTGTGGCAGAACGACTGCATGACCAGGTCAGCAGAGAAGCCCAGGCATGCCA‑GGTCTTTCAGTTCGTCGCGCGTCATCGGGCCAGTGGTGTCCTGGGAGCCGACGGTGGTCATTTTTGGCTCGCAGTAGGCGCCTg                                                                                                                                                   >  2:49463/1‑149 (MQ=255)
        ccGCGGTGTGGCAGAACGACTGCATGACCAGGTCAGCAGAGAAGCCCAGGCATGCCA‑GGTCTTTCAGTTCGTCGCGCGTCATCGGGCCAGTGGTGTCCTGGGAGCCGACGGTGGTCATTTTTGGCTCGCAGTAGGCGCCTGGGCGCACg                                                                                                                                           >  3:307902/1‑149 (MQ=255)
           cgGTGTGGCAGAACGACTGCATGACCAGGTCAGCAGAGAAGCCCAGGCATGCCA‑GGTCTTTCAGTTCGTCGCGCGTCATCGGGCCAGTGGTGTCCTGGGAGCCGACGGTGGTCATTTTTGGCTCGCAGTAGGCGCCTGGGCGCACGccc                                                                                                                                        <  2:303504/149‑1 (MQ=255)
           cgGTGTGGCAGAACGACTGCATGACCAGGTCAGCAGAGAAGCCCAGGCATGCCA‑GGTCTTTCAGTTCGTCGCGCGTCATCGGGCCAGTGGTGTCCTGGGAGCCGACGGTGGTCATTTTTGGCTCGCAGTAGGCGCCTGGGCGCACGccc                                                                                                                                        <  6:158737/149‑1 (MQ=255)
              tgtgGCAGAACGACTGCATGACCAGGTCAGCAGAGAAGCCCAGGCATGCCA‑GGTCTTTCAGTTCGTCGCGCGTCATCGGGCCAGTGGTGTCCTGGGAGCCGACGGTGGTCATTTTTGGCTCGCAGTAGGCGCCTGGGCGCACGCCCTgg                                                                                                                                     <  4:307902/149‑1 (MQ=255)
                tgGCAGAAGGACTGCATGACCAGGTCAGCAGAGAAGGCCAGGCATGCCA‑GGTCTTTCAGTTCGTCGCGCGTCATCGGGCCAGTGGTGTCCTGTGAGCCGACGGTGGTCATTTTTGGCTCGCAGTAGGCGCCTGGGCGCACGCCCTGGcc                                                                                                                                   <  4:317083/149‑1 (MQ=255)
                  gCAGAACGACTGCATGACCAGGTCAGCAGAGAAGCCCAGGCATGCCA‑GGTCTTTCAGTTCGTCGCGCGTCATCGGGCCAGTGGTGTCCTGGGAGCCGACGGTGGTCATTTTTGGCTCGCAGTAGGCGCCTGGGCGCACGCCCTGGCCtt                                                                                                                                 >  6:113632/1‑149 (MQ=255)
                  gCAGAACGACTGCATGACCAGGTCAGCAGAGAAGCCCAGGCATGCCA‑GGTCTTTCAGTTCGTCGCGCGTCATCGGGCCAGTGGTGTCCTGGGAGCCGACGGTGGTCATTTTTGGCTCGCAGTAGGCGCCTGGGCGCACGCCCTGGCCtt                                                                                                                                 <  6:39877/149‑1 (MQ=255)
                  gCAGAACGACTGCATGACCAGGTCAGCAGAGAAGCCCAGGCATGCCA‑GGTCTTTCAGTTCGTCGCGCGTCATCGGGCCAGTGGTGTCCTGGGAGCCGACGGTGGTCATTTTTGGCTCGCAGTAGGCGCCAGGGCGTGCGCCCAGGc                                                                                                                                    <  8:43399/146‑1 (MQ=255)
                  gCAGAACGACTGCATGACCAGGTCAGCAGAGAAGCCCAGGCATGCCA‑GGTCTTTCAGTTCGTCGCGCGTCATCGGGCCAGTGGTGTCCTGGGAGCCGACGGTGGTCATTTTTGGCTCGCAGTAGGCGCCAGGGCGTGCGCCCAGGc                                                                                                                                    >  7:43399/1‑146 (MQ=255)
                              cATGACCAGGTCAGCAGAGAAGCCCAGGCATGCCA‑GGTCTTTCAGTTCGTCGCGCGTCATCGGGCCAGTGGTGTCCTGGGAGCCGACGGTGGTCATTTTTGGCTCGCAGTAGGCGCCTGGGCGCACGCCCTGGCCTTCTGGCAGGCca                                                                                                                      >  4:191309/1‑148 (MQ=255)
                              cATGACCAGGTCAGCAGAGAAGCCCAGGCATGCCA‑GGTCTTTCAGTTCGTCGCGCGTCATCGGGCCAGTGGTGTCCTGGGAGCCGACGGTGGTCATTTTTGGCTCGCAGTAGGCGCCTGGGCGCACGCCCTGGCCTTCTGGCAGGCca                                                                                                                      >  7:64380/1‑148 (MQ=255)
                              cATGACCAGGTCAGCAGAGAAGCCCAGGCATGCCA‑GGTCTTTCAGTTCGTCGCGCGTCATCGGGCCAGTGGTGTCCTGGGAGCCGACGGTGGTCATTTTTGGCTCGCACTAGGCGCCTGGGCGCACGCCCTGGCCTTCTGGCAGGCca                                                                                                                      >  1:236128/1‑148 (MQ=255)
                               aTGACCAGGTCAGCAGAGAAGCCCAGGCATGCCAGGGTCTTTCAGTTCGTCGCGCGTCATCGGGCCAGTGGTGTCCTGGGAGCCGACGGTGGTCATTTTTGGCTCGCAGTAGGCGCCTGGGCGCACGCCCTGGCCTTCTGGCAGGCca                                                                                                                      >  1:231585/1‑148 (MQ=255)
                               aTGACCAGGTCAGCAGAGAAGCCCAGGCATGCCA‑GGTCTTTCAGTTCGTCGCGCGTCATCGGGCCAGTGGTGTCCTGGGAGCCGACGGTGGTCATTTTTGGCTCGCAGTAGGCGCCTGGGCGCACGCCCTGGCCTTCTGGCAGGCca                                                                                                                      >  6:195026/1‑147 (MQ=255)
                                   ccAGGTCAGCAGAGAAGCCCAGGCATGCCA‑GGTCTTTCAGTTCGTCGCGCGTCATCGGGCCAGTGGTGTCCTGGGAGCCGACGGTGGTCATTTTTGGCTCGCAGTAGGCGCCTGGGCGCACGCCCTGGCCTTCTGGCAGGCCACAGGCa                                                                                                                >  1:163854/1‑149 (MQ=255)
                                      ggTCAGCAGAGAAGCCCAGGCATGCCA‑GGTCTTTCAGTTCGTCGCGCGTCATCGGGCCAGTGGTGTCCTGGGAGCCGACGGTGGTCATTTTTGGCTCGCAGTAGGCGCCTGGGCGCACGCCCTGGCCTTCTGGCAGGCCACAGGCACgg                                                                                                             >  8:332725/1‑149 (MQ=255)
                                          agcagAGAAGCCCAGGCATGCCA‑GGTCTTTCAGTTCGTCGCGCGTCATCGGGCCAGTGGTGTCCTGGGAGCCGACGGTGGTCATTTTTGGCTCGCAGTAGGCGCCTGGGCGCACGCCCTGGCCTTCt                                                                                                                               >  7:96251/1‑127 (MQ=255)
                                          agcagAGAAGCCCAGGCATGCCA‑GGTCTTTCAGTTCGTCGCGCGTCATCGGGCCAGTGGTGTCCTGGGAGCCGACGGTGGTCATTTTTGGCTCGCAGTAGGCGCCTGGGCGCACGCCCTGGCCTTCt                                                                                                                               <  8:96251/127‑1 (MQ=255)
                                            cagAGAAGCCCAGGCATGCCA‑GGTCTTTCAGTTCGTCGCGCGTCATCGGGCCAGTGGTGTCCTGGGAGCCGACGGTGGTCATTTTTGGCTCGCAGTAGGCGCCTGGGCGCACGCCCTGGCCTTCTGGCAg                                                                                                                          <  1:212100/130‑1 (MQ=255)
                                            cagAGAAGCCCAGGCATGCCA‑GGTCTTTCAGTTCGTCGCGCGTCATCGGGCCAGTGGTGTCCTGGGAGCCGACGGTGGTCATTTTTGGCTCGCAGTAGGCGCCTGGGCGCACGCCCTGGCCTTCTGGCAg                                                                                                                          >  2:212100/1‑130 (MQ=255)
                                               agaAGCCCAGGCATGCCA‑GGTCTTTCAGTTCGTCGCGCGTCATCGGGCCAGTGGTGTCCTGGGAGCCGACGGTGGTCATTTTTGGCTCGCAGTAGGCGCCTGGGCGCACGCCCTGGCCTTCTGGCAGGCCACAGGCACGGCCGACCAtt                                                                                                    <  1:49463/149‑1 (MQ=255)
                                                      cAGGCATGCCA‑GGTCTTTCAGTTCGTCGCGCGTCATCGGGCCAGTGGTGTCCTGGGAGCCGACGGTGGTCATTTTTGGCTCGCAGTAGGCGCCTGGGCGCACGCCCTGGCCTTCTGGCAGGCCACAGGCACGGCCGACCATTTTCTgcg                                                                                             >  3:338465/1‑149 (MQ=255)
                                                                     cTTTCAGTTCGTCGCGCGTCATCGGGCCAGTGGTGTCCTGGGAGCCGACGGTGGTCATTTTTGGCTCGCAGTAGGCGCCTGGGCGCACGCCCTGGCCTTCTGGCAGGCCACAGGCACGGCCGACCATTTTCTGCGCCAGGGTGAAGccc                                                                               >  4:335537/1‑149 (MQ=255)
                                                                      tttCAGTTCGTCGCGCGTCATCGGGCCAGTGGTGTCCTGGGAGCCGACGGTGGTCATTTTTGGCTCGCAGTAGGCGCCTGGGCGCACGCCCTGGCCTTCTGGCAGGCCACAGGCACGGCCGACCATTTTCTGCGCCAGGGTGAAGCCCt                                                                              >  5:276460/1‑149 (MQ=255)
                                                                           gTTCGTCGCGCGTCATCGGGCCAGTGGTGTCCTGGGAGCCGACGGTGGTCATTTTTGGCTCGCAGTAGGCGCCTGGGCGCACGCCCTGGCCTTCTGGCAGGCCACAGGCACGGCCGACCATTTTCTGCGCCAGGGTGAAGCCCTTGCCg                                                                         <  3:351912/149‑1 (MQ=255)
                                                                                     cgTCATCGGGCCAGTGGTGTCCTGGGAGCCGACGGTGGTCATTTTTGGCTCGCAGTAGGCGCCTGGGCGCACGCCCTGGCCTTCTGGCAGGCCACAGGCACGGCCGACCATTTTCTGCGCCAGGGTGAAGCCCTTGCCGGAATCGGCAg                                                               <  1:30118/149‑1 (MQ=255)
                                                                                     cgTCATCGGGCCAGTGGTGTCCTGGGAGCCGACGGTGGTCATTTTTGGCTCGCAGTAGGCGCCTGGGCGCACGCCCTGGCCTTCTGGCAGGCCACAGGCACGGCCGACCATTTTCTGCGCCAGGGTGAAGCCCTTGCCGGAATCGGCAg                                                               <  2:231585/149‑1 (MQ=255)
                                                                                      gTCATCGGGCCAGTGGTGTCCTGGGAGCCGACGGTGGTCATTTTTGGCTCGCAGTAGGCGCCTGGGCGCACGCCCTGGCCTTCTGGCAGGCCACAGGCACGGCCGA‑CATTTTCTGCGCCAGGGTGAAGCCCTTGCCGGAATCGGCGGgc                                                             >  3:338187/1‑149 (MQ=255)
                                                                                                cAGTGGTGTCCTGGGAGCCGACGGTGGTCATTTTTGGCTCGCAGTAGGCGCCTGGGCGCACGCCCTGGCCTTCTGGCAGGCCACAGGCACGGCCGACCATTTTCTGCGCCAGGGTGAAGCCCTTGCCGGAATCGGCAGGCTGCTCTGGt                                                    <  4:338465/149‑1 (MQ=255)
                                                                                                cAGTGGTGTCCTGGGAGCCGACGGTGGTCATTTTTGGCTCGCAGTAGGCGCCTGGGCGCACGCCCTGGCCTTCTGGCAGGCCACAGGCACGGCCGACCATTTTCTGCGCCAGGGTGAAGCCCTTGCCGGAATCGGCAGGCTGCTCTGGt                                                    <  7:332725/149‑1 (MQ=255)
                                                                                                  gtggtgTCCTGGGAGCCGACGGTGGTCATTTTTGGCTCGCAGTAGGCGCCTGGGCGCACGCCCTGGCCTTCTGGCAGGCCACAGGCACGGCCGACCATTTTCTGCGCCAGGGTGAAGCCCTTGCCGGAATCGGCAGGCTGCTCTGGttt                                                  <  2:163854/149‑1 (MQ=255)
                                                                                                     gtgtCCTGGGAGCCGAGGGTGGTCATTTTTGGCTCGCTGTAGGCGGCTGGGCGGACGCCCTGGCCTTCTGGCAGGCCACAGGCACGGCCGACCATTTTCTGCGCCAGGGTGAAGCCCTTGCCGGAATCGGCAGGCTGCTCTGGTTTCtt                                               <  2:236128/149‑1 (MQ=255)
                                                                                                      tgtCCTGGGAGCCGACGGTGGTCATTTTTGGCTCGCAGTAGGCGCCTGGGCGCACGCCCTGGCCTTCTGGCAGGCCACAGGCACGGCCGACCATTTTCTGCGCCAGGGTGAAGCCCTTGCCGGAATCGGCAGGCTGCTCTGGTTTCTTg                                              <  6:276460/149‑1 (MQ=255)
                                                                                                       gtCCTGGGAGCCGACGGTGGTCATTTTTGGCTCGCAGTAGGCGCCTGGGCGCACGCCCTGGCCTTCTGGCAGGCCACAGGCACGGCCGACCATTTTCTGCGCCAGGGTGAAGCCCTTGCCGGAATCGGCAGGCTGCTCTGGTTTCTTg                                              >  4:123218/1‑148 (MQ=255)
                                                                                                        tCCTGGGAGCCGACGGTGGTCATTTTTGGCTCGCAGTAGGCGCCTGGGCGCACGCCCTGGCCTTCTGGCAGGCCACAGGCACGGCCGACCATTTTCTGCGCCAGGGTGAAGCCCTTGCCGGAATCGGCAGGCTGCTCTGGTTTCTTGaa                                            <  8:64380/149‑1 (MQ=255)
                                                                                                           tGGGAGCCGACGGTGGTCATTTTTGGCTCGCAGTAGGCGCCTGGGCGCACGCCCTGGCCTTCTGGCAGGCCACAGGCACGGCCGACCATTTTCTGCGCCAGGGTGAAGCCCTTGCCGGAATCGGCAGGCTGCTCTGGTTTCTTGAACAg                                         >  1:108839/1‑149 (MQ=255)
                                                                                                              gAGCCGACGGTGGTCATTTTTGGCTCGCAGTAGGCGCCTGGGCGCACGCCCTGGCCTTCTGGCAGGCCACAGGCACGGCCGACCATTTTCTGCGCCAGGGTGAAGCCCTTGCCGGAATCGGCAGGCTGCTCTGGTTTCTTGAACAGGTc                                      <  3:191309/149‑1 (MQ=255)
                                                                                                              gAGCCGACGGTGGTCATTTTTGGCTCGCAGTAGGCGCCTGGGCGCACGCCCTGGCCTTCTGGCAGGCCACAGGCACGGCCGACCATTTTCTGCGCCAGGGTGAAGCCCTTGCCGGAATCGGCAGGCTGCTCTGGTTTCTTGAACAGGTc                                      >  5:196499/1‑149 (MQ=255)
                                                                                                              gAGCCGACGGTGGTCATTTTTGGCTCGCAGTAGGCGCCTGGGCGCACGCCCTGGCCTTCTGGCAGGCCACAGGCACGGCCGACCATTTTCTGCGCCAGGGTGAAGCCCTTGCCGGAATCGGCAGGCTGCTCTGGTTTCTTGAACAGGTc                                      <  5:113632/149‑1 (MQ=255)
                                                                                                                                tttGGCTCGCAGTAGGCGCCTGGGCGCACGCCCTGGCCTTCTGGCAGGCCACAGGCACGGCCGACCATTTTCTGCGCCAGGGTGAAGCCCTTGCCGGAATCGGCAGGCTGCTCTGGTTTCTTGAACAGGTCGGAAGCACCCAGGCCCAg                    >  7:220685/1‑149 (MQ=255)
                                                                                                                                    gCTCGCAGTAGGCGCCTGGGCGCACGCCCTGGCCTTCTGGCAGGCCACAGGCACGGCCGACCATTTTCTGCGCCAGGGTGAAGCCCTTGCCGGAATCGGCAGGCTGCTCTGGTTTCTTGAACAGGTCGGAAGCACCCAGGCCCAGCTCg                <  6:196499/149‑1 (MQ=255)
                                                                                                                                      tCGCAGTAGGCGCCTGGGCGCACGCCCTGGCCTTCTGGCAGGCCACAGGCACGGCCGACCATTTTCTGCGCCAGGGTGAAGCCCTTGCCGGAATCGGCAGGCTGCTCTGGGTTCTTGAACAGGTCGGAAGCACCCAGGCCCAGCTCGGc              >  7:75316/1‑149 (MQ=255)
                                                                                                                                        gCAGTAGGCGCCTGGGCGCACGCCCTGGCCTTCTGGCAGGCCACAGGCACGGCCGACCATTTTCTGCGCCAGGGTGAAGCCCTTGCCGGAATCGGCAGGCTGCTCTGGTTTCTTGAACAGGTCGGAAGCACCCAGGCCCAGCTCGGc              <  2:108839/147‑1 (MQ=255)
                                                                                                                                                 gcCTGGGCGCACGCCCTGGCCTTCTGGCAGGCCACAGGCACGGCCGACCATTTTCTGCGCCAGGGTGAAGCCCTTGCCGGAATCGGCAGGCTGCTCTGGTTTCTTGAACAGGTCGGAAGCACCCAGGCCCAGCTCGGCACGGGCttttt   <  3:335537/149‑1 (MQ=255)
                                                                                                                                                 gcCTGGGCGCACGCCCTGGCCTTCTGGCAGGCCACAGGCACGGCCGACCATTTTCTGCGCCAGGGTGAAGCCCTTGCCGGAATCGGCAGGCTGCTCTGGTTTCTTGAACAGGTCGGAAGCACCCAGGCCCAGCTCGGCACGGGCttttt   >  4:310745/1‑149 (MQ=255)
                                                                                                                                                 gcCAGGGCGTGCGCCCAGGCCTTCTGGCAGGCCACAGGCACGGCCGACCATTTTCTGCGCCAGGGTGAAGCCCTTGCCGGAATCGGc                                                                 >  6:298977/1‑87 (MQ=255)
                                                                                                                                                 gcCAGGGCGTGCGCCCAGGCCTTCTGGCAGGCCACAGGCACGGCCGACCATTTTCTGCGCCAGGGTGAAGCCCTTGCCGGAATCGGc                                                                 <  5:298977/87‑1 (MQ=255)
                                                                                                                                                  ccTGGGCGCTCGCCCTGGCCTTCTGGCAGGCCACAGGCACGGCCGACCATTTTCAGCGCCAGGGTGAAGCCCTTGCCGGAATCGGCAGGCTGCTCTGGTTTCTTGAACAGGTCGGAAGCACCCAGGCCCAGCTCGGCACGGGCTTTTTc  <  4:282294/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CGGATAGGCCGCGGTGTGGCAGAACGACTGCATGACCAGGTCAGCAGAGAAGCCCAGGCATGCCA‑GGTCTTTCAGTTCGTCGCGCGTCATCGGGCCAGTGGTGTCCTGGGAGCCGACGGTGGTCATTTTTGGCTCGCAGTAGGCGCCTGGGCGCACGCCCTGGCCTTCTGGCAGGCCACAGGCACGGCCGACCATTTTCTGCGCCAGGGTGAAGCCCTTGCCGGAATCGGCAGGCTGCTCTGGTTTCTTGAACAGGTCGGAAGCACCCAGGCCCAGCTCGGCACGGGCTTTTTC  >  NC_002947/2670962‑2671255

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