Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	983,256	T→C	8.9%	V89V (GTT→GTC)	fdx →	ferredoxin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	983,256	0	T	C	8.9%	119.7 / 2.8	45	V89V (GTT→GTC)	fdx	ferredoxin
Reads supporting (aligned to +/- strand):  ref base T (21/20);  new base C (2/2);  total (23/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 = 9.33e-01

CGCCTGCGGCGGCGTCAAGGCCTGTACCACTTGCCACTGCATCGTGCGCAAGGGCTTCGATTCGCTCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATGAGTCTGAAATGGGTTGATGTACTTGAGATCGCAATCCAGCTTGCAGAAAGCAAGC  >  NC_002947/983109‑983399                                                                                                                                                    |                                                                                                                                                cgcCTGCGGCGGCGTCAAGGCCTGTACCACTTGCCACTGCATCGTGCGCAAGGGCTTCGATTCGCTCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGttt                                                                                                                                                <  4:223382/149‑1 (MQ=255)   ccTGCGGCTGCGTCAAGGCCTGTACCACTTGCCACTGCATCGTGCGCAAGGGCTTCGATTCGCTCAAAGAGGCTGACGAGCTGGAAGTGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTtc                                                                                                                                              <  3:236512/149‑1 (MQ=255)       cggcggcgTCAAGGCCTGTACCACTTGCCACTGCATCGTGCGCAAGTGCTTCGATTCGCTCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTtcgt                                                                                                                                            <  8:289344/147‑1 (MQ=255)       cggcggcgTCAAGGCCTGTACCACTTGCCACTGCATCGTGCGCAAGGGCTTCGATTCGCTCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTtcgt                                                                                                                                            <  3:228330/147‑1 (MQ=255)               tCAAGGCCTGTACCACTTGCCACTGCATCGTGCGCAAGGGCTTCGATTCGCTCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGatg                                                                                                                                    >  2:306558/1‑145 (MQ=255)                 aaGGCCTGTACCACTTGCCACTGCATCGTGCGCAAGGGCTTCGATTCGCTCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTCTTCGTCAGCGACGaaaa                                                                                                                                >  4:293314/1‑147 (MQ=255)                  aGGCCTGTACCACTTGCCACTGCATCGTGCGCAAGGGCTTCGATTCGCTCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTCTTCGTCAGCGACGAAAAc                                                                                                                               <  3:293314/149‑1 (MQ=255)                    gCCTGTACCACTTGCCACTGCATCGTGCGCAAGGGCTTCGATTCGCTCGAAGAGGCTGACGAGCTGGAAGAGGACAAGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCt                                                                                                                             >  2:136669/1‑149 (MQ=255)                      cTGTACCACTTGCCACTGCATCGTGCGCAAGGGCTTCGATTCGCTCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTCTTCGTCAGCGACGAAAACCTgg                                                                                                                           <  1:45791/149‑2 (MQ=255)                          aCCACTTGCCACTGCATCGTGCGCAAGGGCTTCGATTCGCTCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTtc                                                                                                                                              >  1:83953/1‑126 (MQ=255)                          aCCACTTGCCACTGCATCGTGCGCAAGGGCTTCGATTCGCTCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTtc                                                                                                                                              <  2:83953/126‑1 (MQ=255)                           ccACTTGCCACTGCATCGTGCGCAAGGGCTTCGATTCGCTCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGACCATc                                                                                                                      <  1:316900/149‑1 (MQ=255)                                      tGCATCGTGCGCAAGGGCTTCGATTCGCTCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCt                                                                                                                                        <  5:56972/120‑1 (MQ=255)                                      tGCATCGTGCGCAAGGGCTTCGATTCGCTCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCt                                                                                                                                        >  6:56972/1‑120 (MQ=255)                                      tGCATCGTGCGCAACGGCTTCGATTCGCTCGAAGAGGCTGACGAGCTGCAAGAGGACATGCTGGACAAAGCCTGGGGCCTCGAAGCGCAATCGCGCCTGCGTTGCCAGGATATCCTCGCTGACGAAGACCTGACCAACTAGATcccc                                                                                                             >  2:198863/1‑147 (MQ=255)                                        cATCGTGCGCAAGGGCTTCGATTCGCTCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCaa                                                                                                           <  5:34747/147‑1 (MQ=255)                                        cATCGTGCGCAAGGGCTTCGATTCGCTCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCaa                                                                                                           <  6:116939/147‑1 (MQ=255)                                             tGCGCAAGGGCTTCGATTCGCTCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCAAGTATTCg                                                                                                    <  8:102708/149‑1 (MQ=255)                                                        ttCGATTCGCTCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCAAGTATTCGCTCAACCAtgc                                                                                         >  6:57548/1‑149 (MQ=255)                                                         tCGATTCGCTCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCAAGTATTCGCTCAACCAtgct                                                                                        >  6:11439/1‑149 (MQ=255)                                                                tcTCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAg                                                                                 >  8:298701/2‑149 (MQ=255)                                                                gCTCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAg                                                                                 <  3:57524/149‑1 (MQ=255)                                                                     aaGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCg                                                                            >  3:210858/1‑149 (MQ=255)                                                                           cTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGa                                                                      <  7:298701/149‑1 (MQ=255)                                                                              aCGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAg                                                                     <  5:11439/147‑1 (MQ=255)                                                                                     ggAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTa                                                             >  3:233491/1‑148 (MQ=255)                                                                                        agatGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATg                                                         <  7:205020/149‑1 (MQ=255)                                                                                             acttggtGGAAAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGTGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCGCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATGAGTc                                                     <  5:326046/142‑1 (MQ=255)                                                                                                 gCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATGAGTCTGaaa                                                >  3:251994/1‑149 (MQ=255)                                                                                                  cTGGACATAGGGTGGGGCCTTGAAGCGCAATCGCGGCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGACCTTCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATGAGTCTGAAAt                                               <  1:198863/149‑1 (MQ=255)                                                                                                   tGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATGAGTCTGAAATg                                              >  7:36930/1‑149 (MQ=255)                                                                                                      aCAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATGAGTCTGAAATGGGt                                           <  4:233491/149‑1 (MQ=255)                                                                                                        aaaGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATGAGTCTGAAATGGGtt                                          >  8:27898/1‑148 (MQ=255)                                                                                                                     tGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGACGCTGACGAAGACCTTACCATCCAGATCCCCAAGTATTCGCTCAACAATCATGCCGAAGCACCCCACTGAGGGCAAGTACATGAGTCTGAAATGGGTTGATGTACtt                                 >  4:285981/1‑144 (MQ=255)                                                                                                                          gcgcAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATGAGTCTGAAATGGGTTGATGTACTTGAGATCGCaa                       >  2:93010/1‑149 (MQ=255)                                                                                                                                  gcgcCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATGAGTCTGAAATGGGTTGATGTACTTGAGATCGCAATCCAGCt                <  4:251994/148‑1 (MQ=255)                                                                                                                                        gggTTGCCAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATGAGTCTGAAATGGGTTGATGTACTTGAGATCGCAATCCAGCTTGCAGa          >  5:310358/1‑148 (MQ=255)                                                                                                                                         ggTTGCCAGGTTTTCGTCGCTGATGAAGACCTGACCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATGAGTCTGAAATGGGTTGATGTACTTGAGATCGCAATCCAGCTTGCAGaaa        <  1:306558/149‑1 (MQ=255)                                                                                                                                            ggtcAGGTCTTCGTCAGCGACGAAAACCTGACCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATGAGTCTGAAATGGGTTGATGTACTTGAGATCGCAATCCAGCTTGCAGAaagc      >  4:165765/4‑148 (MQ=255)                                                                                                                                              ccAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATGAGTCTGAAATGGGTTGATGTACTTGAGATCGCAATcc                    <  2:311087/132‑1 (MQ=255)                                                                                                                                              ccAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATGAGTCTGAAATGGGTTGATGTACTTGAGATCGCAATcc                    >  1:311087/1‑132 (MQ=255)                                                                                                                                              ccAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATGAGTCTGAAATGGGTTGATGTACTTGAGATCGCAATCCAGCTTGCAGAaagcaag   >  2:322230/1‑149 (MQ=255)                                                                                                                                              ccAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATCAGTCTGAAATGCGTTGATGTACTTGAGATCGCAATCCAGCTTTCAGAaagcaag   >  7:132686/1‑149 (MQ=255)                                                                                                                                               cAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATGAGTCTGAAATGGGTTGATGTACTTGAGATCGCAATCCAGCTTGCAGAaagcaagc  >  5:4110/1‑149 (MQ=255)                                                                                                                                               cAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATGAGTCTGAAATGGGTTGATGTACTTGAGATCGCAATCCAGCTTGCAGAaagcaagc  >  8:162031/1‑149 (MQ=255)                                                                                                                                               cAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATGAGTCTGAAATGGGTTGATGTACTTGAGATCGCAATCCAGCTTGCAGAaagcaagc  <  8:36930/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                CGCCTGCGGCGGCGTCAAGGCCTGTACCACTTGCCACTGCATCGTGCGCAAGGGCTTCGATTCGCTCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATGAGTCTGAAATGGGTTGATGTACTTGAGATCGCAATCCAGCTTGCAGAAAGCAAGC  >  NC_002947/983109‑983399

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