Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,139,813 A→G 9.3% S123G (AGC→GGC)  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,8130AG 9.3% 114.0 / 2.3 43S123G (AGC→GGC) sfnBSfnB family sulfur acquisition oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (19/20);  new base G (2/2);  total (21/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGTGCCCAAGGCGCACGGCGGCGCCGGGGTATCTTATGTCACCCTCAGCGAGGTGATCGCGACCATCTCCGCCGCCGATTCATCACTCGGTCAGT‑‑TGCCGCAGAACCACTTTGGTGTGCTCAGCAACCTCGGCCTGACCGGCAGCGAGGCGCAGAAGCACCGCTACTACGCCAAGGTGCTGCAAGGCTACCGTTTCGGTAACGCCTTTTCCGAGGCGCGCAGCCAGTACGTCACCACCTTTCAGACGCAGA‑TCCGCTTCGACGGCGACAGCGCCGTGCTCGATGGCGA  >  NC_002947/3139670‑3139958
                                                                                                                                                 |                                                                                                                                                  
gtgtGCCCAAGGCGCACGGCGGCGCCGGGGTATCTTATGTCACCCTCAGCGAGGTGATCGCGACCATCTCCGCCGCCGATTCATCACTCGGTCAGT‑‑TGCCGCAGAACCACTTTGGTGTGCTCAGCAACCTCGGCCTGACCGGCAGCGAg                                                                                                                                               <  2:213055/149‑1 (MQ=255)
   tgCCCAAGGCGCACGGCGGCGCCGGGGTATCTTATGTCACCCTCAGCGAGGTGATCGCGACCATCTCCGCCGCCGATTCATCACTCGGTCAGTTCTATCGGAGAACCACTTTGGTGTGCTCAGCAACCTCGGCCTGACCGGCAGCGAgg                                                                                                                                              >  1:112706/1‑149 (MQ=255)
      ccAAGGCGCACGGCGGCGCCGGGGTATCTTATGTCACCCTCAGCGAGGTGATCGCGACCATCTCCGCCGCCGATTCATCACTCGGTCAGT‑‑TGCCGCAGAACCACTTTGGTGTGCTCAGCAACCTCGGCCTGACCGGCAGCGAGgcg                                                                                                                                            >  3:174649/1‑146 (MQ=255)
         aGGCGCACGGCGGCGCCGGGGTATCTTATGTCACCCTCAGCGAGGTGATCGCGACCATCTCCGCCGCCGATTCATCACTCGGTCAGT‑‑TGCCGCAGAACCACTTTGGTGTGCTCAGCAACCTCGGCCTGACCGGCAGCGAGGCGCAGAAg                                                                                                                                      <  1:131312/149‑1 (MQ=255)
                       gcCGGGGTATCTTATGTCACCCTCAGCGAGGTGATCGCGACCATCTCCGCCGCCGATTCATCACTCGGTCAGT‑‑TGCCGCAGAACCACTTTGGTGTGCTCAGCAACCTCGGCCTGACCGGCAGCGAGGCGCAGAAGCACCGCTACTACGc                                                                                                                        <  6:364741/149‑1 (MQ=255)
                          ggggTATCTTATGTCACCCTCAGCGAGGTGATCGCGACCATCTCCGCCGCCGATTCATCACTCGGTCAGT‑‑TGCCGCAGAACCACTTTGGTGTGCTCAGCAACCTCGGCCTGACCGGCAGCGAGGCGCAGAAGCACCGCTACTACGCCaa                                                                                                                     <  5:195583/149‑1 (MQ=255)
                                     tGTCACCCTCAGCGAGGTGATCGCGACCATCTCCGCCGCCGATTCATCACTCGGTCAGT‑‑TGCCGCAGAACCACTTTGGTGTGCTCAGCAACCTCGGCCTGACCGGCAGCGAGGCGCAGAAGCACCGCTACTACGCCAAGGTGCTGCAAg                                                                                                          >  8:102395/1‑149 (MQ=255)
                                      gTCACCCTCAGCGAGGTGATCGCGACCATCTCCGCCGCCGATTCATCACTCGGTCAGT‑‑TGCCGCAGAACCACTTTGGTGTGCTCAGCAACCTCGGCCTGACCGGCAGCGAGGCGCAGAAGCACCGCTACTACGCCAAGGTGCTGCAAgg                                                                                                         >  4:287687/1‑149 (MQ=255)
                                                   aGGTGATCGCGACCATCTCCGCCGCCGATTCATCACTCGGTCAGT‑‑TGCCGCAGAACCACTTTGGTGTGCTCAGCAACCTCGGCCTGACCGGCAGCGAGGCGCAGAAGCACCGCTACTACGCCAAGGTGCTGCAAGGCTACCGTTTCGGt                                                                                            <  3:287687/149‑1 (MQ=255)
                                                   aGGTGATCGCGACCATCTCCGCCGCCGATTCATCACTCGGTCAGT‑‑TGCCGCAGAACCACTTTGGTGTGCTCAGCAACCTCGGCCTGACCGGCAGCGAGGCGCAGAAGCACCGCTACTACGCCAAGGTGCTGCAAGGCTACCGTTTCGGt                                                                                            <  3:54853/149‑1 (MQ=255)
                                                      tGATCGCGACCATCTCCGCCGCCGATTCATCACTCGGTCAGT‑‑TGCCGCAGAACCACTTTGGTGTGCTCAGCAACCTCGGCCTGACCGGCAGCGAGGCGCAGAAGCACCGCTACTACGCCAAGGTGCTGCAAGGCTACCGTTTCGGTAAc                                                                                         <  6:1903/149‑1 (MQ=255)
                                                           gcgACCATCTCCGCCGCCGATTCATCACTCGGTCAGT‑‑TGCCGCAGAACCACTTTGGTGTGCTCAGCAACCTCGGCCTGACCGGCAGCGAGGCGCAGAAGCACCGCTACTACGCCAAGGTGCTGCAAGGCTACCGTTTCGGTAACGCCtt                                                                                    >  3:180244/1‑149 (MQ=255)
                                                              aCCATCTCCGCCGCCGATTCATCACTCGGTCAGT‑‑TGCCGCAGAACCACTTTGGTGTGCTCAGCAACCTCGGCCTGACCGGCAGCGAGGCGCAGAAGCACCGCTACTACGCCAAGGTGCTGCAAGGCTACCGTTTCGGTAACGCCTTTTc                                                                                 <  8:83829/149‑1 (MQ=255)
                                                                cATCTCCGCCGCCGATTCATCACTCGGTCAGT‑‑TGCCGCAGAACCACTTTGGTGTGCTCAGCAACCTCGGCCTGACCGGCAGCGAGGCGCAGAAGCACCGCTACTACGCCAAGGTGCTGCAAGGCTACCGTTTCGGTAACGCCTTTTcc                                                                                >  5:317834/1‑148 (MQ=255)
                                                                    tcCGCCGCCGATTCATCACTCGGTCAGT‑‑TGCCGCAGAACCACTTTGGTGTGCTCAGCAACCTCGGCCTGACCGGCAGCGAGGCGCAGAAGCACCGCTACTACGCCAAGGTGCTGCAAGGCTACCGTTTCGGTAACGCCTTTTCCGAGgc                                                                           >  1:312571/1‑149 (MQ=255)
                                                                      cgccgccgATTCATCACTCGGTCAGT‑‑TGCCGCAGAACCACTTTGGTGTGCTCAGCAACCTCGGCCTGACCGGCAGCGAGGCGCAGAAGCACCGCTACTACGCCAAGGTGCTGCAAGGCTACCGTTTCGGTAACGCCTTTTCCGAg                                                                             <  4:174649/145‑1 (MQ=255)
                                                                      cgccgccgATTCATCACTCGGTCAGT‑‑TGCCGCAGAACCACTTTGGTGTGCTCAGCAACCTCGGCCTGACCGGCAGCGAGGCGCAGAAGCACCGCTACTACGCCAAGGTGCTGCAAGGCTACCGTTTCGGTAACGCCTTTTCCGAGgcgc                                                                         <  7:102395/149‑1 (MQ=255)
                                                                                  atcaCTCGGTCAGT‑‑TGCCGCAGAACCACTTTGGTGTGCTCAGCAACCTCGGCCTGACCGGCAGCGAGGCGCAGAAGCACCGCTACTACGCCAAGGTGCTGCAAGGCTACCGTTTCGGTAACGCCTTTTCCg                                                                               >  6:94283/1‑131 (MQ=255)
                                                                                  atcaCTCGGTCAGT‑‑TGCCGCAGAACCACTTTGGTGTGCTCAGCAACCTCGGCCTGACCGGCAGCGAGGCGCAGAAGCACCGCTACTACGCCAAGGTGCTGCAAGGCTACCGTTTCGGTAACGCCTTTTCCg                                                                               <  5:94283/131‑1 (MQ=255)
                                                                                   tcaCTCGGTCAGT‑‑TGCCGCAGAACCACTTTGGTGTGCTCAGCAACCTCGGCCTGACCGGCAGCGAGGCGCAGAAGCACCGCTACTACGCCAAGGTGCTGCAAGGCTACCGTTTCGGTAACGCCTTTTCCGAGGCGCGCAGCCAGTACGt                                                            >  3:62274/1‑149 (MQ=255)
                                                                                   tcaCTCGGTCAGT‑‑TGCCGCAGAACCACTTTGGTGTGCTCAGCAACCTCGGCCTGACCGGCAGCGAGGCGCAGAAGCACCGCTACTACGCCAAGGTGCTGCAAGGCTACCGTTTCGGTAACGCCTTTTCCGAGGCGCGCAGCCAGTACGt                                                            >  4:301678/1‑149 (MQ=255)
                                                                                    caCTCGGTCAGT‑‑TGCCGCAGAACCACTTTGGTGTGCTCAGCAACCTCGGCCTGACCGGCAGCGAGGCGCAGAAGCACCGCTACTACGCCAAGGTGCTGCAAGGCTACCGTTTCGGTAACGCCTTTTCCGAGGCGCGCAGCCAGTACGt                                                            >  5:266262/1‑148 (MQ=255)
                                                                                         ggTCAGT‑‑TGCCGCAGAACCACTTTGGTGTGCTCAGCAACCTCGGCCTGACCGGCAGCGAGGCGCAGAAGCACCGCTACTACGCCAAGGTGCTGCAAGGCTACCGTTTCGGTAACGCCTTTTCCGAGGCGCGCAGCCAGTACGTcacca                                                       <  3:301678/148‑1 (MQ=255)
                                                                                             aGT‑‑TGCCGCAGAACCACTTTGGTGTGCTCAGCAACCTCGGCCTGACCGGCAGCGAGGCGCAGAAGCACCGCTACTACGCCAAGGTGCTGCAAGGCTACCGTTTCGGTAACGCCTTTTCCGAGGCGCGCAGCCAGTACGTCACCACCttt                                                  >  4:228129/1‑149 (MQ=255)
                                                                                              gT‑‑TGCCGCAGAACCACTTTGGTGTGCTCAGCAACCTCGGCCTGACCGGCAGCGAGGCGCAGAAACACCGCTACTACGCCAAGGTGCTGCAAGGCTACCGTTTCGGTAACGCCTTTTCCGAGgcgc                                                                         <  3:89959/125‑1 (MQ=255)
                                                                                              gT‑‑TGCCGCAGAACCACTTTGGTGTGCTCAGCAACCTCGGCCTGACCGGCAGCGAGGCGCAGAAACACCGCTACTACGCCAAGGTGCTGCAAGGCTACCGTTTCGGTAACGCCTTTTCCGAGgcgc                                                                         >  4:89959/1‑125 (MQ=255)
                                                                                                   gCCGCAGAACCACTTTGGTGTGCTCAGCAACCTCGGCCTGACCGGCAGCGAGGCGCAGAAGCACCGCTACTACGCCAAGGTGCTGCAAGGCTACCGTTTCGGTAACGCCTTTTCCGAGGCGCGCAGCCAGTACGTCACCACCTTTCAg                                               <  6:266262/148‑1 (MQ=255)
                                                                                                            ccACTTTGGTGTGCTCAGCAACCTCGGCCTGACCGGCAGCGAGGCGCAGAAGCACCGCTACTACGCCAAGGTGCTGCAAGGCTACCGTTTCGGTAACGCCTTTTCCGAGGCGAGCAGCCAGTACGTCACCACCTTTCAGACGCAGA‑Tcc                                    <  4:180244/149‑1 (MQ=255)
                                                                                                             cACTTTGGTGTGCTCAGCAACCTCGGCCTGACCGGCAGCGAGGCGCTGAAGCACCGCTACTACGCCAAGGTGCTGCAAGGCTACCGTTTCGGTAACGCCTTTTCCGAGGCGCGCAGCCAGTACGTCACCACCTTTCAGACGCAGA‑TCCg                                   >  7:163595/1‑149 (MQ=255)
                                                                                                               cTTTGGTGTGCTCAGCAACCTCGGCCTGACCGGCGGCGAGGCGCAGAAGCACCGCTACTACGCCAAGGTGCTGCAAGGCTACCGTTTCGGTAACGCCTTTTCCGAGGCGCGCAGCCAGTACGTCACCACCTTTCAGACGCAGA‑TCCGCt                                 >  7:125546/1‑149 (MQ=255)
                                                                                                                  tGGTGTGCTCAGCAACCTCGGCCTGACCGGCAGCGAGGCGCAGAAGCACCGCTACTACGCCAAGGTGCTGCAAGGCTACCGTTTCGGTAACGCCTTTTCCGAGGCGCGCAGCCAGTACGTCACCACCTTTCAGACGCAGA‑TCCGCtt                                <  2:112706/147‑1 (MQ=255)
                                                                                                                  tGGTGTGCTCAGCAACCTCGGCCTGACCGGCAGCGAGGCGCAGAAGCACCGCTACTACGCCAAGGTGCTGCAAGGCTACCGTTTCGGTAACGCCTTTTCCGAGGCGCGCAGCCAGTACGTCACCACCTTTCAGACGCAGA‑TCCGCTTc                               <  6:317834/148‑1 (MQ=255)
                                                                                                                    gtgtGCTCAGCAACCTCGGCCTGACCGGCGGCGAGGCGCAGAAGCACCGCTACTACGCCAAGGTGCTGCAAGGCTACCGTTTCGGTAACGCCTTTTCCGAGGCGCGCAGCCAGTACGTCACCACCTTTCAGACGCAGA‑TCCGCTTCGAc                            <  8:125546/149‑1 (MQ=255)
                                                                                                                            aGCAACCTCGGCCTGACCGGCAGCGAGGCGCAGAAGCACCGCTACTACGCCAAGGTGCTGCAAGGCTACCGTTTCGGTAACGCCTTTTc                                                                                 >  6:358967/1‑89 (MQ=255)
                                                                                                                            aGCAACCTCGGCCTGACCGGCAGCGAGGCGCAGAAGCACCGCTACTACGCCAAGGTGCTGCAAGGCTACCGTTTCGGTAACGCCTTTTc                                                                                 <  5:358967/89‑1 (MQ=255)
                                                                                                                                 ccTCGGCCTGACCGGCGGCGAGGCGCAGAAGCACCGCTACTACGCCAAGGTGCTGCAAGGCTACCGTTTCGGTAACGCCTTTTCCGAGGCGCGCAGCCAGTACGTCACCACCTTTCAGACGCAGA‑TCCGCTTCGACGGCGACAGCGCCg               >  8:6959/1‑149 (MQ=255)
                                                                                                                                  cTCGGCCTGACCGGCAGCGAGGCGCAGAAGCACCGCTACTACGCCAAGGTGCTGCAAGGCTACCGTTTCGGTAACGCCTTTTCCGAGGCGCGCAGCCAGTACGTCACCACCTTTCAGACGCAGA‑TCCGCTTCGACGGCGACAGCGCCGt              >  2:57042/1‑149 (MQ=255)
                                                                                                                                         tGACCGGCAGCGAGGCGCAGAAGCACCGCTACTACGCCAAGGTGCTGCAAGGCTACCGTTTCGGTAACGCCTTTTCCGAGGCGCGCAGCCAGTACGTCACCACCTTTCAGACGCAGA‑TCCGCTTCGACGGCGACAGCGCCGTGCTCGAt       >  4:244545/1‑149 (MQ=255)
                                                                                                                                          gACCGGCAGCGAGGCGCAGAAGCACCGCTACTACGCCAAGGTGCTGCAAGGCTACCGTTTCGGTAACGCCTTTTCCGAGGCGCGCAGCCAGTACGTCACCACCTTTCAGACGCAGA‑TCCGCTTCGACGGCGACAGCGCCGTGCTCGATg      >  8:111684/1‑149 (MQ=255)
                                                                                                                                           aCCGGCAGCGAGGCGCAGAAGCACCGCTACTACGCCAAGGTGCTGCAAGGCTACCGTTTCGGTAACGCCTTTTCCGAGGCGCGCAGCCAGTACGTCACCACCTTTCAGACGCAGA‑TCCGCTTCGACGGCGACAGCGCCGTGCTCGATgg     <  2:312571/149‑1 (MQ=255)
                                                                                                                                            ccGGCGGCGAGGCGCAGAAGCACCGCTACTACGCCAAGGTGCTGCAAGGCTACCGTTTCGGTAACGCCTTTTCCGAGGCGCGCAGCCAGTACGTCACCACCTTTCAGACGCAGA‑TCCGCTTCGACGGCGACAGCGCCGTGCTCGATGGc    <  7:6959/149‑1 (MQ=255)
                                                                                                                                               gCAGCGAGGCGCAGAAGCACCGCTACTACGCCAAGGTGCTGCAAGGCTACCGTTTCGGTAACGCCTTTTCCGAGGCGCGCAGCCAGTACGTCACCACCTTTCAGACGCAGATTCCGCTTCGACGGCGACAGCGCCGTGCTCGATGGCga  <  1:57042/149‑1 (MQ=255)
                                                                                                                                               gCAGCGAGGCGCAGAAGCACCGCTACTACGCCAAGGTGCTGCAAGGCTACCGTTTCGGTAACGCCTTTTCCGAGGCGCGCAGCCAGTACGTCACCACCTTTCAGACGCAGA‑TCCGCTTCGACGGCGACAGCGCCGTGCTCGATGGc    >  7:14727/1‑146 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
GTGTGCCCAAGGCGCACGGCGGCGCCGGGGTATCTTATGTCACCCTCAGCGAGGTGATCGCGACCATCTCCGCCGCCGATTCATCACTCGGTCAGT‑‑TGCCGCAGAACCACTTTGGTGTGCTCAGCAACCTCGGCCTGACCGGCAGCGAGGCGCAGAAGCACCGCTACTACGCCAAGGTGCTGCAAGGCTACCGTTTCGGTAACGCCTTTTCCGAGGCGCGCAGCCAGTACGTCACCACCTTTCAGACGCAGA‑TCCGCTTCGACGGCGACAGCGCCGTGCTCGATGGCGA  >  NC_002947/3139670‑3139958

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