Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,139,818 G→A 8.7% E124E (GAG→GAA)  sfnB → SfnB family sulfur acquisition oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,139,8180GA 8.7% 125.3 / 2.7 46E124E (GAG→GAA) sfnBSfnB family sulfur acquisition oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (21/21);  new base A (2/2);  total (23/23)
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.78e-01

GTGTGCCCAAGGCGCACGGCGGCGCCGGGGTATCTTATGTCACCCTCAGCGAGGTGATCGCGACCATCTCCGCCGCCGATTCATCACTCGGTCAGTTGCCGCAGAACCACTTTGGTGTGCTCAGCAACCTCGGCCTGACCGGCAGCGAGGCGCAGAAGCACCGCTACTACGCCAAGGTGCTGCAAGGCTACCGTTTCGGTAACGCCTTTTCCGAGGCGCGCAGCCAGTACGTCACCACCTTTCAGACGCAGATCCGCTTCGACGGCGACAGCGCCGTGCTCGATG  >  NC_002947/3139670‑3139954
                                                                                                                                                    |                                                                                                                                        
gtgtGCCCAAGGCGCACGGCGGCGCCGGGGTATCTTATGTCACCCTCAGCGAGGTGATCGCGACCATCTCCGCCGCCGATTCATCACTCGGTCAGTTGCCGCAGAACCACTTTGGTGTGCTCAGCAACCTCGGCCTGACCGGCAGCGAg                                                                                                                                          <  7:246623/149‑1 (MQ=255)
     cccAAGGCGCACGGCGGCGCCGGGGTATCTTATGTCACCCTCAGCGAGGTGATCGCGACCATCTCCGCCGCCGATTCATCACTCGGTCAGTTGCCGGAGAACCACTTTGGTGTGCTCAGCAACCTCGGCCTGACCGGCAGCGAGgcg                                                                                                                                       >  8:204537/1‑147 (MQ=255)
      ccAAGGCGCACGGCGGCGCCGGGGTATCTTATGTCACCCTCAGCGAGGTGATCGCGACCATCTCCGCCGCCGAGTCATCACTCGGTCAGTTGCCGCAGAACCACTTTGGTGTGCTCAGCAACCTCGGCCTGACCGGCAGCGAGgcg                                                                                                                                       >  7:234207/1‑146 (MQ=255)
               aCGGCGGCGCCGGGGTATCTTATGTCACCCTCAGCGAGGTGATCGCGACCATCTCCGCCGCCGATTCATCACTCGGTCAGTTGCCGCAGAACCACTTTGGTGTGCTCAGCAACCTCGGCCTGACCGGCAGCGAGGCGCAGAAGCACCGc                                                                                                                           >  8:207140/1‑149 (MQ=255)
                  gcggcgCCGGGGTATCTTATGTCACCCTCAGCGAGGTGATCGCGACCATCTCCGCCGCCGATTCATCACTCGGTCAGTTGCCGCAGAACCACTTTGGTGTGCTCAGCAACCTCGGCCTGACCGGCAGCGAGGCGCAGAAGCACCGctac                                                                                                                        <  2:284814/149‑1 (MQ=255)
                          ggggTATCTTATGTCACCCTCAGCGAGGTGATCGCGACCATCTCCGCCGCCGATTCATCACTCGGTCAGTTGCCGCAGAACCACTTTGGTGTGCTCAGCAACCTCGGCCTGACCGGCAGCGAGGCGCAGAAGCACCGCTACTACGCCaa                                                                                                                <  3:260940/149‑1 (MQ=255)
                           gggTATCTTATGTCACCCTCAGCGAGGTGATCGCGACCATCTCCGCCGCCGATTCATCACTCGGTCAGTTGCCGCAGAACCACTTTGGTGTGCTCAGCAACCTCGGCCTGACCGGCAGCGAGGCGCAGAAGCACCGCTACTACGCCAAg                                                                                                               >  6:104736/1‑149 (MQ=255)
                            ggTATCTTATGTCACCCTCAGCGAGGTGATCGCGACCATCTCCGCCGCCGATTCATCACTCGGTCAGTTGCCGCAGAACCACTTTGGTGTGCTCAGCAACCTCGGCCTGACCGGCAGCGAGGCGCAGAAGCACCGCTACTACGCCAAgg                                                                                                              >  4:100683/1‑149 (MQ=255)
                                 cTTATGTCACCCTCAGCGAGGTGATCGCGACCATCTCCGCCGCCGATTCATCACTCGGTCAGTTGCCGCAGAACCACTTTGGTGTGCTCAGCAACCTCGGCCTGACCGGCAGCGAGGCGCAGAAGCACCGCTACTACGCCAAGGtgctg                                                                                                         >  2:123198/1‑149 (MQ=255)
                                  ttATGTCACCCTCAGCGAGGTGATCGCGACCATCTCCGCCGCCGATTCATCACTCGGTCAGTTGCCGCAGAACCACTTTGGTGTGCTCAGCAACCTCGGCCTGACCGGCAGCGAGGCGCAGAAGCACCGCTACTACGCCAAGGtgctg                                                                                                         >  6:154694/1‑148 (MQ=255)
                                    aTGTCACCCTCAGCGAGGTGATCGCGACCATCTCCGCCGCCGATTCATCACTCGGTCAGTTGCCGCAGAACCACTTTGGTGTGCTCAGCAACCTCGGCCTGACCGGCAGCGAGGCGCAGAAGCTCCGCTACTACGCCAAGGtgctgc                                                                                                        <  7:204537/147‑1 (MQ=255)
                                    aTGTCACCCTCAGCGAGGTGATCGCGACCATCTCCGCCGCCGATTCATCACTCGGTCAGTTGCCGCAGAACCACTTTGGTGTGCTCAGCAACCTCGGCCTGACCGGCAGCGAGGCGCAGAAGCACCGCTACTACGCCAAGGtgctg                                                                                                         <  1:123198/146‑1 (MQ=255)
                                     tGTCACCCTCAGCGAGGTGATCGCGACCATCTCCGCCGCCGATTCATCACTCGGTCAGTTGCCGCAGAACCACTTTGGTGTGCTCAGCAACCTCGGCCTGACCGGCAGCGAGGCGCAGAAGCACCGCTACTACGCCAAGGTGCTGCaa                                                                                                      <  5:60279/148‑1 (MQ=255)
                                      gTCACCCTCAGCGAGGTGATCGCGACCATCTCCGCCGCCGATTCATCACTCGGTCAGTTGCCGCAGAACCACTTTGGTGTGCTCAGCAACCTCGGCCTGACCGGCAGCGAGGCGCAGAAGCACCGCTACTACGCCAAGGTGCTGCAAgg                                                                                                    >  3:37317/1‑149 (MQ=255)
                                       tCACCCTCAGCGAGGTGATCGCGACCATCTCCGCCGCCGATTCATCACTCGGTCAGTTGCCGCAGAACCACTTTGGTGTGCTCAGCAACCTCGGCCTGACCGGCAGCGAGGCGCAGAAGCACCGCTACTACGCCAAGGTGCTGCAAGGc                                                                                                   >  6:185862/1‑149 (MQ=255)
                                           ccTCAGCGAGGTGATCGCGACCATCTCCGCCGCCGATTCATCACTCGGTCAGTTGCCGCAGAACCACTTTGGTGTGCTCAGCAACCTCGGCCTGACCGGCAGCGAGGCGCAGAAGCACCGCTACTACGCCAAGGTGCTGCAAGGCTAcc                                                                                               >  7:43989/1‑149 (MQ=255)
                                           ccTCAGCGAGGTGATCGCGACCATCTCCGCCGCCGATTCATCACTCGGTCAGTTGCCGCAGAACCACTTTGGTGTGCTCAGCAACCTCGGCCTGACCGGCAGCGAGGCGCAGAAGCACCGCTACTACGCCAAGGTGCTGCAAGGCTAcc                                                                                               >  8:289435/1‑149 (MQ=255)
                                           ccTCAGCGAGGTGATCGCGACCATCTCCGCCGCCGATGCATCACTCGGTCAGTTGCCGCAGAACCCGTTTGGTGTGCTCAGCAAACTCGGCCTGACCGGCAGCGCGGCGCAGAAGCACCGCTACTACGCAAAGGTGCTGCTAGGATAcg                                                                                               >  8:82748/1‑148 (MQ=255)
                                                  gAGGTGATCGCGACCATCTCCGCCGCCGATTCATCACTCGGTCAGTTGCCGCAGAACCACTTTGGTGTGCTCAGCAACCTCGGCCTGACCGGCAGCGAGGCGCAGAAGCACCGCTACTACGCCAAGGTGCTGCAAGGCTACCGTTTCgg                                                                                        <  4:298603/149‑1 (MQ=255)
                                                     gTGATCGCGACCATCTCCGCCGCCGATTCATCACTCGGTCAGTTGCCGCAGAACCACTTTGGTGTGCTCAGCAACCTCGGCCTGACCGGCAGCGAGGCGCAGAAGCACCGCTACTACGCCAAGGTGCTGCAAGGCTACCGTTTCGGTaa                                                                                     >  6:189854/1‑149 (MQ=255)
                                                      tGATCGCGACCATCTCCGCCGCCGATTCATCACTCGGTCAGTTGCCGCAGAACCACTTTGGTGTGCTCAGCAACCTCGGCCTGACCGGCAGCGAGGCGCAGAAGCACCGCTACTACGCCAAGGTGCTGCAAGGCTACCGTTTCGGTAAc                                                                                    <  7:289435/149‑1 (MQ=255)
                                                         tCGCGACCATCTCCGCCGCCGATTCATCACTCGGTCAGTTGCCGCAGAACCACTTTGGTGTGCTCAGCAACCTCGGCCTGACCGGCAGCGAGGCGCAGAAGCACCGCTACTACGCCAAGGTGCTGCAAGGCTACCGTTTCGGTAACGcc                                                                                 >  1:183866/1‑149 (MQ=255)
                                                                      cgccgccgATTCATCACTCGGTCAGTTGCCGCAGAACCACTTTCGGGTGCTCAGCAACCTCGGCCTGACCGGCAGCGAGGCGCAGAAGCACCGCTACTACGCCAAGGTGCTGCAAGGCTACCGTTTCGGTAACGCCTTTTCCGAg                                                                        <  8:75852/145‑1 (MQ=255)
                                                                           ccgATTCATCACTCGGTCAGTTGCCGCAGAACCACTTTGGTGTGCTCAGCAACCTCGGCCTGACCGGCAGCGAGGCGCAGAAGCACCGCTACTACGCCAAGGTGCTGCAAGGCTACCGTTTCGGTAACGCCTTTTc                                                                            >  4:173009/1‑136 (MQ=255)
                                                                           ccgATTCATCACTCGGTCAGTTGCCGCAGAACCACTTTGGTGTGCTCAGCAACCTCGGCCTGACCGGCAGCGAGGCGCAGAAGCACCGCTACTACGCCAAGGTGCTGCAAGGCTACCGTTTCGGTAACGCCTTTTc                                                                            <  3:173009/136‑1 (MQ=255)
                                                                           ccgATTCATCACTCGGTCAGTTGCCGCAGAACCACTTTGGTGTGCTCAGCAACCTCGGCCTGACCGGCAGCGAGGCGCAGAAGCACCGCTACTACGCCAAGGTGCTGCAAGGCTACCGTTTCGGTAACGCCTTTTCCGAGGCGCGCAGc                                                               <  5:189854/149‑1 (MQ=255)
                                                                                 catcaCTCGGTCAGTTGCCGCAGAACCACTTTGGTGTGCTCAGCAACCTCGGCCTGACCGGCAGCGAGGCGCAGAAGCACCGCTACTACGCCAAGGTGCTGCAAGGCTACCGTTTCGGTAACGCCTTTTCCGAGGCGCGCAGCCAGTa                                                          <  3:100683/148‑1 (MQ=255)
                                                                                  atcaCTCGGTCAGTTGCCGCAGAACCACTTTGGTGTGCTCAGCAACCTCGGCCTGACCGGCAGCGAGGCGCAGAAGCACCGCTACTACGCCAAGGTGCTGCAAg                                                                                                     <  6:193156/104‑1 (MQ=255)
                                                                                      cTCGGTCAGTTGCCGCAGAACCACTTTGGTGTGCTCAGCAACCTCGGCCTGACCGGCAGCGAGGCGCAGAAGCACCGCTACTACGCCAAGGTGCTGCAAGGCTACCGTTTCGGTAACGCCTTTTCCGAGGCGCGCAGCCAGTACGTcac                                                    >  3:132097/1‑149 (MQ=255)
                                                                                       tCGGTCAGTTGCCGCAGAGCCACTTTGGTGTGCTCAGCAACCTCGGCCTGACCGGCAGCGAGGCGCAGAAGCACCGCTACTACGCCAAGGTGCTGCAAGGCTACCGTTTCGGTAACGCCTTTTCCGAGGCGCGCAg                                                                >  2:80525/1‑136 (MQ=255)
                                                                                       tCGGTCAGTTGCCGCAGAGCCACTTTGGTGTGCTCAGCAACCTCGGCCTGACCGGCAGCGAGGCGCAGAAGCACCGCTACTACGCCAAGGTGCTGCAAGGCTACCGTTTCGGTAACGCCTTTTCCGAGGCGCGCAg                                                                <  1:80525/136‑1 (MQ=255)
                                                                                          gTCAGTTGCCGCAGAACCACTTTGGTGTGCTCAGCAACCTCGGCCTGACCGGCAGCGAAGCGCAGAAGCACCGCTACTACGCCAAGGTGCTGCAAGGCTACCGTTTCGGTAACGCCTTTTCCGAGGCGCGCAGCCAGTACGTCaccacc                                                >  8:158047/1‑149 (MQ=255)
                                                                                          gTCAGTTGCCGCAGAACCACTTTGGTGTGCTCAGCAACCTCGGCCTGACCGGCAGCGAAGCGCAGAAGCACCGCTACTACGCCAAGGTGCTGCAAGGCTACCGTTTCGGTAACGCCTTTTCCGAGGCGCGCAGCCAGTACGTCaccacc                                                >  8:307187/1‑149 (MQ=255)
                                                                                                   cGCAGAACCACTTTGGTGTGCTCAGCAACCTCGGCCTGACCGGCAGCGAGGCGCAGAAGCACCGCTACTACGCCAAGGTGCTGCAAGGCTACCGTTTCGGTAACGCCTTTTCCGAGGCGCGCAGCCAGTACGTCACCACCTTTCAGACg                                       <  5:185862/149‑1 (MQ=255)
                                                                                                    gCAGAACCACTTTGGTGTGCTCAGCAACCTCGGCCTGACCGGCAGCGAGGCGCAGAAGCACCGCTACTACGCCAAGGTGCTGCAAGGCTACCGTTTCGGTAACGCCTTTTCCGAGGCGCGCAGCCAGTACGTCACCACCTTTCAGACGc                                      <  8:43989/149‑1 (MQ=255)
                                                                                                         aCCACTTTGGTGTGCTCAGCAACCTCGGCCTGACCGGCAGCGAGGCGCAGAAGCACCGCTACTACGCCAAGGTGCTGCAAGGCTACCGTTTCGGTAACGCCTTTTCCGAGGCGCGCAGCCAGTACGTCACCACCTTTCAGACGCAGATc                                 <  5:104736/149‑1 (MQ=255)
                                                                                                             cTTTGGTGTGCTCAGCAACCTCGGCCTGACCGGCAGCGAGGCGCAGAAGCACCGCTACTACGCCAAGGTGCTGCAAGGCTACCGTTTCGGTAACGCCTTTTCCGAGGCGCGCAGCCAGTACGTCACCACCTTTCAGACGCAGATCCGCt                             <  4:132097/149‑1 (MQ=255)
                                                                                                                  gtgtGCTCAGCAACCTCGGCCTGACCGGCAGCGAGGCGCAGAAGCACCGCTACTACGCCAAGGTGCTGCAAGGCTACCGTTTCGGTAACGCCTTTTCCGAGGCGCGCAGCCAGTACGTCACCACCTTTCAGACGCAGATCCGCTTCGAc                        >  2:13435/1‑149 (MQ=255)
                                                                                                                     tgCTCAGCAACCTCGGCCTGACCGGCAGCGAAGCGCAGAAGCACCGCTACTACGCCAAGGTGCTGCAAGGCTACCGTTTCGGTAACGCCTTTTCCGAGGCGCGCAGCCAGTACGTCACCACCTTTCAGACGCAGATCCGCTTCGACGGc                     <  7:158047/149‑1 (MQ=255)
                                                                                                                     tgCTCAGCAACCTCGGCCTGACCGGCAGCGAAGCGCAGAAGCACCGCTACTACGCCAAGGTGCTGCAAGGCTACCGTTTCGGTAACGCCTTTTCCGAGGCGCGCAGCCAGTACGTCACCACCTTTCAGACGCAGATCCGCTTCGACGGc                     <  7:307187/149‑1 (MQ=255)
                                                                                                                        tCAGCAACCTCGGCCTGACCGGCAGCGAGGCGCAGAAGCACCGCTACTACGCCAAGGTGCTGCAAGGCTACCGTTTCGGTAACGCCTTTTCCGAGGCGCGCAGCCAGTACGTCACCACCTTTCAGACGCAGATCCGCTTCGACGGCGa                   >  6:292082/1‑148 (MQ=255)
                                                                                                                        tCAGCAACCTCGGCCTGACCGGCAGCGAGGCGCAGAAGCACCGCTACTACGCCAAGGTGCTGCAAGGCTACCGTTTCGGTAACGCCTTTTCCGAGGCGCGCAGCCAGTACGTCACCACCTTTCAGACGCAGATCCGCTTCGACGGCGa                   >  6:286995/1‑148 (MQ=255)
                                                                                                                          aGCAACCTCGGCCTGACCGGCAGCGAGGCGCAGAAGCACCGCTACTACGCCAAGGTGCTGCAAGGCTACCGTTTCGGTAACGCCTTTTCCGAGGCGCGCAGCCAGTACGTCACCACCTTTCAGACGCAGATCCGCTTCGACGGCGACAg                >  6:123901/1‑149 (MQ=255)
                                                                                                                            cAACCTCGGCCTGACCGGCAGCGAGGCGCAGAAGCACCGCTACTACGCCAAGGTGCTGCAAGGCTACCGTTTCGGTAACGCCTTTTCCGAGGCGCGCAGCCAGTACGTCACCACCTTTCAGACGCAGATCCGCTTCGACGGCGACAgcg              <  4:37317/149‑1 (MQ=255)
                                                                                                                              aCCTCGGCCTGACCGGCAGCGAGGCGCAGAAGCACCGCTACTACGCCAAGGTGCTGCAAGGCTACCGTTTCGGTAACGCCTTTTCCGAGGCGCGCAGCCAGTACGTCACCACCTTTCAGACGCAGATCCGCTTCGACGGCGACAGCGcc            <  5:154694/149‑1 (MQ=255)
                                                                                                                                    gCCTGACCGGCAGCGAGGCGCAGAAGCACCGCTACTACGCCAAGGTGCTGCAAGGCTACCGTTTCGGTAACGCCTTTTCCGAGGCGCGCAGCCAGTACGTCACCACCTTTCAGACGCAGATCCGCTTCGACGGCGACAGCGCCGTGCTc      <  2:183866/149‑1 (MQ=255)
                                                                                                                                         aCCGGCAGCGAGGCGCAGAAGCACCGCTACTACGCCAAGGTGCTGCAAGGCTACCGTTTCGGTAACGCCTTTTCCGAGGCGCGCAGCCAGTACGTCACCACCTTTCAGACGCAGATCCGCTTCGACGGCGACAGCGCCGTGCTCGATg  <  1:13435/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                        
GTGTGCCCAAGGCGCACGGCGGCGCCGGGGTATCTTATGTCACCCTCAGCGAGGTGATCGCGACCATCTCCGCCGCCGATTCATCACTCGGTCAGTTGCCGCAGAACCACTTTGGTGTGCTCAGCAACCTCGGCCTGACCGGCAGCGAGGCGCAGAAGCACCGCTACTACGCCAAGGTGCTGCAAGGCTACCGTTTCGGTAACGCCTTTTCCGAGGCGCGCAGCCAGTACGTCACCACCTTTCAGACGCAGATCCGCTTCGACGGCGACAGCGCCGTGCTCGATG  >  NC_002947/3139670‑3139954

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