Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,608,251	C→G	26.4%	A446G (GCG→GGG)	hyfF →	hydrogenase 4, membrane subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,608,251	0	C	G	26.4%	62.8 / 5.4	43	A446G (GCG→GGG)	hyfF	hydrogenase 4, membrane subunit
Reads supporting (aligned to +/- strand):  ref base C (24/7);  new base G (6/6);  total (30/13)
Fisher's exact test for biased strand distribution p-value = 1.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

ATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCT‑GGTACGGAT‑GGC‑TG‑CGCGGGTGTTAAT‑GGCGAAA‑CCGCCGCAGGCC‑G‑TTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTG‑TCATCAGGATCCTGGCGGGCGCTTCCACT‑ATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAA  >  NC_000913/2608020‑2608489                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                    aTCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCT‑GGGACGGATGGGT‑TC‑CGCGGG‑GTTAAT‑GGCGAAA‑cc                                                                                                                                                                                                                                              >  1:1056006/1‑238 (MQ=255)   ctcAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTAACGCCGGACCGGCACGTAATCAACTGCTGATTATCGTCCTGCTGTTATTGCTTTTAACGCCGGGGCTGGGGGGCCTGGGGAGGGAG‑GGC‑TG‑C‑CGGGGGTTTATAGGGGAAA‑ccc                                                                                                                                                                                                                                             >  2:798301/1‑237 (MQ=255)        cgtcgtCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCT‑GGTACGGAT‑GGC‑TG‑CGCGGGGGTTAAT‑GGGGGAA‑Acgccgcgg                                                                                                                                                                                                                                       >  1:270630/1‑236 (MQ=255)            gtcTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCT‑GGTACGGAT‑GGCTTG‑CCCGGG‑GTTAAT‑GGCGCAA‑CCGCCGCCGGccct                                                                                                                                                                                                                                  >  1:913573/1‑237 (MQ=255)            gtcTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCGGTTAACGCTGGTGCTGGCGGGCCT‑GGTACGGAT‑GGC‑TT‑GGCGGGTGTTAAT‑GTGGCAA‑CCGCCCCCGGCC‑G‑t                                                                                                                                                                                                                                >  1:196156/1‑239 (MQ=255)             tcTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCT‑GGGACGGGT‑GGG‑TG‑CGCGGGGGGTAAT‑GGGGAAA‑CCGCCGCAGGGC‑ggt                                                                                                                                                                                                                                >  1:709097/1‑237 (MQ=255)             tcTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGGGCTGGCGGGCCC‑GGGTAGGGT‑GGC‑TG‑CGCGGGGTTTTAT‑GGCGTAA‑CCGCCgccggcc‑ggt                                                                                                                                                                                                                                >  1:1202557/1‑237 (MQ=255)              cTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCT‑GGCACGGGT‑GGC‑TG‑GCCGGGGGGTAAT‑GCGGAAA‑CCGCCCCCGGCC‑G‑TTa                                                                                                                                                                                                                              >  1:1147743/1‑239 (MQ=255)                                  tGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCT‑GGTACGGAT‑GGC‑TG‑CGCGGGGGTTAAT‑GGCGAAA‑CCGCGCCAGGCG‑G‑TTCACCCGGGTGAGCTCGGcgga                                                                                                                                                                                                          >  1:904821/1‑236 (MQ=255)                                  tGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCT‑GGTACCGAT‑GGC‑TG‑CGCGGGGGGTTAT‑GGCGAAA‑CCGCCGCCGGGC‑GTTTAACCggggggttccgggcgg                                                                                                                                                                                                           >  1:1354001/1‑227 (MQ=255)                                                 tGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCT‑GGTTCGGGT‑GGC‑TT‑CGCGGGGGGTAAT‑GGCGGAA‑CCGCCCCCGGCC‑G‑TTAACCGGGGGGATCCTGGCTGGGTGtccacccccccc                                                                                                                                                                                           >  1:569280/1‑233 (MQ=255)                                                 tGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTAACGCCGGACTGGCACGTAATCACCTTCTGATTTACGTCCTGCTGTTTTTGCTGTTAACGCTGGCGCTGGCGGGGCT‑GGTCCGGAG‑GTC‑GG‑CGCGCGGGTTTGT‑TGCGACA‑CCACCGCCGGCC‑G‑CTATCCGGGGGGGTGTCggtggggggagcccccacgca                                                                                                                                                                                           >  2:200428/1‑220 (MQ=255)                                                      tttGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAAAGCTGGGGCTGGCGGGCTC‑GGTACCGAT‑GGC‑TGCCGCGGG‑GGTAAT‑GGGGGAA‑ACGCCCCCGCCCGG‑TTACCCGGGGGGTCTTCGGCGGGTTGGCCCCCCCCCCAAAgg                                                                                                                                                                                       >  2:1141848/1‑239 (MQ=255)                                                              cggTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCT‑GGGACGGAT‑TGC‑TG‑CGCGGGGGTTAAT‑GGCGAAA‑CCGCCGCAGGCC‑C‑GTTACCGGGGTGGTCTCCGCTGGGTGACCCCACCCCCACTGGGGATTttgc                                                                                                                                                                              >  1:204547/1‑239 (MQ=255)                                                                                              tCAATATTTTCTTAAGCGATTTTAGGCCCATTCCCGCCGGCTGGGCCCTTAATCCCTTGCTGATTACCGTCCTGCTGTTATTGCTGTTAACGTTGGTGCTGGCGGCCCT‑GGTACGGAT‑GGC‑TG‑CGCGGTTTTTAAT‑GGCGAAA‑CCGCCGCAGGCC‑G‑TTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATAtt                                                                                                                                              <  1:798301/239‑1 (MQ=255)                                                                                                                              aaaagagtggggccgcGTAACCCCCTGCTTATTGTGGCCCTGGTTTTTTTGGTGTTACCGGGGGGGTTGGGGGGCGG‑TGCCGGGAG‑GGC‑GG‑GGGGGTTGTTAAG‑GGGTAAA‑CCGCCGCAGGCC‑G‑TTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTG‑TCATCAGGATCCTGGCGGgcgc                                                                                                             <  2:757009/224‑1 (MQ=255)                                                                                                                                    gggcTGGCACGTAATCACCTGGTGATTATCGTCCTGGTGTTTTGGCTGTTACCGCTGGTGCTGGCGGGCGT‑GGTCGGGAT‑GGC‑TG‑CGCGGGTGTTAAG‑GGCGAAA‑CCGCCGCAGGCC‑G‑TTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTG‑TCATCAGGATCCTGGCGGGCGCTTCCa                                                                                                        <  1:460973/235‑1 (MQ=255)                                                                                                                                                    aCCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCT‑GGTACGGAT‑GGC‑TG‑CGCGGGTGTTAAT‑GGCGAAA‑CCGCCGCAGGCC‑G‑TTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGCTGGGAACGCCAATTCCACAACCTT‑TCATCAGGATCCTGGCGGGCGCTTCCACc‑ataactgcgcctttt                                                                                      >  1:757009/1‑227 (MQ=255)                                                                                                                                                    aCCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCT‑GGTACGGAT‑GGC‑TG‑CGCGGGTGTTAAT‑GGCGAAA‑CCGCCGCAGGCC‑G‑TTAACCGGGGGGATCCCGGCTGGTTGACCACCCCGCCAATg                                                                                                                                                                                        >  1:799061/1‑143 (MQ=255)                                                                                                                                                         ctgATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCT‑GGTACGGAT‑GGC‑TG‑CGCGGGGGTTAAT‑GGCGAAA‑CCGCCGCAGGCC‑G‑TTAACCGGGGTGATCTCGGCGGGTTGACCACCTCGCCAATGGGGATTTTCCTGGTCATGATGCGGGGGGTGGGGACGCCAATTCC‑CAACCCGTTCCACAGGATACTGTCGGGGGGTTCCACTCATTGTCCTC‑CAggggagc                                                                                  >  2:147783/1‑234 (MQ=255)                                                                                                                                                                 cggtgttttgTTTTGGCTGTTACCGCTGGTGGTGGCGGGCCG‑GGGGCGGAT‑GGC‑GG‑CGCGGGTGTTAAG‑GGCGAAA‑CCGCCGCAGGCC‑G‑TTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTG‑TCATCAGGATCCTGGCGGGCGCTTCCACT‑ATAGTCCTCTCAGGGACGCACGATctgc                                                                         <  1:1368249/231‑1 (MQ=255)                                                                                                                                                                      gtttttttgtgtgtttagcgtgggtggtGGCGGGCGT‑GGTACGGAT‑GGT‑TG‑CGCGGGTCTAAAG‑GCGGAAA‑CCGGCGCGGGGC‑G‑TTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTG‑TCATCAGGATCCTGGCGGGCGCTTCCACT‑ATAGTCCTCTCAGGGACGCACGATCTGCCTGca                                                                    <  2:1255149/217‑1 (MQ=255)                                                                                                                                                                        ctgTTATTGCTGTTAACGCTGGTGCTGGCGGGCCT‑GGTACGGAT‑GGC‑TG‑CGCGGGTGTTAAT‑GGCGAAA‑CCGCCGCAGGCC‑G‑TTAACCGGGGGGATCTCGGCTGGTTGACCACCCCCCCAAGGGGGATTTTGCCGGGCTGGAGGGTGGGCAGGGGAAAGCAAATTCCCAAACCCG‑TCATCAGGATCCTGGCGGGGGCTTCCACA‑CAAGTCTCCCCAGGGGAGCACGATCTGCCAGcaca                                                                  >  2:1363420/1‑239 (MQ=255)                                                                                                                                                                             ttggtggTTAACGTGGGTGCGGGGGGGCCG‑TGGAGGGTT‑GTC‑TG‑CGGGGTTGTTATG‑GGCGAAC‑CCGCCGCAGGCC‑G‑TTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTG‑TCATCAGGATCCTGGCGGGCGCTTCCACT‑ATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTa                                                             <  1:1314385/233‑1 (MQ=255)                                                                                                                                                                             aTTGCTGTTAACGCTGGTGCTGGCGGGCCT‑GGTACGGAT‑GGC‑TG‑CGCGGGGGTTTAT‑GGCGGAA‑CCGCCGCAGGCC‑G‑GTAACCGGGGTGATCCCGGGTGGTTGACCACCTCGCCAATGGGGATTCTGCTGGTAATGATGCTGGGCGTGGGAACGCATTTTCCAAAACCTG‑GCATCAGGACCCTGGGGGGCGCGTTCACT‑ATAAGCCCCTCCGGGGCCCAcacatggcccgccaaaaggg                                                             >  2:1314385/1‑220 (MQ=255)                                                                                                                                                                                     tAACGCTGGTGCTGGCGGGCCT‑GGTACGGAT‑GGC‑TG‑CGCGGGGGTTAAT‑GGCGAAA‑CCGCCGCAGGCC‑G‑TTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTG‑TCATCAGGATCCTGGCGGGCGCTTCCCCC‑TTAGTTCTCTCAAGGGCGCACGCACTGGCTGCACAACGTTGCAcccg                                                      >  1:866302/1‑236 (MQ=255)                                                                                                                                                                                                 gtggggtgcggggTACGGAT‑GGC‑GG‑CCGGGGTTTTAAG‑GGGGAAC‑CCGCCGCAGGCC‑G‑TTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTG‑TCATCAGGATCCTGGCGGGCGCTTCCACT‑ATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGAttttt                                           <  1:599386/227‑1 (MQ=255)                                                                                                                                                                                                          cgggTACGTAT‑GGC‑GG‑CGGGGGTTTAAAT‑GGCGAAA‑CCGCCGCAGGCC‑G‑TTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTG‑TCATCAGGATCCTGGCGGGCGCTTCCACT‑ATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAgg                                 <  1:964625/237‑1 (MQ=255)                                                                                                                                                                                                                   aT‑GGC‑TG‑CGCGGGTGTTAAT‑GGCGAAA‑CCGCCGCAGGCC‑G‑TTGACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGCTTCTGCTGGTCATGATGCTGGCGCTGGGAACGCATATTCCCCAACCTG‑TCCTCAGGATCCTGGCGGGCGCTTCCCCT‑ATAGACCTTTCAGGGCCGCACCGTCGCCCTGCCCCACGGACCACCTGGCAAGATTTTTTGCCTTCCgggcccccacc                        >  2:938245/1‑230 (MQ=255)                                                                                                                                                                                                                   aT‑GGC‑TG‑CGCGGGTGTTAAT‑GGCGAAA‑CCGCCGCAGGCC‑G‑TTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTG‑TCATCAGGATCCTGGCGGGCGCTTCCACT‑ATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCCCa                          >  1:1103409/1‑237 (MQ=255)                                                                                                                                                                                                                        c‑TG‑CGCGGGTGTTAAT‑GGCGAAA‑CCGCCGCAGGCC‑G‑TTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTG‑TCATCAGGATCCTGGCGGGCGCTTCCACT‑ATAGTCCTCTCAGGGACGCACGATCTGCCTGCCCAACGTAGAACCTGCACTGATTTTTTGCCTTCAGGCACCCCACCAGtt                    >  1:409994/1‑239 (MQ=255)                                                                                                                                                                                                                          tG‑CGCGGGTGTTAAT‑GGCGAAA‑CCGCCGCAGGCC‑G‑TTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCCCAACCTG‑TCATCAGGATCCTGGCGGGCGCTTCCACT‑TTAGTCCCCTCAGGGCCCGACGATCTGCCTGCACACAGAAGCACACGGCTCAGTTTTTTGGCATCAggacacgcaacctttt                   >  2:955766/1‑225 (MQ=255)                                                                                                                                                                                                                          tG‑CGCGGGTGTTAAT‑GGCGAAA‑CCGCCGCAGGCC‑G‑TTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTG‑TCATCAGGATCCTGGCGGGCGCTTCCACT‑ATAGACCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCa                                                          >  1:44817/1‑200 (MQ=255)                                                                                                                                                                                                                          tG‑CGCGGGTGTTAAT‑GGCGAAA‑CCGCCGCAGGCC‑G‑TTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGGGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTG‑TCATCAGGATCCTGGCGGGCGCTTCCACT‑ATAAGCCCCTCCGGGACGCCCGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCCCCCCCTCTtttt                   >  2:1116922/1‑239 (MQ=255)                                                                                                                                                                                                                            gggggggTTTTAAT‑GGCAAAC‑CCGCCGCAGGCC‑G‑TTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTG‑TCATCAGGATCCTGGCGGGCGCTTCCACT‑ATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGttt                   <  1:495751/234‑1 (MQ=255)                                                                                                                                                                                                                             ggggggTTTTTAT‑GGGGACACCCGCCGCAGGCC‑G‑TTAACCGGGGGGATCTCGGGTGGTTGCCCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTG‑TCATCAGGATCCTGGCGGGCGCTTCCACT‑ATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGttt                   <  2:1147743/235‑1 (MQ=255)                                                                                                                                                                                                                             cgcgGGTGTTAAT‑GGCGAAA‑CCGCCGCAGGCC‑G‑TTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTG‑TCATCAGGATCCTGGCGGGCGCTTCCACT‑ATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCCCATCTGTTTCg                 >  1:472458/1‑239 (MQ=255)                                                                                                                                                                                                                              gggggTTTAAAG‑GGCAAAA‑CCGCCGCAGGCC‑G‑TTAACCGGGGTGATCCCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTG‑TCATCAGGATCCTGGCGGGCGCTTCCACT‑ATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCg                 <  1:1116922/236‑1 (MQ=255)                                                                                                                                                                                                                                  gtgtTAAT‑GGCGAAA‑CCGCCGCAGGCC‑G‑TTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTG‑TCATCAGGATCCTGGCGGGCGCTTCCACT‑ATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCACCTGTTTCCGAGaa            >  1:358943/1‑239 (MQ=255)                                                                                                                                                                                                                                    tgaaaT‑GGGTAAA‑CCGCCGCAGGCC‑G‑TTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTG‑TCATCAGGATCCTGGCGGGCGCTTCCACT‑ATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGaaa           <  1:200428/235‑1 (MQ=255)                                                                                                                                                                                                                                    gtTAAT‑GGCGAAA‑CCGCCGCAGGCC‑G‑TTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTG‑TCATCAGGATCCTGGCGGGCGCTTCCACT‑ATAGTCCTCTCAGGGAGCCACGATCTGCCTGCACAACGAAGCACCTGGCATGCTTTTTTGCCTTACGGCACCGCATCTGTTTCGGAGAAac          >  1:990001/1‑239 (MQ=255)                                                                                                                                                                                                                                     ttAAT‑GGCGAAA‑CCGCCGCAGGCC‑G‑TTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTG‑TCATCAGGATCCTGGCGGGCGCTTCCACT‑ATAGTCCTCTCAGGGACGCACGCACTGCCCGCACAACGTAACACCCGGCATGAGTTTTTGCCCTCACGGACCGCACCTTGTTCGGAGaaaaa         >  2:522198/1‑237 (MQ=255)                                                                                                                                                                                                                                     ttAAG‑GGGGAAA‑CCGCCGCAGGCC‑G‑TTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTG‑TCATCAGGATCCTGGCGGGCGCTTCCACT‑ATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAaca         <  1:1131280/239‑1 (MQ=255)                                                                                                                                                                                                                                             cGAAA‑CCGCCGCAGGCC‑G‑TTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTG‑TCATCAGGATCCTGGCGGGCGCTTCCCCT‑ATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCCAATCTGTTTCGGGGaaaaccaaggaa  >  2:1229709/1‑230 (MQ=255)                                                                                                                                                                                                                                             cGAAA‑CCGCCGCAGGCC‑G‑TTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTG‑TCATCAGGATCCTGGCGGGCGCTTCCACT‑ATAGTCCTTTCAGGGACGCACGAACGGCCGGCCGAACGTACCACCTGGCATGATTTTTTGCCTTTGGGCACCCCATCTTTTTCGGCAAAAAACAGGGaa  >  2:189007/1‑239 (MQ=255)                                                                                                                                                                                                                                             cGAAA‑CCGCCGCAGGCC‑G‑TTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTG‑TCATCAGGATCCTGGCGGGCGCTTCCACT‑ATAGTCCTCTCAGGGACGACCGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCCCCGCATCTGTTTCGGCAAAAAACAGTGa   >  1:506843/1‑238 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                    ATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCT‑GGTACGGAT‑GGC‑TG‑CGCGGGTGTTAAT‑GGCGAAA‑CCGCCGCAGGCC‑G‑TTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTG‑TCATCAGGATCCTGGCGGGCGCTTCCACT‑ATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAA  >  NC_000913/2608020‑2608489

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