Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,594,558	G→A	20.0%	L56F (CTC→TTC)	lgoT ←	putative L‑galactonate:H+ symporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,594,558	0	G	A	20.0%	89.6 / 4.7	58	L56F (CTC→TTC)	lgoT	putative L‑galactonate:H+ symporter
Reads supporting (aligned to +/- strand):  ref base G (31/13);  new base A (2/10);  total (33/25)
Fisher's exact test for biased strand distribution p-value = 1.80e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.84e-01

ACACGAACTGCGTAAAGTTGTGCACCATGCCAGACATTGCCTGGAACAGTGACCAGAAGAACATCCCCAGTCCCAGCATCAGGCGTGGGCCTTTACGATCCAATAGTGGGCCGCAAGGAAGTTGCGCAATCCCGTAAGCGAGTGAAAACAC‑GGAGAGCAAAGCGCC‑GATTTCGGTGGCAC‑TTAATCCCAATTCT‑TC‑AC‑G‑AATCGTTAAA‑TTT‑GCTACCG‑ACAGCGAACT‑GCGGT‑CGAG‑ATAATTGATTACCGCCGCAAAAAATAATAACAACATGGC‑GGTGGTTTGAATGCGTTTAATTCGGGTACTGCGTTGAACTAA‑TCC‑ATCTGGTGGCACGGGAATATCTGCCGACGAAGATGGAGTAAATGAAGAACGCGGATCGATGGTGATATTTTCTTTTTCCACGCCGGACTCCTGTATTTATAATTTTAATTGTTTACGTAATT  >  NC_000913/4594321‑4594764                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                   acccGAACTGCGTAAAGTTGTGCACCATGCCAGACATTGCCTGGAACAGTGACCAGAAGAACATCCCCAGTCCCAGCATCAGGCGTGGGCCTTTACGATCCAATAGTGGGCCGCAAGGAAGTTGCGCAATCCCGTAAGCGAGTGAAAACAC‑GGAGAGCAAAGCGCC‑GATTTCGGTGGCAC‑TTAATCCCAATTCT‑TC‑AC‑G‑AATCGTTAAA‑TTT‑GCTACCG‑ACAGCGAACT‑GCGGT‑Cgag‑a                                                                                                                                                                                                                  <  1:56733/236‑1 (MQ=255)  cacGAACTGCGTAAAGTTGTGCACCATGCCAGACATTGCCTGGAACAGTGACCAGAAGAACATCCCCAGTCCCAGCATCAGGCGTGGGCCTTTACGATCCAATAGTGGGCCGCAAGGAAGTTGCGCAATCCCGTAAGCGAGTGAAAACAC‑GGAGAGCAAAGCGCC‑GATTTCGGTGGCAC‑TTAATCCCAATTCT‑TC‑AC‑G‑AATCGTTAAA‑TTT‑GCTACCG‑ACAGCGAACT‑GCGGT‑CGAG‑At                                                                                                                                                                                                                 >  1:490171/1‑239 (MQ=255)       aCTGCGTAAAGTTGTGCACCATGCCAGACATTGCCTGGAACAGTGACCAGAAGAACATCCCCAGTCCCAGCATCAGGCGTGGGCCTTTACGATCCAATAGTGGGCCGCAAGGAAGTTGCGCAATCCCGTAAGCGAGTGAAAACAC‑GGAGAGCAAAGCGCC‑GATTTCGGTGGCAC‑TTAATCCCAATTCT‑TC‑AC‑G‑AATCGTTAAA‑TTT‑GCTACCG‑ACAGCGAACT‑GCGGT‑CGAG‑AAAATTg                                                                                                                                                                                                            >  1:176244/1‑239 (MQ=255)           cGTAAAGTTGTGCACCATGCCAGACATTGCCTGGAACAGTGACCAGAAGAACATCCCCAGTCCCAGCATCAGGCGTGGGCCTTTACGATCCAATAGTGGGCCGCAAGGAAGTTGCGCAATCCCGTAAGCGAGTGAAAACAC‑GGAGAGCAAAGCGCC‑GATTTCGGTGGCAC‑TTAATCCCAATTCT‑TC‑AC‑G‑AATCGTTAAA‑TTT‑GCTACCG‑ACAGCGAACT‑GCGGT‑CGAG‑ATAATTGAtt                                                                                                                                                                                                         >  1:1083042/1‑238 (MQ=255)                     tgCACCATGCCAGACATTGCCTGGAACAGTGACCAGAAGAACATCCCCAGTCCCAGCATCAGGCGTGGGCCTTTACGATCCAATAGAGGGCCGCAAGGAAGTTGCGCAATCCCGTAAGCGAGTGAAAACAC‑GGAGAGCAAAGCGCC‑GATTTCGGTGGCAC‑TTAATCCCAATTCT‑TC‑AC‑G‑AATCGTTAAA‑TTT‑GCTACCG‑ACACCGAACA‑GCGGT‑CAAA‑ATAATTGATTACCGCAGCaaa                                                                                                                                                                                              >  2:245758/1‑239 (MQ=255)                            tGCCAGACATTGCCTGGAACAGTGACCAGAAGAACATCCCCAGTCCCAGCATCAGGCGTGGGCCTTTACGATCCAATAGTGGGCCGCAAGGAAGTTGCGCAATCCCGTAAGCGAGTGAAACCAC‑GGAGAGCAAAGCGCC‑GATTTCGGTGGCAC‑TTAATCCCAATTCT‑TC‑AC‑G‑AATCGTTAAA‑TTT‑GCTACCG‑ACACCAAAAT‑CCGGT‑CGAG‑TAAATTAGTTACCACCGaaaaaaaaaaa                                                                                                                                                                                       >  2:273947/1‑238 (MQ=255)                            tGCCAGACATTGCCTGGAACAGTGACCAGAAGAACATCCCCAGTCCCAGCATCAGGCGTGGGCCTTTACGATCCAATAGTGGGCCGCAAGGAAGTTGCGCAATCCCGTAAGCGAGTGAAAACAC‑GGAGAGCAAAGCGCC‑GATTTCGGTGGCAC‑TTAATCCCAATTCT‑TC‑AC‑G‑AATCGTTAAA‑TTT‑GCTACCG‑ACAGCGAACT‑GCGGT‑CGAG‑ATAATTGATTACCGCCGCAAAAaataat                                                                                                                                                                                       >  1:356530/1‑239 (MQ=255)                            tGCCAGACATTGCCTGGAACAGTGACCAGAAGAACATCCCCAGTCCCAGCATCAGGCGTGGGCCTTTACGATCCAATAGTGGGCCGCAAGGAAGTTGCGCAATCCCGTAAGCGAGTGAAAACAC‑GGAGAGCAAAGCGCC‑GATTTCGGTGGCAC‑TTAATCCCAATTCT‑TC‑AC‑G‑AATCGTTAAA‑TTT‑GCTACCG‑ACAGCGAACT‑GCGGT‑CGAG‑ATAATTGATTACCGCCGCAAAAaataat                                                                                                                                                                                       >  1:928048/1‑239 (MQ=255)                                      tGCCTGGAACAGTGACCAGAAGAACATCCCCAGTCCCAGCATCAGGCGTGGACCTTTACGATCCAATAGTGGGCCGCCAGGAAGTTGCGCAATCCCGTAAGCGAGTGAAAACTC‑GGAGAGCAAACCGCC‑GGTTTCGGTGGCAC‑TTAATCCCAATTCT‑TC‑AC‑G‑CATCGTTAAA‑TTT‑GCACCCG‑ACACCGGACT‑CCGGT‑AGAG‑ATATTTGATTAACACCGCCAAAACAAAAAacaacaatg                                                                                                                                                                             >  2:114880/1‑236 (MQ=255)                                                  tGACCAGAAGAACATCCCCAGTCCCAGCATCAGGCGTGGGCCTTTACGATCCAATAGTGGGCCGCAAGGAAGTTGCGCAATCCCGTAAGCGAGTGAAAACAC‑GGAGAGCAAAGCGCC‑GATTTCGGTGGCAC‑TTAATCCCAATTCT‑TC‑AC‑G‑AATCGTTAAA‑TTT‑GCTACCG‑ACAGCGAACT‑GCGGT‑CGAG‑ATAATTGATTACCGCCGCAaaaaaaaaaaaaaaaaagggggggggggtaa                                                                                                                                                                 >  2:375415/1‑222 (MQ=255)                                                       agaagaACATCCCCAGTCCCAGCATCAGGCGTGGGCCTTTACGATCCAATAGTGGGCCGCAAGGAAGTTGCGCAATCCCGTAAGCGAGTGAAAACAC‑GGAGAGCAAAGCGCC‑GATTTCGGTGGCAC‑TTAATCCCAATTCT‑TC‑AC‑G‑AATCGTTAAA‑TTT‑GCTACCG‑AAGGCGCACC‑GGGGT‑CGAA‑GAAATTGTTTACCACCGCAAAAAAAAAAACAAACATGGC‑GGGGGTTTGGAttggt                                                                                                                                                           >  2:1051655/1‑235 (MQ=255)                                                                       tCCCAGCATCAGGCGTGGGCCTTTACGATCCAATAGTGGGCCGCAAGGAAGTTGCGCAATCCCGTAAGCGAGTGAAAACAC‑GGAGAGCAAAGCGCC‑GATTTCGGTGGCAC‑TTAATCCCAATTCT‑TC‑AC‑G‑AATCGTTAAA‑TTT‑GCTACCG‑ACAGCGAACT‑GCGGT‑CGAG‑ATAATTGATTACCGCCGCAAAAAATAATAACAACATGGC‑GGAGGTTTGAATGCGTTTAATTTGGGTAATGc                                                                                                                                           >  1:1068871/1‑239 (MQ=255)                                                                            gCATCAGGCGTGGGCCTTTACGATCCAATAGTGGGCCGCAAGGAAGTTGCGCAATCCCGTAAGCGAGTGAAAACAC‑GGAGAGCAAAGCGCC‑GATTTCGGTGGCAC‑TTAATCCCAATTCT‑TC‑AC‑G‑AATCGTTAAA‑TTT‑GCTACCG‑ACAGCGAACT‑GCGGT‑CGAG‑ATAATTGATTACCGCCGCAAAAAATAATAACAACATGGC‑GGCGGTTTGAATTCGTTTTATTCGGGTAcctgtctct                                                                                                                                      >  1:174892/1‑231 (MQ=255)                                                                            acATGAGGCGGGGGTCTTTACGACCAAAAGGGGGGCCAAAAGGAAGTCGCGCAATCCGGTAAGAGAGTAAAAACGC‑GGAAAGAAAGGCCCC‑GATTTGGGTGACAT‑TTAACCCAAATTCT‑TC‑AC‑A‑AATCGTAAAA‑TTT‑GCTACCG‑ACAGCGAACT‑GCGGT‑CGAA‑ATATTTGTTTACCGCCGCAAAAAATAATAACAACATGGC‑GGTGGTTTGAATGCGTTTAATTCGGGTACTGCGTTGa                                                                                                                                      <  1:738970/238‑1 (MQ=255)                                                                                     ggggtcttttgagatcaaaaggggggCCGAAAGAAAGTGGCCAAACCCCGTAAGAGA‑TAAAAACCCGGGAGACAAAAGCGCG‑TATGGCGGGGGCAC‑TTAA‑CCCTATTCC‑TCAAC‑A‑ATTCGTTAAT‑TTC‑GCACCCA‑ACAGCGAACT‑GCGGT‑AGAG‑ATATTTGATTCCCGCCGCAAAAAATAATAACAACATGGC‑GGTGGTTTGAATGCGTTTAATTCGGGTACTGCGTTGAACTAA‑Tcc                                                                                                                             <  2:928048/215‑1 (MQ=255)                                                                                       gggCCTTTACGATCCAATAGTGGGCCGCAAGGAAGTTGCGCAATCCCGTAAGCGAGTGAAAACAC‑GGAGAGCAAAGCGCC‑GATTTCGGTGGCAC‑TTAATCCCAATTCT‑TC‑AC‑G‑AATCGTTAAA‑TTT‑GCTACCG‑ACAGCGAACT‑GCGGT‑CGAG‑ATAATTGATTACCGCCGCAAAAAATAATAACAACATGGC‑GGTGGTTTTAATGCGTTTAATTCCGGGACTGCGTTTAACTAA‑Tcaa                                                                                                                            >  1:1006962/1‑235 (MQ=255)                                                                                       gggCCTTTACGATCCAATAGTGGGCCGCAAGGAAGTTGCGCAATCCCGTAAGCGAGTGAAAACAC‑GGAGAGCAAAGCGCC‑GATTTCGGTGGCAC‑TTAATCCCAATTCT‑TC‑AC‑G‑AATCGTTAAA‑TTT‑GCTACCG‑ACAGCGAACT‑GCGGT‑CGAG‑ATAATTGATTACCGCCGCAAAAAATAATAAAAACACGGC‑GGGGGTTTTAAGCGGTTTAATTTGGGGACTCAGTTGtaataaacatct                                                                                                                           >  2:633134/1‑233 (MQ=255)                                                                                        ggCCTTTACGATCCAATAGTGGGCCGCAAGGAAGTTGCGCAATCCCGTAAGCGAGTGAAAACAC‑GGAGAGCAAAGCGCC‑GATTTCGGTGGCAC‑TTAATCCCAATTCT‑TC‑AC‑G‑AATCGTTAAA‑TTT‑GCTACCG‑ACAGCGAACT‑GCGGT‑CGAG‑ATAATTGATTACCGCCGCAAAAAATAATAACAACATGGC‑GGTGGTTTGAATGCGTTTAATTCGGGGACTGCGTTTAACTAA‑ACC‑a                                                                                                                           >  1:241228/1‑236 (MQ=255)                                                                                              ttagagaccatatgggggCGCAAGGGAGTTGCGCAATCCCGTAAGAGAGAGAAAAAAC‑AGAGAAAAAAGCGGC‑GGTTGTGGGGGCACATTAACCCCAATTCT‑CT‑AG‑A‑AATTGTAAAA‑TTT‑GCCACCA‑GCAGGAAAGG‑GGGGT‑GGAA‑AAAATTGATTACCGCCGCAAAAAATAATAACAACATGGC‑GGTGGTTTGAATGCGTTTAATTCGGGTACTGCGTTGAACTAA‑TCC‑ATCtggtgg                                                                                                                   <  2:241228/225‑1 (MQ=255)                                                                                                   tCCAATAGTGGGCCGCAAGGAAGTTGCGCAATCCCGTAAGCGAGTGAAAACAC‑GGAGAGCAAAGCGCC‑GAATTCGGCGGCAC‑TTAATCCCAATTCT‑TC‑AC‑G‑AATCGTTAAA‑TTT‑GCTACCG‑ACAGCGAACT‑CCGGT‑CCAG‑ATAATTGATTACCGCAGCCAAAAATACTAACAAACTGGC‑GGAGTCTTTGATGCGTTTAAATACGGGACTGCGATGAACTaaaaaacatggggggcgggg                                                                                                               >  2:725659/1‑221 (MQ=255)                                                                                                             ggCCGCAAGGAAGTTGCGCAATCCCGTAAGCGAGTGAAAACAC‑GGAGAGCAAAGCGCC‑GATTTCGGTGGCAC‑TTAATCCCAATTCT‑TC‑AC‑G‑AATCGTTAAA‑TTT‑GCTACCG‑ACAGCGAACT‑GCGGT‑CGAG‑ATaa                                                                                                                                                                                                               >  2:12780/1‑134 (MQ=255)                                                                                                             ggCCGCAAGGAAGTTGCGCAATCCCGTAAGCGAGTGAAAACAC‑GGAGAGCAAAGCGCC‑GATTTCGGTGGCAC‑TTAATCCCAATTCT‑TC‑AC‑G‑AATCGTTAAA‑TTT‑GCTACCG‑ACAGCGAACT‑GCGGT‑CGAG‑ATAATTGATTACCGCCGCAAAAAATAATAACAACATGGCGGGGGGTTTGATTGGGGTTATTTCGGGCAAACGGTTGAA‑TCA‑TCCTCTCTGGTGGCCCGGGaaatattcc                                                                                                    >  2:953672/1‑232 (MQ=255)                                                                                                             ggCCGCAAGGAAGTTGCGCAATCCCGTAAGCGAGTGAAAACAC‑GGAGAGCAAAGCGCC‑GATTTCGGTGGCAC‑TTAATCCCAATTCT‑TC‑AC‑G‑AATCGTTAAA‑TTT‑GCTACCG‑ACAGCGAACT‑GCGGT‑CGAG‑ATAATTGATTACCGCCGCAAAAAATAATAAAAACATGGC‑GGGGGTTTGAGTGCGTTTTATTCGGGGACTGCGTTTGACAAA‑TAC‑CACTGGTGGGACggggaaatattgc                                                                                                   >  2:828666/1‑229 (MQ=255)                                                                                                             ggCCGCAAGGAAGTTGCGCAATCCCGTAAGCGAGTGAAAACAC‑GGAGAGCAAAGCGCC‑GATTTCGGTGGCAC‑TTAATCCAAATTCT‑TC‑AC‑G‑AATCGTTAAA‑TTT‑GCTACCG‑ACAGCGAACT‑GCGGT‑CGAG‑ATaa                                                                                                                                                                                                               <  1:12780/134‑1 (MQ=255)                                                                                                                 aaaaGGAAGTGGCGCAACCCCGGAAGGGAGAGAAAACAG‑GGAGAAAAAAGCGCC‑TATTTCGGTGGAAC‑TTAACCCAAATTC‑‑TC‑ACAG‑AATCTTAAAA‑TTT‑GCTACCG‑CGGGCGAACT‑GGGGT‑AGAA‑ATAATTAATTACCGCCGCAAAAAATAATAACAACATGGC‑GGTGGTTTGAATGCGTTTAATTCGGGTACTGCGTTGAACTAA‑TCC‑ATCTGGTGGCACGGGAATATCTGCcgacg                                                                                               <  1:1051655/237‑1 (MQ=255)                                                                                                                  cAAGGAAGTTGCGCAATCCCGTAAGCGAGTGAAAACAC‑GGAGAGCAAAGCGCC‑GATTTCGGTGGCAC‑TTAATCCCAATTCT‑TC‑AC‑G‑AATCGTTAAA‑TTT‑GCTACCG‑ACAGCGAACT‑GCGGT‑CGAG‑CTAATTGATTACCGCCGCAAAAAATAATAACAACATGGC‑GGTGGGTTTGATTCGGTTTATTCTGGACATCCGTTGAATCAA‑TCC‑CACTGGTGGCAAGGGagaatcggcgccgga                                                                                              >  2:223142/1‑225 (MQ=255)                                                                                                                  cAAGGAAGTTGCGCAATCCCGTAAGCGAGTGAAAACAC‑GGAGAGCAAAGCGCC‑GATTTCGGTGGCAC‑TTAATCCCAATTCT‑TC‑AC‑G‑AATCGTTAAA‑TTT‑GCTACCG‑ACAGCGAACT‑GCGGT‑CGAG‑ATAATTGATTACCGCCGCAAAAAATAATAACAACATGGC‑GGGGGGTTTGATGCTGTTAATTCGGGGTCCTCGTTTGACTAA‑ACC‑ATTCGGGGGGACggggtaatttcggcaccg                                                                                              >  2:955049/1‑224 (MQ=255)                                                                                                                           tgccaaaccccGTAGTCGAGACAAAGCAC‑GGGAAACAAAGCGCC‑TATTTCGGAGTCAC‑TTAATCCCAATTCT‑GT‑AC‑G‑AATTCTTAAA‑TTT‑TCTTCCA‑ACAGCAAATT‑GCTGT‑‑GAGAATAATTGATAACCGCCGCAAAAAATAATAACAACATGGC‑GGTGGTTTGAATGCGTTTAATTCGGGTACTGCGTTGAACTAA‑TCC‑ATCTGGTGGCACGGGAATATCTGCCGACGAAGATGGAg                                                                                      <  1:114880/230‑1 (MQ=255)                                                                                                                                aaTCCGGAAAGCGAGGAAAAACGG‑GTAGAGC‑AAGCGCCAGTTTCCGGTGCCAC‑TTAACCCCATTTCT‑TC‑AC‑G‑AACCGTAAAT‑TTT‑GCTACCA‑ACAGCAAACT‑GGGGT‑CGAG‑ATAATTGATTACCGCCGCAAAAAATAATAACAACATGGC‑GGTGGTTTGAATGCGTTTAATTCGGGTACTGCGTTGAACTAA‑TCC‑ATCTGGTGGCACGGGAATATCTGCCGACGAAGATGGAGTAAATg                                                                                <  1:160700/239‑1 (MQ=255)                                                                                                                                    ccGTAAGCGAGTGAAAACAC‑GGAGAGCAAAGCGCC‑GATTTCGGTGGCAC‑TTAATCCCAATTCT‑TC‑AC‑G‑AATCGTTAAA‑TTT‑GCTACCG‑ACAGCGAACT‑GCGGT‑CGAG‑ATAATTGATTACCGCCGCAAAAAATAATAACAACATGGC‑GGTGGTTTTAATGCGTTTAATTCGGGGACTGCGTTGAACTAA‑TCC‑ATCTTGTGGCACGGGAAAATCTGCCCACGAAGATGGGGGAAATgagaa                                                                            >  1:392177/1‑236 (MQ=255)                                                                                                                                      ctaaaggAGGGAAAACAC‑GGAAAACGAAGCCCC‑TGTTTCGGTGGCTC‑TTAACCCCAATTTC‑TC‑AA‑A‑AATTGTTAAATTTT‑GCCAACG‑ACAGCGAACG‑GGGGT‑CGAA‑ATATTATATTCCCGCCGCAAAAAATAATAACAACATGGC‑GGTGGTTTGAATGCGTTTAATTCGGGTACTGCGTTGAACTAA‑TCC‑ATCTGGTGGCACGGGAATATCTGCCGACGAAGATGGAGTAAATgaagaa                                                                           <  2:392177/233‑1 (MQ=255)                                                                                                                                       tAAGCGAGTGAAAACAC‑GGAGAGCAAAGCGCC‑GATTTCGGTGGCAC‑TTAATCCCAATTCT‑TC‑AC‑G‑AATCGTTAAA‑TTT‑GCTACCG‑ACAGCGAACT‑GCGGT‑CGAG‑ATAATTGATTACCGCCGCAAAAAATAATAACAACATGGC‑GGCGGTTTTAATGGGTTTAATTCGGGTACTGCGTTGAACTAA‑TCC‑ATCTGGTGGCACGGGAATATCTGCCGACGAAGCTGGGGTAAATgaagacag                                                                         >  1:12736/1‑236 (MQ=255)                                                                                                                                              gTGAAAACAC‑GGAGAGCAAAGCGCC‑GATTTCGGTGGCAC‑TTAATCCCAATTCT‑TC‑AC‑G‑AATCGTTAAA‑TTT‑GCTACCG‑ACAGCGAACT‑GCGGT‑CGAG‑ATAATTGATTGCCGCCGCAAAAAATAATAACAACATGGG‑GGGGGGTTGCATGCGTTTTATTTCGGAAATGCGTTGAAATGA‑TCC‑CTCTGGGGGGAGGGGAAAATATGCCGACGAGAATGGAGGAAAAGAAGAAAGCgggggg                                                                  >  2:191090/1‑235 (MQ=255)                                                                                                                                               tGAAAACAC‑GGAGAGCAAAGCGCC‑GATTTCGGTGGCAC‑TTAATCCCAATTCT‑TC‑AC‑G‑AATCGTTAAA‑TTT‑GCTACCG‑ACAGCGAACT‑GCGGT‑CGAG‑ATAATTGATTACCGCCGCAAAAAATAATAACAACATGGC‑GGTGGTTTGAATGGGTTTAATTCGGGGACTGGGTTGAACTAA‑TCC‑ATCTGTGGGCACGGGAATATTTGCCGCCGAAGATGGGGTAAATgaagaaacgggaccgt                                                                 >  1:283174/1‑229 (MQ=255)                                                                                                                                                            gagCAAAGCCCC‑GATTTCGGAGGCAC‑TTAATCCCAATTCT‑TC‑AC‑G‑AATCGTTAAA‑TTT‑GCTACCG‑ACACCGAACT‑GCGGT‑CGAG‑ATAAATGATTACCGCCGCAAAAAATAATAACAACATGGC‑GGTAGTTTGTATGAGTTTTATTTAGGAACTGCCTTGAACTAACTCC‑TCCCGGTGG‑ACGGGAAGAGCTTGCGAAGCAGATGGATTAAATtaaaaaacgggttcacggggtagttttt                                                     >  2:1025872/1‑217 (MQ=255)                                                                                                                                                             agCAAAGCGCC‑GATTTCGGTGGCAC‑TTAATCCCAATTCT‑TC‑AC‑G‑AATCGTTAAA‑TTT‑GCTACCG‑ACAGCGAACT‑GCGGT‑CGAG‑ATAATTGATTACCGCCGCAAAAAATAATAACAACATGGC‑GGTGGGTTGGATGCGTTTAATTTGGGTACTGCGTTGAACTAA‑TCC‑ATCTTGTGGACCGGGAATATCTGCCGACCAAGATGGGGGAAATGAGAAACGCGGATCGATgggggtacgggt                                                    >  1:588963/1‑229 (MQ=255)                                                                                                                                                                 aaGCGCC‑GATTTCGGTGGCAC‑TTAATCCCAATTCT‑TC‑AC‑G‑AATCGTTAAA‑TTT‑GCTACCG‑ACAGCGAACT‑GCGGT‑CGAG‑ATAATTGATTACCGCCGCAAAAAATAATAACAACATGGC‑GGTGGTTTGAATGCGTTTAATTCGGGGAATGCGTTGAACTAA‑TCC‑ATCTGGTGGCACGGGAATATTTGCCGCCGAAGATGGAGTAAATGAAGAACGCGGGTCGATGGTGATATtttttttt                                                >  2:681649/1‑239 (MQ=255)                                                                                                                                                                       cagattcggggcgCAT‑TTAATCCAAAGTCT‑TC‑AC‑G‑AATCGTTAAA‑TTT‑GCTACCA‑ACAGCGAACG‑GGGGT‑CGAA‑ATAATGGATGACCGCCGCAAAAAATAATAACAACAGGGC‑GGTGGTTTGAATGCGTTTAATTCGGGTACTGCGTTGAACTAA‑TCC‑ATCTGGTGGCACGGGAATATCTGCCGACGAAGATGGAGTAAATGAAGAACGCGGATCGATGGTGATATTTTCTTTTTCCa                                            <  1:135814/226‑1 (MQ=255)                                                                                                                                                                           tttCGGGGCCAC‑TTAACCCAAATTCG‑TC‑AA‑A‑AATCGTAAAA‑TTT‑ATAACAA‑AAGGAGAACG‑GCTGT‑GGAG‑TAAGTTTATTACCGCCGCAAAAAATAATAACAACATGGC‑GGTGGTTTGAATGCGTTTAATTCGGGTACTGCGTTGAACTAA‑TCC‑ATCTGGTGGCACGGGAATATCTGCCGACGAAGATGGAGTAAATGAAGAACGCGGATCGATGGTGATATTTTCTTTTTCCACGCCg                                       <  2:1068871/239‑1 (MQ=255)                                                                                                                                                                           ttgggggggCAT‑TAAATCCCAATTCA‑TC‑AC‑G‑AATCTATAAT‑TTG‑CCTCCCG‑ACAGCGAACT‑GCGGT‑GGAG‑ATAATTGATTACCGCCGCAAAAAATAATAACAACATGGC‑GGTGGTTTGAATGCGTTTAATTCGGGTACTGCGTTGAACTAA‑TCC‑ATCTGGTGGCACGGGAATATCTGCCGACGAAGATGGAGTAAATGAAGAACGCGGATCGATGGTGATATTTTCTTTTTCCACGCCg                                       <  1:681649/232‑1 (MQ=255)                                                                                                                                                                           ttgcGGGGGCAT‑TTAACCCCATTTCT‑TC‑AA‑A‑AATTGTTAAT‑TTT‑CCTACGC‑ACGGCGAACT‑GGGGT‑CGAG‑ATAAGTGATAACCGCCGCAAAAAATAATAACAACATGGC‑GGTGGTTTGAATGCGTTTAATTCGGGTACTGCGTTGAACTAA‑TCC‑ATCTGGTGGCACGGGAATATCTGCCGACGAAGATGGAGTAAATGAAGAACGCGGATCGATGGTGATATTTTCTTTTTCCACGCCg                                       <  1:273947/236‑1 (MQ=255)                                                                                                                                                                           tgtggggtgCAT‑TTAACCGCAATTTC‑AC‑AC‑G‑AATGGTTAAT‑TTT‑GCTACCG‑GAAGAGAACT‑GGGGG‑CGAG‑TAAATTGATTACCGCCGCAAAAAATAATAACAACATGGC‑GGTGGTTTGAATGCGTTTAATTCGGGTACTGCGTTGAACTAA‑TCC‑ATCTGGTGGCACGGGAATATCTGCCGACGAAGATGGAGTAAATGAAGAACGCGGATCGATGGTGATATTTTCTTTTTCCACGCCg                                       <  1:1069075/231‑1 (MQ=255)                                                                                                                                                                            tttttaaaagacacagccccATTTCT‑TC‑AC‑G‑AATCGTTAAA‑TTT‑GCTACCG‑ACAGCGACCT‑GCGGT‑CAAG‑ATAATTTATTACCGCCGCAAAAAATAATAACAACATGGC‑GGTGGTTTGAATGCGTTTAATTCGGGTACTGCGTTGAACTAA‑TCC‑ATCTGGTGGCACGGGAATATCTGCCGACGAAGATGGAGTAAATGAAGAACGCGGATCGATGGTGATATTTTCTTTTTCCACGCCg                                       <  1:269285/222‑1 (MQ=255)                                                                                                                                                                            ttggggtgCAC‑TTAAACCCAATTCTCAC‑AA‑G‑AATCGTAAAA‑TTT‑GCCACCG‑AAAGAGAACG‑GCGGT‑GGAA‑ATAATTGATTACCGCCGCAAAAAATAATAACAACATGGC‑GGTGGTTTGAATGCGTTTAATTCGGGTACTGCGTTGAACTAA‑TCC‑ATCTGGTGGCACGGGAATATCTGCCGACGAAGATGGAGTAAATGAAGAACGCGGATCGATGGTGATATTTTCTTTTTCCACGCCg                                       <  1:223142/232‑1 (MQ=255)                                                                                                                                                                                   tgcttttAACCCCTATTCT‑TA‑AA‑G‑AATTTTTAAA‑TTT‑GC‑ACCGCAAAGAGAGCG‑GGGGT‑GGAG‑ATAATTGATTACCGCCGCAAAAAATAATAACAACATGGC‑GGTGGTTTGAATGCGTTTAATTCGGGTACTGCGTTGAACTAA‑TCC‑ATCTGGTGGCACGGGAATATCTGCCGACGAAGATGGAGTAAATGAAGAACGCGGATCGATGGTGATATTTTCTTTTTCCACGCCGGACTcc                                 <  2:1052143/233‑1 (MQ=255)                                                                                                                                                                                    ccctttatcaccATTTTT‑TA‑AC‑A‑AATCGTAAAT‑TTT‑CCTACCC‑ATAGCAACGT‑GGGGG‑CGAG‑ATATTTGAGTCCCGCCGCAAAAAATAATAACAACATGGC‑GGTGGTTTGAATGCGTTTAATTCGGGTACTGCGTTGAACTAA‑TCC‑ATCTGGTGGCACGGGAATATCTGCCGACGAAGATGGAGTAAATGAAGAACGCGGATCGATGGTGATATTTTCTTTTTCCACGCCGGACTCCt                                <  2:1006962/228‑1 (MQ=255)                                                                                                                                                                                          aaTCCCAATTCT‑TC‑AC‑G‑AATCGTTAAA‑TTT‑GCTACCG‑ACAGCGAACT‑GCGGT‑CGAG‑ATAATTGATTACCGCCGCAAAAAATAATAACAACATGGC‑GGTGGTTTGAATGCGTTTAATTCGGGTACTGCGTTGAACTAC‑TCC‑ATCTTGTGGCCCGGGAATATTTGCCGACGAAAGTGGAGTAAATGAAGAACGCGGATCGATGGTGGTATTTTTTTTTTTCACCCCGGACTCCTGGATTta                         >  1:276893/1‑239 (MQ=255)                                                                                                                                                                                           acccaaaTTCT‑GC‑AC‑A‑AATCGTAAAT‑TTT‑GCTCCCG‑‑CAGCGGATTGGCGGT‑GGAA‑AAATTTGTTTACCGCCGCAAAAAATAATAACAACATGGC‑GGTGGTTTGAATGCGTTTAATTCGGGTACTGCGTTGAACTAA‑TCC‑ATCTGGTGGCACGGGAATATCTGCCGACGAAGATGGAGTAAATGAAGAACGCGGATCGATGGTGATATTTTCTTTTTCCACGCCGGACTCCTGTATTtat                        <  1:948852/234‑1 (MQ=255)                                                                                                                                                                                           accacAATTCC‑TC‑AC‑G‑AATCTTAAAA‑TTG‑GATACGG‑ACAGGAAAGG‑GGGGT‑GGAA‑ATATTTGTTTACCGCCGCAAAAAATAATAACAACATGGC‑GGTGGTTTGAATGCGTTTAATTCGGGTACTGCGTTGAACTAA‑TCC‑ATCTGGTGGCACGGGAATATCTGCCGACGAAGATGGAGTAAATGAAGAACGCGGATCGATGGTGATATTTTCTTTTTCCACGCCGGACTCCTGTATTtat                        <  1:315158/235‑1 (MQ=255)                                                                                                                                                                                           aacccAATTCT‑TC‑AA‑GAAATCTTAAAA‑TTT‑CCT‑CCC‑AAAGGGAAAC‑GGGGTAGGAA‑ATAATTATTTACCGCCGCAAAAAATAATAACAACATGGC‑GGTGGTTTGAATGCGTTTAATTCGGGTACTGCGTTGAACTAA‑TCC‑ATCTGGTGGCACGGGAATATCTGCCGACGAAGATGGAGTAAATGAAGAACGCGGATCGATGGTGATATTTTCTTTTTCCACGCCGGACTCCTGTATTta                         <  1:728662/237‑1 (MQ=255)                                                                                                                                                                                                aaTTCT‑TC‑AC‑G‑AATCGTTAAA‑TTT‑GCTACCG‑ACAGCGAACT‑GCGGT‑CGAG‑ATAATTGATTACCGCCGC‑AAAAATAATAACAACATGGC‑GGTGGGTTGAATGCGTTTAAATCCGGTACTGCGTTGaa                                                                                                                                     >  1:992642/1‑126 (MQ=255)                                                                                                                                                                                                  ttct‑tc‑AC‑G‑AATCGTTAAA‑TTT‑GCTACCG‑ACAGCGAACT‑GCGGT‑CGAG‑ATAATTGATTACCGCCGCAAAAAATAATAACAACATGGC‑GGGGGTTTGAATGCGTTTAATTCGGGTACTGCGTTGAACTAA‑t                                                                                                                               >  1:512638/1‑130 (MQ=255)                                                                                                                                                                                                    ttgtc‑AC‑G‑AATTGTTAAA‑TTT‑ACACACA‑AGAGCGAATG‑GCGGT‑CGAG‑ATAATTGATTACCGCCGCAAAAAATAATAACAACATGGC‑GGTGGTTTGAATGCGTTTAATTCGGGTACTGCGTTGAACTAA‑TCC‑ATCTGGTGGCACGGGAATATCTGCCGACGAAGATGGAGTAAATGAAGAACGCGGATCGATGGTGATATTTTCTTTTTCCACGCCGGACTCCTGTATTTATAATTTTaa                <  2:356530/236‑1 (MQ=255)                                                                                                                                                                                                        c‑AC‑G‑AATCGTTAAA‑TTT‑GCTACCG‑ACAGCGAACT‑GCGGT‑CGAG‑ATAATTGATTACCGCCGCAAAAAATAATAACAACATGGC‑GGTGGTTTGAATGCGTTTAATTCGGGGACTGCGTTGAACTAA‑TCC‑ATCTGGGGGCACGGGAATATCTGCCGACGAAGATGGAGTAAATGAAGAAAGCGGGTCGATGGGGATATTTTCTTTTTCCACGCCGGGACCCTGTATTTTTAATTTTAATTGt            >  1:653175/1‑239 (MQ=255)                                                                                                                                                                                                        c‑AC‑G‑AATCGTTAAA‑TTT‑GCTACCG‑ACAGCGAACT‑GCGGT‑CGAG‑ATAATTGATTACCGCCGCAAAAAATAATAACAACATGGC‑GGTGGGTTTAATGAGTTTAATTCGGGTAATGCATTTAACTAA‑TCC‑CTCTGGGGGCACGGGAATATCTGGCGGCGAAGATGGAGAAAAAGAAAGAAGGGGAGGGGAGGGGAATTTTTTTTTTTCCCCGACGGACACCCGtgtcggtctcttataaaca            >  2:269285/1‑221 (MQ=255)                                                                                                                                                                                                             g‑AATCGTAAAA‑TTT‑CCTCCCA‑GCGGAGAACT‑GCGGG‑CGAG‑ATAATGGATTACCGCCGCAAAAAATAATAACAACATGGC‑GGTGGTTTGAATGCGTTTAATTCGGGTACTGCGTTGAACTAA‑TCC‑ATCTGGTGGCACGGGAATATCTGCCGACGAAGATGGAGTAAATGAAGAACGCGGATCGATGGTGATATTTTCTTTTTCCACGCCGGACTCCTGTATTTATAATTTTAATTGTTTa         <  1:245758/239‑1 (MQ=255)                                                                                                                                                                                                             cgcgaagctttaagtgtccaaCCG‑ACAGCGATGG‑GCGAT‑GAAT‑ATAATTACTTCCCGCCGTAAAAAATAATAACAACATGGC‑GGTGGTTTGAATGCGTTTAATTCGGGTACTGCGTTGAACTAA‑TCC‑ATCTGGTGGCACGGGAATATCTGCCGACGAAGATGGAGTAAATGAAGAACGCGGATCGATGGTGATATTTTCTTTTTCCACGCCGGACTCCTGTATTTATAATTTTAATTGt            <  1:514817/219‑1 (MQ=255)                                                                                                                                                                                                               aaagtgtAAA‑TTTCGCAACCA‑CCAGCGAACG‑GGGGT‑CAAT‑ATATTTTATTACCGCCGCAAAAAATAATAACAACATGGC‑GGTGGTTTGAATGCGTTTAATTCGGGTACTGCGTTGAACTAA‑TCC‑ATCTGGTGGCACGGGAATATCTGCCGACGAAGATGGAGTAAATGAAGAACGCGGATCGATGGTGATATTTTCTTTTTCCACGCCGGACTCCTGTATTTATAATTTTAATTGTTTa         <  1:725659/233‑1 (MQ=255)                                                                                                                                                                                                                      aaatttt‑GCCACCG‑ACAGCAACCG‑GGGGT‑GAAA‑AAATTTGATACCCGCCGCAAAAAATAATAACAACATGGC‑GGTGGTTTGAATGCGTTTAATTCGGGTACTGCGTTGAACTAA‑TCC‑ATCTGGTGGCACGGGAATATCTGCCGACGAAGATGGAGTAAATGAAGAACGCGGATCGATGGTGATATTTTCTTTTTCCACGCCGGACTCCTGTATTTATAATTTTAATTGTTTACGTAAtt  <  1:375415/235‑1 (MQ=255)                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                   ACACGAACTGCGTAAAGTTGTGCACCATGCCAGACATTGCCTGGAACAGTGACCAGAAGAACATCCCCAGTCCCAGCATCAGGCGTGGGCCTTTACGATCCAATAGTGGGCCGCAAGGAAGTTGCGCAATCCCGTAAGCGAGTGAAAACAC‑GGAGAGCAAAGCGCC‑GATTTCGGTGGCAC‑TTAATCCCAATTCT‑TC‑AC‑G‑AATCGTTAAA‑TTT‑GCTACCG‑ACAGCGAACT‑GCGGT‑CGAG‑ATAATTGATTACCGCCGCAAAAAATAATAACAACATGGC‑GGTGGTTTGAATGCGTTTAATTCGGGTACTGCGTTGAACTAA‑TCC‑ATCTGGTGGCACGGGAATATCTGCCGACGAAGATGGAGTAAATGAAGAACGCGGATCGATGGTGATATTTTCTTTTTCCACGCCGGACTCCTGTATTTATAATTTTAATTGTTTACGTAATT  >  NC_000913/4594321‑4594764

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