Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_012967	3,927,576	Δ1 bp	coding (16/330 nt)	trxA →	thioredoxin TrxA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_012967	3,927,576	0	A	.	100.0%	146.3 / NA	29	coding (16/330 nt)	trxA	thioredoxin TrxA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (17/12);  total (17/12)

GTTAATGTTACACCAACAACGAAACCAACACGCCAGGCTTATTCCTGTGGAGTTATATATGAGCGATAAAATTATTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCGGACGGGGCGATCCTCGTCGATTTCTGGG  >  NC_012967/3927503‑3927648                                                                          |                                                                         gTTAATGTTACACCAACAACGAAACCAACACGCCAGGCTTATTCCTGTGGAGTTATATATGAGCGATAAAATT‑TTCACctga                                                                 >  1:410167/1‑82 (MQ=255)      tGTTACACCAACAACGAAACCAACACGCCAGGCTTATTCCTGTGGAGTTATATATGAGCGATAAAATT‑TTCACCtgactgac                                                            >  1:2711717/1‑82 (MQ=255)        ttACACCAACAACGAAACCAACACGCCAGGCTTATTCCTGTGGAGTTATATATGAGCGATAAAATT‑TTCACCTGACTgacga                                                          >  1:1430443/1‑82 (MQ=255)        ttACACCAACAACGAAACCAACACGCCAGGCTTATTCCTGTGGAGTTATATATGAGCGATAAAATT‑TTCACCTGACTgacga                                                          >  1:2667076/1‑82 (MQ=255)         tACACCAACAACGAAACCAACACGCCAGGCTTATTCCTGTGGAGTTATATATGAGCGATAAAATT‑TTCACCTGACTgacgac                                                         <  1:4092547/82‑1 (MQ=255)           cacCAACAACGAAACCAACACGCCAGGCTTATTCCTGTGGAGTTATATATGAGCGATAAAATT‑TTCACCTGACTGACGACAg                                                       <  1:2595821/82‑1 (MQ=255)               aacaacGAAACCAACACGCCAGGCTTATTCCTGTGGAGTTATATATGAGCGATAAAATT‑TTCACCTGACTGACGACAGtttt                                                   >  1:85184/1‑82 (MQ=255)                    cGAAACCAACACGCCAGGCTTATTCCTGTGGAGTTATATATGAGCGATAAAATT‑TTCACCTGACTGACGACAGTTTTGacac                                              >  1:3647712/1‑82 (MQ=255)                        aCCAACACGCCAGGCTTATTCTTGTGGAGTTATATATGAGCGATAAAATT‑TTCACCTGACTGACGACAGTTTTGACACGGAt                                          >  1:3415443/1‑82 (MQ=255)                        aCCAACACGCCAGGCTTATTCCTGTGGAGTTATATATGAGCGATAAAATT‑TTCACCTGACTGACGACAGTTTTGACACGGAt                                          <  1:2335633/82‑1 (MQ=255)                         ccAACACGCCAGGCTTATTCCTGTGGAGTTATATATGAGCGATAAAATT‑TTCACCTGACTGACGACAGTTTTGACACGGATg                                         <  1:1337350/82‑1 (MQ=255)                         ccAACACGCCAGGCTTATTCCTGTGGAGTTATATATGAGCGATAAAATT‑TTCACCTGACTGACGACAGTTTTGACACGGATg                                         <  1:3094865/82‑1 (MQ=255)                           aaCACGCCAGGCTTATTCCTGTGGAGTTATATATGAGCGATAAAATT‑TTCACCTGACTGACGACAGTTTTGACACGGATGTa                                       >  1:154997/1‑82 (MQ=255)                               cGCCAGGCTTATTCCTGTGGAGTTATATATGAGCGATAAAATT‑TTCACCTGACTGACGACAGTTTTGACACGGATGTACTCa                                   <  1:3398537/82‑1 (MQ=255)                                   aGGCTTATTCCTGTGGAGTTATATATGAGCGATAAAATT‑TTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGc                               >  1:1672714/1‑82 (MQ=255)                                   aGGCTTATTCCTGTGGAGTTATATATGAGCGATAAAATT‑TTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGc                               <  1:2490399/82‑1 (MQ=255)                                   aGGCTTATTCCTGTGGAGTTATATATGAGCGATAAAATT‑TTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGc                               <  1:4302522/82‑1 (MQ=255)                                        tATTCCTGTGGCGTTATATATGAGCGATAAAATT‑TTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCGGACg                          >  1:743333/1‑82 (MQ=255)                                        tATTCCTGTGGAGTTATATATGAGCGATAAAATT‑TTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCGGACg                          >  1:1139156/1‑82 (MQ=255)                                        tATTCCTGTGGAGTTATATATGAGCGATAAAATT‑TTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCGGACg                          >  1:3669408/1‑82 (MQ=255)                                              tgtgGAGTTATATATGAGCGATAAAATT‑TTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCGGACGGGGCGa                    <  1:143084/82‑1 (MQ=255)                                                tgGAGTTATATATGAGCGATAAAATT‑TTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCGGACGGGGCGATc                  >  1:1768983/1‑82 (MQ=255)                                                       atatatGAGCGATAAAATT‑TTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCGGACGGGGCGATCCtcgtcg           >  1:2541347/1‑82 (MQ=255)                                                          tatGAGCGATAAAATT‑TTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCGGACGGGGCGATCCTCGTCGAtt        <  1:77138/82‑1 (MQ=255)                                                            tGAGCGATAAAATT‑TTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCGGACGGGGCGATCCTCGTCGATTTc      >  1:3800067/1‑82 (MQ=255)                                                             gAGCGATAAAATT‑TTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCGGACGGGGCGATCCTCGTCGATTTCt     >  1:1142699/1‑82 (MQ=255)                                                             gAGCGATAAAATT‑TTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCGGACGGGGCGATCCTCGTCGATTTCt     <  1:4425578/82‑1 (MQ=255)                                                             gAGCGATAAAATT‑TTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCGGACGGGGCGATCCTCGTCGATTTCt     >  1:775107/1‑82 (MQ=255)                                                                cGATAAAATT‑TTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCGGACGGGGCGATCCTCGTCGATTTCTggg  <  1:1541603/82‑1 (MQ=255)                                                                          |                                                                         GTTAATGTTACACCAACAACGAAACCAACACGCCAGGCTTATTCCTGTGGAGTTATATATGAGCGATAAAATTATTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCGGACGGGGCGATCCTCGTCGATTTCTGGG  >  NC_012967/3927503‑3927648

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑