Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 664,889 C→T 15.4% A584A (GCC→GCT yfhM → lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947664,8890CT15.4% 113.7 / 12.3 51A584A (GCC→GCTyfhMlipoprotein
Reads supporting (aligned to +/- strand):  ref base C (19/24);  new base T (4/4);  total (23/28)
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.97e-01

GCTACCAATTCGGCTCGGTCACCGAAACCGATCTGAACCAGGACCTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTGGAATTCGAGTTCCTGGTGGCCGA  >  NC_002947/664741‑665038
                                                                                                                                                    |                                                                                                                                                     
gccaCCAATTCGGCTCCGTCGCCGAAACCGATCTGAACCAGGACCTGGAACTGGACGAACTCCCTTTATACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCa                                                                                                                                                      <  1:449422/147‑1 (MQ=255)
  tACCAATTCGGCTCGGTCACCGAAACCGATCTGAACCAGGACCTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATAGCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGc                                                                                                                                                    <  2:564133/150‑1 (MQ=255)
           ggCTCGGTCACCGAAACCGATCTGAACCAGGACCTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAg                                                                                                                                           >  1:1121626/1‑150 (MQ=255)
             cTCGGTCACCTAAACCGATCTGAACCATGACGGGGAACTGGACGAAGTCACCCTTGACCAGGCCGGGAAGGCGCTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCCTGCAGGCTGGCCTGCacggtca                                                                                                                                         <  2:220545/150‑7 (MQ=255)
             cTCGGTCACCGAAACCGATCTGAACCAGGACCTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTc                                                                                                                                         >  2:940841/1‑150 (MQ=255)
                      gtaaaCCGATCTGAACCAGGACCTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTggccg                                                                                                                                >  2:794811/3‑150 (MQ=255)
                          gccGATCTGAACCAGGACCTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGgccggccg                                                                                                                            >  2:1209937/2‑150 (MQ=255)
                                 tGAACCAGGACCTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCAccc                                                                                                                     >  1:514516/1‑150 (MQ=255)
                                    aCCAGGTCCGGGAATTGGACGACGCCACCCGTCTTCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGATCGGCCGGCCACCCgac                                                                                                                  <  1:191799/150‑3 (MQ=255)
                                      cAGGACCTGTAACTGGACGAAGTCACCCTTTTCCGGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGgcgc                                                                                                                <  1:483275/150‑1 (MQ=255)
                                           ccTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGa                                                                                                           <  2:1426570/150‑1 (MQ=255)
                                                 aCTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCgg                                                                                                                                       >  2:1273682/1‑116 (MQ=255)
                                                 aCTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCgg                                                                                                                                       <  1:1273682/116‑1 (MQ=255)
                                                     gACGAAGTCCCCCTTGCCCCGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATc                                                                                                 <  1:265610/150‑1 (MQ=255)
                                                      aCGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCt                                                                                                <  1:1209937/150‑1 (MQ=255)
                                                        gAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCgg                                                                                                                                       >  1:801330/1‑109 (MQ=255)
                                                        gAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCgg                                                                                                                                       <  2:801330/109‑1 (MQ=255)
                                                         aaGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGc                                                                                             <  1:794811/150‑1 (MQ=255)
                                                           gTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCg                                                                                           >  2:917138/1‑150 (MQ=255)
                                                                ccTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGa                                                                                      >  1:296383/1‑150 (MQ=255)
                                                                 cTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAg                                                                                     >  2:123101/1‑150 (MQ=255)
                                                                     aCCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGcc                                                                                 >  2:1412850/1‑150 (MQ=255)
                                                                           cccgaagggCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCgg                                                                           <  1:566679/143‑1 (MQ=255)
                                                                           ccGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCgg                                                                           <  2:146623/150‑1 (MQ=255)
                                                                           ccGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCgg                                                                           <  1:123101/150‑1 (MQ=255)
                                                                                     ggtggtCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAgg                                                                 <  1:917138/150‑1 (MQ=255)
                                                                                      gtggtCGATATCGAAACCCGTTGGGCAGCACACCTTTCGCCGCGGGAAGTTACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGc                                                                <  2:514516/150‑1 (MQ=255)
                                                                                           cGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGtt                                                           <  2:579770/150‑1 (MQ=255)
                                                                                                  gAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCCACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGa                                                                                                                           >  2:520335/1‑79 (MQ=255)
                                                                                                  gAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGa                                                                                                                           <  1:520335/79‑1 (MQ=255)
                                                                                                  gAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGt                                                    <  1:124426/150‑1 (MQ=255)
                                                                                                       ccGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTgagga                                               >  2:292216/1‑150 (MQ=255)
                                                                                                        cGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTgaggag                                              >  2:577548/1‑150 (MQ=255)
                                                                                                                cGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATa                                       >  2:332555/1‑149 (MQ=255)
                                                                                                                cGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATa                                       <  1:332555/149‑1 (MQ=255)
                                                                                                                      ccGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACgg                                <  1:940841/150‑1 (MQ=255)
                                                                                                                               gCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTGGa                       >  1:922598/1‑150 (MQ=255)
                                                                                                                                      cTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATa                                       <  1:881732/127‑1 (MQ=255)
                                                                                                                                      cTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATa                                       >  2:881732/1‑127 (MQ=255)
                                                                                                                                      cTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTGGAATTCGAg                >  1:171222/1‑150 (MQ=255)
                                                                                                                                         ctcctGCAGGCTGGCCTGCACGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTGGAATTCGa                 <  2:246641/142‑1 (MQ=255)
                                                                                                                                         ctcctGCAGGCTGGCCTGCACGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTGGAATTCGa                 >  1:246641/5‑146 (MQ=255)
                                                                                                                                         aCCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTgag                                                 >  2:1250950/1‑114 (MQ=255)
                                                                                                                                         aCCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTgag                                                 <  1:1250950/114‑1 (MQ=255)
                                                                                                                                           cctGCAGGCTGGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTGGAATTCGAGTTCCt           >  1:288356/3‑150 (MQ=255)
                                                                                                                                            ctGCAGGCTGGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTc                                                          >  1:965694/2‑102 (MQ=255)
                                                                                                                                            ctGCAGGCTGGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTc                                                          <  2:965694/101‑1 (MQ=255)
                                                                                                                                             tGCAGGCTGGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATa                                       <  2:1268834/120‑1 (MQ=255)
                                                                                                                                             tGCAGGCTGGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATa                                       >  1:1268834/1‑120 (MQ=255)
                                                                                                                                             tGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTGGAATTCGAGTTCCtgg         <  1:577548/150‑1 (MQ=255)
                                                                                                                                               cAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTGGAATTCGAGTTCCtggtg       <  2:1121626/150‑1 (MQ=255)
                                                                                                                                                    cAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTGGAATTCGAGTTCCTGGTGGCCGa  >  1:485104/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCTACCAATTCGGCTCGGTCACCGAAACCGATCTGAACCAGGACCTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTGGAATTCGAGTTCCTGGTGGCCGA  >  NC_002947/664741‑665038

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