Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,608,253	C→G	26.2%	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,253	0	C	G	26.2%	52.4 / 3.9	36	A446G (GCG→GGG)	hyfF	hydrogenase 4, membrane subunit
Reads supporting (aligned to +/- strand):  ref base C (19/7);  new base G (6/4);  total (25/11)
Fisher's exact test for biased strand distribution p-value = 4.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCT‑GGTACGGAT‑GGC‑TG‑CGCGGGTGTTAAT‑GGCGAAA‑CCGCCGCAGGCC‑GTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTG‑TCATCAGGATCCTGGCGGGCGCTTCCACT‑ATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAA  >  NC_000913/2608022‑2608491                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                   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‑Gt                                                                                                                                                                                                                                >  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‑GTTa                                                                                                                                                                                                                              >  1:1147743/1‑239 (MQ=255)                                  tGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCT‑GGTACGGAT‑GGC‑TG‑CGCGGGGGTTAAT‑GGCGAAA‑CCGCGCCAGGCG‑GTTCACCCGGGTGAGCTCGGcgga                                                                                                                                                                                                          >  1:904821/1‑236 (MQ=255)                                                 tGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCT‑GGTTCGGGT‑GGC‑TT‑CGCGGGGGGTAAT‑GGCGGAA‑CCGCCCCCGGCC‑GTTAACCGGGGGGATCCTGGCTGGGTGtccacccccccc                                                                                                                                                                                           >  1:569280/1‑233 (MQ=255)                                                      tttGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAAAGCTGGGGCTGGCGGGCTC‑GGTACCGAT‑GGC‑TGCCGCGGG‑GGTAAT‑GGGGGAA‑ACGCCCCCGCCCGGTTACCCGGGGGGTCTTCGGCGGGTTGGCCCCCCCCCCAAAgg                                                                                                                                                                                       >  2:1141848/1‑239 (MQ=255)                                                              cggTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCT‑GGGACGGAT‑TGC‑TG‑CGCGGGGGTTAAT‑GGCGAAA‑CCGCCGCAGGCC‑CGTTACCGGGGTGGTCTCCGCTGGGTGACCCCACCCCCACTGGGGATTttgc                                                                                                                                                                              >  1:204547/1‑239 (MQ=255)                                                                                              tCAATATTTTCTTAAGCGATTTTAGGCCCATTCCCGCCGGCTGGGCCCTTAATCCCTTGCTGATTACCGTCCTGCTGTTATTGCTGTTAACGTTGGTGCTGGCGGCCCT‑GGTACGGAT‑GGC‑TG‑CGCGGTTTTTAAT‑GGCGAAA‑CCGCCGCAGGCC‑GTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATAtt                                                                                                                                              <  1:798301/239‑1 (MQ=255)                                                                                                                                    gggcTGGCACGTAATCACCTGGTGATTATCGTCCTGGTGTTTTGGCTGTTACCGCTGGTGCTGGCGGGCGT‑GGTCGGGAT‑GGC‑TG‑CGCGGGTGTTAAG‑GGCGAAA‑CCGCCGCAGGCC‑GTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTG‑TCATCAGGATCCTGGCGGGCGCTTCCa                                                                                                        <  1:460973/235‑1 (MQ=255)                                                                                                                                                    aCCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCT‑GGTACGGAT‑GGC‑TG‑CGCGGGTGTTAAT‑GGCGAAA‑CCGCCGCAGGCC‑GTTAACCGGGGGGATCCCGGCTGGTTGACCACCCCGCCAATg                                                                                                                                                                                        >  1:799061/1‑143 (MQ=255)                                                                                                                                                         ctgATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCT‑GGTACGGAT‑GGC‑TG‑CGCGGGGGTTAAT‑GGCGAAA‑CCGCCGCAGGCC‑GTTAACCGGGGTGATCTCGGCGGGTTGACCACCTCGCCAATGGGGATTTTCCTGGTCATGATGCGGGGGGTGGGGACGCCAATTCC‑CAACCCGTTCCACAGGATACTGTCGGGGGGTTCCACTCATTGTCCTC‑CAggggagc                                                                                  >  2:147783/1‑234 (MQ=255)                                                                                                                                                                 cggtgttttgTTTTGGCTGTTACCGCTGGTGGTGGCGGGCCG‑GGGGCGGAT‑GGC‑GG‑CGCGGGTGTTAAG‑GGCGAAA‑CCGCCGCAGGCC‑GTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTG‑TCATCAGGATCCTGGCGGGCGCTTCCACT‑ATAGTCCTCTCAGGGACGCACGATctgc                                                                         <  1:1368249/231‑1 (MQ=255)                                                                                                                                                                        ctgTTATTGCTGTTAACGCTGGTGCTGGCGGGCCT‑GGTACGGAT‑GGC‑TG‑CGCGGGTGTTAAT‑GGCGAAA‑CCGCCGCAGGCC‑GTTAACCGGGGGGATCTCGGCTGGTTGACCACCCCCCCAAGGGGGATTTTGCCGGGCTGGAGGGTGGGCAGGGGAAAGCAAATTCCCAAACCCG‑TCATCAGGATCCTGGCGGGGGCTTCCACA‑CAAGTCTCCCCAGGGGAGCACGATCTGCCAGcaca                                                                  >  2:1363420/1‑239 (MQ=255)                                                                                                                                                                             ttggtggTTAACGTGGGTGCGGGGGGGCCG‑TGGAGGGTT‑GTC‑TG‑CGGGGTTGTTATG‑GGCGAAC‑CCGCCGCAGGCC‑GTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTG‑TCATCAGGATCCTGGCGGGCGCTTCCACT‑ATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTa                                                             <  1:1314385/233‑1 (MQ=255)                                                                                                                                                                                     tAACGCTGGTGCTGGCGGGCCT‑GGTACGGAT‑GGC‑TG‑CGCGGGGGTTAAT‑GGCGAAA‑CCGCCGCAGGCC‑GTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTG‑TCATCAGGATCCTGGCGGGCGCTTCCCCC‑TTAGTTCTCTCAAGGGCGCACGCACTGGCTGCACAACGTTGCAcccg                                                      >  1:866302/1‑236 (MQ=255)                                                                                                                                                                                                 gtggggtgcggggTACGGAT‑GGC‑GG‑CCGGGGTTTTAAG‑GGGGAAC‑CCGCCGCAGGCC‑GTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTG‑TCATCAGGATCCTGGCGGGCGCTTCCACT‑ATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGAttttt                                           <  1:599386/227‑1 (MQ=255)                                                                                                                                                                                                          cgggTACGTAT‑GGC‑GG‑CGGGGGTTTAAAT‑GGCGAAA‑CCGCCGCAGGCC‑GTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTG‑TCATCAGGATCCTGGCGGGCGCTTCCACT‑ATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAgg                                 <  1:964625/237‑1 (MQ=255)                                                                                                                                                                                                                   aT‑GGC‑TG‑CGCGGGTGTTAAT‑GGCGAAA‑CCGCCGCAGGCC‑GTTGACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGCTTCTGCTGGTCATGATGCTGGCGCTGGGAACGCATATTCCCCAACCTG‑TCCTCAGGATCCTGGCGGGCGCTTCCCCT‑ATAGACCTTTCAGGGCCGCACCGTCGCCCTGCCCCACGGACCACCTGGCAAGATTTTTTGCCTTCCgggcccccacc                        >  2:938245/1‑230 (MQ=255)                                                                                                                                                                                                                   aT‑GGC‑TG‑CGCGGGTGTTAAT‑GGCGAAA‑CCGCCGCAGGCC‑GTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTG‑TCATCAGGATCCTGGCGGGCGCTTCCACT‑ATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCCCa                          >  1:1103409/1‑237 (MQ=255)                                                                                                                                                                                                                        c‑TG‑CGCGGGTGTTAAT‑GGCGAAA‑CCGCCGCAGGCC‑GTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTG‑TCATCAGGATCCTGGCGGGCGCTTCCACT‑ATAGTCCTCTCAGGGACGCACGATCTGCCTGCCCAACGTAGAACCTGCACTGATTTTTTGCCTTCAGGCACCCCACCAGtt                    >  1:409994/1‑239 (MQ=255)                                                                                                                                                                                                                          tG‑CGCGGGTGTTAAT‑GGCGAAA‑CCGCCGCAGGCC‑GTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTG‑TCATCAGGATCCTGGCGGGCGCTTCCACT‑ATAGACCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCa                                                          >  1:44817/1‑200 (MQ=255)                                                                                                                                                                                                                          tG‑CGCGGGTGTTAAT‑GGCGAAA‑CCGCCGCAGGCC‑GTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGGGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTG‑TCATCAGGATCCTGGCGGGCGCTTCCACT‑ATAAGCCCCTCCGGGACGCCCGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCCCCCCCTCTtttt                   >  2:1116922/1‑239 (MQ=255)                                                                                                                                                                                                                            gggggggTTTTAAT‑GGCAAAC‑CCGCCGCAGGCC‑GTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTG‑TCATCAGGATCCTGGCGGGCGCTTCCACT‑ATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGttt                   <  1:495751/234‑1 (MQ=255)                                                                                                                                                                                                                             ggggggTTTTTAT‑GGGGACACCCGCCGCAGGCC‑GTTAACCGGGGGGATCTCGGGTGGTTGCCCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTG‑TCATCAGGATCCTGGCGGGCGCTTCCACT‑ATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGttt                   <  2:1147743/235‑1 (MQ=255)                                                                                                                                                                                                                             cgcgGGTGTTAAT‑GGCGAAA‑CCGCCGCAGGCC‑GTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTG‑TCATCAGGATCCTGGCGGGCGCTTCCACT‑ATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCCCATCTGTTTCg                 >  1:472458/1‑239 (MQ=255)                                                                                                                                                                                                                              gggggTTTAAAG‑GGCAAAA‑CCGCCGCAGGCC‑GTTAACCGGGGTGATCCCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTG‑TCATCAGGATCCTGGCGGGCGCTTCCACT‑ATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCg                 <  1:1116922/236‑1 (MQ=255)                                                                                                                                                                                                                                  gtgtTAAT‑GGCGAAA‑CCGCCGCAGGCC‑GTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTG‑TCATCAGGATCCTGGCGGGCGCTTCCACT‑ATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCACCTGTTTCCGAGaa            >  1:358943/1‑239 (MQ=255)                                                                                                                                                                                                                                    tgaaaT‑GGGTAAA‑CCGCCGCAGGCC‑GTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTG‑TCATCAGGATCCTGGCGGGCGCTTCCACT‑ATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGaaa           <  1:200428/235‑1 (MQ=255)                                                                                                                                                                                                                                    gtTAAT‑GGCGAAA‑CCGCCGCAGGCC‑GTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTG‑TCATCAGGATCCTGGCGGGCGCTTCCACT‑ATAGTCCTCTCAGGGAGCCACGATCTGCCTGCACAACGAAGCACCTGGCATGCTTTTTTGCCTTACGGCACCGCATCTGTTTCGGAGAAac          >  1:990001/1‑239 (MQ=255)                                                                                                                                                                                                                                     ttAAT‑GGCGAAA‑CCGCCGCAGGCC‑GTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTG‑TCATCAGGATCCTGGCGGGCGCTTCCACT‑ATAGTCCTCTCAGGGACGCACGCACTGCCCGCACAACGTAACACCCGGCATGAGTTTTTGCCCTCACGGACCGCACCTTGTTCGGAGaaaaa         >  2:522198/1‑237 (MQ=255)                                                                                                                                                                                                                                     ttAAG‑GGGGAAA‑CCGCCGCAGGCC‑GTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTG‑TCATCAGGATCCTGGCGGGCGCTTCCACT‑ATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAaca         <  1:1131280/239‑1 (MQ=255)                                                                                                                                                                                                                                             cGAAA‑CCGCCGCAGGCC‑GTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTG‑TCATCAGGATCCTGGCGGGCGCTTCCCCT‑ATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCCAATCTGTTTCGGGGaaaaccaaggaa  >  2:1229709/1‑230 (MQ=255)                                                                                                                                                                                                                                             cGAAA‑CCGCCGCAGGCC‑GTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTG‑TCATCAGGATCCTGGCGGGCGCTTCCACT‑ATAGTCCTTTCAGGGACGCACGAACGGCCGGCCGAACGTACCACCTGGCATGATTTTTTGCCTTTGGGCACCCCATCTTTTTCGGCAAAAAACAGGGaa  >  2:189007/1‑239 (MQ=255)                                                                                                                                                                                                                                             cGAAA‑CCGCCGCAGGCC‑GTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTG‑TCATCAGGATCCTGGCGGGCGCTTCCACT‑ATAGTCCTCTCAGGGACGACCGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCCCCGCATCTGTTTCGGCAAAAAACAGTGa   >  1:506843/1‑238 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                   ATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCT‑GGTACGGAT‑GGC‑TG‑CGCGGGTGTTAAT‑GGCGAAA‑CCGCCGCAGGCC‑GTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTG‑TCATCAGGATCCTGGCGGGCGCTTCCACT‑ATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAA  >  NC_000913/2608022‑2608491

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