Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,093,146	G→T	17.1%	Q98K (CAA→AAA)	rhaM ←	L‑rhamnose mutarotase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,093,146	0	G	T	17.1%	134.4 / 5.2	64	Q98K (CAA→AAA)	rhaM	L‑rhamnose mutarotase
Reads supporting (aligned to +/- strand):  ref base G (22/27);  new base T (10/2);  total (35/29)
Fisher's exact test for biased strand distribution p-value = 2.37e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

TGAGACAGAGTGAAAGGTCAGATTTGCGGAGTAATGCACATAATGGTTATTTAAATAAACCACATGAATCATTAATGGGTTATTCATTATTTTTGTGATTTACTCACAATGTATCA‑GGGAAATATAACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGC‑AGC‑TCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATA‑ACATCGGTCAT‑ATA‑TTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGC‑ATTCCAGCGTTCTTCAGATTCAATCTCTACCATGGCAAACAG‑CAGATTACGCGCTTTGTCGAGATAGATGGCGTAGTTATGCGCACCGTGA‑GATTTCAGCACTGCTTCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGA  >  NC_000913/4092913‑4093384                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                       tGAGACAGAGTGAAAGGTCAGATTTGCGGAATAATGCACATAATGGTTATTTAAATAAACCACATGAATCATTAATGGGTTATTCATTATTTTTGTGATTTACTCACAATGTATCA‑GGGAAATATAACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGC‑AGC‑t                                                                                                                                                                                                                                                 <  1:1292905/239‑1 (MQ=255)         agTGAAAGGTCAGATTTGCGGAGTAATGCACATAATGGTTATTTAAATAAACCACATGAATCATTAATGGGTTATTCATTATTTTTGTGATTTACTCACAATGTATCA‑GGGAAATATAACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGGAAAAAAATCCTTTC‑GGC‑CCGCGACta                                                                                                                                                                                                                                         >  2:633356/1‑238 (MQ=255)                  tCAGATTTGCTGAGTAATGCACTTAATGGTTAGTTAAATAAACCACATGAATCATTAATGGGTTATTCATTATTTTTGTGATTTACTCACAATGTATCA‑GGGAAATATAACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGC‑AGC‑TCGCTACTCACCGGGc                                                                                                                                                                                                                                  <  2:51016/237‑1 (MQ=255)                   cAGATTTGCGGAGTAATGCACATAATGGTTATTTAAATAAACCACATGAATCATTAATGGGTTATTCATTATTTTTGTGATTTACTCACAATGTATCA‑GGGAAATATAACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAAACATTCTTTC‑AGcccccaacctccccggccgt                                                                                                                                                                                                                                >  2:617448/1‑220 (MQ=255)                   cAGATTTGCGGAGTAATGCACATAATGGTTATTTAAATAAACCACATGAATCATTAATGGGTTATTCATTATTTTTGTGATTTACTCACAATGTATCA‑GGGAAATATAACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGGAAAAAACTTCTTTC‑AAC‑TCGCTACTCCCCCGGgtgt                                                                                                                                                                                                                               >  2:1047193/1‑239 (MQ=255)                    aGATTTGCGGAGTAATGCACATAATGGTTATTTAAATAAACCACATGAATCATTAATGGGTTATTCATTATTTTTGTGATTTACTCACAATGTATCA‑GGGAAATATAACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAAAAACTCTTGC‑AGC‑TCGCCACTCCCCCGGCTggt                                                                                                                                                                                                                              >  2:749247/1‑237 (MQ=255)                       tttgggggatAATGGACATAATGGTTATTTAAATAAACCACATGAATCATTAAGTGGTTATTTATTATTTTGTTGATTTACTCACAATGTATCA‑GGGAAATATAACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGC‑AGC‑TCGCTACTCACCGGGCTGTTATc                                                                                                                                                                                                                           <  1:186277/230‑1 (MQ=255)                                          aaTGGTTATTTAAATAAACCACATTAATCATTAATGGGTTATTCATTATTTTTGTGATTTACTCACAATGTATCA‑GGGAAATATAACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGC‑AGC‑TCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATA‑ACa                                                                                                                                                                                                       <  2:1232572/239‑1 (MQ=255)                                          aaTGGTTATTTAAATAAACCACATGAATCATTAATGGGTTATTCATTATTTTTGTGATTTACTCACAATGTATCA‑GCGAAATATAACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAGATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTTCGTCCGGGTAAACACGTATTCC‑AGC‑GCGCGACTACCCCGGCTGTTATCCCGGGTCCCGGGCATA‑Ac                                                                                                                                                                                                        >  2:776966/1‑238 (MQ=255)                                           tttgTTTTTTAATTTACCCACATTAATCATTAATTGTTTTTTCATTATTTTTGTGATTTACTCACAATGTATCA‑GGGAAATATAACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGC‑AGC‑TCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATA‑ACAt                                                                                                                                                                                                      <  1:1252135/236‑1 (MQ=255)                                           aTGGTTATTTAAATAAACCACATGAATCATTAATGGGTTATTCATTAGTTTTGTGATTTACCCACAATGTATCA‑GGGAAATATAACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGC‑AGC‑TCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATA‑ACAt                                                                                                                                                                                                      <  2:574771/239‑1 (MQ=255)                                             ggTTATTTAAATAAACCACATGAATCATTAATGGGTTATTCATTATTTTTGTGATTTACTCACAATGTATCA‑GGGAAATATAACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGGAAAACACTTCTTTC‑CAC‑TCCCTACTCACCCGGCTGTTATTCCGGGTAGCGGGGATa‑aaaacg                                                                                                                                                                                                    >  2:1183169/1‑234 (MQ=255)                                                 aTTTAAATAAACCACATGAATCATTAATGGGTTATTCATTATTTTTGTGATTTACTCACAATGTATCA‑GGGAAATATAACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAAACATTCTTGC‑AGC‑TCGCTACTCAACGGGCCGGTCTCCCGGGTCACGTGCATA‑AAATcaggcc                                                                                                                                                                                                >  2:1144297/1‑234 (MQ=255)                                                  tttAAATTATCCATTGGAATTATTAGTGTGTTTTTCATTCTTTTTGTGATTTACTAATAATGTAACA‑GGGAAATAGAACTTACCGTCTGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCTTCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACGTCTTTC‑AGG‑TCGCTACTCACCGGTCTGTTATCCGAGTTATCGGGCATA‑ACATCGGTCat                                                                                                                                                                                               <  1:40273/239‑1 (MQ=255)                                                  tttAAATAAACCCAAGTAATCATTAATGGGTTATTTATTATTTTTTGGATTTACTTACAAAGTTTCA‑TGGAAATATAACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCCCATTGGGCAATTACGGCAGGTAAAACACTTCTTGC‑AGC‑TCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATA‑ACATCGGTCat                                                                                                                                                                                               <  1:1047193/239‑1 (MQ=255)                                                  tttAAATAAACCACATGAATCATTAATGGGTTATTCATTATTTTTGTGATTTACTCACAATGTATCA‑GGGAAATATAACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTTAAAAACTTCTTTC‑AGC‑CCGCCACTCACCCGGCTGGTTTCCGGGTTAGCGGGCCTa‑aaaatgggcca                                                                                                                                                                                               >  2:734921/1‑229 (MQ=255)                                                  tttAAATAAACCACATGAATCATTAATGGGTTATTCATTATTTTTGTGATTTACTCACAATGTATCA‑GGGAAATATAACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTAC‑AGC‑TCCCTACTCAACCGGCTGTTTTCCGGGGTAGCGGGCCTC‑AAATAGGTcct                                                                                                                                                                                               >  2:76247/1‑237 (MQ=255)                                                  tttAAATAAACCACATGAATCATTAATGGGTTATTCATTATTTTTGTGATTTACTCACAATGTATCA‑GGGAAATATAACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAAAACTTCTTGG‑AGC‑TCGCTACTCACCGGGCTGTTTTCCCGGGTAAGGGGGATA‑AAATCGGGCat                                                                                                                                                                                               >  2:854905/1‑239 (MQ=255)                                                    taaataaaCCACATGAATCATTAATGGGTTATTCATTATTTTTGTGATTTACTCACAATGTATCA‑GGGAAATATAACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGACAGGTAAAAAACATCTTAC‑AGC‑ACGTTACTACCCGGCTTGTTATCCAGGTAAggaggggaaaaaacgggacaaa                                                                                                                                                                                              >  1:574771/1‑218 (MQ=255)                                                                aTGAATCATTAATGGGTTATTCATTATTTTTGTGATTTACTCACAATGTATCA‑GGGAAATATAACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGGAAAAAACCTCTTGC‑AGC‑CCGCCACTTACCCGGCTGTTTTCCCGGTTTCGGGGCATA‑AAATCCGGCCA‑ATA‑TTTCCCCCAAc                                                                                                                                                                               >  2:618752/1‑239 (MQ=255)                                                                       aTTAATGGGTTATTCATTATTTTTGTGATTTACTCACAATGTATCA‑GGGAAATATAACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGC‑AGC‑TCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATA‑ACATCGGTCAT‑ATA‑TTTCCACCAACGTTGGCa                                                                                                                                                                        <  2:125437/239‑1 (MQ=255)                                                                       aTTAATGGGTTATTCATTATTTTTGTGATTTACTCACAATGTATCA‑GGGAAATATAACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGGAAAAAACTTCTTGC‑AAA‑TCGCTACCCCCCCGGCTGTTTTCCGGGTTTACCGGCAAA‑AACTCGGGCAA‑TAA‑TTTTCACCCACCCTGGCa                                                                                                                                                                        >  2:486171/1‑239 (MQ=255)                                                                              ggTTATTCATTATTTTTGTGATTTACTCACAATGTATCA‑GGGAAATATAACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGTCAGGGAAAACACTTTTTTC‑CCG‑TCCCCATCTACCCGGCTGTTATCCTGGTTAGCGCGCATACACA‑CCGTCCT‑ATT‑TTTTCACCACCGTTGGCaaaactca                                                                                                                                                                 >  2:651974/1‑234 (MQ=255)                                                                                          tttttGTGATTTACTCACAATGTATCA‑GGGAAATATAACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAAAAATTCTTGC‑AGA‑TAGCTAATCACCGGGCTGGTTTACGGGTTAACGGGCATG‑AAATACGGCAT‑AAA‑TTTTCACCCACCTTgggacaaaaacagttctgcaacg                                                                                                                                                     >  1:23192/1‑218 (MQ=255)                                                                                                  aTTTACTCATAATGT‑TCATGGGAATAATAAATTACCGGGAGATGTAAATTGTTTAATTTTTCAAATTGAAATTTAAAATTTTGTGCCGGAGGCATCTCTGGCACATTGGGCACTGTCGGCAGGTAAAACACTTCTTGC‑AGC‑TCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATA‑ACATCGGTCAT‑ATA‑TTTCCGACAACGTTGGCAAACATCAGTGCTGGCTACCGC‑ATTCCa                                                                                                                                            <  1:76247/239‑1 (MQ=255)                                                                                                     tataaaagacagaTCA‑GGGAAATATAACTTACCGGGAGATGCAATTGGTTTAAATTTTCAAATCGAAATTTAAAATCTTGTGCCGGAGTCATCTCTGGCCAATTGGGCAATTACGGCAGGTAAAACTCTTCTTGC‑AGC‑TCGCTACTCACCGGGCTGGTATCCGGGTTAGCGGGCATA‑ACATCGGTCAT‑ATA‑TTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGC‑ATTCCAGCg                                                                                                                                         <  1:820152/227‑1 (MQ=255)                                                                                                         caaaaTGTATCA‑GGGAAATATAACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGC‑AGC‑TCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATA‑ACATCGGTCAT‑ATA‑TTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGC‑ATTCCAGCGttct                                                                                                                                     <  2:23192/236‑1 (MQ=255)                                                                                                          acaATGTATCA‑GGGAAATATAACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGC‑AAC‑TCGCTACTCACCGGGCTGTTATCCGGGCTAGCGGGCATT‑ACATCGGTCATAATA‑TTT‑CACCAACGTTGGCAAAAATCAGGTATGGCACCCGC‑ATTCCACCGttct                                                                                                                                     >  2:683183/1‑238 (MQ=255)                                                                                                            aaTGTATCA‑GGGAAATATAACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAACAAACTCTTTC‑AGT‑GTCCCAGTAACCGGGCTATTCAACGAGGTAGCGGGAATA‑ACATCAGGCAT‑AAATTTTACCACAACG‑TGGCCAACACAAATACAGGCAACCGA‑ACTCAAGCGttcgtac                                                                                                                                  >  1:1006009/1‑235 (MQ=255)                                                                                                               gTATCA‑GGGAAATATAACTTACCGGGAAATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGC‑AGC‑TCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATA‑ACATCGGTCAT‑ATA‑TTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGC‑ATTCCAGCGTTCTTCAGAt                                                                                                                               <  2:547735/239‑1 (MQ=255)                                                                                                                  tCA‑GGGAAATATAACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTGAAAAACTTCTTGC‑CGC‑TCCCTACTCTCCTGGCTGTTTTCCCGGTTTGCCGGCCTG‑ACATCGGTCCT‑ATA‑TTTTCACCAACGTTGGCAACCATCAACGCTGGCAAACGC‑ATTCCAGCGTTTTTTCGATTCa                                                                                                                            >  2:78555/1‑239 (MQ=255)                                                                                                                      gggAAATATAACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAAAACTTCTTGC‑AGCACCGCTA‑TCACCGGGCTTTTATTCGGGTTTGAGCGGAGA‑ACATCGGTCAT‑ATA‑TTTCCCCCCACGTTGGGCAACATAAAGATTGGTAGCCGC‑AATCCAGCggtatttagtatagaaa                                                                                                                         >  1:271894/1‑223 (MQ=255)                                                                                                                       ggAAATATAACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTTC‑AGC‑TCCCTACTCACCCGGCTGTTATCCGGGTTAGCGGGCCTA‑ACATCGGGCAC‑ATA‑TTTCCACCCACGTTGGCCAACATCAGTGCTGGCAACCGC‑CTTCCAGCGttctttcgagttcatc                                                                                                                         >  2:728933/1‑227 (MQ=255)                                                                                                                             ataACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGC‑AGC‑TCGCTACTCACCGGGCTGTTATCCGGGTTTGCGGGCATA‑ACATCGGTCAT‑TTA‑TTTCCACCAACGTTGGCAAACATCCATGCTGGGCACCCC‑CTTCCCACGTTCTTCACATTCTGTctctttt                                                                                                                   >  2:820152/1‑235 (MQ=255)                                                                                                                                 cTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGGTAAACACTTCCTGC‑AGC‑TCCCTACTTCCCGCGCTGTTATCCGGGTTAGAGGGCATA‑AAAACGGGCAT‑ATA‑TTTTCACCCACGTTTGGAAACATCAGTACTGGCCACCGCAATTCCA‑CGGTCTTCAGTATCAAAATCTACCAAggg                                                                                                              >  1:621168/1‑238 (MQ=255)                                                                                                                                        ggAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGC‑AGC‑TCGCTACTCACCGGGCTGTTATTCGGGGTAGCGGGCCTA‑AAATCGGTCCT‑ATA‑TTTCCACCAACGTTTGCAAACATCAGTGCTGGCAAACGC‑ATTCCAGCCTTCTTCAGATTCAATCTCTACCATTGGAAAcaac                                                                                                       >  2:717785/1‑237 (MQ=255)                                                                                                                                               tAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGC‑AGC‑TCGCTACTCCCCCGGCTGTTATCCGGGTTAGCGGGCATA‑ACATCGGTCAT‑ATA‑TTTCCACCAACGTTTGCAAACATCCGTGCTGGCAAccc                                                                                                                                                    >  2:171069/1‑187 (MQ=255)                                                                                                                                                    tGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGC‑AGC‑TCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATA‑ACATCGGTCAT‑ATA‑TTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGC‑ATTCCAGCGTTCTTCAGATTCAATCTCTACCATGGCAAACAG‑CAGATTACGCGCt                                                                                          <  2:923418/239‑1 (MQ=255)                                                                                                                                                      tttAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGC‑AGC‑TCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATA‑ACAt                                                                                                                                                                                                      <  1:955625/133‑1 (MQ=255)                                                                                                                                                      tttAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGC‑AGC‑TCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATA‑ACAt                                                                                                                                                                                                      >  2:955625/1‑133 (MQ=255)                                                                                                                                                      tttAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACATCTTTC‑AGA‑CCCCTACTCACCAGGCTGTTTTGCGGGTTAGCGGGCTTA‑ACATCGGTGAT‑ATC‑TTTCCACCACCGTTGGCGAACATCAGTGCTGGCAACCGC‑ATTCCAACGGTCTTCAAATTTAAACTCAACCCCGGCaaaaacacaatatcccgcttt                                                                                         >  1:441147/1‑221 (MQ=255)                                                                                                                                                       ttAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGC‑AGC‑TCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATA‑ACATCGGTCAT‑ATA‑TTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGC‑ATTCCAGCGTTCTTCAGATTCAATCTCTACCATGGCAAACAG‑CAGATTACGCGCTTTg                                                                                       <  2:291074/239‑1 (MQ=255)                                                                                                                                                       ttAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGC‑AGC‑TCGCTACTCACCCGGCTGTTATCCGGGGTAGCGGGCATA‑ACATCGGTCAT‑ATA‑TTTCCACCAACGTTTGCAAACATCAGTGCTGGCAAACGC‑ATTCCAGCGTTCTTCCGATTCCATTTCTACCATGGCAAACAG‑CAGATTACACGCTTTg                                                                                       >  2:307659/1‑239 (MQ=255)                                                                                                                                                       ttAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGGAAAAAAACTCTTTCAACC‑TCCTCATTCCCCGGCATGTT‑TCCCGGGTAGCGGGGCTA‑ACATAGGGCAT‑ATA‑TTTTCAACAACCTTGGGAACCATCCGGGCTGGAAAACGC‑ATTCCAGCGTTCTTTAGATTAAAATCCTACCAAGGCCAAAAGACACATAACGCGtttg                                                                                        >  1:923418/1‑238 (MQ=255)                                                                                                                                                         aaTTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGC‑AGC‑TCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATA‑ACATCGGTCAT‑ATA‑TTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGC‑ATTCCAGCGTTCTTCAGATTCAATCTCTACCATGGCAAACAG‑CAGATTACGCGCTTTGTc                                                                                     <  2:441147/239‑1 (MQ=255)                                                                                                                                                               gcAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACCCTTCTTGC‑AGC‑TCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATA‑ACATCGGTCAT‑ATA‑TTTCCACCAACGTTGGAAAACATCAGTGCTGGCAACCGC‑ATTCCAGCGTTCTTCAGATTCAATCTCTACCATGGCAAACAG‑CAGATTACGCGCTTTGTCGagata                                                                               <  2:45994/238‑1 (MQ=255)                                                                                                                                                                                    tgtgCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGC‑AGC‑TCGCTACTCACCGGGGTGTTATCCGGGTTAGCGGGCATA‑ACATCGGTCAT‑ATA‑TTTCCACCAACGTTGGCAAACATCAGTGCTGGAAACCGC‑ATTCCAGCGTTCTTCAGATTCAATTTCTACCCTGGCAAACAC‑CAGATTACGCGCTTTGTCGGGTTAACTGGGGTACTTATTCGCAcc                                                          >  2:1231344/1‑239 (MQ=255)                                                                                                                                                                                    tgtgCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGGAAAACACTTCTTGC‑AGC‑TCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATA‑ACATCGGGCAT‑ATA‑TTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGC‑ATTCCAGCGGTCTTCAGATTCAATCTCTACCCTGGCAAACAG‑CAGATTACGCGCTTTGTCGAGATAGATGGCGTAGTTATGCGCAcc                                                          >  2:245102/1‑239 (MQ=255)                                                                                                                                                                                    tgtgCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGGAAAACACTTCTTGC‑AGC‑TCGCTACTCACCGGGCTGTTATCCGGGGTAACGGGCATA‑ACATCGGTCCT‑ATA‑TTTCCACCCACGTTGGCAAACATCAGTGCTGGCAACCGC‑ATTCCAGCGTTCTTCAGATTCAATCTCTACCATGGCAAACAG‑CCGATTACGCGCTTTGTTGGGGTAGATTGCGTGGTTATTCGCAcc                                                          >  2:591624/1‑239 (MQ=255)                                                                                                                                                                                                  ttttGGCTCTTTGGGCAATTACGGGAGGTAAACCTCTTTTTGG‑AGC‑CCGCTACGCGCCGGGCTGTTATCCGGGTTAGCCGGCATA‑ACATCGGTCAT‑ATA‑TTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGC‑ATTCCAGCGTTCTTCAGATTCAATCTCTACCATGGCAAACAG‑CAGATTACGCGCTTTGTCGAGATAGATGGCGTAGTTATGCGCACCGTGA‑GATTTCAGCa                                           <  1:617448/236‑1 (MQ=255)                                                                                                                                                                                                  ctctTGGACCTTGGGGAATGACGGCAGGTAAAACCCTTCCTGC‑AGC‑TCCCTTCTCACCGGGCTGTTATCCGGGTTAGCGGGCTTA‑ACATCGGTCAT‑ATA‑TGTCCACCAACGTTGGCAACCATCAGTGCTGGCAACCGC‑ATTCCAGCGTTCTTCAGATTCAATCTCTACCATGGCAAACAG‑CAGATTACGCGCTTTGTCGAGATAGATGGCGTAGTTATGCGCACCGTGA‑GATTTCAGCa                                           <  1:854905/239‑1 (MQ=255)                                                                                                                                                                                                  ctctGGCACATTGGGCAATTACGGCAGGTAAAACCCTTCTTGC‑AGC‑TCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATA‑ACATCGGTCAT‑ATA‑TTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGC‑ATTCCAGCGTTCTTCAGATTCAATCTCTACCATGGCAAACAG‑CAGATTACGCGCTTTGTCGAGATAGATGGCGTAGTTATGCGCACCGTGA‑GATTTCAGCa                                           <  2:271894/239‑1 (MQ=255)                                                                                                                                                                                                    ctgtgcccttGGGCTATTCCGGCCGGTAAAAACCTTCTTTC‑AAC‑TCCATACCCTCCGGGCTGTTATCCGGGTTAGCGGGCATA‑ACATCGGGCAT‑ATA‑TTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGC‑ATTCCAGCGTTCTTCAGATTCAATCTCTACCATGGCAAACAG‑CAGATTACGCGCTTTGTCGAGATAGATGGCGTAGTTATGCGCACCGTGA‑GATTTCAGCACt                                         <  1:734921/231‑1 (MQ=255)                                                                                                                                                                                                             tGGGCAATTACGGCAGGGAAAACACTTCTTGC‑AGC‑CCGCTACCCACCCGGCTGTTATCCGGGGTAGCGGGCATA‑ACATTGGTCAT‑ATA‑TTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGC‑ATTCCAGCGTTCTTCAGATTCAATCTCCACCATGGCAAACAG‑CAGATTTAGAGCTTTGTCGAAGTTATTGGCGGAGTTATCCGCACCGCGG‑GATTTCAGCCCTGCTTCCaag                                >  2:313091/1‑237 (MQ=255)                                                                                                                                                                                                                    gataagagacagaaaCACTTCTTGC‑AAC‑TCGCTACTCACCGGGCTTTTATCCGGGTTAGCGGGCATA‑ACATCGGTCAT‑ATA‑TTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGC‑ATTCCAGCGTTCTTCAGATTCAATCTCTACCATGGCAAACAG‑CAGATTACGCGCTTTGTCGAGATAGATGGCGTAGTTATGCGCACCGTGA‑GATTTCAGCACTGCTTCCAGTTCTGGc                          <  2:1061561/226‑1 (MQ=255)                                                                                                                                                                                                                       cgggaagTAAAACCCTTCTTGC‑AGC‑TCGCTTCGCACCGGGCTGTTATCCGGGTTAGCGGGCATA‑ACATCGGTCAT‑ATA‑TTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGC‑ATTCCAGCGTTCTTCAGATTCAATCTCTACCATGGCAAACAG‑CAGATTACGCGCTTTGTCGAGATAGATGGCGTAGTTATGCGCACCGTGA‑GATTTCAGCACTGCTTCCAGTTCTGGCCAGa                      <  2:284205/233‑1 (MQ=255)                                                                                                                                                                                                                       cGGCAGGTAAAACACTTCTTGC‑AGC‑TCGGTACTCCCCGGGCTGTTAACCGGGTTAGCGGGCATA‑ACATCGGTCAT‑ATA‑TTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGC‑ATTCCAGCGTTCTTCAGATTCAATCTCTACCATGGCAAACAG‑CAGATTACGCGCTTTGTCGAGATAGATGGCGTAGTTATGCGCACCGTGA‑GATTTCAGCACTGCTTCCAGTTCTGGCCAGa                      <  2:593531/239‑1 (MQ=255)                                                                                                                                                                                                                       cGGCAGGGAAAACTCTTTCTGC‑AGC‑TCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATA‑ACATCGGTCAT‑ATA‑TTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGC‑ATTCCAGCGTTCTTCAGATTCAATCTCTACCATGGCAAACAG‑CAGATTACGCGCTTTGTCGAGATAGATGGCGTAGTTATGCGCACCGTGA‑GATTTCAGCACTGCTTCCAGTTCTGGCCAGa                      <  2:1006009/239‑1 (MQ=255)                                                                                                                                                                                                                              tAAAACCCTTCTTGC‑AGC‑TCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATA‑ACATCGGTCAT‑ATA‑TTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGC‑ATTCCAGCGTTCTTCAGATTCAATCTCTACCATGGCAAACAG‑CAGATTACGCGCTTTGTCGAGATAGATGGCGTAGTTATGCGCACCGTGA‑GATTTCAGCACTGCTTCCAGTTCTGGCCAGATGGGa                 <  2:1231224/237‑1 (MQ=255)                                                                                                                                                                                                                                     cttcttGC‑AGC‑TCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATA‑ACATCGGTCAT‑ATA‑TTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGC‑ATTCCAGCGTTCTTCAGATTCAATCTCTACCATGGCAAACAG‑CAGATTACGCGCTTTGTCGAGATAGATGGCGTAGTTATGCGCACCGTGA‑GATTTCAGCACTGCTTCCAGTTCTGGCCAGATGGGATTATGCCGa        <  2:976261/239‑1 (MQ=255)                                                                                                                                                                                                                                       tcttGC‑AGC‑TCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATA‑ACATCGGTCAT‑ATA‑TTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGC‑ATTCCAGCGTTCTTCAGATTCAATCTCTACCATGGCAAACAG‑CAGATTACGCGCTTTGTCGAGATAGATGGCGTAGTTATGCGCACCGTGA‑GATTTCAGCACTGCTTCCAGTTCTGGCCAGATGGGATTATGCCGACg      <  2:770158/239‑1 (MQ=255)                                                                                                                                                                                                                                       tcttGC‑AGC‑TCGCTACTCACCGGGCTGTTATCCGGGATAGCGGGCATA‑ACATCGGTCAT‑ATA‑TTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACACC‑ATTCCAGCGTTCTTCAGATTCAAACTCTACCATGGCAAAAAG‑CAGCTTACGCGCTTTGACGACATAGATGGCCTCAGTACGCGC‑CCCTGACGATTTCCGCACCCTCTCCAGTTCTGACCAGACGGGAATATGcgaag       >  1:605278/1‑234 (MQ=255)                                                                                                                                                                                                                                        cttGC‑AGC‑TCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATA‑ACATCGGTCAT‑ATA‑TTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGC‑ATTCCAGCGTTCTTCAGATTCAATCTCTACCATGGCAAACAG‑CAGATTACGCGCTTTGTCGAGATAGATGGCGTAGTTATGCGCACCGTGA‑GATTTCAGCACTGCTTCCAGTTCTGGCCAGATGGGATTATGCCGAc       >  2:1174527/1‑237 (MQ=255)                                                                                                                                                                                                                                        cttGC‑AGC‑TCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATA‑ACATCGGTCAT‑ATA‑‑TTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGC‑ATTCCAGCCTTCTTCAGATTAAATCTCTACCATGGCAAACAG‑CAGATTACCCGCTTTGTCGAGATAGTTGGCTTAGTTCTGCGCACCGTGT‑GATTTCAGACCTGCTTCCCGTTCTGCCACGATggactttgtccacacat    >  1:10717/1‑224 (MQ=255)                                                                                                                                                                                                                                           gC‑AGC‑TCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATA‑ACATCGGTCAT‑ATA‑TTTCCACCAACGTTGGCAAACCTCAGTGCTGGCAACCGC‑ATTCCAGCGTTCTTCAGATTCAATCTCTACCATGGCAAACAG‑CAGCTTACGCCCTTTTTCGAGATAGATGGCTTCGTACGTCGCCACGTCA‑GATTTCAGCACTGCTTCCATTACTGCCCAGATGGGaattgttctcttaata  >  1:1061561/1‑224 (MQ=255)                                                                                                                                                                                                                                           gC‑AGC‑TCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATA‑ACATCGGTCAT‑ATA‑TTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGC‑ATTCCAGCGTTCTTCAGATTCAATCTCTACCATGGCAAACAG‑CAGATTACGCGCTTTGTCGAGATAGATGGCGTAGTTATGCGCACCGTGA‑GATTTCAGCACTGCTTCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGa  <  2:420225/239‑1 (MQ=255)                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                       TGAGACAGAGTGAAAGGTCAGATTTGCGGAGTAATGCACATAATGGTTATTTAAATAAACCACATGAATCATTAATGGGTTATTCATTATTTTTGTGATTTACTCACAATGTATCA‑GGGAAATATAACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGC‑AGC‑TCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATA‑ACATCGGTCAT‑ATA‑TTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGC‑ATTCCAGCGTTCTTCAGATTCAATCTCTACCATGGCAAACAG‑CAGATTACGCGCTTTGTCGAGATAGATGGCGTAGTTATGCGCACCGTGA‑GATTTCAGCACTGCTTCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGA  >  NC_000913/4092913‑4093384

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