Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	983,026	T→C	12.9%	F12L (TTC→CTC)	fdx →	ferredoxin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	983,026	0	T	C	12.9%	75.4 / 5.2	31	F12L (TTC→CTC)	fdx	ferredoxin
Reads supporting (aligned to +/- strand):  ref base T (12/15);  new base C (2/2);  total (14/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.58e-01

GTCTGTCGCAGGTGACCGATGCATTTGCCGCCCGTCGCCTTGATTCGACGGTCAAAGCCGCACTGGCCGGGCGCAACCTGAATGAGATCGAGGAGTAATCGATGCCGCTGGTGACATTCCTGCCGCATGAGAAGTTCTGCCCCGAGGGGCTGACCGTGGAGGTCAAGCCCGGGACCAACATCCTGGAGCTGGCCCACGACCATCACATCGAGATGGAAAGCGCCTGCGGCGGCGTCAAGGCCTGTACCACTTGCCACTGCATCG  >  NC_002947/982892‑983155                                                                                                                                       |                                                                                                                                  gTCTGTCGCAGGTGACCGATGCATTTGCCGCCCGTCGCCTTGATTCGACGGTCAAAGCCGCACTGGCCGGGCGCAACCTGAATGAGATCGAGGAGTAATCGATGCCGCTGGTGACATTCCTGCCGCATGAGAAGtt                                                                                                                                  <  2:411610/136‑1 (MQ=255)                    tGCATTTGCCGCCCGTCGCCTTGATTCGACGGTCAAAGCCGCACTGGCCGGGCGCAACCTGAATGAGATCGAGGAGTAATCGATGCCGCTGGTGACATTCCTGCCGCATGAGAAGTTCTGCCCCGAGGGGCTGAcc                                                                                                               <  1:215037/136‑1 (MQ=255)                    tGCATTTGCCGCCCGTCGCCTTGATTCGACGGTCAAAGCCGCACTGGCCGGGCGCAACCTGAATGAGATCGAGGAGTAATCGATGCCGCTGGTGACATTCCTGCCGCATGAGAAGTTCTGCCCCGAGGGGCTGAcc                                                                                                               <  1:28392/136‑1 (MQ=255)                       aTTTGCCGCCCGTCGCCGTGATTCCACGGTCAAAGCCGCACTGGCCGGGCGCAACCTGAATGAGATCGAGGAGTAATCGATGCCGCTGGTGACATTCCTGCCGCATGAGAAGTTCAGCCCCGAGGGGCTGACCGTg                                                                                                            >  2:404477/1‑136 (MQ=255)                                 cGTCGCCTTGATTCGACGGTCAAAGCCGCACTGGCCGGGCGCAACCTGAATGAGATCGAGGAGTAATCGATGCCGCTGGTGACATTCCTGCCGCATGAGAAGTTCTGCCCCGAGGGGCTGACCGTGGAGGTCAAGc                                                                                                  >  1:329102/1‑136 (MQ=255)                                     gCCTTGATTCGACGGTCAAAGCCGCACTGGCCGGGCGCAACCTGAATGAGATCGAGGAGTAATCGATGCCGCTGGTGACATTCCTGCCGCATGAGAAGTTCTGCCCCGAGGGGCTGACCGTGGAGGTCAAGCCCgg                                                                                              <  2:243994/136‑1 (MQ=255)                                     gCCTTGATTCGACGGTCAAAGCCGCACTGGCCGGGCGCAACCTGAATGAGATCGAGGAGTAATCGATGCCGCTGGTGACATTCCTGCCGCATGAGAAGTTCTGCCCCGAGGGGCTGACCGTGGAGGTCAAGCCCgg                                                                                              <  2:342306/136‑1 (MQ=255)                                      ccTTGATTCGACGGTCAAAGCCGCACTGGCCGGGCGCAACCTGAATGAGATCGAGGAGTAATCGATGCCGCTGGTGACATTCCTGCCGCATGAGAAGTTCTGCCCCGAGGGGCTGACCGTGGAGGTCAAGCCCggg                                                                                             >  2:131717/1‑136 (MQ=255)                                        ttGATTCGACGGTCAAAGCCGCACTGGCCGGGCGCAACCTGAATGAGATCGAGGAGTAATCGATGCCGCTGGTGACATTCCTGCCGCATGAGAAGTTCTGCCCCGAGGGGCTGACCGTGGAGGTCAAGCCCGGGAc                                                                                           <  1:92123/136‑1 (MQ=255)                                        ttGATTCGACGGTCAAAGCCGCACTGGCCGGGCGCAACCTGAATGAGATCGAGGAGTAATCGATGCCGCTGGTGACATTCCTGCCGCATGAGAAGTTCTGCCCCGAGGGGCTGACCGTGGAGGTCAAGCCCGGGAc                                                                                           <  2:8477/136‑1 (MQ=255)                                                   gTCAAAGCCGCACTGGCCGGGCGCAACCTGAATGAGATCGAGGAGTAATCGATGCCGCTGGTGACATTCCTGCCGCATGAGAAGTTCTGCCCCGAGGGGCTGACCGTGGAGGTCACGCCCGGGACCAACATCCTgg                                                                                >  2:256401/1‑136 (MQ=255)                                                   gTCAAAGCCGCACTGGCCGGGCGCAACCTGAATGAGATCGAGGAGTAATCGATGCCGCTGGTGACATTCCTGCCGCATGAGAAGTTCTGCCCCGAGGGGCTGACCGTGGAGGTCAAGCCCGGGACCAACATCCTgg                                                                                >  2:431164/1‑136 (MQ=255)                                                   gTCAAAGCCGCACTGGCCGGGCGCAACCTGAATGAGATCGAGGAGTAATCGATGCCGCTGGTGACATTCCTGCCGCATGAGAAGTTCTGCCCCGAGGGGCTGACCGTGGAGGTCAAGCCCGGGACCAACATCCTgg                                                                                <  1:459682/136‑1 (MQ=255)                                                    ccAAAGCCGCACTGGCCGGGCGCAACCTGAATGAGATCGAGGAGTAATCGATGCCGCTGGTGACATTCCTGCCGCATGAGAAGTTCTGCCCCGAGGGGCTGACCGTGGAGGTCAAGCCCGGGACCAACATCCTGGa                                                                               <  1:465090/135‑1 (MQ=255)                                                            gCACTGGCCGGGCGCAACCTGAATGAGATCGAGGAGTAATCGATGCCGCTGGTGACATTCCTGCCGCATGAGAAGTTCTGCCCCGAGGGGCTGACCGTGGAGGTCAAGCCCGGGACCAACATCCTGGAGCTGGccc                                                                       <  1:50212/136‑1 (MQ=255)                                                                tGGCCGGGCGCAACCTGAATGAGATCGAGGAGTAATCGATGCCGCTGGTGACATTCCTGCCGCATGAGAAGCTCTGCCCCGAGGGGCTGACCGTGGAGGTCAAGCCCGGGACCAACATCCTGGAGCTGGCCCACGa                                                                   >  2:502479/1‑136 (MQ=255)                                                                tGGCCGGGCGCAACCTGAATGAGATCGAGGAGTAATCGATGCCGCTGGTGACATTCCTGCCGCATGAGAAGCTCTGCCCCGAGGGGCTGACCGTGGAGGTCAAGCCCGGGACCAACATCCTGGAGCTGGCCCACGa                                                                   >  2:504511/1‑136 (MQ=255)                                                                     gggCGCAACCTGAATGATATCGAGGAGTAATCGATGCCGCTGGTGACATTCCTGCCGCATGAGAAGTTCTGCCCCGAGGGGCTGACCGTGGAGGTCAAGCCCGGTACCAACATCCTGGAGCTGGCCCACGACCATc                                                              >  2:117247/1‑136 (MQ=255)                                                                     gggCGCAACCTGAATGAGATCGAGGAGTAATCGATGCCGCTGGTGACATTCCTGCCGCATGAGAAGTTCTGCCCCGAGGGGCTGACCGTGGAGGTCAAGCCCGGGACCAACATCCTGGAGCTGGCCCACGACCATc                                                              >  1:390538/1‑136 (MQ=255)                                                                       gcgcAACCTGAATGAGATCGAGGAGTAATCGATGCCGCTGGTGACATTCCTGCCGCATGAGAAGCTCTGCCCCGAGGGGCTGACCGTGGAGGTCAAGCCCGGGACCAACATCCTGGAGCTGGCCCACGACCATcac                                                            <  1:502479/136‑1 (MQ=255)                                                                       gcgcAACCTGAATGAGATCGAGGAGTAATCGATGCCGCTGGTGACATTCCTGCCGCATGAGAAGCTCTGCCCCGAGGGGCTGACCGTGGAGGTCAAGCCCGGGACCAACATCCTGGAGCTGGCCCACGACCATcac                                                            <  1:504511/136‑1 (MQ=255)                                                                                                  ttCGATGCCGCTGGTGACAATCCTGCCGCATTAGAAGTTATGCCCCGAGGGGCTGACCGTGGAGGTCAAGCCCGTGACCAACATCCTGGAGCTGGCCCACGACCATCACATCGAGATTGAAAGCGCCTgcggcggc                                 <  1:300738/135‑1 (MQ=255)                                                                                                     gATGCCGCTGGTGACATTCCTGCCGCATGAGAAGTTCTGCCCCGAGGGGCTGACCGTGGAGGTCAAGCCCGGGACCAACATCCTGGAGCTGGCCCACGACCATCACATCGAGATGGAAAGCGCCTGCGGCGGCGTc                              <  1:387444/136‑1 (MQ=255)                                                                                                           gCTGGTGACATTCCTGCCGCATGAGAAGTTCTGCCCCGAGGGGCTGACCGTGGAGGTCAAGCCCGGGACCAACATCCTGGAGCTGGCCCACGACCATCACATCGAGATGGAAAGCGCCTGCGGCGGCGTCAAGGcc                        >  1:476481/1‑136 (MQ=255)                                                                                                           gCTGGTGACATTCCTGCCGCATGAGAAGTTCTGCCCCGAGGGGCTGACCGTGGAGGTCAAGCCCGGGACCAACATCCTGGAGCTGGCCCACGACCATCACATCGAGATGGAAAGCGCCTGCGGCGGCGTCAAGGcc                        >  2:490987/1‑136 (MQ=255)                                                                                                           gCTGGTGACATTCCTGCCGCATGAGAAGTTCTGCCCCGAGGGGCTGACCGTGGAGGTCAAGCCCGGGACCAACATCCTGGAGCTGGCCCACGACCATCACATCGAGATGGAAAGCGCCTGCGGCGGCGTCAAGGcc                        >  1:213755/1‑136 (MQ=255)                                                                                                             tGGTGACATTCCTGCCGCATGAGAAGTTCTGCCCCGAGGGGCTGACCGTGGAGGTCAAGCCCGGGACCAACATCCTGGAGCTGGCCCACGACCATCACATCGAGATGGAAAGCGCCTGCGGCGGCGTCAAGGCCTg                      >  2:465696/1‑136 (MQ=255)                                                                                                                 ttcATTCGTGGCGCATGAGAAGTTCTGCCCCGAGGGGCTGACGGTGGAGGGCAAGCCCGGGACCTACATCCTGGAGCTGGCCGACGACCATCACATCGAGATGGAAAGCGCCTGCGGCGGGGTCAAGGCCTGTAcc                  <  1:171416/134‑1 (MQ=255)                                                                                                                 ttcATTCCTGCCGCATGAGAAGTTCTGCCCCGAGGGGCTGACCGTGGAGGTCAAGCCCGTGACCAACATCCTGGAGCTGGCCCACGACCATCACATCGAGATGGAAAGCGCCTGCGGCGGCGTCAAGGCCTGTAcc                  <  1:306069/134‑1 (MQ=255)                                                                                                                 gACATTCCTGCCGCATGAGAAGTTCTGCCCCGAGGGGCTGACCGTGGAGGTCAAGCCCGGGACCAACATCCTGGAGCTGGCCCACGACCATCACATCGAGATGGAAAGCGCCTGCGGCGGCGTCAAGGCCTGTAcc                  <  1:287430/136‑1 (MQ=255)                                                                                                                   cATTCCTGCCGCATGAGAAGTTCTGCCCCGAGGGGCTGACCGTGGAGGTCAAGCCCGGTACCAACATCCTGGAGCTGGCCCACGACCATCACATCGAGATGGAAAGCGCCTGCGGCGGCGTCAAGGCCTGTACCAc                <  1:117247/136‑1 (MQ=255)                                                                                                                                 gagaAGTTCTGCCCCGAGGGGCTGACCGTGGAGGTCAAGCCCGGGACCAACATCCTGGAGCTGGCCCACGACCATCACATCGAGATGGAAAGCGCCTGCGGCGGCGTCAAGGCCTGTACCACTTGCCACTGCATCg  >  1:391900/1‑136 (MQ=255)                                                                                                                                       |                                                                                                                                  GTCTGTCGCAGGTGACCGATGCATTTGCCGCCCGTCGCCTTGATTCGACGGTCAAAGCCGCACTGGCCGGGCGCAACCTGAATGAGATCGAGGAGTAATCGATGCCGCTGGTGACATTCCTGCCGCATGAGAAGTTCTGCCCCGAGGGGCTGACCGTGGAGGTCAAGCCCGGGACCAACATCCTGGAGCTGGCCCACGACCATCACATCGAGATGGAAAGCGCCTGCGGCGGCGTCAAGGCCTGTACCACTTGCCACTGCATCG  >  NC_002947/982892‑983155

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