Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,128,949 C→G 8.0% L66V (CTG→GTG)  dmoA‑I → dimethyl‑sulfide monooxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,128,9490CG 8.0% 154.6 / 3.6 50L66V (CTG→GTG) dmoA‑Idimethyl‑sulfide monooxygenase
Reads supporting (aligned to +/- strand):  ref base C (23/23);  new base G (2/2);  total (25/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.86e-01

CATCAGTCGCCCGGACTCTGGAGGCACCCCAAGGATCGCGCCTGGCAATACAAGGATCTGGAATACTGGACGGACCTGGCCAGGTTGCTTGAGCGCGGTAAGTTTGACGGCATCTTCATCGCCGACGTGATCGGCATCTACGACGTGCTGGGCGGCAATGGCGATGCGGCCATCCGTCAGGCGACCCAGGTTCCGGTCAATGATCCGCTGGCGCTGATCACCCCGATGGCGCTGGTTACCGAGCACCTCGGTTTTGGCCTCACCGCCTCGCTGACTTTCGAGCACCCGTAT  >  NC_002947/3128802‑3129092
                                                                                                                                                   |                                                                                                                                               
cATCAGTCGCCCGGACTCTGGAGGCACCCCAAGGATCGCGCCTGGCAATACAAGGATCTGGAATACTGGACGGACCTGGCCAGGTTGCTTGAGCGCGGTAAGTTTGACGGCATCTTCATCGCCGACGTGATCGGCATCTACGACGTGCt                                                                                                                                                <  1:165178/149‑1 (MQ=255)
cATCAGTCGCCCGGACTCTGGAGGCACCCCAAGGATCGCGCCTGGCAATACAAGGATCTGGAATACTGGAAGGACCTGGCCAGGTTGCATGTGCGAGGTAAGTTTGACGGCATCTTCATCGCGGACGTGATAGGCAAGTACGACGAGCt                                                                                                                                                >  7:283596/1‑149 (MQ=255)
        gCCCGGACTCTGGAGGCACCCCAAGGATCGCGCCTGGCAATACAAGGATCTGGAATACTGGACGGACCTGGCCAGGTTGCTTGAGCGCGGTAAGTTTGACGGCATCTTCATCGCCGACGTGATCGGCATCTACGACGTGCTGggcggc                                                                                                                                         <  3:114398/148‑1 (MQ=255)
                  tGGAGGCACCCCAAGGATCGCGCCTGGCAATACAAGGATCTGGAATACTGGACGGACCTGGCCAGGTTGCTTGAGCGCGGTAAGTTTGACGGCATCTTCATCGCCGACGTGATCGGCATCTACGACGTGCTGGGCGGCAATgg                                                                                                                                    <  8:286791/143‑1 (MQ=255)
                  tGGAGGCACCCCAAGGATCGCGCCTGGCAATACAAGGATCTGGAATACTGGACGGACCTGGCCAGGTTGCTTGAGCGCGGTAAGTTTGACGGCATCTTCATCGCCGACGTGATCGGCATCTACGACGTGCTGGGCGGCAATgg                                                                                                                                    >  7:286791/1‑143 (MQ=255)
                    gAGGCACCGCAAGGATCGCGCCTGGCAATACAAGGATCTGGAATACTGGACGGACCTGGCCAGGTTGCTTGAGCGCGGTAAGTTTGACGGCATCTTCATCGCCGACGTGATCGGCATCTACGACGTGCTGGGCGGCAATGGCGATGCg                                                                                                                             <  6:307982/148‑1 (MQ=255)
                         aCCCCAAGGATCGCGCCTGGCAATACAAGGATCTGGAATACTGGACGGACCTGGCCAGGTTGCTTGAGCGCGGTAAGTTTGACGGCATCTTCATCGCCGACGTGATCGGCATCTACGACGTGCTGGGCGGCAATGGCGATGCGGCCATc                                                                                                                       >  5:11692/1‑149 (MQ=255)
                          ccccAACGATCGCGCCTGGCAATACAAGGATCTGGAATACTGGACTGACCTGGCCAGGTTGCTTGAGCGCGGTAAGTTTGACGGCATCTTCATCGCCGACGTGATCGGCATCTACGACGTGCTGGGCGGCAATGGCGATGCGGCCATcc                                                                                                                      <  6:154903/149‑1 (MQ=255)
                           cccAAGGATCGCGCCTGGCAATACAAGGATCTGGAATACTGGACGGACCTGGCCAGGTTGCTTGAGCGCGGTAAGTTTGACGGCATCTTCATCGCCGACGTGATCGGCATCTACGACGTGCTGGGCGGCAATGGCGATGCGGCCATCCg                                                                                                                     >  1:273368/1‑149 (MQ=255)
                            ccAAGGATCGCGCCTGGCAATACAAGGATCTGGAATACTGGACGGACCTGGCCAGGTTGCTTGAGCGCGGTAAGTTTGACGGCATCTTCATCGCCGACGTGATCGGCATCTACGACGTGCTGGGCGGCAATGGCGATGCg                                                                                                                             <  1:131996/140‑1 (MQ=255)
                            ccAAGGATCGCGCCTGGCAATACAAGGATCTGGAATACTGGACGGACCTGGCCAGGTTGCTTGAGCGCGGTAAGTTTGACGGCATCTTCATCGCCGACGTGATCGGCATCTACGACGTGCTGGGCGGCAATGGCGATGCg                                                                                                                             >  2:131996/1‑140 (MQ=255)
                             cAAGGATCGCGCCTGGCAATACAAGGATGTGGAATACTGGACGGACCTGGCCAGGTTGCTTGAGCGCGGTAAGTTTGACGGCATCTTCATCGCCGACGTGATCGGCATCTACGACGTGCTGGGCGGCAATGGCGATGCGGCCATCCGt                                                                                                                    >  1:177491/1‑148 (MQ=255)
                              aaGGATCGCGCCTGGCAATACAAGGATCTGGAATACTGGACGGACCTGGCCGGGTTGCTTGAGCGCGGTAAGTTTGACGGCATCTTCATCGCCGACGTGATCGGCATCTACGACGTGCTGGGCGGCAATGGCGATGCGGCCATCCGTCa                                                                                                                  <  5:286443/149‑1 (MQ=255)
                                ggATCGCGCCTGGCAATACAAGGATCTGGAATACTGGACGGACCTGGCCAGGTTGCTTGAGCGCGGTAAGTTTGACGGCATCTTCATCGCCGACGTGATCGGCATCTACGACGTGCTGGGCGGCAATGGCGATGCGGCCATCCGTCAgg                                                                                                                <  1:103578/149‑1 (MQ=255)
                                       gcCTGGCAATACAAGGATCTGGAATACTGGACGGACCTGGCCAGGTTGCTTGAGCGCGGTAAGTTTGACGGCATCTTCATCGCCGACGTGATCGGCATCTACGACGTGCTGGGCGGCAATGGCGATGCGGCCATCCGTCAGGCGACCCa                                                                                                         <  8:41616/149‑1 (MQ=255)
                                              aaTACAAGGATCTGGAATACTGGACGGACCTGGCCAGGTTGCTTGAGCGCGGTAAGTTTGACGGCATCTTCATCGCCGACGTGATCGGCATCTACGACGTGCTGGGCGGCAATGGCGATGCGGCCATCCGTCAGGCGACCCAGGTTCCg                                                                                                  >  3:140747/1‑149 (MQ=255)
                                              aaTACAAGGATCTGGAATACTGGACGGACCTGGCCAGGTTGCTTGAGCGCGGTAAGTTTGACGGCATCTTCATCGCCGACGTGATCGGCATCTACGACGTGCTGGGCGGCAATGGCGATGCGGCCATCCGTCAGGCGACCCAGGTTCCg                                                                                                  >  2:248269/1‑149 (MQ=255)
                                               aTACAAGGATCTGGAATACTGGACGGACCTGGCCAGGTTGCTTGAGCGCGGTAAGTTTGACGGCATCTTCATCGCCGACGTGATCGGCATCTACGACGTGCTGGGCGGCAATGGCGATGCGGCCATCCGTCAGGCGACCCAGGTTCCgg                                                                                                 <  4:310440/149‑1 (MQ=255)
                                                      gATCTGGAATACTGGACGGACCTGGCCAGGTTGCTTGAGCGCGGTAAGTTTGACGGCATCTTCAGCGCCGACGTGATCGGCATCTACGACGTGCTGGGCGGCAATGGCGATGCGGCCATCCGTCAGGCGACCCAGGTTCCGGTCAATGa                                                                                          <  8:12122/149‑1 (MQ=255)
                                                              aTACTGGACGGACCTGGCCAGGTTGCTTGAGCGCGGTAAGTTTGACGGCATCTTCATCGCCGACGTGATCGGCATCTACGACGTGCTGGGCGGCAATGGCGATGCGGCCATCCGTCAGGCGACCCAGGTTCCGGTCAATGATCCGCTgg                                                                                  <  1:116098/149‑1 (MQ=255)
                                                                  tGGACGGACCTGGCCAGGTTGCTTGAGCGCGGTAAGTTTGACGGCATCTTCATCGCCGACGTGATCGGCATCTACGACGTGCTGGGCGGCAATGGCGATGCGGCCATCCGTCAGGCGACCCAGGTTCCGGTCAATGATCCGCTGGCGCt                                                                              >  4:138013/1‑149 (MQ=255)
                                                                    gacggacCTGGCCAGGTTGCTTGAGCGCGGTAAGTTTGACGGCATCTTCATCGCCGACGTGATCGGCATCTACGACGTGCTGGGCGGCAATGGCGATGCGGCCATCCGTCAGGCGACCCAGGTTCCGGTCAATGATCCGCTGGCGCTGa                                                                            >  7:116823/1‑149 (MQ=255)
                                                                         acCTGGCCAGGTTGCTTGAGCGCGGTAAGTTTGACGGCATCTTCATCGCCGACGTGATCGGCATCTACGACGTGCTGGGCGGCAATGGCGATGCGGCCATCCGTCAGGCGACCCAGGTTCCGGTCAATGATCCGCTGGCGCTGATCAcc                                                                       >  7:195738/1‑149 (MQ=255)
                                                                         acCTGGCCAGGTTGCTTGAGCGCGGTAAGTTTGACGGCATCTTCATCGCCGACGTGATCGGCATCTACGACGTGCTGGGCGGCAATGGCGATGCGGCCATCCGTCAGGCGACCCAGGTTCCGGTCAATGATCCGCTGGCGCTGATCAcc                                                                       <  6:11692/149‑1 (MQ=255)
                                                                            tGGCCAGGTTGCTTGAGCGCGGTAAGTTTGACGGCATCTTCATCGCCGACGTGATCGGCATCTACGACGTGCTGGGCGGCAATGGCGATGCGGCCATCCGTCAGGCGACCCAGGTTCCGGTCAATGATCCGCTGGCGCTGATCa                                                                         >  6:220158/1‑144 (MQ=255)
                                                                            tGGCCAGGTTGCTTGAGCGCGGTAAGTTTGACGGCATCTTCATCGCCGACGTGATCGGCATCTACGACGTGCTGGGCGGCAATGGCGATGCGGCCATCCGTCAGGCGACCCAGGTTCCGGTCAATGATCCGCTGGCGCTGATCa                                                                         <  5:220158/144‑1 (MQ=255)
                                                                                      gCTTGAGCGCGGTAAGTTTGACGGCATCTTCATCGCCGACGTGATCGGCATCTACGACGTGCTGGGCGGCAATGGCGATGCGGCCATCCGTCAGGCGACCCAGGTTCCGGTCAATGATCCGCTGGCGCTGATCACCCCGATGGCGCTg                                                           <  3:206487/148‑1 (MQ=255)
                                                                                        ttGAGCGCGGTAAGTTTGACGGCATCTTCATCGCCGACGTGATCGGCATCTACGACGTGGTGGGCGGCAATGGCGATGCGGCCATCCGTCAGGCGACCCAGGTTCCGGTCAATGATCCGCTGGCGCTGATCACCCCGATGGCGCTGGt                                                         >  5:200297/1‑148 (MQ=255)
                                                                                        ttGAGCGCGGTAAGTTTGACGGCATCTTCATCGCCGACGTGATCGGCATCTACGACGTGGTGGGCGGCAATGGCGATGCGGCCATCCGTCAGGCGACCCAGGTTCCGGTCAATGATCCGCTGGCGCTGATCACCCCGATGGCGCTGGt                                                         >  5:204933/1‑148 (MQ=255)
                                                                                        ttGAGCGCGGTAAGTTTGACGGCATCTTCATCGCCGACGTGATCGGCATCTACGACGTGCTGGGCGGCAATGGCGATGCGGCCATCCGTCAGGCGACCCAGGTTCCGGTCAATGATCCGCTGGCGCTGATCACCCCGATGGCGCTGGt                                                         >  6:169628/1‑148 (MQ=255)
                                                                                            gcgcGGTAAGTTTGACGGCATCTTCATCGCCGACGTGATCGGCATCTACGACGTGCTGGGCGGCAATGGCGATGCGGCCATCCGTCAGGCGACCCAGGTTCCGGTCAATGATCCGCTGGCGCTGATCACCCCGATGGCGCTGGTTAcc                                                     <  5:169628/148‑1 (MQ=255)
                                                                                               cgGTAAGTTTGACGGCATCTTCATCGCCGACGTGATCGGCATCTACGACGTGCTGGGCGGCAATGGCGATGCGGCCATCCGTCAGGCGACCCAGGTTCCGGTCAATGATCCGCTGGCGCTGATCACCCCGATGGCGCTGGTTACCGAg                                                  >  7:90773/1‑148 (MQ=255)
                                                                                                ggTAAGTTTGACGGCATCTTCATCGCCGACGTGATCGGCATCTACGACGTGCTGGGCGGCAATGGCGATGCGGCCATCCGTCAGGCGACCCAGGTTCCGGTCAATGATCCGCTGGCGCTGATCACCCCGATGGCGCTGGTTACCGAGCa                                                <  8:116823/149‑1 (MQ=255)
                                                                                                     gTTTGACGGCATCTTCATCGCCGACGTGATCGGCATCTACGACGTGCTGGGCGGCAATGGCGATGCGGCCATCCGTCAGGCGACCCAGGTTCCGGTCAATGATCCGCTGGCGCTGATCACCCCGATGGCGCTGGTTACCGAGCACCTCg                                           >  5:4451/1‑149 (MQ=255)
                                                                                                            ggCATCTTCATCGCCGACGTGATCGGCATCTACGACGTGGTGGGCGGCAATGGCGATGCGGCCATCCGTCAGGCGACCCAGGTTCCGGTCAATGATCCGCTGGCGCTGATCACCCCGATGGCGCTGGTTACCGAGCACCTCGGTTTTgg                                    <  6:204933/149‑1 (MQ=255)
                                                                                                            ggCATCTTCATCGCCGACGTGATCGGCATCTACGACGTGGTGGGCGGCAATGGCGATGCGGCCATCCGTCAGGCGACCCAGGTTCCGGTCAATGATCCGCTGGCGCTGATCACCCCGATGGCGCTGGTTACCGAGCACCTCGGTTTTgg                                    <  6:200297/149‑1 (MQ=255)
                                                                                                             gCATCTTCATCGCCGACGTGATCGGCATCTACGACGTGCTGGGCGGCAATGGCGATGCGGCCATCCGTCAGGCGACCCAGGTTCCGGTCAATGATCCGCTGGCGCTGATCACCCCGATGGCGCTGGTTACCGAGCACCTCGGTTTTGGc                                   <  1:248269/149‑1 (MQ=255)
                                                                                                                tCTTCATCGCCGACGTGATCGGCATCTACGACGTGCTGGGCGGCAATGGCGATGCGGCCATCCGTCAGGCGACCCAGGTTCCGGTCAATGATCCGCTGGCGCTGATCACCCCGATGGCGCTGGTTACCGAGCACCTCGGTTTTGGCCt                                 >  5:116690/1‑148 (MQ=255)
                                                                                                                   tCATCGCCGACGTGATCGGCATCTACGACGTGCTGGGCGGCAATGGCGATGCGGCCATCCGTCAGGCGACCCAGGTTCCGGTCAATGATCCGCTGGCGCTGATCACCCCGATGGCGCTGGTTACCGAGCACCTCGGTTTTGGCCTCAcc                             >  8:140272/1‑149 (MQ=255)
                                                                                                                         ccGACGTGATCGGCATCTACGACGTGCTGGGCGGCAATGGCGATGCGGCCATCCGTCAGGCGACCCAGGTTCCGGTCAATGATCCGCTGGCGCTGATCACCCCGATGGCGCTGGTTACCGAGCACCTCGGTTTTGGCCTCACCGCCTCg                       >  1:26592/1‑149 (MQ=255)
                                                                                                                          cGACGTGATCGGCATCTACGACGTGCTGGGCGGCAATGGCGATGCGGCCATCCGTCAGGCGACCCAGGTTCCGGTCAATGATCCGCTGGCGCTGATCACCCCGATGGCGCTGGTTACCGAGCAc                                               >  4:61129/1‑124 (MQ=255)
                                                                                                                          cGACGTGATCGGCATCTACGACGTGCTGGGCGGCAATGGCGATGCGGCCATCCGTCAGGCGACCCAGGTTCCGGTCAATGATCCGCTGGCGCTGATCACCCCGATGGCGCTGGTTACCGAGCAc                                               <  3:61129/124‑1 (MQ=255)
                                                                                                                            aCGTGATCGGCATCTACGACGTGCTGGGCGGCAATGGCGATGCGGCCATCCGTCAGGCGACCCAGGTTCCGGTCAATGATCCGCTGGCGCTGATCACCCCGATGGCGATGGTTACCGAGCACCTCGGTTTTGGCCTCACCGCCTCGCt                     >  8:225180/1‑148 (MQ=255)
                                                                                                                                  tCGGCATCTACGACGTGCTGGGCGGCAATGGCGATGCGGCCATCCGTCAGGCGACCCAGGTTCCGGTCAATGATCCGCTGGCGCTGATCACCCCGATGGCGCTGGTTACCGAGCACCTCGGTTTTGGCCTCACCGCCTCGCTGACTTTc              <  8:90773/149‑1 (MQ=255)
                                                                                                                                  tCGGCATCTACGACGTGCTGGGCGGCAATGGCGATGCGGCCATCCGTCAGGCGACCCAGGTTCCGGTCAATGATCCGCTGGCGCTGATCACCCCGATGGCGCTGGTTACCGAGCACCTCGGTTTTGGCCTCACCGCCTCGCTGACTTTc              >  7:149812/1‑149 (MQ=255)
                                                                                                                                    ggCATCTACGACGTGCTGGGCGGCAATGGCGATGCGGCCATCCGTCAGGCGACCCAGGTTCCGGTCAATGATCCGCTGGCGCTGATCACCCCGATGGCGCTGGTTACCGAGCACCTCGGTTTTGGCCTCACCGCCTCGCTGACTTTCg             <  4:140747/148‑1 (MQ=255)
                                                                                                                                     gCATCTACGACGTGCTGGGCGGCAATGGCGATGCGGCCATCCGTCAGGCGACCCAGGTTCCGGTCAATGATCCGCTGGCGCTGATCACCCCGATGGCGCTGGTTACCGAGCAACTCGGTTTTGGCCTCACCGCCTCGCTGACTTTCGAg           >  3:235240/1‑149 (MQ=255)
                                                                                                                                               cgTGCTTGGCGGAAATTGCGATGCGGCCATCCGTCAGGCGACCCAGGTTCCGGTCAATGATCCGCTGGCGCTGATCACCCCGATGGCGCTGGTTACCGAGCACCTCGGTTTTGGCCTCACCGCCTCGCTGACTTTCGAGCACCCGTAt  <  2:177491/148‑1 (MQ=255)
                                                                                                                                               cgTGCTGGGCGGCAATGGCGATGCGGCCATCCGTCAGGCGACCCAGGTTCCGGTCAATGATCCGCTGGCGCTGATCACCCCGATGGCGCTGGTTACCGAGCACCTCGGTTTTGGCCTCACCGCCTCGCTGACTTTCGAGCAc        >  7:50992/1‑142 (MQ=255)
                                                                                                                                               cgTGCTGGGCGGCAATGGCGATGCGGCCATCCGTCAGGCGACCCAGGTTCCGGTCAATGATCCGCTGGCGCTGATCACCCCGATGGCGCTGGTTACCGAGCACCTCGGTTTTGGCCTCACCGCCTCGCTGACTTTCGAGCAc        <  8:50992/142‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
CATCAGTCGCCCGGACTCTGGAGGCACCCCAAGGATCGCGCCTGGCAATACAAGGATCTGGAATACTGGACGGACCTGGCCAGGTTGCTTGAGCGCGGTAAGTTTGACGGCATCTTCATCGCCGACGTGATCGGCATCTACGACGTGCTGGGCGGCAATGGCGATGCGGCCATCCGTCAGGCGACCCAGGTTCCGGTCAATGATCCGCTGGCGCTGATCACCCCGATGGCGCTGGTTACCGAGCACCTCGGTTTTGGCCTCACCGCCTCGCTGACTTTCGAGCACCCGTAT  >  NC_002947/3128802‑3129092

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