Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,139,457	A→G	8.1%	Q4R (CAG→CGG)	sfnB →	SfnB family sulfur acquisition oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,139,457	0	A	G	8.1%	135.5 / 3.3	49	Q4R (CAG→CGG)	sfnB	SfnB family sulfur acquisition oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (21/24);  new base G (2/2);  total (23/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.96e-01

GCCCGGCAACAGCCAATCCACAGTTTGCCGCCACAGGTGCTTTTCCGATTATCTGAATTGCACGTCAAGTCATTGAATAAGAAGGATTTCCTGAGTTGGTACGGGCTGTGCACTAGCCCTGGGGAATCCATCACGCT‑TGGAGCAACAGATGACCCAGTCACAACCCGCCGATCTGCGTGAGCAGGTCAGCTCGGCGCCACGGCGCCCGGCACAGATCATCCGTAGCGACGCCGAAGCCATCGAGGTGGCGCAGCGCTTGGCCGAGGATTTTGCCCGTGAGGCGGCCTGGCGT  >  NC_002947/3139311‑3139602                                                                                                                                                    |                                                                                                                                                  gCCCGGCAACAGCCAATCCACAGTTTGCCGCCACAGGTGCTTTTCCGATTATCTGAATTGCACGTCAAGTCATTGAATAAGAAGGATTTCCTGAGTTGGTACGGGCTGTGCACTAGCCCTGGGGAATCCATCACGCT‑TGGAGCAACAGa                                                                                                                                                 <  3:67717/149‑1 (MQ=255)    cGGCAACAGCCAATCCACAGTTTCCCGACAGAGGTGCTTTTCCGATTATCTGAATTGCACGTCAAGTCATTGAATAAGAAGGATTTCCTGAGGTGGTACGGGCTGTGCACTAGCCCTGGGGAATCCATCACGCT‑TGGAGCAACAGATGa                                                                                                                                              <  3:275401/149‑1 (MQ=255)            gCCAATCCACAGTTTGCCGCCACAGGTGCTTTTCCGATTATCTGAATTGCACGTCAAGTCATTGAATAAGAAGGATTTCCTGAGTTGGTACGGGCTGTGCACTAGCCCTGGGGAATCCATCACGCT‑TGGAGCAACAGATGACCCAGTca                                                                                                                                      >  7:95411/1‑149 (MQ=255)              cAATCCACAGTTTGCCGCCACAGGTGCTTTTCCGATTATCTGAATTGCACGTCAAGTCATTGAATAAGAAGGATTTCCTGAGTTGGTACGGGCTGTGCACTAGCCCTGGGGAATCCATCACGCT‑TGGAGCAACAGATGACCCAGTcaca                                                                                                                                    <  5:222518/149‑1 (MQ=255)                  ccACAGTTTGCCGCCACAGGTGCTTTTCCGATTATCTGAATTGCACGTCAAGTCATTGAATAAGAAGGATTTCCTGAGTTGGTACGGGCTGTGCACTAGCCCTGGGGAATCCATCACGCT‑TGGAGCAACAGATGACCCAGTCACAAccc                                                                                                                                >  5:314307/1‑149 (MQ=255)                    acaGTTTGCCGCCACAGGTGCTTTTCCGATTATCTGAATTGCACGTCAAGTCATTGAATAAGAAGGATTTCCTGAGTTGGTACGGGCTGTGCACTAGCCCTGGGGAATCCATCACGCT‑TGGAGCAACAGATGACCCAGTCACAACccgc                                                                                                                              >  2:194475/1‑149 (MQ=255)                       gTTTGCCGCCACAGGTGCTTTTCCGATTATCTGAATTGCACGTCAAGTCATTGAATAAGAAGGATTTCCTGAGTTGGTACGGGCTGTGCACTAGCCCTGGGGAATCCATCACGCT‑TGGAGCAACAGATGACCCAGTCACAACCCGCCGa                                                                                                                           <  7:113478/149‑1 (MQ=255)                                 acaGGTGCTTTTCCGATTATCTGAATTGCACGTCAAGTCATTGAATAAGAAGGATTTCCTGAGTTGGTACGGGCTGTGCACTAGCCCTGGGGAATCCATCACGCT‑TGGAGCAACAGATGACCCAGTCACAACCCGCCGATCTGCGTGAg                                                                                                                 >  6:234752/1‑149 (MQ=255)                                    ggTGCTTTTCCGATTATCTGAATTGCACGTCAAGTCATTGAATAAGAAGGATTTCCTGAGTTGGTACGGGCTGTGCACTAGCCCTGGGGAATCCATCACGCT‑TGGAGCAACAGATGACCCAGTCACAACCCGCCGATCTGCGTGAGCAg                                                                                                              <  3:70588/149‑1 (MQ=255)                                               gATTATCTGAATTGCACGTCAAGTCATTGAATAAGAAGGATTTCCTGAGTTGGTACGGGCTGTGCACTAGCCCTGGGGAATCCATCACGCT‑TGGAGCAACAGATGACCCAGTCACAACCCGCCGATCTGCGTGAGCAGGTCAGCTCGgc                                                                                                   <  5:234752/149‑1 (MQ=255)                                                      tGAATTGCACGTCAAGTCATTGAATAAGAAGGATTTCCTGAGTTGGTACGGGCTGTGCACTAGCCCTGGGGAATCCATCACGCT‑TGGAGCAACAGATGACCCAGTCACAACCCGCCGATCTGCGTGAGCAGGTCAGCTCGGCGCCACgg                                                                                            >  3:9421/1‑149 (MQ=255)                                                      tGAATTGCACGTCAAGTCATTGAATAAGAAGGATTTCCTGAGTTGGTACGGGCTGTGCACTAGCCCTGGGGAATCCATCACGCT‑TGGAGCAACAGATGACCCAGTCACAACCCGCCGATCTGCGTGAGCAGGTCAGCTCGGCGCCACgg                                                                                            >  5:131061/1‑149 (MQ=255)                                                                    aGTCATTGAATAAGAAGGATTTCCTGAGTTGGTACGGGCTGTGCACTAGCCCTGGGGAATCCATCACGCT‑TGGAGCAACAGATGACCCAGTCACAACCCGCCGATCTGCGTGAGCAGGTCAGCTCGGCGCCACGGCGCCCGGca                                                                                   >  8:36919/1‑144 (MQ=255)                                                                     gTCATTGAATAAGAAGGATTTCCTGAGTTGGTACGGGCTGTGCACTAGCCCTGGGGAATCCATCACGCT‑TGGAGCAACAGATGACCCAGTCACAACCCGCCGATCTGCGTGAGCAGGTCAGCTCGGCGCCACGGCGCCCGGCACAGatc                                                                             <  8:95411/149‑1 (MQ=255)                                                                      tCATTGAATAAGAAGGATTTCCTGAGTTGGTACGGGCTGTGCACTAGCCCTGGGGAATCCATCACGCTCTGGAGCAACAGATGACCCAGTCACAACCCGCCGATCTGCGTGAGc                                                                                                                <  1:245295/114‑1 (MQ=255)                                                                        aTTGAATAAGAAGGATTTCCTGAGTTGGTACGGGCTGTGCACTAGCCCTGGGGAATCCATCACGCT‑TGGAGCAACAGATGACCCAGTCACAACCCGCCGATCTGCGTGAGCAGGTCAGCTCGGCGCCACGGCGCCCGGCACAGatcatc                                                                          >  6:91558/1‑149 (MQ=255)                                                                        aTTGAATAAGAAGGATTTCCTGAGTTGGTACGGGCTGTGCACTAGCCCTGGGGAATCCATCACGCT‑TGGAGCAACAGATGACCCAGTCACAACCCGCCGATCTGCGTGAGCAGGTCAGCTCGGCGCCACGGCGCCCGGCACAGatcatc                                                                          >  4:7030/1‑149 (MQ=255)                                                                        aTTGAATAAGAAGGATTTCCTGAGTTGGTACGGGCTGTGCACTAGCCCTGGGGAATCCATCACGCT‑TGGAGCAACAGATGACCCAGTCACAACCCGCCGATCTGCGTGAGCAGGTCAGCTCGGCGCCACGGCGCCCGGCACAGat                                                                              >  6:107624/1‑145 (MQ=255)                                                                        aTTGAATAAGAAGGATTTCCTGAGTTGGTACGGGCTGTGCACTAGCCCTGGGGAATCCATCACGCT‑TGGAGCAACAGATGACCCAGTCACAACCCGCCGATCTGCGTGAGCAGGTCAGCTCGGCGCCACGGCGCCCGGCACAGat                                                                              <  5:107624/145‑1 (MQ=255)                                                                         ttGAATAAGAAGGATTTCCTGAGTTGGTACGGGCTGTGCACTAGCCCTGGGGAATCCATCACGCT‑TGGAGCAACAGATGACCCAGTCACAACCCGCCGATCTTCGTGAGCAGGTCAGCTCGGCGCCACGGCGCCCGGCACAGATCATcc                                                                         >  5:284860/1‑149 (MQ=255)                                                                         ttGAATAAGAAGGATTTCCTGAGTTGGTACGGGCTGTGCACTAGCCCTGGGGAATCCATCACGCT‑TGGAGCAACAGATGACCCAGTCACAACCCGCCGATCTGCGTGAGCAGGTCAGCTCGGCGCCAc                                                                                              <  1:133781/128‑1 (MQ=255)                                                                         ttGAATAAGAAGGATTTCCTGAGTTGGTACGGGCTGTGCACTAGCCCTGGGGAATCCATCACGCT‑TGGAGCAACAGATGACCCAGTCACAACCCGCCGATCTGCGTGAGCAGGTCAGCTCGGCGCCAc                                                                                              >  2:133781/1‑128 (MQ=255)                                                                         ttGAATAAGAAGGATTTCCTGAGTTGGTACGGGCTGTGCACTAGCCCTGGGGAATCCATCACGCT‑TGGAGCAACAGATGACCCAGTCACAACCCGCCGATCTGCGTGAGCAGGTCAGCTCGGCGCCACGGCGCCCGGCACAGATCATcc                                                                         <  6:131061/149‑1 (MQ=255)                                                                                agaagGATTTCCTGAGTTGGTACGGGCTGTGCACTAGCCCTGGGGAATCCATCACGCT‑TGGAGCAACAGATGACCCAGTCACAACCCGCCGATCTGCGTGAGCAGGTCAGCTCGGCGCCACGGCGCCCGGCACAGATCATCCGTAGCGa                                                                  >  7:207582/1‑149 (MQ=255)                                                                                 gaagGATTTCCTGAGTTGGTACGGGCTGTGCACTAGCCCTGGGGAATCCATCACGCT‑TGGAGCAACAGATGACCCAGTCACAACCCGCCGATCTGCGTGAGCAGGTCAGCTCGGCGCCACGGCGCCCGGCACAGATCATCCGTAGCGAc                                                                 >  8:182981/1‑149 (MQ=255)                                                                                    ggATTTCCTGAGTTGGTACGGGCTGTGCACTAGCCCTGGGGAATCCATCACGCT‑TGGAGCAACAGATGACCCAGTCACAACCCGCCGATCTGCGTGAGCAGGTCAGCTCGGCGCCACGGCGCCCGGCACAGATCATCCGTAGCGACGcc                                                              >  4:250166/1‑149 (MQ=255)                                                                                       tttCCTGAGTTGGTACGGGCTGTGCACTAGCCCTGGGGAATCCATCACGCT‑TGGAGCAACAGATGACCCAGTCACAACCCGCCGATCTGCGTGAGCAGGTCAGCTCGGCGCCACGGCGCCCGGCACAGATCATCCGTAGCGACGCCGaa                                                           <  7:182981/149‑1 (MQ=255)                                                                                         tCCTGAGTTGGTACGGGCTGTGCACTAGCCCTGGGGAATCCATCACGCT‑TGGAGCAACAGATGACCCAGTCACAACCCGCCGATCTGCGTGAGCAGGTCAGCTCGGCGCCACGGCGCCCGGCACAGATCATCCGTAGCGACGCCGAAg                                                          <  8:207582/148‑1 (MQ=255)                                                                                          ccTGAGTTGGTACGGGCTGTGCACTAGCCCTGGGGAATCCATCACGCT‑TGGAGCAACAGATGACCCAGTCACAACCCGCCGATCTGCGTGAGCAGGTCAGCTCGGCGCCACGGCGCCCGGCACAGat                                                                              >  5:10564/1‑127 (MQ=255)                                                                                          ccTGAGTTGGTACGGGCTGTGCACTAGCCCTGGGGAATCCATCACGCT‑TGGAGCAACAGATGACCCAGTCACAACCCGCCGATCTGCGTGAGCAGGTCAGCTCGGCGCCACGGCGCCCGGCACAGat                                                                              <  6:10564/127‑1 (MQ=255)                                                                                            tGAGTTGGTACGGGCTGTGCACTAGCCCTGGGGAATCCATCACGCT‑TGGAGCAACGGATGACCCAGTCACAACCCGCCGATCTGCGTGAGCAGGTCAGCTCGGCGCCACGGCGCCCGGCACAGatca                                                                            >  4:105650/1‑127 (MQ=255)                                                                                            tGAGTTGGTACGGGCTGTGCACTAGCCCTGGGGAATCCATCACGCT‑TGGAGCAACGGATGACCCAGTCACAACCCGCCGATCTGCGTGAGCAGGTCAGCTCGGCGCCACGGCGCCCGGCACAGatca                                                                            >  4:100808/1‑127 (MQ=255)                                                                                            tGAGTTGGTACGGGCTGTGCACTAGCCCTGGGGAATCCATCACGCT‑TGGAGCAACGGATGACCCAGTCACAACCCGCCGATCTGCGTGAGCAGGTCAGCTCGGCGCCACGGCGCCCGGCACAGatca                                                                            <  3:105650/127‑1 (MQ=255)                                                                                            tGAGTTGGTACGGGCTGTGCACTAGCCCTGGGGAATCCATCACGCT‑TGGAGCAACGGATGACCCAGTCACAACCCGCCGATCTGCGTGAGCAGGTCAGCTCGGCGCCACGGCGCCCGGCACAGatca                                                                            <  3:100808/127‑1 (MQ=255)                                                                                               gTTGGTACGGGCTGTGCACTAGCCCTGGGGAATCCATCACGCT‑TGGAGCAACAGATGACCCAGTCACAACCCGCCGATCTGCGTGAGCAGGTCAGCTCGGCGCCACGGCGCCCGGCACAGTTCATCCGTAGCGACGCCGAAGCCATCGa                                                   <  1:194475/149‑1 (MQ=255)                                                                                                ttGGTACGGGCTGTGCACTAGCCCTGGGGAATCCATCACGCT‑TGGAGCAACAGATGACCCAGTCACAACCCGCCGATCTGCGTGAGCAGGTCAGCTCGGCGCCACGGCGCCCGGCACAGATCATCCGTAGCGACGCCGAAGCCATCGAg                                                  <  6:284860/149‑1 (MQ=255)                                                                                                   gTACGGGCTGTGCACTAGCCCTGGGGAATCCATCACGCT‑TGGAGCAACAGATGACCCAGTCACAACCCGCCGATCTGCGTGAGCAGGTCAGCTCGGCGCCACGGCGCCCGGCACAGATCATCCGTAGCGACGCCGAAGCCATCGAGGTg                                               <  5:91558/149‑1 (MQ=255)                                                                                                                aCTAGCCCTGGGGAATCCATCACGCT‑TGGAGCAACAGATGACCCAGTCACAACCCGCCGATCTGCGTGAGCAGGTCAGCTCGGCGCCACGGCGCCCGGCACAGATCATCCGTAGCGACGCCGAAGCCATCGAGGTGGCGCAGCGCTTgg                                  <  3:7030/149‑1 (MQ=255)                                                                                                                 cTAGCCCTGGGGAATCCATCACGCT‑TGGAGCAACAGATGACCCAGTCACAACCCGCCGATCTGCGTGAGCAGGTCAGCTCGGCGCCACGGCGCCCGGCACAGATCATCCGTAGCGACGCCGAAGCCATCGAGGTGGCGCAGCGCTTgg                                  <  6:169644/148‑1 (MQ=255)                                                                                                                 cTAGCCCTGGGGAATCCATCACGCT‑TGGAGCAACAGATGACCCAGTCACAACCCGCCGATCTGCGTGAGCAGGTCAGCTCGGCGCCACGGCGCCCGGCACAGATCATCCGTAGCGACGCCGAAGCCATCGAGGTGGCGCAGCGCTTgg                                  >  5:169644/1‑148 (MQ=255)                                                                                                                     cccTGGGGAATCCATCACGCT‑TGGAGCAACAGATGACCCAGTCACAACCCGCCGATCTGCGTGAGCAGGTCAGCTCGGCGCCACGGCGCCCGGCACAGATCATCCGTAGCGACGCCGAAGCCATCGAGGTGGCGCAGCGCTTGGCCGAg                             >  4:298990/1‑149 (MQ=255)                                                                                                                      ccTGGGGAATCCATCACGCT‑TGGAGCAACAGATGACCCAGTCACAACCCGCCGATCTGCGTGAGCAGGTCAGCTCGGCGCCACGGCGCCCGGCACAGATCATCCGTAGCGACGCCGAAGCCATCGAGGTGGCGCAGCGCTTGGCCGAgg                            <  4:9421/149‑1 (MQ=255)                                                                                                                        tGGGGAATCCATCACGCT‑TGGAGCAACAGATGACCCAGTCACAACCCGCCGATCTGCGTGAGCAGGTCAGCTCGGCGCCACGGCGCCCgg                                                                                     <  7:52318/90‑1 (MQ=255)                                                                                                                        tGGGGAATCCATCACGCT‑TGGAGCAACAGATGACCCAGTCACAACCCGCCGATCTGCGTGAGCAGGTCAGCTCGGCGCCACGGCGCCCgg                                                                                     >  8:52318/1‑90 (MQ=255)                                                                                                                            gAATCCATCACGCT‑TGGAGCAACAGATGACCCAGTCACAACCCGCCGATCTGCGTGAGCAGGTCAGCTCAGCGCCACg                                                                                             <  4:14492/78‑1 (MQ=255)                                                                                                                            gAATCCATCACGCT‑TGGAGCAACAGATGACCCAGTCACAACCCGCCGATCTGCGTGAGCAGGTCAGCTCAGCGCCACg                                                                                             >  3:14492/1‑78 (MQ=255)                                                                                                                                  aTCACGCT‑TGGAGCAACAGATGACCCAGTCACAACCCGCCGATCTGCGTGAGCAGGTCAGCTCGGCGCCACGGCGCCCGGCACAGATCATCCGTAGCGACGCCGAAGCCATCGAGGTGGCGCAGCGCTTGGCCGAGGATTTTGCCCGt                 >  6:247974/1‑148 (MQ=255)                                                                                                                                               gCAACAGATGACCCAGTCACAACCCGCCGATCTGCGTGAGCAGGTCAGCTCGGCGCCACGGCGCCCGGCACAGATCATCCGTAGCGACGCCGAAGCCATCGAGGTGGCGCAGCGCTTGGCCGAGGATTTTGCCCGTGAGGCGGCCTGGc    <  6:314307/149‑1 (MQ=255)                                                                                                                                                 aaCAGATGACCCAGTCACAACCCGCCGATCTGCGTGAGCAGGTCAGCTCGGCGCCACGGCGCCCGGCACAGATCATCCGTAGCGACGCCGAAGCCATCGAGGTGGCGCAGCGCTTGGCCGAGGATTTTGCCCGTGAGGCGGCCTGGCGt  <  3:298990/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                  GCCCGGCAACAGCCAATCCACAGTTTGCCGCCACAGGTGCTTTTCCGATTATCTGAATTGCACGTCAAGTCATTGAATAAGAAGGATTTCCTGAGTTGGTACGGGCTGTGCACTAGCCCTGGGGAATCCATCACGCT‑TGGAGCAACAGATGACCCAGTCACAACCCGCCGATCTGCGTGAGCAGGTCAGCTCGGCGCCACGGCGCCCGGCACAGATCATCCGTAGCGACGCCGAAGCCATCGAGGTGGCGCAGCGCTTGGCCGAGGATTTTGCCCGTGAGGCGGCCTGGCGT  >  NC_002947/3139311‑3139602

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