Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,840,541 C→G 7.3% E14Q (GAG→CAG)  sseA ← 3‑mercaptopyruvate sulfurtransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,840,5410CG 7.3% 168.9 / 3.0 55E14Q (GAG→CAG) sseA3‑mercaptopyruvate sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base C (26/25);  new base G (2/2);  total (28/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.33e-01

TCACCGGACCGCTGAGATCCCGGTCCAGGTCGGCAAAATGCGCCCCCTCAATATGCCCCAGGGCATAGCTGCGCTGGCCATAGTCCACATCCTCGAGGGCAAAGCGGCAGTCGAGGATCACCAACTTGGGCGAGCCCAGACGCTCGGCCAACTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTC‑CTGATCGATGGCACAAGAGCCCTGGGCTAGAGCTCCAGGGCCTGGTTGAACGGTACATGGAATTCCTGACACAGCGCTTCTACGGCAGTACGGGCAT  >  NC_002947/5840397‑5840686
                                                                                                                                                |                                                                                                                                                  
tCACCGGACCGCTGAGATCCCGGTCCAGGTCGGCAAAATGCGCCCCCTCAATATGCCCCAGGGCATAGCTGCGCTGGCCATAGTCCACATCCTCGAGGGCAAAGCGGCAGTCGAGGATCACCAACTTGGGCGAGCCCAGACGCTCGGc                                                                                                                                                 >  5:162456/1‑148 (MQ=255)
  aCCGGACCGCTGAGATCCCGGTCCAGGTCGGCAAAATGCGCCCCCTCAATATGCCCCAGGGCATACCTGCGCTGGCCATAGTCCACATCCTCGAGGGCAAAGCGGCAGTCGAGGATCACCAACTTGGCCGAGCCCAGACGCTCGCCCaa                                                                                                                                              >  3:130699/1‑149 (MQ=255)
    cGGACCGCTGAGATCCCGGTCCAGGTCGGCAAAATGCGCCCCCTCAATATGCCCCAGGGCATAGCTGCGCTGGCCATAGTCCACATCCTCGAGGGCAAAGCGGCAGTCGAGGATCACCAACTTGGGCGAGCCCAGACGCTCGGCCAAct                                                                                                                                            >  8:180613/1‑149 (MQ=255)
          gCTGAGATCCCGGTCCAGGTCGGCAAAATGCGCCCCCTCAATATGCCCCAGGGCATAGCTGCGCTGGCCATAGTCCACATCCTCGAGGGCAAAGCGGCAGTCGAGGATCACCAACTTGGGCGATCCCAGACGCTCGGCCAACTGCTGCg                                                                                                                                      >  8:136904/1‑149 (MQ=255)
           cTGAGATCCCGGTCCAGGTCGGCAAAATGCGCCCCCTCAATATGCCCCAGGGCATAGCTGCGCTGGCCATAGTCCACATCCTCGAGGGCAAAGCGGCAGTCGAGGATCACCAACTTGGGCGAGCCCAGACGCTCGGCCAACTGCTGCgg                                                                                                                                     >  3:18118/1‑149 (MQ=255)
            tGAGATCCCGGTCCAGGTCGGCAAAATGCGCCCCCTCAATATGCCCCAGGGCATAGCTGCGCTGGCCATAGTCCACATCCTCGAGGGCAAAGCGGCAGTCGAGGATCACCAACTTGGGCGAGCCCAGACGCTCGGCCAACTGCTGCggg                                                                                                                                    <  4:103562/149‑1 (MQ=255)
               gaTCCCGGTCCAGGTCGGCAAAATGCGCCCCCTCAATATGCCCCAGGGCATAGCTGCGCTGGCCATAGTCCACATCCTCGAGGGCAAAGCGGCAGTCGAGGATCACCAACTTGGGCGAGCCCAGACGCTCGGCCAACTGCTGCGGGGTg                                                                                                                                 <  6:162456/149‑1 (MQ=255)
                 tCCCGGTCCAGGTCGGCAAAATGCGCCCCCTCAATATGCCCCAGGGCATAGCTGCGCTGGCCATAGTCCACATCCTCGAGGGCAAAGCGGCAGTCGAGGATCACCAACTTGGGCGAGCCCAGACGCTCGGCCAACTGCTGCGGGGTGAt                                                                                                                               <  4:204348/149‑1 (MQ=255)
                  cccGGTCCAGGTCGGCAAAATGCGCCCCCTCAATATGCCCCAGGGCATAGCTGCGCTGGCCATAGTCCACATCCTCGAGGGCAAAGCGGCAGTCGAGTATCACCAACTTGGGCGAGCCCAGACGCTCGGCCAACTGCTGCGGGGTGAt                                                                                                                               >  6:301663/1‑148 (MQ=255)
                               ggCAAAATGCGCCCCCTCAATATGCCCCAGGGCATAGCTGCGCTGGCCATAGTCCACATCCTCGAGGGCAAAGCGGCAGTCGAGGATCACCAACTTGGTCGAGCCCAGACGCTCGGCCAACTGCTGCGGGGTGATCAATTGCGCAAgg                                                                                                                  <  1:76299/148‑1 (MQ=255)
                                 cAAAATGCGCCCCCTCAATATGCCCCAGGGCATAGCTGCGCTGGCCATAGTCCACATCCTCGAGGGCAAAGCGGCAGTCGAGGATCACCAACTTGGGCGAGCCCAGACGCTCGGCCAACTGCTGCGGGGTGAt                                                                                                                               >  3:185505/1‑133 (MQ=255)
                                 cAAAATGCGCCCCCTCAATATGCCCCAGGGCATAGCTGCGCTGGCCATAGTCCACATCCTCGAGGGCAAAGCGGCAGTCGAGGATCACCAACTTGGGCGAGCCCAGACGCTCGGCCAACTGCTGCGGGGTGAt                                                                                                                               <  4:185505/133‑1 (MQ=255)
                                   aaaTGCGCCCCCTCAATATGCCCCAGGGCATAGCTGCGCTGGCCATAGTCCACATCCTCGAGGGCAAAGCGGCAGTCGAGGATCACCAACTTGGGCGAGCCCAGACGCTCGGCCAACTGCTGCGGGGTGATCAATTGCGCAAGGGGCAt                                                                                                             <  4:18118/149‑1 (MQ=255)
                                             ccTCAATATGCCCCAGGGCATAGCTGCGCTGGCCATAGTCCACATCCTCGAGGGCAAAGCGGCAGTCGAGGATCACCAACTTGGGCGAGCCCAGACGCTCGGCCAACTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTc‑c                                                                                                  <  2:234512/149‑1 (MQ=255)
                                                agctatGCCCCAGGGCATAGCTGCGCTGGCCATAGTCCACATCCTCGAGGGCAAAGCGGCAGTCGAGGATCACCAACTTGGGCGAGCCCAGACGCTCGGCCAACTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTC‑CTGa                                                                                               >  2:38360/4‑149 (MQ=255)
                                                        cccAGGGCATAGCTGCGCTGGCCATAGTCCACATCCTCGAGGGCAAAGCGGCAGTCGAGGATCACCAACTTGGGCGAGCCCAGACGCTCGGTCAACTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTC‑CTGATCGATgg                                                                                        >  6:34551/1‑148 (MQ=255)
                                                        cccAGGGCATAGCTGCGCTGGCCATAGTCCACATCCTCGAGGGCAAAGCGGCAGTCGAGGATCACCAACTTGGGCGAGCCCAGACGCTCGGCCAACTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTC‑CTGATCGATgg                                                                                        >  7:17956/1‑148 (MQ=255)
                                                        cccAGGGCATAGCTGCGCTGGCCATAGTCCACATCCTCGAGGGCAAAGCGGCAGTCGAGGATCACCAACTTGGGCGAGCCCAGACGCTCGGCCAACTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTC‑CTGATCGATgg                                                                                        >  1:274089/1‑148 (MQ=255)
                                                             ggCATAGCTGCGCTGGCCATAGTCCACATCCTCGAGGGCAAAGCGGCAGTCGAGGATCACCAACTTGGGCGAGCCCAGACGCTCGGCCAACTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTC‑CTGATCGATGGcac                                                                                     <  6:152657/146‑1 (MQ=255)
                                                             ggCATAGCTGCGCTGGCCATAGTCCACATCCTCGAGGGCAAAGCGGCAGTCGAGGATCACCAACTTGGGCGAGCCCAGACGCTCGGCCAACTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTC‑CTGATCGATGGcac                                                                                     >  5:152657/1‑146 (MQ=255)
                                                                    cTGCGCTGGCCATAGTCCACATCCTCGAGGGCAAAGCGGCAGTCGAGGATCACCAACTTGGGCGAGCCCAGACGCTCGGCCAACTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTC‑CTGATCGATGGCACAAGAGCCCTg                                                                           <  5:301663/149‑1 (MQ=255)
                                                                       ggcTGGGCATAGTCCAGATCCTCGAGGGCAATGCGGCAGTCGAGGATCACCAACTTGGGCGAGCCCAGAGGCTCGGCGAAGTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTC‑CTGATCGATGGCACAAGAGCCCTgg                                                                          <  4:130699/146‑1 (MQ=255)
                                                                         cTGGCCATAGTCCACATCCTCGAGGGCAAAGCGGCAGTCGAGGATCACCAACTTGGGCGAGCCCAGACGCTCGGCCAACTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTC‑CTGATCGATGGCACAAGAGCCCTGGGCTa                                                                      <  7:180613/149‑1 (MQ=255)
                                                                          tGGCCATAGTCCACATCCTCGAGGGCAAAGCGGCAGTCGAGGATCACCAACTTGGGCGAGCCCAGACGCTGGGCCAACTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTC‑CTGATCGATGGCACAAGAGCCCTGGGc                                                                        >  2:250976/1‑146 (MQ=255)
                                                                          tGGCCATAGTCCACATCCTCGAGGGCAAAGCGGCAGTCGAGGATCACCAACTTGGGCGAGCCCAGACGCTGGGCCAACTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTC‑CTGATCGATGGCACAAGAGCCCTGGGc                                                                        >  2:247390/1‑146 (MQ=255)
                                                                             ccATAGTCCACATCCTCGAGGGCAAAGCGGCAGTCGAGGATCACCAACTTGGGCGAGCCCAGACGCTCGGCCAACTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTC‑CTGATCGATGGCACAAGAGCCCTGGGCTAGAGc                                                                  <  3:97894/149‑1 (MQ=255)
                                                                              cATAGTCCACATCCTCGAGGGCAAAGCGGCAGTCGAGGATCACCAACTTGGGCGAGCCCAGACGCTCGGCCAACTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTC‑CTGATCGATGGCACAAGAGCCCTGGGCTAGAGCt                                                                 >  7:4683/1‑149 (MQ=255)
                                                                                  gTCCACATCCTCGAGGGCAAAGCGGCAGTCGAGGATCACCAACTTGGGCGAGCCCAGACGCTCGGCCAACTGCTGCGGGGTGATCAATTGCGCAAgg                                                                                                                  <  6:252660/97‑1 (MQ=255)
                                                                                  gTCCACATCCTCGAGGGCAAAGCGGCAGTCGAGGATCACCAACTTGGGCGAGCCCAGACGCTCGGCCAACTGCTGCGGGGTGATCAATTGCGCAAgg                                                                                                                  >  5:252660/1‑97 (MQ=255)
                                                                                  gTCCACATCCTCGAGGGCAAAGCGGCAGTCGAGGATCACCAACTTGGGCGAGCCCAGACGCTCGGCCAACTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTC‑CTGATCGATGGCACAAGAGCCCTGGGCTAGAGCTCCAg                                                             >  7:12638/1‑149 (MQ=255)
                                                                                   tCCACATCCTCGAGGGCAAAGCGGCAGTCGAGGATCACCAACTTGGGCGAGCCCAGACGCTCGGCCAACTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTC‑CTGATCGATg                                                                                         <  2:147013/120‑1 (MQ=255)
                                                                                   tCCACATCCTCGAGGGCAAAGCGGCAGTCGAGGATCACCAACTTGGGCGAGCCCAGACGCTCGGCCAACTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTC‑CTGATCGATg                                                                                         >  1:147013/1‑120 (MQ=255)
                                                                                       caTCCTCGAGGGCAAAGCGGCAGTCGAGGATCACCAACTTGGGCGAGCCCAGACGCTCGGCCAACTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTC‑CTGATCGATGGCACAAGAGCCCTGGGCTAGAGCTCCAGGGCCt                                                        >  5:59599/1‑149 (MQ=255)
                                                                                            tCGAGGGCAAAGCGGCAGTCGAGGATCACCAACTTGGGCGAGCCCAGACGCTCGGCCAACTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTCTCTGATCGATGGCACAAGAGCCCTGGGCTAGAGCTCCAGGGCCTg                                                       <  1:218751/146‑1 (MQ=255)
                                                                                            tCGAGGGCAAAGCGGCAGTCGAGGATCACCAACTTGGGCGAGCCCAGACGCTCGGCCAACTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTC‑CTGATCGATGGCACAAGAGCCCTGGGCTAGAGCTCCAGGGCCTgg                                                      >  5:269867/1‑146 (MQ=255)
                                                                                            tCGAGGGCAAAGCGGCAGTCGAGGATCACCAACTTGGGCGAGCCCAGACGCTCGGCCAACTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTC‑CTGATCGATGGCACAAGAGCCCTGGGCTAGAGCTCCAGGGCCTgg                                                      <  6:269867/146‑1 (MQ=255)
                                                                                            tCGAGGGCAAAGCGGCAGTCGAGGATCACCAACTTGGGCGAGCCCAGACGCTCGGCCAACTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTC‑CTGATCGATGGCACAAGAGCCCTGGGCTAGAGCTCCAGGGCCTg                                                       >  2:218751/1‑145 (MQ=255)
                                                                                                  gCAAAGCGGCAGTCGAGGATCACCAACTTGGGCGAGCCCAGACGCTCGGTCAACTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTC‑CTGATCGATGGCACAAGAGCCCTGGGCTAGAGCTCCAGGGCCTGGTTGAACGGt                                             <  5:34551/149‑1 (MQ=255)
                                                                                                      aGCGGCAGTCGAGGATCACCAACTTGGGCGAGCCCAGACGCTCGGCCAACTGCTGCGGGGTGATAAATTGCGCAAGGGCCATGACAATCTC‑CTGATCGATGGCACAAGAGCCCTGCGCTAGAGCTCCAGGGCCTGGTTGAACGGTACAt                                         >  5:301327/1‑149 (MQ=255)
                                                                                                               cGAGGATCACCAACTTGGGCGAGCCCAGACGCTCGGCCAACTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTC‑CTGATCGATGGCACAAGAGCCCTGGGCTAGAGCTCCAGGGCCTGGTTGAACGGTACATGGAATTCCt                                <  8:12638/149‑1 (MQ=255)
                                                                                                                 aGGATCACCAACTTGGGCGAGCCCAGACGCTCGGCCAACTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTC‑CTGATCGATGGCACAAGAGCCCTGGGCTAGAGCTCCAGGGCCTGGTTGAACGGTACATGGAATTCCTGa                              <  2:274089/149‑1 (MQ=255)
                                                                                                                  ggATCACCAAGTTGGGCTAGCACAGACGCTCGGCCAACTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTC‑CTGATCGATGGCACAAGAGCCCTGGGCTAGAGCTCCAGGGCCTGGTTGAACGGTACATGGAATTCCTGa                              <  6:301327/148‑1 (MQ=255)
                                                                                                                  ggATCACCAACTTGGGCGAGCCCAGACGCTCGGCCAACTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTC‑CTGATCGATGGCACAAGAGCCCTGGGCTAGAGCTCCAGGGCCTGGTTGAACGGTACATGGAATTCCTGa                              >  8:53350/1‑148 (MQ=255)
                                                                                                                        ccAACTTGGGCGAGCCCAGACGCTCGGCCAACTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTC‑CTGATCGATGGCACAAGAGCCCTGGGCTAGAGCTCCAGGGCCTGGTTGAACGGTACATGGAATTCCTGACACAgcg                       >  2:97974/1‑149 (MQ=255)
                                                                                                                         cAACTTGGGCGAGCCCAGACGCTCGGCCAACTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTC‑CTGATCGATGGCACAAGAGCCCTGGGCTAGAGCTCCAGGGCCTGGTTGAACGGTACATGGAATTCCTGACACAgcgc                      >  4:98571/1‑149 (MQ=255)
                                                                                                                                ggCGAGCCCAGACGCTCGGCCAACTGCTGCGGGTTGATCAATTGCGCAAGGGGCATGACAATCTC‑CTGATCGATGGCACAAGAGCACTGGGCTAGAGCTCCAGGGCCTGGTTGAACGGTACATGGAATTCCTGACACAGCGCTTCTAc                <  7:136904/148‑1 (MQ=255)
                                                                                                                                ggCGAGCCCAGACGCTCGGCCAACTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTC‑CTGATCGATGGCACAAGAGCCCTGGGCTAGAGCTCCAGGGCCTGGTTGAACGGTACATGGAATTCCTGACACAGCGCTTCTACg               <  3:98571/149‑1 (MQ=255)
                                                                                                                                 gCGAGCCCAGACGCTGGGCCAACTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTC‑CTGATCGATGGCACAAGAGCCCTGGGCTAGAGCTCCAGGGCCTGGTTGAACGGTACATGGAATTCCTGACACAGCGCTTCTACgg              <  1:250976/149‑1 (MQ=255)
                                                                                                                                 gCGAGCCCAGACGCTGGGCCAACTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTC‑CTGATCGATGGCACAAGAGCCCTGGGCTAGAGCTCCAGGGCCTGGTTGAACGGTACATGGAATTCCTGACACAGCGCTTCTACgg              <  1:247390/149‑1 (MQ=255)
                                                                                                                                   gAGCCCAGACGCTCGGCCAACTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTC‑CTGATCGATGGCACAAGAGCCCTGGGCTAGAGCTCCAGGGCCTGGTTGAACGGTACATGGAATTCCTGACACAGCGCTTCTACgg              >  4:269855/1‑147 (MQ=255)
                                                                                                                                        cAGACGCTCGGCCAACTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTC‑CTGATCGATGGCACAAGAGCCCTGGGCTAGAGCTcc                                                               <  4:215332/93‑1 (MQ=255)
                                                                                                                                        cAGACGCTCGGCCAACTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTC‑CTGATCGATGGCACAAGAGCCCTGGGCTAGAGCTcc                                                               >  3:215332/1‑93 (MQ=255)
                                                                                                                                         aGACGCTCGGCCAACTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTC‑CTGATCGATGGCACAAGAGCCCTGGGCTAGAGCTCCAGGGCCTGGTTGAACGGTACATGGAATTCCTGACACAGCGCTTCTACGGCAGTACgg      <  7:53350/149‑1 (MQ=255)
                                                                                                                                          gACGCTCGGCCAACTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTC‑CTGATCGATGGCACAAGAGCCCTGGGCTAGAGCTCCAGGGCCTGGTTGAACGGTACATGGAATTCCTGACACAGCGCTTCTACGGCAGTACggg     <  1:38360/149‑1 (MQ=255)
                                                                                                                                             gCTCGGCCAACTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTC‑CTGATCGATGGCACAAGAGCCCTGGGCTAGAGCTCCAGGGCCTGGTTGAACGGTACATGGAATTCCTGACACAGCGCTTCTACGGCAGTACGGGCAt  >  4:291839/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
TCACCGGACCGCTGAGATCCCGGTCCAGGTCGGCAAAATGCGCCCCCTCAATATGCCCCAGGGCATAGCTGCGCTGGCCATAGTCCACATCCTCGAGGGCAAAGCGGCAGTCGAGGATCACCAACTTGGGCGAGCCCAGACGCTCGGCCAACTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTC‑CTGATCGATGGCACAAGAGCCCTGGGCTAGAGCTCCAGGGCCTGGTTGAACGGTACATGGAATTCCTGACACAGCGCTTCTACGGCAGTACGGGCAT  >  NC_002947/5840397‑5840686

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑