Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_007779	2,606,882	T→G	23.6%	V437G (GTA→GGA)	hyfF →	hydrogenase 4 subunit F

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007779	2,606,882	0	T	G	23.6%	57.2 / 4.2	48	V437G (GTA→GGA)	hyfF	hydrogenase 4 subunit F
Reads supporting (aligned to +/- strand):  ref base T (19/17);  new base G (10/2);  total (29/19)
Fisher's exact test for biased strand distribution p-value = 9.05e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-01

CCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAAC‑GCTGGTGCTGGCGGGCCTGGTA‑CGGATGGCTGCGCGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCA‑AT‑GGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTG  >  NC_007779/2606645‑2607095                                                                                                                                                                                                                                               |                                                                                                                                                                                                                         ccGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAAC‑GCTGGTGCTGGCGGGGCTGGTa                                                                                                                                                                                                                         >  1:1203703/1‑239 (MQ=255)         gtccTGCTCAGTTACGGCCGGGGCGATCTCATGGTGGTGGGTGGGCTTGTCAAAATCAGCCCATTTACCGCCGTGCGGTTTGGCGGTGGTGCGCTGGCGCTGGCAGGGAGGCCGCCGTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTACTGTTAAC‑GCTGGTGCTGGCGGGCCTGGTA‑CGGAt                                                                                                                                                                                                                   <  2:589318/233‑1 (MQ=255)         gTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGCATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAAC‑GCTGGTGCTGGCGGGACTGGTA‑CGGATggg                                                                                                                                                                                                                >  1:1517292/1‑238 (MQ=255)             tGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTGCCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAAC‑GCTGGTGCTGGCGGGCCTGGTA‑CGGATGGCTgcg                                                                                                                                                                                                            <  1:325776/239‑1 (MQ=255)             ggataaaGTCCGGCACGCGGGATCTCAACGTCGTCTGTGGGATGCCAAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAAC‑GCTGGTGCTGGCGGGCCTGGTA‑CGGATGGCTgcg                                                                                                                                                                                                            <  2:406146/234‑1 (MQ=255)              gtccaatgccGGCACGCGCGATCTCACCGTCGTCTGTGGGAGGCTCAAAATCGTGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAAC‑GCTGGTGCTGGCGGGCCTGGTA‑CGGATGGCTgcg                                                                                                                                                                                                            <  2:159321/229‑1 (MQ=255)              gCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAAC‑GCTGGTGCTGGCGGGGCTGGTA‑CCGGTGGCTgggc                                                                                                                                                                                                           >  1:2266437/1‑239 (MQ=255)               cTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGAGCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAAC‑GCTGGTGCTGGCGGGGCAGGGACCGGATGGCTccgc                                                                                                                                                                                                           >  1:9166/1‑239 (MQ=255)               cTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAAC‑GCTGGTGCTGGCGGGGCTGGTA‑CGGGTGGCTGcgcg                                                                                                                                                                                                          >  1:2103411/1‑239 (MQ=255)               cTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAAC‑GCTGGTGCTGGCGGGCCTGGTA‑CGGGTGGCTGcgcg                                                                                                                                                                                                          >  1:823515/1‑239 (MQ=255)                    gTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACCTTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAAC‑GCTGGTGCTGGCGGGCCCGGTA‑Cggg                                                                                                                                                                                                                    >  1:577133/1‑223 (MQ=255)                    gTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAAC‑GCTGGTGCTGGCGGGCCTGGGA‑AGGATGGCTGCGCgggggg                                                                                                                                                                                                     >  1:794172/1‑236 (MQ=255)                       cGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCCCGTAATCACCTGCTTATTATCGTCCTGCTGTTATTGCTGTTAAC‑GCTGGTGCTGGCGGGCCCGGTT‑CGGATGGCTGCGCGGGGGTTat                                                                                                                                                                                                  >  1:1035337/1‑238 (MQ=255)                         gCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAAC‑GCTGGTGCTGGCGGGCCTGGTA‑CGGATGGGTGCGCggggggttatg                                                                                                                                                                                                >  1:360154/1‑231 (MQ=255)                         gCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAAC‑GCTGGTGCTGGCGGGCCTGGGA‑AGGGTGGCTGGGGGGGGGGTaa                                                                                                                                                                                                  >  1:159321/1‑237 (MQ=255)                              cgcgATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAAC‑GCTGGTGCTGGCGGGCCCTGTA‑CGGGTGGCTTCGCGGGGGGTTATggggga                                                                                                                                                                                           >  1:1573646/1‑235 (MQ=255)                                             gtcTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCCTCAACATTTTTCTTAGCGAATTTATTACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAAC‑GCTGGTGCTGGCGGGGCTGGTACCGGGGGGCGGGCGGGGGGTTAAGGGGGAAACCCCCCCAGGccc                                                                                                                                                                             >  2:1641909/1‑238 (MQ=255)                                                                   tGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAAC‑GCTGGTGCTGGCGGGCCTGGGA‑AGGATGGGTTCGCGGGGGTTTATGGCGAAACCGCCGCCGGCCGGTAACCGGGGTGGTTGCGGCTgg                                                                                                                                                      >  1:460270/1‑239 (MQ=255)                                                                        tttACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAAC‑GCTGGTGCTGGCGGGCCTGGGA‑CGGATGGCTGCGCGGGGGTTAATGGCGAAACCGCCCCCGGCCGTTCACCGGGGGGGTCTTGGCCggggga                                                                                                                                                  >  1:1688654/1‑234 (MQ=255)                                                                          tACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAAC‑GCTGGTGCTGGCGGGCCTGGGA‑CGGGTGGCTGCGGGGGGGTTAATGGCGAAACCGCCGCCGGCCGGTTACCGGGGGGAGCTCGggtgggtggac                                                                                                                                                >  1:545388/1‑235 (MQ=255)                                                                               ccgTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAAC‑GCTGGTGCTGGCGGGCCCGGTA‑CGGATGGCTGGGCGGGTGGTTATGGCGAAACCGCCGCACGGCGTTAACCCGGGGGAGCTCGGCTGGTTGaccccccc                                                                                                                                           >  1:531999/1‑236 (MQ=255)                                                                                       tttGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAAC‑GCTGGTGCTGGCGGGGCTGGGA‑CGGATGGCTGGGCGGGTGGTAATGGCGGAACCGCCGCCGGGCCGTAAACGGGGGGGTGTCGGGTGGGTGGACACCCCGCCC‑AT‑gg                                                                                                                                 >  1:208407/1‑238 (MQ=255)                                                                                                                           cgagtgaccgtttttCTTGGCTATTTTATGCCCTTGCCCGCGGGCTGGGCAGGTCTTCCCCTGTGGTTTATCGTCCTGCTGTTTTTGCTGTTAAG‑GTGGTTGCTGGGGGGCCGGTTG‑CGGAGGGCGGCGCGGTTGTGAATGGCAAACCCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCA‑AT‑GGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCa                                                                                            <  2:1963302/229‑1 (MQ=255)                                                                                                                           cccTTCATATTTTTTCTTGGCGTTTTTATGACCTTTCCCGCGGGGTGGGCCCGTATTCCCCTGCTGATTTTCGTCCTGTTGTTGTTGTTGTTAAC‑GGGGGTGCTGGCGGCCCTGGTG‑GGGAGGGCTGCGGGGGTGTTATGGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCA‑AT‑GGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCa                                                                                            <  2:431732/239‑1 (MQ=255)                                                                                                                             cTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAAC‑GCTGGTGCTGGCGGGCCTGGGA‑CGGATGGCTGCGCGGGTGTTAATGGGGCAACCGCCGCCGGGCGTTAACCGGGGGGGTCTCGGCTGGGTTGCCACCCCCCCC‑CT‑GGGGAGTCTGTCGGGCATGGTGCTGGGGAGGGGGACGCaca                                                                                          >  2:833955/1‑237 (MQ=255)                                                                                                                                          cTTAGCGAATTTATGACCTTTACCGCGGGCCGGGCGCGTATTCACCTGTTGATTATCGTCCGGTTGTTATTGCTGTTAAC‑TCTGGTGCTGGCGGGCCTGGTA‑CGGATGGCTGCGCGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCA‑AT‑GGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGt                                                                             <  1:486081/239‑1 (MQ=255)                                                                                                                                                            aTTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAAC‑GCTGGTGCTGGCGGGCCTGGTA‑CGGATGGCTGCGCGGGTGGTAATGGCGAAAACGCCGCCGGCCGGTAAACGGGGGGATCTCGGGGGGGTGACCCCCCCCCCC‑AA‑CGGGGTTTTGGCGGGCATGGTTCTTGCGATGGGGACCGATATTACCCAAACCGTCATCAAGACACTgggggg                                                           >  2:1471804/1‑239 (MQ=255)                                                                                                                                                              tACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAAC‑GCTGGTGCTGGCGGGCCTGGTA‑CGGATGGCTGCGCGGGTGTTAATGGCGAAACCGCCGCCGGCCGTTAACCGGGGGGGTCGCGGGTGGGTGGGCACCACGCCC‑AC‑GGGGATTCTGGCGGGCATGATGGCGGGGGGGGGGACCCACATATCCCAAACTGTCCTCAGGATCCCggcgggcg                                                         >  2:54818/1‑239 (MQ=255)                                                                                                                                                              tACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAAC‑GCTGGTGCTGGCGGGCCTCGTT‑CGGATGGCTTCGCGGGTGTTAATGGCGAAAACGCCGCAGGCCGGTAACCGGGGGGATCTCGGCTGGGTGACCACCCACCCA‑ATGGGTGATCTTGC‑GGGCATGATGCTGGCGGTGGGGACGCAAATTTCAAAACCTTTCCCCAGGAGCACggcggggg                                                         >  2:2044672/1‑237 (MQ=255)                                                                                                                                                                  gccgGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAAC‑GCTGGTGCTGGCGGGCCTGGTA‑CGGATGGCTGCGCGGGTGTTAATGGCGAAACCGCCGCAGGCCGGTAACCGGGGGGGTCTCGGCTGGGTGACCACCCCCCCC‑AT‑GGTGGATCTGCTGGGCATGATGCTGGGGGGGGGGAAGCCTATTCCCCAACCTGTCCTCCGGGTCCTgggccggtctct                                                     >  2:547853/1‑229 (MQ=255)                                                                                                                                                                    cgGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAAC‑GCTGGTGCTGGCGGGCCTGGTA‑CGGGTGGCTGCGCGGGGGTTAATGGCGGAAACGCCGCAGGGCGCTAACCGGGGGGGGCTCGGCTGGGTGGCCACCCCGCCC‑AA‑GGGGGTTCTGCTGGGCATGATGCTGGGGGTGGGGACACATATTCCACAACCCGGCATCAGGATCCTCGCGGgggccttc                                                    >  2:757974/1‑234 (MQ=255)                                                                                                                                                                      gACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAAC‑GCTGGTGCTGGCGGGCCTGGTA‑CGGGTGGCTGCGCGGGGGTTTATGGCGAAACCGCCCCCGGCCGTTGACCGGGGGGGGCCCGGCTGGGTGACCCCCCCCCCC‑AC‑GGGGGATCTGCTGGGGCTGATTCTGGCGGTGGGGAACCCTCTTCCACAACCCTTCCTCAAGGTCCTcgcggggcggtccc                                                   >  2:1608897/1‑229 (MQ=255)                                                                                                                                                                      gACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAAC‑GCTGGTGCTGGCGGGCCTGGGA‑AGGAAGGCTGCGCGGGGGGTAATGGCGCAACCGCCGCAGGCCGTTAAACGGGGGGATCTCGGCTGGGTTCCCAACCCCCCACAT‑GG‑GTGTCTGCTGGGCATGATGCTGGTGTAGGGGAAGCATAATTCACAACCTGTACTCAGGAGTCTGGGGGGGGCCTCCACt                                                 >  2:1608907/1‑239 (MQ=255)                                                                                                                                                                           gCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAAC‑GCTGGTGCTGGCGGGCCTGGGA‑CGGGTGGCTGCGCGGGGGGTAATGGCGGAACCGCCGCCGGGCGCTAACCCGGGGGATCTCGGCTGGGTGGCCCCCCCCCCC‑AA‑GGTGGATCTTCCGGGTCTGATGGTGGGGGGGGGGACGCCTATTTCACAACACCTCCTCACGAACCTGGCGGGGGCGTCCCCcataat                                            >  2:1369248/1‑237 (MQ=255)                                                                                                                                                                            cgcttatcccccTGCTGATTTTCGTCTTG‑TGTTTTTGGGGTTTACGGGGGGTGTGGGCGGGCCTGGTC‑GGGAGGGCGGCGCGGGTGTAAAGGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCA‑AT‑GGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGt                                            <  2:1573646/228‑1 (MQ=255)                                                                                                                                                                            cACGTAATCCCGTGCTGTTTATCGTCTTGCTGTTATTGCTGTTAGC‑GCGGGTGCTGGCGGGCGTGGTG‑CGGATGGGTGGGCGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCA‑AT‑GGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGt                                            <  1:948575/238‑1 (MQ=255)                                                                                                                                                                            cACGTAACCGCCTGCTGATTATCGGCCTGCTTTTTTTGCTGTTACC‑GCGGGTGGGGGGGGGCCTGGGA‑GGGAGGGCTGCGCGGGTTTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGaccac                                                                                                                                              <  1:1750191/142‑1 (MQ=255)                                                                                                                                                                                          ctgTTTATCGTCCTGTTTTTTTGGCTGTTAAC‑GGTGGTGGTGGCGGGCCTGGTA‑CGGATGGCGGCGCGGGTTTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCA‑AT‑GGTGATTCTGCTGGTCATGATGCTGGCGATggg                                                                                                  <  1:431491/170‑1 (MQ=255)                                                                                                                                                                                             atatgaaacgggcttTTATTGTTGTTAAC‑GCTGGTGCTGGCGGCCGGGGTG‑CGGATGCCTGCGGGGTTGTTAAGGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCA‑AT‑GGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGa                          <  1:1550082/228‑1 (MQ=255)                                                                                                                                                                                                  cgcccTGCTGTTTTTGCTGTTAAC‑GCTGGTGCGGGCGGGCCTGGTC‑GGGAGGGCGGGGCGGGTGTTAAGGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCA‑AT‑GGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATctgc                     <  1:461512/236‑1 (MQ=255)                                                                                                                                                                                                   gTCCTGGTGTTTTGGCGGTTAGC‑GGGGGGGCTGGCGGGCCGGGGA‑GGGGGGGTTGGGCGGGGGTTAATGGGGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCA‑AT‑GGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATctgc                     <  1:1221911/238‑1 (MQ=255)                                                                                                                                                                                                       tgctgTTCTTGCTGTTACC‑GCTGGTGCTGGCGGCCCGGGTC‑GGGATGGCGGCGCGGTTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCA‑AT‑GGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGca                <  1:2106930/239‑1 (MQ=255)                                                                                                                                                                                                        gctgTTGTTGTGGTTACC‑GTGGGTGGGGGCGGGCTTGGTG‑GGGATGGCGGCGCGGGTGTGAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCA‑AT‑GGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGca                <  2:753619/238‑1 (MQ=255)                                                                                                                                                                                                          tgTTTTTGCTGTTAAC‑GCTGGTGCGGGCGGGCCTGGTA‑CGGATGGCTGCGCGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCA‑AT‑GGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACta                                                <  1:1333999/204‑1 (MQ=255)                                                                                                                                                                                                           tttATTGCGGTTAAC‑GCTGGGGCGGGCGGGCCGGGTC‑CGGATGGCTGCGTGGGTGTAAATGGCAAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCA‑AT‑GGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAAc            <  1:1265765/238‑1 (MQ=255)                                                                                                                                                                                                                      tAAC‑GCTGGTGCTGGCGGGCCTGGTA‑CGGATGGCTGCGCGGGGGTTAATGGCGAAACCGCCCCCGGCCGGTAACCCGGGGGATCTCGGCTGGGTGGCCACCCCCCCC‑AA‑GGGGATTCTGCTTGTCATGATGCTTGCGATGGGAACGCATATTCCCCAACCTGTCCTCAGGATCCTGGCGGGGGCTTCCCCCACAGTCCCCCCCGGGACGCCCGACCAGCCTGCACaaactataccctc  >  2:2145910/1‑227 (MQ=255)                                                                                                                                                                                                                           gcgggTGCTGGCGGGCCGGGTA‑CGGGTGGCTGCGGGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCA‑AT‑GGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTc                                                     <  1:749919/180‑1 (MQ=255)                                                                                                                                                                                                                               gTGCTGGCGGGCCTGGTA‑CGGATGGCTGCGCGGGGGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCC‑AT‑gggggtt                                                                                                                            >  1:1285232/1‑103 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                         CCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAAC‑GCTGGTGCTGGCGGGCCTGGTA‑CGGATGGCTGCGCGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCA‑AT‑GGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTG  >  NC_007779/2606645‑2607095

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