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%	113.1 / NA	35	E232G (GAA→GGA)	torC	pentaheme c‑type cytochrome TorC
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (26/9);  total (26/9)

GACACGCTGTACTCTATTGATATTAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGAAGTGAAAGTCCTTAAACGTGACGGCGACTGGCTGCAAATTGAAATTACCGGCTGGACGGAAAGCGCCGGACGTCAGCGTGTACTCA  >  NZ_CP009273/1054149‑1054320                                                                                      |                                                                                       gACACGCTGTACTCTATTGATATTAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTg                                                                                    >  1:1206585/1‑90 (MQ=255)     cGCTGTACTCTATTGATATTAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAg                                                                                <  2:167751/90‑1 (MQ=255)      gCTGTACTCTATTGATATTAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGt                                                                               >  1:232167/1‑90 (MQ=255)      gCTGTACTCTATTGATATTAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGt                                                                               >  2:564462/1‑90 (MQ=255)       cTGTACTCTATTGATATTAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTTCCTGCTTCGGGAGTGAAAGTc                                                                              >  2:471060/1‑90 (MQ=255)              ctATTGATATTAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGTCCTTAAAc                                                                       >  1:1424299/1‑90 (MQ=255)                aTTGATATTAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGTCCTTAAACGt                                                                     >  1:838173/1‑90 (MQ=255)                   gATATTAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGAc                                                                  >  2:514673/1‑90 (MQ=255)                     tatTAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGACgg                                                                >  2:1081752/1‑90 (MQ=255)                     tatTAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGACgg                                                                >  2:134296/1‑90 (MQ=255)                     tatTAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGACgg                                                                >  2:27340/1‑90 (MQ=255)                        tAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGACGGCGa                                                             >  2:415181/1‑90 (MQ=255)                        tAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGACGGCGa                                                             >  1:655631/1‑90 (MQ=255)                        tAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGACGGCGa                                                             >  1:689664/1‑90 (MQ=255)                        tAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGACGGCGa                                                             >  2:832095/1‑90 (MQ=255)                        tAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGACGGCGa                                                             >  1:1028270/1‑90 (MQ=255)                          aGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGACGGCGACt                                                           >  2:614365/1‑90 (MQ=255)                          aGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGACGGCGACt                                                           >  2:989909/1‑90 (MQ=255)                          aGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGACGGCGACt                                                           >  1:281023/1‑90 (MQ=255)                               aTTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGACGGCGACTGGCTg                                                      >  1:91116/1‑90 (MQ=255)                                   aTGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGACGGCGACTGGCTGCaaa                                                  <  1:436597/90‑1 (MQ=255)                                         cgAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGACGGCGACTGGCTGCAAATTGaaa                                            >  1:409860/1‑90 (MQ=255)                                                  aTAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGACGGCGACTGGCTGCAAATTGAAATTACCGGCt                                   >  2:1249240/1‑90 (MQ=255)                                                  aTAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGACGGCGACTGGCTGCAAATTGAAATTACCGGCt                                   >  2:1501914/1‑90 (MQ=255)                                                  aTAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGACGGCGACTGGCTGCAAATTGAAATTACCGGCt                                   <  2:1541388/90‑1 (MQ=255)                                                  aTAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGACGGCGACTGGCTGCAAATTGAAATTACCGGCt                                   >  1:806564/1‑90 (MQ=255)                                                                   tctcTGCTGCCTGCTTCGGGAGTGAAAGACCTTAAACGTGACGGCGACTGGCTGCAAATTGAAATTACCGGCTGGACGGAAAGCCCCggg                  >  2:1430166/1‑89 (MQ=255)                                                                      ctgctgCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGACGGCGACTGGCTGCAAATTGAAATTACCGGCTGGACGGAAAGCGCCGGACGt               >  2:1190643/1‑90 (MQ=255)                                                                      ctgctgCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGACGGCGACTGGCTGCAAATTGAAATTACCGGCTGGACGGAAAGCGCCGGACGt               <  2:1415612/90‑1 (MQ=255)                                                                      ctgctgCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGACGGCGACTGGCTGCAAATTGAAATTACCGGCTGGACGGAAAGCGCCGGACGt               <  1:832095/90‑1 (MQ=255)                                                                       tgctgcCTGCTTCGGGAGTGAAAGTCCTTAAACGTGACGGCGACTGGCTGCAAATTGAAATTACCGGCTGGACGGAAAGCGCCGGACGTc              <  1:102051/90‑1 (MQ=255)                                                                           gcctgcTTCGGGAGTGAAAGTCCTTAAACGTGACGGCGACTGGCTGCAAATTGAAATTACCGGCTGGACGGAAAGCGCCGGACGTCAGCg          >  1:1386133/1‑90 (MQ=255)                                                                                cTTCGGGAGTGAAAGTCCTTAAACGTGACGGCGACTGGCTGCAAATTGAAATTACCGGCTGGACGGAAAGCGCCGGACGTCAGCGTGTAc     <  1:27340/90‑1 (MQ=255)                                                                                 ttCGGGAGTGAAAGTCCTTAAACGTGACGGCGACTGGCTGCAAATTGAAATTACCGGCTGGACGGAAAGCGCCGGACGTCAGCGTGTACt    <  1:292244/90‑1 (MQ=255)                                                                                   cGGGAGTGAAAGTCCTTAAACGTGACGGCGACTGGCTGCAAATTGAAATTACCGGCTGGACGGAAAGCGCCGGACGTCAGCGTGTACTCa  <  1:1430166/90‑1 (MQ=255)                                                                                      |                                                                                       GACACGCTGTACTCTATTGATATTAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGAAGTGAAAGTCCTTAAACGTGACGGCGACTGGCTGCAAATTGAAATTACCGGCTGGACGGAAAGCGCCGGACGTCAGCGTGTACTCA  >  NZ_CP009273/1054149‑1054320

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