Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,493,580 A→C 8.3% M42L (ATG→CTG) ‡ PP_2186 → formate dehydrogenase subunit delta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,493,5800AC 8.3% 142.9 / 4.0 48M42L (ATG→CTG) ‡PP_2186formate dehydrogenase subunit delta
Reads supporting (aligned to +/- strand):  ref base A (21/23);  new base C (2/2);  total (23/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 = 9.66e-01

CGGAGAAAGCCGAGGTTCGCCGATGAGCAGTGAAAACCTGGTCAAGATGGCCAACCAGATTGCCCATTACTTCGACAGCGAGCCCGACCGGGCGCTGGCGGTGCAAGGGGTGCGGCAGCATTTGCAGAGCTTCTGGACGCCGGCGATGCGCCGGCAGCTGGCGGAATGGAGCGAGGCCCATGCGG‑‑GGGAGGGGCTTGATCC‑‑GAAGGTGCGGGAGGCGTTGGTCTAGACTCCCGGGCGGCAACGGAACCTTGGGGCGGGTTTGCCGGCGAATGTGGCGTTGAAATCAA  >  NC_002947/2493435‑2493721
                                                                                                                                                 |                                                                                                                                                 
cGGAGAAAGCCGAGGTTCGCCGATGAGCAGTGAAAACCTGGTCAAGATGGCCAACCAGATTGCCCATTACTTCGACAGCGAGCCCGACCGGGCGCTGGCGGTGCAAGGGGTGCGGCAGCATTTGCAGAGCTTCTGGACGCCGGCGATgc                                                                                                                                                <  3:163576/149‑1 (MQ=255)
   agaAAGCCGAGGTTCGCCGATGAGCAGTGAAAACCTGGTCAAGATGGCCAACCAGATTGCCCATTACTTCGACAGCGAGCCCGACCGGGCGCTGGCGGTGCAAGGGGTGCGGCAGCATTTGCAGAGCTTCTGGACGCCGGCGATGCGcc                                                                                                                                             >  1:50150/1‑149 (MQ=255)
         ccGAGGTTCGCCGATGAGCAGTGAAAACCTGGTCAAGATGGCCAACCAGATTGCCCATTACTTCGACAGCGAGCCCGACCGGGCGCTGGCGGTGCAAGGGGTGCGGCAGCATTTGCAGAGCTTCTGGACGCCGGCGATGCGCCGGCAGc                                                                                                                                       >  1:218538/1‑149 (MQ=255)
                  gCCGATGAGCAGTGAAAACCTGGTCAAGATGGCCAACCAGATTGCCCATTACTTCGACAGTGAGCCCGACCGGGCGCTGGCGGTGCAAGGGGTGCGGCAGCATTTGCAGAGCTTCTGGACGCCGGCGATGCGCCGGCAGCTGGCGGaa                                                                                                                               <  1:4346/148‑1 (MQ=255)
                       tGAGCAGTGAAAACCTGGTCAAGATGGCCAACCAGATTGCCCATTACTTCGACAGCGAGCCCGACCGGGCGCTGGCGGTGCAAGGGGTGCGGCAGCATTTGCAGAGCTTCTGGACGCCGGCGATGCGCCGGCAGCTGGCGGAATGGAGc                                                                                                                         <  2:50150/149‑1 (MQ=255)
                            aGTGAAAACCTGGTCAAGATGGCCAACCAGATTGCCCATTACTTCGACAGCGAGCCCGACCGGGCGCTGGCGGTGCAAGGGGTGCGGCAGCATTTGCAGAGCTTCTGGACGCCGGCGATGCGCCGGCAGCTGGCGGAATGGAGCGAGGc                                                                                                                    >  6:266785/1‑149 (MQ=255)
                                     cTGGTCAAGATGGCCAACCAGATTGCCCATTACTTCGACAGCGAGCCCGACCGGGCGCTGGCGGTGCAAGGGGTGCGGCAGCATTTGCAGAGCTTCTGGACGCCGGCGATGCGCCGGCAGCTGGCGGAATGGAGCGAGGCCCATGCgg‑‑g                                                                                                         >  1:174293/1‑149 (MQ=255)
                                     cTGGTCAAGATGGCCAACCAGATTGCCCATTACTTCGACAGCGAGCCCGACCGGGCGCTGGCGGTGCAAGGGGTGCGGCAGCATTTGCAGAGCTTCTGGACGCCGGCGATGCGCCGGCAGCTGGCGGAATGGAGCGAGGCCCATGCgg‑‑g                                                                                                         <  2:174293/149‑1 (MQ=255)
                                        gTCAAGATGGCCAACCAGATTGCCCATTACTTCGACAGCGAGCCCGACCGGGCGCTGGCGGTGCAAGGGGTGCGGCAGCATTTGCAGAGCTTCTGGACGCCGGCGATGCGCCGGCAGCTGGCGGAATGGAGCGAGGCCCATGCgg‑‑gg                                                                                                        <  5:266785/147‑1 (MQ=255)
                                          agtaaTGGGCAATCTGGTTGGCCATTACTTCGACAGCGAGCCCGACCGGGCGCTGGCGGTGCAAGGGGTGCGGCAGCATTTGCAGAGCTTCTGGACGCCGGCGATGCGCCGGCAGCTGGCGGAATGGAGCGAGGCCCATGCGG‑‑GGGAgg                                                                                                    >  6:18498/5‑149 (MQ=255)
                                              aTGGCCAACCAGATTGCCCATTACTTCGACAGCGAGCCCGACCGGGCGCTGGCGGTGCAAGGGGTGCGGCAGCATTTGCAGAGCTTCTGGACGCCGGCGATGCGCCGGCAGCTGGCGGAATGGAGCGAGGCCCATGCGG‑‑GGGAGGGGCt                                                                                                <  3:31901/149‑1 (MQ=255)
                                                 gCCAACCAGATTGCCCATTACTTCGACAGCGAGCCCGACCGGGCGCTGGCGGTGCAAGGGGTGCGGCAGCATTTGCAGAGCTTCTGGACGCCGGCGATGCGCCGGCAGCTGGCGGAATGGAGCGAGGCCCATGCGG‑‑GGGAGGGGCTTg                                                                                              <  5:10416/148‑1 (MQ=255)
                                                 gCCAACCAGATGGCCGATTACTTCGACAGCGAGCCCGACCGGGCGCTGGCGGTGCAAGGGGTGCGGCAGCATTTGCAGAGCTTCTGGACGCCGGCGATGCGCCGGCAGCTGGCGGAATGGAGCGAGGCCCATGCGG‑‑GGGAGGGGCTTGa                                                                                             <  4:145014/149‑1 (MQ=255)
                                                    aaCCAGATTGCCCATTACTTCGACAGCGAGCCCGACCGGGCGCTGGCGGTGCAAGGGGTGCGGCAGCATTTGCAGAGCTTCTGGACGCCGGCGATGCGCCg                                                                                                                                            <  6:102015/101‑1 (MQ=255)
                                                    aaCCAGATTGCCCATTACTTCGACAGCGAGCCCGACCGGGCGCTGGCGGTGCAAGGGGTGCGGCAGCATTTGCAGAGCTTCTGGACGCCGGCGATGCGCCg                                                                                                                                            >  5:102015/1‑101 (MQ=255)
                                                               ccATTACTTCGACAGCGAGCCCGACCGGGCGCTGGCGGTGCAAGGGGTGCGGCAGCATTTGCAGAGCTTCTGGACGCCGGCGATGCGCCGGCAGCTGGCGGAATGGAGCGAGGCCCATGCGG‑‑GGGAGGGGCTTGATCC‑‑GAAGGTGCgg                                                                              >  2:197441/1‑148 (MQ=255)
                                                                 aTTACTTCGACAGCGAGCCCGACCGGGCGCTGGCGGTGCAAGGGGTGCGGCAGCATTTGCAGTGCTTCTGGACGCCGGCGATGCGCCGGCAGCTGGCGGAATGGAGCGAGGCCCATGCGG‑‑GGGAGGGGCTTGATCC‑‑GAAgg                                                                                   >  5:20648/1‑141 (MQ=255)
                                                                 aTTACTTCGACAGCGAGCCCGACCGGGCGCTGGCGGTGCAAGGGGTGCGGCAGCATTTGCAGTGCTTCTGGACGCCGGCGATGCGCCGGCAGCTGGCGGAATGGAGCGAGGCCCATGCGG‑‑GGGAGGGGCTTGATCC‑‑GAAgg                                                                                   <  6:20648/141‑1 (MQ=255)
                                                                     cTTCGACAGCGAGCCCGACCGGGCGCTGGCGGTGCAAGGGGTGCGGCAGCATTTGCAGAGCTTCTGGACGCCGGCGATGCGCCGGCAGCTGGCGGAATGGAGCGAGGCCCATGCGG‑‑GGGAGGGGCTTGATCC‑‑GAAGGTGCGGGAGGCGt                                                                       >  5:100219/1‑149 (MQ=255)
                                                                      ttCGACAGCGAGCCCGACCGGGCGCTGGCGGTGCAAGGGGTGCGGCAGCATTTGCAGAGCTTCTGGACGCCGGCGATGCGCCGGCAGCTGGCGGAATGGAGCGAGGCCCATGCGG‑‑GGGAGGGGCTTGATCC‑‑GAAGGTGCGGGAGGCGtt                                                                      >  1:296568/1‑149 (MQ=255)
                                                                      ttCGACAGCGAGCCCGACCGGGCGCTGGCGGTGCAAGGGGTGCGGCAGCATTTGCAGAGCTTCTGGACGCCGGCGATGCGCCGGCAGCTGGCGGAATGGAGCGAGGCCCATGCGG‑‑GGGAGGGGCTTGATCC‑‑GAAGGTGCGGGAGGCGtt                                                                      >  3:183791/1‑149 (MQ=255)
                                                                        cGACCGCGACGCCGACCGGGCGCTTGCGGTGCAAGGGGGGCGGCAGCATTTGCACAGCTTCTGGACGCCTGCGATGCGCCGGCAGCTGGCGGAATGGAGCGAGGCCCATGCGG‑‑GGGAGGGGCTTGATCC‑‑GAAGGTGCGGGAGGCGt                                                                       <  4:184134/146‑1 (MQ=255)
                                                                            aGCGAGCCCGACCGGGCGCTGGCGGTGCAAGGGGTGCGGCAGCATTTGCAGAGCTTCTGGACGCCGGCGATGCGCCGGCAGCTGGCGGAATGGAGCGAGGCCCATGCGG‑‑GGGAGGGGCTTGATCC‑‑GAAGGTGCGGGAGGCGTTGGTCTa                                                                <  4:183791/149‑1 (MQ=255)
                                                                            aGCGAGCCCGACCGGGCGCTGGCGGTGCAAGGGGTGCGGCAGCATTTGCAGAGCTTCTGGACGCCGGCGATGCGCCGGCAGCTGGCGGAATGGAGCGAGGCCCATGCGG‑‑GGGAGGGGCTTGATCC‑‑GAAGGTGCGGGAGGCGTTGGTCTa                                                                <  5:18498/149‑1 (MQ=255)
                                                                            aGCGAGCCCGACCGGGCGCTGGCGGTGCAAGGGGTGCGGCAGCATTTGCAGAGCTTCTGGACGCCGGCGATGCGCCGGCAGCTGGCGGAATGGAGCGAGGCCCATGCGG‑‑GGGAGGGGCTTGATCC‑‑GAAGGTGCGGGAGGCGTTGGTCTa                                                                <  5:141546/149‑1 (MQ=255)
                                                                                    cGACCGGGCGCTGGCGGTGCAAGGGGTGCGGCAGCATTTGCAGAGCTTCTGGACGCCGGCGATGCGCCGGCAGCTGGCGGAATGGAGCGAGGCCCATGCGG‑‑GGGAGGGGCTTGATc                                                                                           <  4:48573/116‑1 (MQ=255)
                                                                                    cGACCGGGCGCTGGCGGTGCAAGGGGTGCGGCAGCATTTGCAGAGCTTCTGGACGCCGGCGATGCGCCGGCAGCTGGCGGAATGGAGCGAGGCCCATGCGG‑‑GGGAGGGGCTTGATc                                                                                           >  3:48573/1‑116 (MQ=255)
                                                                                        cGGGCGCTGGCGGTGCAAGGGGTGCGGCAGCATTTGCAGAGCTTCTGGACGCCGGCGATGCGCCGGCAGCTGGCGGAATGGAGCGAGGCCCATGCGG‑‑GGGAGGGGCTTGATCC‑‑GAATGTGCGGGAGGCGTTTGTCTAGACTCCCGggcg                                                    <  2:218538/149‑1 (MQ=255)
                                                                                                 gCGGTGCAAGGGGTGCGGCAGCATTTGCAGAGCTTCTGGACGCCGGCGATGCGCCGGCAGCTGGCGGAATGGAGCGAGGCCCATGCGG‑‑GGGAGGGGCTTGATCC‑‑GAAGGTGCGGGAGGCGTTGGTCTAGACTCCCGGGCGGCAACGGaa                                           >  6:97954/1‑149 (MQ=255)
                                                                                                       cAAGGGGTGCGGCAGCATTTGCAGAGCTTCTGGACGCCGGCGATGCGCCGGCAGCTGGCGGAATGGAGCGAGGCCCATGCGG‑‑GGGAGGGGCTTGATCC‑‑GAAGGTGCGGGAGGCGTTGGTCTAGACTCCCGGGCGGCAACGGAACCTTgg                                     <  6:100219/149‑1 (MQ=255)
                                                                                                        aaGGGGTGCGGCAGCATTTGCAGAGCTTCTGGACGCCGGCGATGCGCCGGCAGCTGGCGGAATGGAGCGAGGCCCATGCGG‑‑GGGAGGGGCTTGATCC‑‑GAAGGTGCGGGAGGCGTTGGTCTAGACTCCCGGGCGGCAACGGAACCTTggg                                    <  2:296568/149‑1 (MQ=255)
                                                                                                               gCGGCAGCATTTGCAGAGCTTCTGGACGCCGGCGATGCGCCGGCAGCTGGCGGAATGGAGCGAGGCCCATGCGG‑‑GGGAGGGGCTTGATCC‑‑GAAGGTGCGGGAGGCGTTGGTCTAGACTCCCGGGCGGCAACGGAACCTTGGGGCGGGtt                             >  4:175165/1‑149 (MQ=255)
                                                                                                                  gcagcaTTTGCAGAGCTTCTGGACGCCGGCGATGCGCCGGCAGCTGGCGGAATGGAGCGAGGCCCATGCGGATGAGGGGGGCTTGATCCTGGAAGGTGGGGGAGGCGTTGGTCTAGAATCCCGGGGGGCAACGTACCCTTGGGGCGGGt                              >  1:257004/1‑149 (MQ=255)
                                                                                                                  gcagcaTTTGCAGAGCTTCTGGACGCCGGCGATGCGCCGGCAGCTGGCGGAATGGAGCGAGGCCCATGCGG‑‑GGGAGGGGCTTGATCC‑‑GAAGGTGCGGGAGGCGTTGGTCTAGACTCCCGGGCGGCAACGGAACCTTGGGGCggg                               <  4:75239/144‑1 (MQ=255)
                                                                                                                  gcagcaTTTGCAGAGCTTCTGGACGCCGGCGATGCGCCGGCAGCTGGCGGAATGGAGCGAGGCCCATGCGG‑‑GGGAGGGGCTTGATCC‑‑GAAGGTGCGGGAGGCGTTGGTCTAGACTCCCGGGCGGCAACGGAACCTTGGGGCggg                               >  3:75239/1‑144 (MQ=255)
                                                                                                                  gcagcaTTTGCAGAGCTTCTGGACGCCGGCGATGCGCCGGCAGCTGGCGGAATGGAGCGAGGCCCATGCGG‑‑GGGAGGGGCTTGATCC‑‑GAAGGTGCGGGAGGCGTTGGTCTAGACTCCCGGGCGGCAACGGAACCTTGGGGCGGGTTTGc                          >  3:120625/1‑149 (MQ=255)
                                                                                                                   cagcaTTTGCAGAGCTTCTGGACGCCGGCGATGCGCCGGCAGCTGGCGGAATGGAGCGAGGCCCATGCGG‑‑GGGAGGGGCTTGATCC‑‑GAAGGTGCGGGAGGCGTTGGTCTAGACTCCCGGGCGGCAACGGAACCTTGGGGCGGGTTTGcc                         <  3:175165/149‑1 (MQ=255)
                                                                                                                       aTTTGCAGAGCTTCTGGACGCCGGCGATGCGCCGGCAGCTGGCGGAATGGAGCGAGGCCCATGCGG‑‑GGGAGGGGCTTGATCC‑‑GAAGGTGCGGGAGGCGTTGGTCTAGACTCCCGGGCGGCAACGGAACCTTGGGGCGGGTTTGCCGGc                      >  4:131075/1‑148 (MQ=255)
                                                                                                                                 cTTCTGGACGCCGGCGATGCGCCGGCAGCTGGCGGAATGGAGCGAGGCCCATGCGG‑‑GGGAGGGGCTTGATCC‑‑GAAGGTGCGGGAGGCGTTGGTCTAGACTCCCGGGCGGCAACGGAACCTTGGGGCGGGTTTGCCGGCGAATGTGGCGt           <  1:197441/149‑1 (MQ=255)
                                                                                                                                      ggACGCCGGCGATGCGCCGGCAGCTGGCGGAATGGAGCGAGGCCCATGCGG‑‑GGGAGGGGCTTGATCC‑‑GAAGGTGCGGGAGGCGTTGGTCTAGACTCCCGGGCGGCAACGGAACCTTGGGGCGGGTTTGCCGGCGAATGTGGCGTTGaaa      >  8:100243/1‑149 (MQ=255)
                                                                                                                                      ggACGCCGGCGATGCGCCGGCAGCTGGCGGAATGGAGCGAGGCCCATGCGG‑‑GGGAGGGGCTTGATCC‑‑GAAGGTGCGGGAGGCGTTGGTCTAGACTCCCGGGCGGCAACGGAACCTTGGGGCGGGTTTGCCGGCGAATGTGGCGTTGaaa      >  2:168973/1‑149 (MQ=255)
                                                                                                                                      agctgCCGGCGCATCGCCGGCAGCTGGCGGAATGGAGCGAGGCCCATGCGG‑‑GGGAGGGGCTTGATCC‑‑GAAGGTGCGGGAGGCGTTGGTCTAGACTCCCGGGCGGCAACGGAACCTTGGGGCGGGTTTGCCGGCGAATGTGGCGTTGaaa      <  7:275535/145‑1 (MQ=255)
                                                                                                                                       gctgCCGGCGCATCGCCGGCAGCTGGCGGAATGGAGCGAGGCCCATGCGG‑‑GGGAGGGGCTTGATCC‑‑GAAGGTGCGGGAGGCGTTGGTCTAGACTCCCGGGCGGCAACGGAACCTTGGGGCGGGTTTGCCGGCGAATGTGGCGTTGaaa      >  8:275535/4‑148 (MQ=255)
                                                                                                                                        ctgCCGGCGCATCGCCGGCAGCTGGCGGAATGGAGCGAGGCCCATGCGG‑‑GGGAGGGGCTTGATCC‑‑GAAGGTGCGGGAGGCGTTGGTCTAGACTCCCGGGCGGCAACGGAACCtt                                       >  6:87180/3‑114 (MQ=255)
                                                                                                                                        ctgCCGGCGCATCGCCGGCAGCTGGCGGAATGGAGCGAGGCCCATGCGG‑‑GGGAGGGGCTTGATCC‑‑GAAGGTGCGGGAGGCGTTGGTCTAGACTCCCGGGCGGCAACGGAACCtt                                       <  5:87180/112‑1 (MQ=255)
                                                                                                                                         cGCCGGCGATGCGCCGGCAGCTGGCGGAATGGAGCGAGGCCCATGCGG‑‑GGGAGGGGCTTGATCC‑‑GAAGGTGCGGGAGGCGTTg                                                                     <  6:120506/83‑1 (MQ=255)
                                                                                                                                         cGCCGGCGATGCGCCGGCAGCTGGCGGAATGGAGCGAGGCCCATGCGG‑‑GGGAGGGGCTTGATCC‑‑GAAGGTGCGGGAGGCGTTg                                                                     >  5:120506/1‑83 (MQ=255)
                                                                                                                                          gCCGGCGATGCGCCGGCAGCTGGCGGAATGGAGCGAGGCCCATGCGG‑‑GGGAGGGGCTTGATCC‑‑GAAGGTGCGGGAGGCGTTGGTCTAGACTCCCGGGCGGCAACGGAACCTTGGGGCGGGTTTGCCGGCGAATGTGGCGTTGAAATCaa  <  1:168973/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
CGGAGAAAGCCGAGGTTCGCCGATGAGCAGTGAAAACCTGGTCAAGATGGCCAACCAGATTGCCCATTACTTCGACAGCGAGCCCGACCGGGCGCTGGCGGTGCAAGGGGTGCGGCAGCATTTGCAGAGCTTCTGGACGCCGGCGATGCGCCGGCAGCTGGCGGAATGGAGCGAGGCCCATGCGG‑‑GGGAGGGGCTTGATCC‑‑GAAGGTGCGGGAGGCGTTGGTCTAGACTCCCGGGCGGCAACGGAACCTTGGGGCGGGTTTGCCGGCGAATGTGGCGTTGAAATCAA  >  NC_002947/2493435‑2493721

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