Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	910,811	T→G	23.6%	H80P (CAC→CCC)	poxB ←	pyruvate dehydrogenase, thiamine triphosphate‑binding, FAD‑binding

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	910,811	0	T	G	23.6%	94.0 / 6.3	66	H80P (CAC→CCC)	poxB	pyruvate dehydrogenase, thiamine triphosphate‑binding, FAD‑binding
Reads supporting (aligned to +/- strand):  ref base T (22/27);  new base G (11/6);  total (33/33)
Fisher's exact test for biased strand distribution p-value = 2.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAACCGAAACGCCACGGTTAAGCACCGCTTTGCGCATGGCAATCGCCAGTACTTGTGGGATCTGCTCCGGGCTGGAAACCAGCTCGCAATAGTGACTACATTCGCGGAATAGCTCTTGTGGGTGGGTTTCCTGGAAA‑TAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAA‑T‑CGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGAT‑TAAGTGCAGGTT‑GCCGGGGCCGCACGATCCG‑GCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGC‑CAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACTCGATGGTGCCCATGCGATTAAGACTGTCACTAAGACCGTTCAGAGAGTCGCCTGTGACTCCCCAGATGCGTTTCACCCCTGCCGATTCGAGTGTTTTGGCGATATAAGCTGCAACCGTTTGTTT  >  NC_000913/910581‑911046                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                               cAACCGAAACGCCACGGTTAAGCACCGCTTTGCGCATGGCAATCGCCAGTACTTGTGGGATCTGCTCCGGGCTGGAAACCAGCTCGCAATAGTGACTACATTCGCGGAATAGCTCTTGTGGGTGGGTTTCCTGTAAA‑TAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAA‑T‑CGCCAGTACCGGAACTCGATTCAGGTGGCAATCGAACAGGCCGTCGAT‑TAAGTGCAGGtg                                                                                                                                                                                                                                         >  1:828214/1‑237 (MQ=255) cAACCGAAACGCCACGGTTAAGCACCGCTTTGCGCATGGCAATCGCCAGTACTTGTGGGATCTGCTCCGGGCTGGAAACCAGCTCGCAATAGTGACTACATTCGCGGAATAGCTCTTGTGGGTGGGTTTCCTGGAAA‑TAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAA‑T‑CGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAAAAGGCCGTTGAT‑TAAGTGAAGGt                                                                                                                                                                                                                                          >  1:1201703/1‑237 (MQ=255)            ccgcGGTTAAGCCCCGCTTTGCCATTGGCAATCGCCAGTACTTGTGGGATCTGCTCCGGGCTGGAACCCAGCTCGCAATAGTGACTACATTCGCGGAATAGCTCTTGTGGGTGGGTTTCCTGGAAA‑TAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAA‑T‑CGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGAT‑TAAGTGCAGGTT‑GCCGGGGCCGCa                                                                                                                                                                                                                            <  1:305634/236‑1 (MQ=255)            ccccGGTTAAGCACCGCTTTGGGAATGGCAATCGCCAGTACTTGTGGGATCTGCTCCGGGCTGGAAACCAGCTCGCAATAGTGCCTACATTCGCGGAATAGCTCTTGTGGGTGGGTTTCCTGGAAA‑TAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAA‑T‑CGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGAT‑TAAGTGCAGGTT‑GCCGGGGCCGCa                                                                                                                                                                                                                            <  1:570158/236‑1 (MQ=255)            ccACGGTTAAGCACCGCTTTGCGCATGGCAATCGCCAGTACTTGTGGGATCTGCTCCGGGCTGGAAACCAGCTCGCAATAGTGACTACATTCGCGGAATAGCTCTTGTGGGTGGGTTTCCTGGAAA‑TAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAA‑T‑CGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGAT‑TAAGTGCAGGTT‑GCCGGGGCCGCa                                                                                                                                                                                                                            <  1:1122712/239‑1 (MQ=255)                        acccttttGCGCATGGCAATCGCCAGTCTTTGTGGGATTTCCCCTGGGCTGGAAACCAGTTCGTAATAGGGACTACATTCGCGGAATAGCTCTTGTGGGTGGGTTTCCTGGAAA‑TAGCCGCTGCCAATTTCGCTGGGGGGAATATGAGCGGCAA‑T‑CGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGAT‑TAAGTGCAGGTT‑GCCGGGGCCGCACGATCCG‑gcg                                                                                                                                                                                                                 <  2:1178228/232‑1 (MQ=255)                           gtttggccattGGGAACCGCCATTACTTGGGGGATCTGCCCCGGGCTGGAAACCAGCTCGCAATAGTGACTACATTCGCGGAATAGCTCTTGTGGGTGGGTTTCCTGGAAA‑TAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAA‑T‑CGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGAT‑TAAGTGCAGGTT‑GCCGGGGCCGCACGATCCG‑GCGCAGAc                                                                                                                                                                                                            <  1:916631/229‑1 (MQ=255)                            cTTTGCGCATGGCAATCGCCAGTACTTGTGGGATCTGCTCCGGGCTGGAAACCAGCTCGCAATAGTGACTACATTCGCGGAATAGCTCTTGTGGGTGGGTTTCCTGGAAA‑TAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAA‑T‑CGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGAT‑TAAGTGCAGGTT‑TCCGGGGCCGCACGATCCG‑GCGCAGAc                                                                                                                                                                                                            >  1:195886/1‑238 (MQ=255)                            cTTTGCGCATGGCAATCGCCAGTACTTGTGGGATCTGCTCCGGGCTGGAAACCAGCTCGCAATAGTGACTACATTCGCGGAATAGCTCTTGTGGGTGGGTTTCCTGGAAA‑TAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAA‑T‑CGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGAT‑TAAGTGCAGGTT‑GCCGGGGCCGCCCGATCCG‑GCCCAGAc                                                                                                                                                                                                            >  1:1025636/1‑238 (MQ=255)                            cTTTGCGCATGGCAATCGCCAGTACTTGTGGGATCTGCTCCGGGCTGGAAACCAGCTCGCAATAGTGACTACATTCGCGGAATAGCTCTTGTGGGTGGGTTTCCTGGAAA‑TAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAA‑T‑CGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGAT‑TAAGTGCAGGTT‑GCCGGGGCCGCACGATCCG‑GGGCAGAc                                                                                                                                                                                                            >  1:68069/1‑238 (MQ=255)                            cTTTGCGCATGGCAATCGCCAGTACTTGTGGGATCTGCTCCGGGCTGGAAACCAGCTCGCAATAGTGACTACATTCGCGGAATAGCTCTTGTGGGTGGGTTTCCTGGAAA‑TAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAA‑T‑CGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGAT‑TAAGTGCAGGTT‑GCCGGGGCCGCACGATCCG‑GCGCAGAc                                                                                                                                                                                                            >  1:489825/1‑238 (MQ=255)                            cTTTGCGCATGGCAATCGCCAGTACTTGTGGGATCTGCTCCGGGCTGGAAACCAGCTCGCAATAGTGACTACATTCGCGGAATAGCTCTTGTGGGTGGGTTTCCTGGAAA‑TAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAA‑T‑CGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGAT‑TAAGTGCAGGTT‑GCCGGGGCCGCACGATCCG‑GCGCAGAc                                                                                                                                                                                                            >  1:607300/1‑238 (MQ=255)                            cTTTGCGCATGGCAATCGCCAGTACTTGTGGGATCTGCTCCGGGCTGGAAACCAGCTCGCAATAGTGACTACATTCGCGGAATAGCTCTTGTGGGTGGGTTTCCTGGAAA‑TAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAA‑T‑CGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGAT‑TAAGGGCAGGTT‑GCCGGGGCCGCACGATCCG‑GGCCAGAc                                                                                                                                                                                                            >  1:578938/1‑238 (MQ=255)                            cTTTGCGCATGGCAATCGCCAGTACTTGTGGGATCTGCTCCGGGCTGGAAACCAGCTCGCAATAGTGACTACATTCGCGGAATAGCTCTTGTGGGTGGGTTTCCTGGAAA‑TAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAA‑T‑CGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGAT‑ACAGGGCAGGTG‑GCCGGGGCCGCACGATCAG‑GCGCAAAc                                                                                                                                                                                                            >  1:626204/1‑238 (MQ=255)                               tGCGCATGGCAATCGCCAGTACTTGTGGGATCTGCTCCGGGCTGGAAACCAGCTCGCAATAGTGACTACATTCGCGGAATAGCTCTTGTGGGTGGGTTTCCTGGAAA‑TAGCCGCTGCCAATTTCGCTGGAGGGAATCTGAGCGGCAA‑T‑CGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGAT‑TAAGGGCAGGTT‑GCCGGGGCCGCACGTACCG‑GCGCAGACcccc                                                                                                                                                                                                        >  2:707770/1‑239 (MQ=255)                                 cgcATGGCAATCGCCAGTACTTGTGGGATCTGCTCCGGGCTGGAAACCAGCTCGCAATAGTGACTACATTCGCGGAATAGCTCTTGTGGGTGGGTTTCCTGGAAA‑TAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAA‑T‑CGCCAGTACCGGAACGTGCTTGCGGCGGCAATCGCACAGGCCGTTGAT‑TAATGGCAGGTT‑GCCGGGGCAGCACGATCCG‑GCGCAGACCGCCAg                                                                                                                                                                                                      >  1:83642/1‑239 (MQ=255)                                   gttGGCAGTCGCGTGTATTGGTGGGTTCTGCCCGGGGCGGGAAACCAGCTAGAAATAGTGACTACATTCGCGGAATAGCTTTTGTGGGTGGGTTTCCTGGAAA‑TAGCGGAGGCCAATTTCGCTGGAGGGAATATGAGCGGCAA‑T‑CGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGAT‑TAAGTGCAGGTT‑GCCGGGGCCGCACGATCCG‑GCGCAGACCGCCAGtt                                                                                                                                                                                                    <  2:437758/237‑1 (MQ=255)                                   cATGGCAATCGCTAGTCCTTGTGGGATTTGCTCCGGGCTGGAACCCAGCTCGCAATAGTGACTACATTCGCGGAATAGCTCTTGTGGGTGGGTTTCCTGGAAA‑TAGCCGCTGCCAATTGCGCTGGAGGGAATATGAGCGGCAA‑T‑CGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGAT‑TAAGTGCAGGTT‑GCCGGGGCCGCACGATCCG‑GCGCAGACCGCCAGtt                                                                                                                                                                                                    <  1:372520/239‑1 (MQ=255)                                    ctGGCAATCGCCAGTACTTGGGGGATCTGCCCCGGGCTGGAACCCAGCTCGCAATAGTGACTACATTCGCGGAATAGCTCTTGGGGGTGGGTTTCCTGGAAA‑TAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAA‑T‑CGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGAT‑TAAGTGCAGGTT‑GCCGGGGCCGCACGATCCG‑GCGCAGACCGCCAGtt                                                                                                                                                                                                    <  1:175489/237‑1 (MQ=255)                                    atggacaTCGCCAGTACTTGGGGGATCTGCTCGGGGCTGGAAACCAGCTCGCAATAGTGACTACATTCGCGGAATAGCTCTTGGGGGTGGGTTTCCTGGAAA‑TAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAA‑T‑CGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGAT‑TAAGTGCAGGTT‑GCCGGGGCCGCACGATCCG‑GCGCAGACCGCCAGtt                                                                                                                                                                                                    <  1:1371652/232‑1 (MQ=255)                                    agggCAATCGCCGGGTCTTGTGGGAACTGCCCCGGGCTGGAAACCAGCTCGCAATAGTGACTACAGTCGCGGAATAGCTCTTGTGGGTGGGTTTCCTGGAAA‑TAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAA‑T‑CGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGAT‑TAAGTGCAGGTT‑GCCGGGGCCGCACGATCCG‑GCGCAGACCGCCAGtt                                                                                                                                                                                                    <  1:837906/236‑1 (MQ=255)                                    aTGGCAATCGCCAGTACTTGTGGGATCTGCCCCGGGCTGGAACCCAGCTCGCAATAGTGACTACATTCGCGGAATAGCTCTTGTGGGTGGGTTTCCTGGAAA‑TAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAA‑T‑CGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGAT‑TAAGTGCAGGTT‑GCCGGGGCCGCACGATCCG‑GCGCAGACCGCCAGtt                                                                                                                                                                                                    <  1:707770/238‑1 (MQ=255)                                      ggCAATCGCCAGTACTTGTGGGATCTGCTCCGGGCTGGAAACCAGCTCGCAATAGTGACTACATTCGCGGAATAGCTCTTGTGGGTGGGTTTCCTGGAAA‑TAGCCGCTGCCAATTTCCCTGGAGGGAATATGAGCGGAAA‑T‑CCCCCGTCCCGAACCTGGGTTGGGGGGACAACGGAAC‑GGCGCGTGATGGAAATGCAGGTT‑GCCCGGGCCCGCCGACCCG‑GCGCAGAACGCCaatcggc                                                                                                                                                                                                 >  2:837906/1‑233 (MQ=255)                                      ggCAATCGCCAGTACTTGTGGGATCGGCTCCGGGCTGGAAACCAGCTCGCAATAGTGACTACATTCGCGGAATAGCTCTTGTGGGTGGGTTTCCTGGAAA‑TAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAA‑T‑CGCCAGTACCGGACCGTGCTTGCGATGGCAATCGAAAAGGCCGTTGAT‑TAATTCAAGGCT‑GCCGGGGCCGCACgcaccgccaccagaccccagtgac                                                                                                                                                                                                   >  2:207562/1‑215 (MQ=255)                                        caaccccttgttcTTTTGGGTTTTGCTCGGGGTTGGAAACCAGCTCGCAATAGTGACTACATTCGCGGAATAGCTCTTGTGGGTGGTTTTCCTGGAAA‑TAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAA‑T‑CGCCAGTACCGGAACGTGTTTGCGGTGGCAATCGAACAGGCCGTTGAT‑TAAGTGCAGGTT‑GCCGGGGCCGCACGATCCG‑GCGCAGACCGCCAGTTCTCCg                                                                                                                                                                                               <  2:489825/227‑1 (MQ=255)                                           tCGCCAGTACTTGTGGGATCTGCTCCGGGCTGGAACCCAGCTGGCAATAGTGACTACATTCGCGGAATAGCTCTTGGGGGTGGGTTTCCTGGAAA‑TAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAA‑T‑CGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGAT‑TAAGTGCAGGTT‑GCCGGGGCCGCACGATCCG‑GCGCAGACCGCCAGTTCTCCGCTa                                                                                                                                                                                            <  1:541941/239‑1 (MQ=255)                                           tCGCCAGTACTTGTGGGATCTGCTCCGGGCTGGAAACCAGCTCGCAATAGTGACTACATTCGCGGAATAGCTCTTGTGGGTGGGTTTCCTGGAAA‑TAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAA‑T‑CGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGAT‑TAATGGCAGGTT‑GCCGGGGCCGCACCATCCg‑ggccaaccccccaggtccccggta                                                                                                                                                                                            >  2:1031685/1‑216 (MQ=255)                                                              ttGCCCCGGGGTGGAACCACGCTCGAAATGGGGATTGCATTCGCGAAATAGCTCTTGTGGGTTGGTTTCGGGGAAA‑TGGCCGCTGCCAATGTCGCGGAAGGGAATAGGAGCGGCAA‑T‑CGCCAGTCCCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGAT‑TAAGTGCAGGTT‑GCCGGGGCCGCACGATCCG‑GCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGC‑Ca                                                                                                                                                                         <  2:644671/237‑1 (MQ=255)                                                               tGCTCCGGGCTGGAAACCAGCTCGCAATAGTGACTACATTCGCGGAATAGCTCTTGTGGGTGGGTTTCCTGGAAA‑TAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAA‑T‑CGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGAT‑TAAGTGCAGGTT‑GCCGGGGCCGCACGATCCG‑GCGCCGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGC‑Cacc                                                                                                                                                                       >  1:550534/1‑237 (MQ=255)                                                                 cTCCGGGCTGGAAACCAGCTCGCAATAGTGACTACATTCGCGGAATAGCTCTTGTGGGTGGGTTTCCTGGAAA‑TAGCCGCTGCCAATTTCGCTGGAGGGAATATGCGCGGCAA‑T‑CGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGAT‑TAAGTGCAGGTT‑GCCGGGGCCGCACGATCCG‑GCGCAGACCGCCAGTTCTCCGCTCATTTTGGCTTCGGCGc‑cgacg                                                                                                                                                                      >  1:925698/1‑234 (MQ=255)                                                                 cTCAGGGCTGGAAACCAGCTCGCAATAGTGACTACATTCGCGGAATAGCTCTTGTGGGTGGGTTTCCTGGAAA‑TAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAA‑T‑CGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACTGGCCGTTGAT‑TAAGGGCAGGTT‑GCCGGGGCCGCACGATCCG‑GCGCCGCCAGCACGTATTCAGCTAGGTTGGGCTTCAGGGc‑cgacg                                                                                                                                                                      >  2:41758/1‑234 (MQ=255)                                                                    cGGGCTGGAAACCAGCTCGCAATAGTGACTACATTCGCGGAATAGCTCTTGTGGGTGGGTTTCCTGGAAA‑TAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAA‑T‑CGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGGT‑TAAGGGCAGGTT‑GCCGGGGCCGCACGACCCG‑GGGCAGCCCGCCCGTTCTCCCCGAAGtgggggttcgtggcccagggcaa                                                                                                                                                                   >  2:1092345/1‑217 (MQ=255)                                                                    cGGGCTGGAAACCAGCTCGCAATAGCGACTACATTCGCGGAATAGCTCTTGTGGGTGGGTTTCCTGGAAA‑TAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAA‑T‑CGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGAT‑TAAGGGCAGGTT‑GCCGGGGCCGCACGATCCG‑GGGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGC‑CAGCGGCaa                                                                                                                                                                  >  1:701241/1‑239 (MQ=255)                                                                         gggaacccggctcgatatgtttaCTACATTGGGGGAATAGCTGTTGGGGGGGGTTGTCCTGAAAA‑GACCCGCTGCCAATTCCGGCGGGGGAAATTGGGGGGGCAA‑C‑CGCCAGTCCCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGAT‑TAAGTGCAGGTT‑GCCGGGGCCGCACGATCCG‑GCGCAGCCCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGC‑CAGCGGCAAAGgc                                                                                                                                                              <  2:1025636/216‑1 (MQ=255)                                                                           gAAACCAGCTCGCAATAGTGACTACATTCGCGGAATAGCTCTTGTGGGTGGGTTTCCTGGAAA‑TAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAA‑T‑CGCCAGTACCGGAACGTGATTTCGGTGGCAATCGAACAGGCCGTTGAT‑TAAGTGCAGGTT‑GCCGGGGCCGCACGATCCGAGCACA‑ACCGCCAATTCTCCCCGAACTTTGGCTTCAACGCACAG‑GGCAAAGggcgc                                                                                                                                                           >  2:140565/1‑239 (MQ=255)                                                                                cAGCTCGCAATAGTGACTACATTCGCGGAATAGCTCTTGTGGGTGGGTTTCCTGGAAA‑TAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAA‑T‑CGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGGT‑TAAGGGCCGGTT‑GCCGGGGCCGCCCGATCCG‑GCGCCGACCTCCGGTTCTCCCGTAAGTTGTGTTCCAGCCC‑CACCGGCAAAAGCCCCCAcc                                                                                                                                                       >  2:222733/1‑237 (MQ=255)                                                                                   cTCGCAATAGTGACTACATTCGCGGAATAGCTCTTGTGGGTGGGTTTCCTGGAAA‑TAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAA‑T‑CGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGAT‑TAAGTGCAGGTT‑GCCGGGGCCGCACGATCCG‑GCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGC‑CAGCGGCAAAggaccccccgtct                                                                                                                                                    >  1:240051/1‑227 (MQ=255)                                                                                               aCTACATTCGCGGAATAGCTCTTGTGGGTGGGTTTCCTGGAAA‑TAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAA‑T‑CGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGAT‑TAAGTGCAGGTT‑GCCGGGGCCGCa                                                                                                                                                                                                                            <  1:948694/156‑1 (MQ=255)                                                                                               aCTACATTCGCGGAATAGCTCTTGTGGGTGGGTTTCCTGGAAA‑TAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAA‑T‑CGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGAT‑TAAGTGCAGGTT‑GCCGGGGCCGCa                                                                                                                                                                                                                            >  2:948694/1‑156 (MQ=255)                                                                                               aCTACATTCGCGGAATAGCTCTTGTGGGTGGGTTTCCTGGAAA‑TAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAA‑T‑CGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGAT‑TAAGTGCAGGTT‑GCCGGGGCCGCACGATCCG‑GCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGC‑CAGCGGCAAAGGCCGCCACTTCTTCGGGGCgggggg                                                                                                                                       >  1:73694/1‑239 (MQ=255)                                                                                                cTACATTCGCGGAATAGCTCTTGTGGGTGGGTTTCCTGGAAA‑TAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAA‑T‑CGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGAT‑TAGGTGCAGGTT‑GCCGGGGCCGCACGATCCG‑GCGCAGACCGCCCGTTCTCCGCTAAGTTGGGCTTCACCGC‑AAGGGGCAAAGGCCGCCACTCCTCCGggggggggggg                                                                                                                                      >  2:523538/1‑238 (MQ=255)                                                                                                 ttcctttgccggatcgctgtgggggggggggTTCCTGGAAAGTACCGCCTGCAAATTTCGCGGGGGGAAAAATGGACGGAAA‑T‑GCCAAGTACGGGACGGTGATGGGGGGGGCAAGCAAACAGGCCGTTGAT‑TAAGGGCAGGTT‑GCGGGGGCCGCACGATCCG‑GCGCAGACCGCCGGTTCTCCGCTAAGTTGTGCTTCAGCGC‑CAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGa                                                                                                                                      <  2:1112667/216‑1 (MQ=255)                                                                                                 tacatcgggggattggCTCTTGGGGGGGGTTTTCCGGGAAA‑TGCCCGCTGCCAATTTCGGTGGAGGGAATTTGAGCGGCAA‑T‑CGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGAT‑TAAGTGCAGGTT‑GCCGGGGCCGCACGATCCG‑GCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGC‑CAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGAc                                                                                                                                     <  1:867004/224‑1 (MQ=255)                                                                                                    aTTCGCGGAATAGCTCTTGTGGGTGGGTTTCCTGGAAA‑TAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAA‑T‑CGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGAT‑TAAGTGCAGGTT‑GCCGGGGCCGCACGATCCG‑GCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCACCGC‑CAGCGGCAAAGGCCGCCACTTCTTCGTGG‑GGGTgtgggct                                                                                                                                  >  1:1112517/1‑232 (MQ=255)                                                                                                    aTTCGCGGAATAGCTCTTGTGGGTGGGTTTCCTGGAAA‑TAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAA‑T‑CGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGAT‑TAAGGGCAGGTT‑GCCGGGGCCGCACGATCCG‑GCGCAGACCGCCAGTTCTCGGCGAAGTTGGGCTTCAGCGC‑CGGCGGCAAAAGCCCCCACTTCTTCGgggggggggggagc                                                                                                                                   >  2:242008/1‑234 (MQ=255)                                                                                                      tCGCGGAATAGCTCTTGTGGGTGGGTTTCCTGGAAA‑TAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAA‑T‑CGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGAT‑TAAGGGCAGGTT‑GCCGGGGCCGCACGATCCG‑GCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGC‑CAGCGGCAAAGGCCGCCACTTCTTCTGGGCggggggggacgc                                                                                                                                 >  2:663232/1‑232 (MQ=255)                                                                                                                 ctctTGTGGGTGGGTTTCCTGGAAA‑TAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAA‑T‑CGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACCGGCCGTTGAT‑TAAGTGCAGGTT‑GCCGGGGCCGCACGATCCG‑GCGCAGACCGCCAGATCTCCGCTAAGTTGTGCTTCAGCGC‑CAGCGGCAAAGGCCGCCACTTCTTCGTGACGGGGTGACATCCCCTCCgtggtg                                                                                                                      >  1:486058/1‑238 (MQ=255)                                                                                                                     tgggggtgggtTCCCGGGAAA‑TTGCCGCTGCCATTTCCGTTGGAGGGAATATGGGCGGCAA‑T‑CGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGAT‑TAAGTGCAGGTT‑GCCGGGGCCGCACGATCCG‑GCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGC‑CAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACTCGATGGTGCCCAt                                                                                                                 <  1:66576/236‑1 (MQ=255)                                                                                                                      gtgggtggGTTTCCTGGAAA‑TAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAA‑T‑CGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGAT‑TAAGTGCAGGTT‑GCCGGGGCCGCACGATCCG‑GCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGC‑CAGCGGCAAAGGCCGCCACTTCTTCGTGGCggggggggcttccccccatgggggccat                                                                                                                 >  1:1001820/1‑216 (MQ=255)                                                                                                                      gtgggtggGCTTCCTGGAAA‑TAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAA‑T‑CGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGAT‑TAAGTGCAGGTT‑GCCGGGGCCGCACGATCCG‑GCGCAGACCGCCAGTTCTCCGCCAAGGTGTGCTTCAGCGC‑AAACAGCAAAGACCTCCACTTCTTCGGGGCGGGTggggccgcaccactcagtgcccccg                                                                                                                >  2:455742/1‑216 (MQ=255)                                                                                                                      gtggggggGTTTTCTGGAAA‑TAGCCCGTGCCAATTGTGGGGGGGGGAAATTGGGGGGGAATT‑CGCACTGACGGGACGGTGGTTGGGGGGGTAATCGAACAGGCCGTTGAT‑TAAGGCCCGGTT‑GCCGGGGCCGCAGGATCCG‑GCGAAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGC‑AAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACTCGATGGTGCCCAt                                                                                                                 <  2:1014587/239‑1 (MQ=255)                                                                                                                           tgggtTTCCTGGAAA‑TAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAA‑T‑CGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGAT‑TAAGTGCAGGTT‑GCCGGGGCCGCCCGATCCG‑GCGCAGACCGCCAGTTATCCGCTAAGTTGGGCTTCAGCGC‑CAGCGGCAAAGGCCCCCACTTCTTCGTGGGGGGGGGGCATcacccccgaggggcccccgcggtt                                                                                                           >  2:414982/1‑216 (MQ=255)                                                                                                                                 tCCTGGAAA‑TGGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAA‑T‑CGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGAT‑TAAGTGCAGGTT‑GCCGGGGCCGCACGATCCG‑GCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGC‑CAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACTCGATGGTGCCCATGCGATTAAGACt                                                                                                     <  1:242008/239‑1 (MQ=255)                                                                                                                                 cccggaaaa‑TGGCGGCTGCCAATTTCGCTGGAGGGAATATGGGCGGCAA‑T‑CGCCAGTACCGGACGGTGATTGCGGGGGCAATCGAACAGGCCGTTGAT‑TAAGTGCAGGTT‑GCCGGGGCCGCACGATCCG‑GCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGC‑CAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACTCGATGGTGCCCAt                                                                                                                 <  1:1368767/221‑1 (MQ=255)                                                                                                                                  gcggtgaaattcccgtggCCATTTTCGGTGGGGGTAATAGGGGGGGCAA‑C‑CGCCGGGACGGGAGGGTGTTGGCGGTGCCAATGGATCAGGCCGTTGAT‑GAAGT‑CGGGTTGGCCGGGGCCGCACGATCCG‑GCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGC‑CAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACTCGATGGTGCCCATGCGATTAAGACt                                                                                                     <  2:607300/222‑1 (MQ=255)                                                                                                                                   tccgtgaaagaccGGTGCCAATTTCCCTGGGGGGAAATTGGGGGGGAA‑TCCGCACGGCCCGGAACTGGATGGGGTGGGCAATGGAACAGGCCGTTGAT‑TAAGTGCAGGTG‑GCCGGGGCCGCCCGTTCCG‑GCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGT‑CAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACTCGATGGTGCCCATGCGATTAAGACt                                                                                                     <  2:550534/228‑1 (MQ=255)                                                                                                                                               cggtGCAAATTTCGTTGGAGGGAATATGAGCGGCAA‑T‑CGCCAGGACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGAT‑TAAGTGCAGGTT‑GCCGGGGCCGCACGATCCG‑GCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGC‑CAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACTCGATGGTGCCCATGCGATTAAGACTGTCACTAAGAc                                                                                          <  1:1007576/234‑1 (MQ=255)                                                                                                                                                   cgcattttCGATGGGGGGAATAGGGGGGGAAA‑T‑CGACGGGCCGGGAAGGGGATGGGGGTGCCAAGCGAAAGGGCCGTTGAT‑TAAGGGCAGGTT‑GCGGGGCCCGCACGACCCG‑GCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGC‑CAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACTCGATGGTGCCCATGCGATTAAGACTGTCACTAAGACCGTTc                                                                                     <  2:291038/233‑1 (MQ=255)                                                                                                                                                                gAGGGAATATGAGCGGCAA‑T‑CGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGAT‑TAAGTGCAGGTT‑GCCGGGGCCGCACGATCCG‑GCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGC‑CAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGGGGACATCCACTCGATGGGGCCCCTGCGCTTAGAAATGTCCCTAAAACCGCTctcgaagaggcccc                                                                        >  1:291038/1‑225 (MQ=255)                                                                                                                                                                   ggaaatggAGCGGCAA‑T‑CGCAAGTACCGGAACGTGATGGGGGTGGCAATCGAAAAGGCCGTTGAT‑TAAGTGCAGGTT‑GCCGGGGCCGCACGATCCG‑GCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGC‑CAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACTCGATGGTGCCCATGCGATTAAGACTGTCACTAAGACCGTTCAGAGAGTCGCCTGTGa                                                                     <  1:918274/231‑1 (MQ=255)                                                                                                                                                                   gaaaTAGGGGCGCCAC‑T‑CGCAAGTACCGGAACGGTATTGCGGTGGCAATCGAACAGGCCGTTGAT‑TAAGTGCAGGTT‑GCCGGGGCCGCACGATCCG‑GCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGC‑CAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACTCGATGGTGCCCATGCGATTAAGACTGTCACTAAGACCGTTCAGAGAGTCGCCTGTGa                                                                     <  1:455742/236‑1 (MQ=255)                                                                                                                                                                     aaTAGGAGCGGCAA‑T‑CGCCATTACCGGACCGTGATGGCGGTGCAAATCGAGCAGGCCGTTGAT‑TAAGTGCAGGTG‑TCCGGGCCCGCACGAGCCG‑GCGCAGACCGCCAGTTCTCCGCTAAGTTGGGCTTCAGCGC‑CAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACTCGa                                                                                                                           <  1:668625/182‑1 (MQ=255)                                                                                                                                                                              ctgtaaatcgcaatgtcgggAACGTGATTGCGGTGGAAACCGACCAGGCAGTTGCT‑TAAGGGCAGGTT‑GCCGGGGCCGCACGATCCG‑GCCAAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGC‑CAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACTCGATGGTGCCCATGCGATTAAGACTGTCACTAAGACCGTTCAGAGAGTCGCCTGTGACTCCCCAGa                                                            <  1:1031685/220‑1 (MQ=255)                                                                                                                                                                                              ccGGAATGTTGTTGGGGGGGCAATCGACCCGGCCGTTGAT‑TAGGTGCAGGTT‑GCCGGGGCCGCACGATCCG‑GCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGC‑CAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACTCGATGGTGCCCATGCGATTAAGACTGTCACTAAGACCGTTCAGAGAGTCGCCTGTGACTCCCCAGATGCGTTTCACCCCTGcc                                           <  1:523538/239‑1 (MQ=255)                                                                                                                                                                                                     gttgttGCGGTGCCAATGAACCAGGCGTTGGAT‑TGAGGGGCGGTC‑GCGGGGCCCGCACTGCCCG‑GCGCAGCCCGCCAGTCCTCCGCTAGGTTGTGCTTCAGCGC‑CAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACTCGATGGTGCCCATGCGATTAAGACTGTCACTAAGACCGTTCAGAGAGTCGCCTGTGACTCCCCAGATGCGTTTCACCCCTGCCGATTCGa                                    <  2:195886/235‑1 (MQ=255)                                                                                                                                                                                                                                        aaGGGCAGGTT‑GCGGGGGCCGCACGACCCG‑GCCCAGCCCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGC‑CAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACTCGATGGTGCCCATGCGATTAAGACTGTCACTAAGACCGTTCAGAGAGTCGCCTGTGACTCCCCAGATGCGTTTCACCCCTGCCGATTCGAGTGTTTTGGCGATATAAGCTGCAACCgtttgttt  <  1:663232/239‑1 (MQ=255)                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                               CAACCGAAACGCCACGGTTAAGCACCGCTTTGCGCATGGCAATCGCCAGTACTTGTGGGATCTGCTCCGGGCTGGAAACCAGCTCGCAATAGTGACTACATTCGCGGAATAGCTCTTGTGGGTGGGTTTCCTGGAAA‑TAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAA‑T‑CGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGAT‑TAAGTGCAGGTT‑GCCGGGGCCGCACGATCCG‑GCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGC‑CAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACTCGATGGTGCCCATGCGATTAAGACTGTCACTAAGACCGTTCAGAGAGTCGCCTGTGACTCCCCAGATGCGTTTCACCCCTGCCGATTCGAGTGTTTTGGCGATATAAGCTGCAACCGTTTGTTT  >  NC_000913/910581‑911046

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