Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	1,054,234	A→G	E232G (GAA→GGA)	torC →	pentaheme c‑type cytochrome TorC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	1,054,234	0	A	G	100.0%	96.8 / NA	29	E232G (GAA→GGA)	torC	pentaheme c‑type cytochrome TorC
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (16/13);  total (16/13)

CACGCTGTACTCTATTGATATTAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGAAGTGAAAGTCCTTAAACGTGACGGCGACTGGCTGCAAATTGAAATTACCGGCTGGACGGAAAGCGCCGGACGTCAGCG  >  NZ_CP009273/1054151‑1054312                                                                                    |                                                                               cacGTTGTACTCTATTGATATTAAGCCAGTTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGaa                                                                          <  1:362867/90‑1 (MQ=255) cacGCTGTACTCTATTGATATTAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGaa                                                                          <  1:725836/90‑1 (MQ=255)    gCTGTACTCTATTGATATTAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGt                                                                       >  1:727626/1‑90 (MQ=255)            ctATTGATATTAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGTCCTTAAAc                                                               >  2:154530/1‑90 (MQ=255)              aTTGATATTAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGTCCTTAAACGt                                                             >  1:714321/1‑90 (MQ=255)                 gATATTAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGa                                                           >  2:249854/1‑89 (MQ=255)                 gATATTAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGa                                                           <  1:249854/89‑1 (MQ=255)                 gATATTAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGAc                                                          >  2:99109/1‑90 (MQ=255)                 gATATTAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGAc                                                          >  2:527783/1‑90 (MQ=255)                 gATATTAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGAc                                                          >  1:373805/1‑90 (MQ=255)                 gATATTAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGAc                                                          >  1:232916/1‑90 (MQ=255)                   tatTAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGACgg                                                        >  1:287178/1‑90 (MQ=255)                   tatTAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGACgg                                                        >  2:938328/1‑90 (MQ=255)                      tAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGACGGCGa                                                     >  1:598910/1‑90 (MQ=255)                      tAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGACGGCGa                                                     >  1:290570/1‑90 (MQ=255)                      tAAGCCAATTTATGCGGCGAAAGGCGATAAAAAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGGGAAAGTCCTTAAACGTGACGGCGa                                                     >  1:983398/1‑90 (MQ=255)                        aGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGACGGCGACt                                                   >  2:244383/1‑90 (MQ=255)                            aaTTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGACGGCGACTGGCt                                               >  1:315925/1‑90 (MQ=255)                            aaTTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGACGGCGACTGGCt                                               >  1:369834/1‑90 (MQ=255)                             aTTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGACGGCGACTGGCTg                                              <  2:650842/90‑1 (MQ=255)                                                        gCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGACGGCGACTGGCTGCAAATTGAAATTACCGGCTGGACGGaa                   <  1:99109/90‑1 (MQ=255)                                                             tGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGACGGCGACTGGCTGCAAATTGAAATTACCGGCTGGACGGAAAgcgc              <  1:938328/90‑1 (MQ=255)                                                             tGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGACGGCGACTGGCTGCAAATTGAAATTACCGGCTGGACGGAAAgcgc              <  2:287178/90‑1 (MQ=255)                                                                    ctgctgCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGACGGCGACTGGCTGCAAATTGAAATTACCGGCTGGACGGAAAGCGCCGGACGt       <  1:155619/90‑1 (MQ=255)                                                                    ctgctgCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGACGGCGACTGGCTGCAAATTGAAATTACCGGCTGGACGGAAAGCGCCGGACGt       <  1:230476/90‑1 (MQ=255)                                                                    ctgctgCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGACGGCGACTGGCTGCAAATTGAAATTACCGGCTGGACGGAAAGCGCCGGACGt       <  1:244383/90‑1 (MQ=255)                                                                     tgctgcCTGCTTCGGGAGTGAAAGTCCTTAAACGTGACGGCGACTGGCTGCAAATTGAAATTACCGGCTGGACGGAAAGCGCCGGACGTc      <  2:764488/90‑1 (MQ=255)                                                                      gctgcCTGCTTCGGGAGTGAAAGTCCTTAAACGTGACGGCGACTGGCTGCAAATTGAAATTACCGGCTGGACGGAAAGCGCCGGACGTCa     <  2:899731/90‑1 (MQ=255)                                                                         gcctgcTTCGGGAGTGAAAGTCCTTAAACGTGACGGCGACTGGCTGCAAATTGAAATTACCGGCTGGACGGAAAGCGCCGGACGTCAGCg  <  2:112295/90‑1 (MQ=255)                                                                                    |                                                                               CACGCTGTACTCTATTGATATTAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGAAGTGAAAGTCCTTAAACGTGACGGCGACTGGCTGCAAATTGAAATTACCGGCTGGACGGAAAGCGCCGGACGTCAGCG  >  NZ_CP009273/1054151‑1054312

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