Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,671,116 A→G 8.2% V401A (GTG→GCG) ‡ acnB ← aconitate hydratase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,671,1160AG 8.2% 132.2 / 3.0 50V401A (GTG→GCG) ‡acnBaconitate hydratase B
Reads supporting (aligned to +/- strand):  ref base A (22/23);  new base G (2/2);  total (24/26)
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.42e-01

CCGCGGTGTGGCAGAACGACTGCATGACCAGGTCAGCAGAGAAGCCCAGGCATGCCA‑GGTCTTTCAGTTCGTCGCGCGTCATCGGGCCAGTGGTGTCCTGGGAGCCGACGGTGGTCATTTTTGGCTCGCAGTAGGCGCCTGGGCGCACGCCCTGGCCTTCTGGCAGGCCACAGGCACGGCCGACCATTTTCTGCGCCAGGGTGAAGCCCTTGCCGGAATCGGCAGGCTGCTCTGGTTTCTTGAACAGGTCGGAAGCACCCAGGCCCAGCTCGGCACGGGCTTTTTCGGT  >  NC_002947/2670970‑2671258
                                                                                                                                                   |                                                                                                                                              
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                                                                                                                             >  2:212100/1‑130 (MQ=255)
                                    cagAGAAGCCCAGGCATGCCA‑GGTCTTTCAGTTCGTCGCGCGTCATCGGGCCAGTGGTGTCCTGGGAGCCGACGGTGGTCATTTTTGGCTCGCAGTAGGCGCCTGGGCGCACGCCCTGGCCTTCTGGCAg                                                                                                                             <  1:212100/130‑1 (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                                         >  5:196499/1‑149 (MQ=255)
                                                                                                      gAGCCGACGGTGGTCATTTTTGGCTCGCAGTAGGCGCCTGGGCGCACGCCCTGGCCTTCTGGCAGGCCACAGGCACGGCCGACCATTTTCTGCGCCAGGGTGAAGCCCTTGCCGGAATCGGCAGGCTGCTCTGGTTTCTTGAACAGGTc                                         <  5:113632/149‑1 (MQ=255)
                                                                                                      gAGCCGACGGTGGTCATTTTTGGCTCGCAGTAGGCGCCTGGGCGCACGCCCTGGCCTTCTGGCAGGCCACAGGCACGGCCGACCATTTTCTGCGCCAGGGTGAAGCCCTTGCCGGAATCGGCAGGCTGCTCTGGTTTCTTGAACAGGTc                                         <  3:191309/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      >  4:310745/1‑149 (MQ=255)
                                                                                                                                         gcCTGGGCGCACGCCCTGGCCTTCTGGCAGGCCACAGGCACGGCCGACCATTTTCTGCGCCAGGGTGAAGCCCTTGCCGGAATCGGCAGGCTGCTCTGGTTTCTTGAACAGGTCGGAAGCACCCAGGCCCAGCTCGGCACGGGCttttt      <  3:335537/149‑1 (MQ=255)
                                                                                                                                         gcCAGGGCGTGCGCCCAGGCCTTCTGGCAGGCCACAGGCACGGCCGACCATTTTCTGCGCCAGGGTGAAGCCCTTGCCGGAATCGGc                                                                    >  6:298977/1‑87 (MQ=255)
                                                                                                                                         gcCAGGGCGTGCGCCCAGGCCTTCTGGCAGGCCACAGGCACGGCCGACCATTTTCTGCGCCAGGGTGAAGCCCTTGCCGGAATCGGc                                                                    <  5:298977/87‑1 (MQ=255)
                                                                                                                                          ggatctcgcACGCCCTGGCCTTCTGGCAGGCCACAGGCACGGCCGACCATTTTCTGCGCCAGGGTGAAGCCCTTGCCGGAATCGGCAGGCTGCTCTGGTTTCTTGAACAGGTCGGAAGCACCCAGGCCCAGCTCGGCACGGGCTTTTTc     <  3:343403/143‑1 (MQ=255)
                                                                                                                                          ccTGGGCGCTCGCCCTGGCCTTCTGGCAGGCCACAGGCACGGCCGACCATTTTCAGCGCCAGGGTGAAGCCCTTGCCGGAATCGGCAGGCTGCTCTGGTTTCTTGAACAGGTCGGAAGCACCCAGGCCCAGCTCGGCACGGGCTTTTTc     <  4:282294/149‑1 (MQ=255)
                                                                                                                                              ggCGCACGCCCTGGCCTTCTGGCAGGCCACAGGCACGGCCGACCATTTTCTGCGCCAGGGTGAAGCCCTTGCCGGAATCGGCAGGCTGCTCTGGTTTCTTGAACAGGTCGGAAGCACCCAGGCCCAGCTCGGCACGGGCTTTTTCGGt  >  2:311820/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
CCGCGGTGTGGCAGAACGACTGCATGACCAGGTCAGCAGAGAAGCCCAGGCATGCCA‑GGTCTTTCAGTTCGTCGCGCGTCATCGGGCCAGTGGTGTCCTGGGAGCCGACGGTGGTCATTTTTGGCTCGCAGTAGGCGCCTGGGCGCACGCCCTGGCCTTCTGGCAGGCCACAGGCACGGCCGACCATTTTCTGCGCCAGGGTGAAGCCCTTGCCGGAATCGGCAGGCTGCTCTGGTTTCTTGAACAGGTCGGAAGCACCCAGGCCCAGCTCGGCACGGGCTTTTTCGGT  >  NC_002947/2670970‑2671258

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