Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,728,950	T→C	15.1%	E141G (GAG→GGG)	xylB ←	xylulokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,728,950	0	T	C	15.1%	189.3 / 3.4	78	E141G (GAG→GGG)	xylB	xylulokinase
Reads supporting (aligned to +/- strand):  ref base T (45/15);  new base C (2/13);  total (48/30)
Fisher's exact test for biased strand distribution p-value = 1.86e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAACAAAGCACCAGTAATTTCGCTGCCTT‑CGTATAATGCGGGCATCTGGTCACGAGATAAGTCGCAAGCCTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGTGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGTAATAATA‑CTTTGTCGATTTGACGGAATATCTCCGGCTCA‑TGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCTTGCGCACAGCGCCCGTCGTTCCACAAAATGGCAGGGCGTAACACCCGTTGCTGAGCATCCAGCAAGGTTGCTCCGTGCATCTGGCCGGCAATACCCAATGCTTTAACGTCCT  >  NC_000913/3728733‑3729175                                                                                                                                                                                                                            |                                                                                                                                                                                                                                   aaaaaaaGCCCCAGGAAGTGCGCTGCCTTGTGAATAGGGCGGAC‑TCTGGTCACGAGGAAGGTCGCAAGCCCGAAGTCGGACGTCACTCCAGTCACGCTTTGCAACATCCAGCCAAATGGTGCCAGCTGCGTCAGACATAGCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGTAATAATA‑CTTTGTCGATTTGACGGAATATCTCCGGCTCA‑Tg                                                                                                                                                                                                                         <  2:504821/224‑1 (MQ=255) aaaCAAAGCACCACTCATTCCCCTGCATA‑CTTATAATGCGGGAAACTGGTCACGAGATGAGCCGCAACCCTGCAGCTTCCGGTCACTCCAATCACTCTTTGCGACATCCAGCCACCTGGTGCCAGTTGCGTCAGACATATCGCTGGCATACTCCCCCGTCATACGCAGACGCAAGTAGTCTTTCGGTAATAATC‑CTTTGTCGATTTGACGGAATATCTCCGGCTCA‑Tg                                                                                                                                                                                                                         <  2:204755/227‑1 (MQ=255)               tAATTTCGCTGCCTT‑CGTATAATGCGGGCATCTGGTCACGAGATAAGTCGCAAGCCTGCAGCATGACGTCCCTCTAGTCACGCTTTGCGACATCAAGCCACATGGTGCCAGCTGCGTCCGACATCTCGCCGGCAAACTCCCCCGTCATACGCAGACGCAAGAAATCTTTCGGTAATAATA‑CTTTGCCGCTTTGCCGACACACGCCCGGCTCA‑GGCCACAGAACCCATa                                                                                                                                                                                                           >  2:297605/1‑228 (MQ=255)                         gCCTT‑CGTATAATGCGGGCATCTGGTCACGAGATAAGTCGCAAGCCTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCACCCACATGGTGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCACACGCCAGCAATCTTTCGGTAATAATA‑CTTTGTCGATTTGACGGAATATCTCCGGCTCC‑TGCCCCTAAACCCATAGCAATTTAga                                                                                                                                                                                                 >  2:137259/1‑227 (MQ=255)                         gCCTT‑CGTATAATGCGGGAATCTGGTCACGAGATAAGTCGCAAGCCTGCAGCATGTCGTAACTCCAGTCACGCATTGCGACATCCAGACACATGGTGCCAGCTGCGTCAGACATACAGCTGGCAAAATCCCCATTCATACGCAGACGCAAGTAATCCTTCGGTACTACTA‑CTTTGTGTGTTTGACGTAATAGCTCCGACTCA‑GTCCGCTTGACCCTTAGCAACTTAg                                                                                                                                                                                                  >  1:319176/1‑227 (MQ=255)                         gCCTT‑AGTATAATGCGGGCATCTGGTCACGAAATAAGTCGCAAGCCTGCAGCATGCCGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGTGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGTAATAATA‑CTTTGTCGATTTGACGGAATATCTCCGGCTCA‑TGCCGCTGAACCCATAGCAATTTAgg                                                                                                                                                                                                 <  1:75147/228‑1 (MQ=255)                              ttggataACGCGGGCATCTGGTACCGAGATAAGTCGCAAGCCTGCAGCATGACGTCCCCCCAGTCACGCTTTGCGACATCCAGCCACAGGGTGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCGTTCGGTAATAATA‑CTTTGTCGATTTGACGGAATATCTCCGGCTCA‑TGCCGCTGAACCCATAGCAATTTAGGCGCa                                                                                                                                                                                             <  1:592557/224‑1 (MQ=255)                                      tGCGGGCATCTGGTCACGAGATAAGTCGCAAGCCTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGTGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGAAGACGCAAGGAATCTTTCGGTAATAATA‑CTTTGTCGATTGGACGGAATATCCCCGG‑TCACTGCCCCTGAACCCTTACCAATTACGGCGCAGAAAATcc                                                                                                                                                                                     >  2:75147/1‑228 (MQ=255)                                          ggCATCTGGTCACGAGATAAGTCGCAAGCCTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGCTACATGGTGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATCCGCAGACGCAAGTAATCTTTCGGTAATAATA‑CTTTGTCGATTTGACGGAACATCTCCGGCACA‑TGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCggg                                                                                                                                                                                  >  1:90089/1‑227 (MQ=255)                                          ggCATCTGGTCACGAGATAAGTCGCAAGCCTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGTGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGTAATAATA‑CTTTGTCGATTTGACGGAGTATCTCCGGCTCA‑TGCCGCTGAACCCATCACAATTTAGGCGCAGTAAATCCggg                                                                                                                                                                                  >  2:613599/1‑227 (MQ=255)                                          ggCATCTGGTCACGAGATAAGTCGCAAGCCTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGTGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGTAATAATA‑CTTTGTCGATTTGACGGAATATCTCCGGCTCA‑TGCCGCTGAACCCATAGCAATTTAGGCGCAGCAAATCCGGGc                                                                                                                                                                                 >  1:513278/1‑228 (MQ=255)                                          ggCATCTGGTCACGAGATAAGTCGCAAGCCTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGTGCCACCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATCCGCAGACGCAAGGAATCTTTCGGGAATAATA‑CTTTGTCGATTTGACGGAAAAACACCGACTCA‑TGCCCCTGAACCCTCAGCACTTTAGGGCCAGTAAACCCcgg                                                                                                                                                                                  >  2:62226/1‑227 (MQ=255)                                           gCATCTGGTCACGAGATAAGTCGCAAGCCTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGTGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGTAATAATA‑CTTTGTCGCTTTGACGGAATATCTCCGGCTCA‑TGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGc                                                                                                                                                                                 >  1:304624/1‑227 (MQ=255)                                            cATCTGGTAACGAGATAAGTCGCAAGCCTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGTGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGCAATCTTTCGGTAATAATA‑CTTTGTCGATTTGACGGAATATCTCCGGCTCA‑TGCGGCTGAACCCATAGCAATTTAGGCGCAGTACATCCGGGca                                                                                                                                                                                >  1:203898/1‑227 (MQ=255)                                                    cccGAGATAAGCCGCAAGCCTGCAGCTGTACGTCCTCCCAGTCACCTTTTGCGACTACCAGCCCCATGGTGCCAGCTGCGTTCGACATAGCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGTAATAATA‑CTTTGTCGATTTGACGTAATATCTCCGGCTCA‑TGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGt                                                                                                                                                                       <  2:501604/226‑1 (MQ=255)                                                        agaTAAGTCGCAAGCCTGAAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGTGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGTAATAATA‑CTTTGTCGATTTGACGGAATATCTCCGGCTCA‑TGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGc                                                                                                                                                                    <  1:652792/227‑1 (MQ=255)                                                                        tGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGTGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGTACTAATA‑CTTTGTCGATTTGACGGAATATCTCCGGCTCA‑TGCCGCTGAACCCATAGCAATTTAGGGGCAGTAAATCCTGGCACCATCCGGACGCCGGTAAACACCCGCGaa                                                                                                                                                   >  2:115180/1‑227 (MQ=255)                                                                        tGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGTGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGTAATAATA‑CTTTGTCGATTTGACGGAATATCTCCGGCCCA‑TGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGa                                                                                                                                                    >  1:124658/1‑227 (MQ=255)                                                                         gcagcaTGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGTGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGTAATAATA‑CTTTGTCGATTTGACGGAATATCTCCGGCTCA‑TGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGAt                                                                                                                                                   >  1:240657/1‑227 (MQ=255)                                                                          cagcaTGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGTGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGGAATAATA‑CTTTGTCGATTTGACGGAATATCTCCGGCTCA‑TGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCcgggctt                                                                                                                                                  >  1:231076/1‑222 (MQ=255)                                                                                       cTCCAGTCACGCTTTGCGACATCCAGCCACATGGTGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGTAATAATA‑CTTTGTCGATTTGACGGAATATCTCCGGCTCA‑TGCCGCTGAACCCTTACCAATTTAGCCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCt                                                                                                                                    >  2:8710/1‑228 (MQ=255)                                                                                        tCCAGTCACGCTTTGCGACATCCAGCCACATGGTGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGTAATAATA‑CTTTGTCGATTTGCCGGAATATCTCCGGCTCA‑TGCCGCTGAACCCATAGCAATTTAgg                                                                                                                                                                                                 >  2:173076/1‑166 (MQ=255)                                                                                        tCCAGTCACGCTTTGCGACATCCAGCCACATGGTGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGTAATAATA‑CTTTGTCGATTTGACGGAATATCTCCGGCTCA‑TGCCGCTGAACCCATAGCAATTTAgg                                                                                                                                                                                                 <  1:173076/166‑1 (MQ=255)                                                                                                  cTTTGCGACATCCAGCCACATGGTGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGTAATAATA‑CTTTGTCGATTTGACGGAATATCTCCGGCTCA‑TGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTCGGGAACTCGCGCTTCCAGCAAAg                                                                                                                          >  2:462190/1‑227 (MQ=255)                                                                                                    ttGCGACATCCAGCCACATGGTGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGTAATAATA‑CTTTGTCGATTTGACGGAATATCTCCGGCTCA‑TGCCGCTGAACCCATAGCAATTTAGGGGCCGTAAATCCGGGCATCCTCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTTCCAGCaaaagga                                                                                                                       >  2:45911/1‑224 (MQ=255)                                                                                                     tGCCACACCCCGCCACATGGTGCCAGGTGCGTCAGACATTTCGCTGGAAAACTCCCCCGTCATTCGCAGACGCAAGTAATCTTCCGGTAATAATA‑CTTTGTCGATTTGACGGAATATCCCCGGCTCA‑TGCCGCTGATCCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCa                                                                                                                      <  2:240657/228‑1 (MQ=255)                                                                                                       cGACATCCAGCCACATGGTGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGTAATAATA‑CTTTGTCGATTTGACGGAATATCTCCGGCTCA‑TGCCGCTGAACCCATAGCAATTTAGGCGCATTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGGGCAcc                                                                                                                    >  2:549594/1‑227 (MQ=255)                                                                                                               cgCCACATGGTGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATGCGCACACTCAAGTAATCTTTCGGTAATAATA‑CTTTGTCGATTTGACGGAATATCTCCGGCTCA‑TGCCGCTGAACCCAAAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCTTGCGca                                                                                                            <  1:401020/227‑1 (MQ=255)                                                                                                               aGCCACCTGGTTCCAGCTGCGTCAGACCTATTGCTGGCAAACTCCCCCGTCATCAGCAGACGAAAGTAATCTTTCGGTAATAATA‑CTTTGGCGATTTGACGGAATATCTCCGGCTCA‑TGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCTTGCGca                                                                                                            <  1:337387/228‑1 (MQ=255)                                                                                                                 ccACATGGTGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGTAATAATA‑CTTTGCCGATTTGACGGAATATCTCCGGCTCA‑TGTCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATACTCCGGTGTCCGGTAATAACCCGGGATTGCGGAACTCGCACTTCCAGCAAAAGGCACCCTTgtgc                                                                                                             >  2:404453/1‑225 (MQ=255)                                                                                                                         tGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGTAATAATA‑CTTTGTCGATTTGACGGAATATCTCCGGCTCA‑TGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCTTCATCAGATTGCCGGTAATCACCCGCGATTGCAGAACTCGCACTTCCAGCAAAAGGCACTCTTGCCCCCAGCCCCcgc                                                                                                  >  2:679765/1‑227 (MQ=255)                                                                                                                           ccAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGTAATAATA‑CTTTGTCGATTTGACGGAATATCTCCGGCTCA‑TGCCGCTGAACCCATAGCTATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTTCCAGCAAAATGCACTCTTGCTCCCAGCGCCcgtcg                                                                                                >  1:19374/1‑228 (MQ=255)                                                                                                                              gCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGTAATAATA‑CTTTGTCGATTTGACGGAATATCTCCGGCTCA‑TGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCTGGCATCATCAGGTTGCCGGAAATAACCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGGGGCCTCTGCCGCCCAGCGCCCGTCGTTc                                                                                             >  2:13841/1‑228 (MQ=255)                                                                                                                                  cGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGGCGCAAGTAATCTTTCGGTAATAATA‑CTTTGTCGATTTGACGGAATATCTCCGGCTCA‑TGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCTTGCGCACAGCGCCCGTCGTTCcaca                                                                                         <  1:137259/228‑1 (MQ=255)                                                                                                                                   gTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGTAATAATA‑CTTTGTCGATTTGACGGAATATCTCCGGCTCA‑TGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCTTCCGCACAGCGCCCGTCGTTCCACaa                                                                                        >  1:111247/1‑228 (MQ=255)                                                                                                                                         cATACCGCTGGCAAATCACCCCGTCATACGGAGACGCAAGTAATCTTTCGGTAATAATA‑CTTTGTCGATTTGACGGAATATCTCCGGCTCA‑TGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCTTGCGCACAGCGCCCGTCGTTCCACAAAATGGc                                                                                  <  1:9361/228‑1 (MQ=255)                                                                                                                                          ttttctctgcaaaaCTCCCCCGTCATCCGCAGACGCAATTAAGCTTTCGGTAATAATA‑CTTTGGCGATTTGACGGAATATCCCCGGCTCA‑TGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCTTGCGCACAGCGCCCGTCGTTCCACAAAATGGCa                                                                                 <  2:111247/217‑1 (MQ=255)                                                                                                                                          atatCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGTAATAATC‑CTTTGTCGATTTGACGGAATATCTCCGGCTCA‑TGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCTTGCGCACAGCGCCCGTCGTTCCACAAAATGGCa                                                                                 <  1:449129/228‑1 (MQ=255)                                                                                                                                           tatggagtgCAAACCCCCCCGTTAAACGCCAACGCAAGTAATCTTTCGGTAATAATA‑CTTTGCCGATTTGGCGGAATATCGCCCGCTCA‑TGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCTTGCGCACAGCGCCCGTCGTTCCACAAAATGGCAg                                                                                <  2:513278/220‑1 (MQ=255)                                                                                                                                            atCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGTAATAATA‑CTTTGTCGATTTGACGGAATATCTCCGGCTCA‑TGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCTTGCGCACACCGCCCGTCGTTCCCCAAAATGGCAgg                                                                               >  1:420884/1‑228 (MQ=255)                                                                                                                                            atCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGTAATAATA‑CTTTGTCGATTTGACGGAATATCTCCGGCTCA‑TGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCCTCCCGCAAAGTGCACTCTTGCGCACAGCCCCCGTCGTTCCCCAAAATGCCAgg                                                                               >  1:327457/1‑228 (MQ=255)                                                                                                                                                  ggCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGTAATAATA‑CTTTGTCGATTTGACGGAATATCTCCGGCTCA‑TGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCTTGCGCACAGCGCCCGTCGTTCCACAAAATGGCAGGGCg                                                                            >  1:173574/1‑225 (MQ=255)                                                                                                                                                      taCTCCCCCGTCATACGCAGACGCAAGGAAGCGTCCGGTAATAATA‑CGTTGGCGATGTGACGGAATATCTCCGGCTCA‑TGCCGCTGAACCCATAGCAATTTAGGCGCAGGAAATCCGGGCATCATCAGGTTGCCGGTAATCCCCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCTTGCGCACAGCGCCCGTCGTTCCACAAAATGGCAGGGCGTAACAcc                                                                     <  1:418317/227‑1 (MQ=255)                                                                                                                                                                tCATACGCAGACGCAAGTAATCTTTCGGTAATAATA‑CTTTGTCGATTTGACGGAATATCTCCGGCTCA‑TGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTTCCACCAAAGTGACCTCTTCCGCCCAGCGCCCGCCGTTCCACAAAATGGCAGGCCGAAACAcctctcccttaa                                                           >  2:451470/1‑218 (MQ=255)                                                                                                                                                                 cATACGCAGACGCAAGTAATCTTTCGGTAATAATA‑CTTTGTCGATTTGACGGAATATCTCCGGCTCA‑TGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCTTGCGCACAGCGCCCGTCGTTCCACAAAATGGCa                                                                                 <  1:451470/205‑1 (MQ=255)                                                                                                                                                                       cAGACGGAAGTCATCTTTCCGTTA‑AATACCTTTGTCGATTTGACGGAATATCCCCGGCTCA‑TGCCCCTGGACCCATAGCAATTTAGGCGCAGTAAATCCTGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCTTGCGCACAGCGCCCGTCGTTCCACAAAATGGCAGGGCGTAACACCCGTTGCTGAGCATCCAg                                                    <  2:203898/228‑1 (MQ=255)                                                                                                                                                                       aaagcGCAAGTAATCTGTCGGTAATAATA‑TTTTGTCGATTTGACGGAATATCCCCGGCTCA‑TGCCGCTGAACCCAGAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCCCCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCTTGCGCACAGCGCCCGTCGTTCCACAAAATGGCAGGGCGTAACACCCGTTGCTGAGCATCCAg                                                    <  2:173574/224‑1 (MQ=255)                                                                                                                                                                         gaagCAAGTTATTTTTTGGTAATAATT‑CTTTTTCAATTTTACGGAATTTCCCCGGCTCA‑TTCCGCTGGACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCTTGCGCACAGCGCCCGTCGTTCCACAAAATGGCAGGGCGTAACACCCGTTGCTGAGCATCCAGCa                                                  <  2:420884/225‑1 (MQ=255)                                                                                                                                                                          aCGCAAGTAATCTTTCGGTAATAATA‑CTTTGTCGATTTGACGGAATATCTCCGGCTCA‑TGCCGCTGAACCCATAGCACTTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCTTGCGCACAGCGCCCGCCGTTCCACAAAATGGCAGGGCGTAACACCCGTTGCTGAGAATCCAGcca                                                 >  1:92515/1‑226 (MQ=255)                                                                                                                                                                          aCGCAAGTAATCTTTCGGTAATAATA‑CTTTGTCGATTTGACGGAATATCTCCGGCTCA‑TGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCTTGCGCACAGCGCCCGTCGTTCCCCAAAATGGCAGGGCGTACACCCCGTTGCTGGCCCTCCACCaa                                                 >  2:123368/1‑228 (MQ=255)                                                                                                                                                                          aCGCAAGTAATCTTTCGGTAATAATA‑CTTTGTCGATTTGACGGAATATCTCCGGCTCA‑TGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCTTCCGCACAGCGCCCGACGTTCCACAAAATGGACGGGCGTAACACCGCTTTCTGacataccaccaa                                                 >  1:330447/1‑217 (MQ=255)                                                                                                                                                                                        tCGGTAATAATA‑CTTTGTCGATTTGACGGAATATCTCCGGCTCA‑TGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCTTGCGCACAGCGCCCGTCGTTCCCCAAAATGGCAGGGCGTAACACCCGTTGCTGCGCATCCCGCAAGGTTGCTCCGTGCa                                   >  1:342521/1‑228 (MQ=255)                                                                                                                                                                                        tCGGTAATAATA‑CTTTGTCGATTTGACGGAATATCTCCGGCTCA‑TGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCTTGCGCACAGCGCCCGTCGTTCCACAAAATGGCAGGGCGTAACACCCGTTGCTGAGCATCCAGCAAGGTTGCTcctggca                                   >  1:514009/1‑223 (MQ=255)                                                                                                                                                                                        tCGGTAATAATA‑CTTTGTCGATTTGACGGAATATCTCCGGCTCA‑TGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCTTGCGCACAGCGCCCGTCGTTCCACAAAATGGCAGGGCGTAAAACCCGTTGCTTAGCATCCAGCAAGGTTGCTCCGTGCa                                   >  1:621837/1‑228 (MQ=255)                                                                                                                                                                                           cttattaaatcTTTGTCGATTTTACGGAATTACGCCCGGTCA‑AGCCGCTGGACCCATAGCCATTTAGGCGCAGAAAATCCGGGCATCGTCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCCTGCCCACAGTGCCCGTCGTTCCACAAAATGGCAGGGCGTAACCCCCGTTGCTGAGCATCCAGCAAGGTTGCTCCGTGCAt                                  <  2:305983/217‑1 (MQ=255)                                                                                                                                                                                            taataaTA‑CTTTGTCGATTTGACGGAATATCTCCGGCTCA‑TGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCTTGCGCACAGCGCCCGTCGTTCCACAAAATGGCAGGGCGTACCACCCGTTTCTGAGCATCCAGCAAGGTTGCTCCGCGCa                                   >  1:305983/1‑224 (MQ=255)                                                                                                                                                                                              acaatt‑CTTTGGCGATTTGACGGAAATCCCCCGGCTAA‑GGCCGCGGAACCCATAGCAATTTAGGCGCAGTAAAACCCGGCATCATCAGGTTGCCGGTAATCCCCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCTTGCGCACAGCGCCCGTCGTTCCACAAAATGGCAGGGCGTAACCCCCGTTGCTGAGCATCCAGCAAGGTTGCTCCGTGCATCTGGc                             <  2:514009/226‑1 (MQ=255)                                                                                                                                                                                               ataaccatttGTCGATTTTGCGGAATTTTCCCGGGCAA‑CGCCGCTGAACCCCAAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCTTGCGCACAGCGCCCGTCGTTCCACAAAATGGCAGGGCGTAACACCCGTTGCTGAGCATCCAGCAAGGTTGCTCCGTGCATCTGGc                             <  2:265356/221‑1 (MQ=255)                                                                                                                                                                                                   a‑CTTTTTCGTTTTGGCGGAATATCCCCGGCTCA‑TGCCGTTGAACCCAGAGCAATTTGGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAACCACCCGCGATTGCGGAACCCGCGCTTCCAGCAAAGTGCACTCTTGCGCACAGCGCCCGTCGTTCCACAAAATGGCAGGGCGTAACCCCCGTTGCTGAGCATCCAGCAAGGTTGCTCCGTGCATCTGGCCGGCa                        <  2:327457/228‑1 (MQ=255)                                                                                                                                                                                                    ttttttTCGATTTGCAGGAATTTCCCCGGCACA‑TGCCGCTGAACCCAAAGCAATTGGGGCGCAGTAAATCCGGGCATCACCAGGTTGCCCGTAATCACCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCTTGCGCACAGCGCCCGTCGTTCCACAAAATGGCAGGGCGTAACGCCCGTTGCTGAGCATCCAGCAAGGTTGCTCCGTGCATCTGGCCGGCa                        <  2:342521/226‑1 (MQ=255)                                                                                                                                                                                                       ttGTCGATTTGACGGAATATCTCCGGCTCA‑TGCCGCTGAACTCATAGCAATTTAGGCTCAGTAAATCCAGGACACATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTACCAGCAAAATGCACTCTTGCGCCAACTGCCCGACGCTCCCGAAAATAGCCAGGCGACAGACCCGTTCCGTGGACTCCAGCAAAGATGCTCCAGCCACATGGCCGGCacaa                     >  2:671067/1‑225 (MQ=255)                                                                                                                                                                                                       ttGTCGATTTGACGGAATATCTCCGGCTCA‑TGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGCTTGCGGAACTCGCGCTTCCCGCAACGAGCACCCTTGCGCCCCGCGCCCGTCGTTCCGCCAAATGGCAGGGCGTACAACCCGTGGCTCAGCATCCCGCAACGTTGCTCCTGGGCTCTGGCCGCCatac                     >  2:132321/1‑225 (MQ=255)                                                                                                                                                                                                       ttGTCGATTTGACGGAATATCTCCGGCTCA‑TGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTTCCCGCAAAGTGCACTCTTGCGCACAGCGCCcgtc                                                                                                 >  2:268590/1‑152 (MQ=255)                                                                                                                                                                                                       ttGTCGATTTGACGGAATATCTCCGGCTCA‑TGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTTCCAGCCAAGTGCACTCTTCCGCACAACGCCCGTCGTTCACCAAAATGCCAGGGCTCACCACCCGTGCCCGAGCATCCAGACAGGCTGCTCCGCGAAACTGGCCGCCaaaa                     >  2:323558/1‑226 (MQ=255)                                                                                                                                                                                                       ttGTCGATTTGACGGAATATCTCCGGCTCA‑TGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCTTGCGCACAGCGCCCGTCGTTCCACAAAATGGCAGGGCGTAACACCCGTTGCTGAGCATCCAGCAAGGCTGTGCCGTGCATCTGCCCGGCAATa                     >  1:560574/1‑228 (MQ=255)                                                                                                                                                                                                       ttGTCGATTTGACGGAATATCTCCGGCTCA‑TGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCTTGCGCACAGCGCCCGTCGTTCCACAAAATGGCAGGGCGTAACACCCGTTGCTGAGCATCAAGCAAGGTTGCTCCGTGCACCTGACCGGCCATa                     >  1:352608/1‑228 (MQ=255)                                                                                                                                                                                                       ttGTCGATTTGACGGAATATCTCCGGCTCA‑TGCCGCTGAACCCATAGCAATTTAGGCGAAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGATTCCAGCAAAGTTCACTCTTGCGCCCAGCGCCCGCCGTTCCACAAAATGGACGGGCGAACCACCCGTTGCTGACCACCAGGCAAGGATCCTCCGAGCATGTGCCCAGCaaac                     >  2:199800/1‑226 (MQ=255)                                                                                                                                                                                                       gtgccGATTTGACGGAATACGTCTGGCTCA‑AGACGCTGAAACCTTAACAATTTAGGGGCAGTAAATCAGGGCATCATCAGGTTGCCGGTAATCAACCGCGATTGCGGAACTCGCGCTTCCAGCAAAGAGCACACTTGCGCACAGCGCACGTCGTTCCCCAAAATGGCAGTGCGTACCACCCGTTGCTCGGAACCCAGCACGACGGTCCCGGGCTTCTGGCCGGCaaaa                     >  2:643560/5‑226 (MQ=255)                                                                                                                                                                                                        tGTCGATTTGACGGAATATCTCCGGCTCA‑TGCCGCTGAACCCCTAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGCTTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCTTGCGCACAGCGCCCGTCGTTCCACCAACTGGCAGGGCGTAACACCCGTCGCTGAGCCTCCAGCACGCTGGCTCCCTGCCCCAGTCCGGCaacaa                    >  2:535020/1‑225 (MQ=255)                                                                                                                                                                                                        tGTCGATTTGACGGAATATCTCCGGCTCA‑TGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCTTGCGCACAGCGCCcgtcg                                                                                                <  1:268590/152‑1 (MQ=255)                                                                                                                                                                                                          tCGATTTGACGGAATATCTCCGGCTCA‑TGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCTTGCGCACAGCGCCCGTCGTTCCCCAAAATGGCAGGGCGTAACACCCGTTCCTGCGCATCCACCAAGGTAGCTCCGGGCATCTGCCCGCCAAAAcca                  >  2:110993/1‑227 (MQ=255)                                                                                                                                                                                                           ggttttaccGGGATACCCCCGGCACG‑TGCCGCTGCACCCACAGCAATTTGGGCGCAGTACATCCGGGACTCCTCAGGTTGCCGGTAATCCCCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCTTGCGAACAGCGCCCGTCGTTCCACCAAATGGCAGGGCGTAACACCCGTTGCTGAGCATCCAGCAAGGTTGCTCCGTGCATCTGGCCGGCAATACCCa                 <  2:525962/220‑1 (MQ=255)                                                                                                                                                                                                           cGATTGGACGGAAATTCCCCGGCTCC‑TGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCTTGCGCACAGCGCCCGTCGTTCCACAAAATGGCAGGGCTTAACACCCGTTGCTGAGCATCCAGCAAGGTTGCTCCGTGCATCTGGCCGGCAATACCCa                 <  1:671067/228‑1 (MQ=255)                                                                                                                                                                                                            gATTTGGCGGAATTTCCCCGGCCCA‑TTCCGCTGAACCCATACAAATTTAGGCGACGTAAATCCGGGCACCATCAGGTTGCGCGTAATCCCCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCTTGCGCACAGCGCCCGTCGTTCCACAAAATGGCAGGGCGTAACACCCGTTGCTGAGCATCCAGCAAGGTTGCTCCGTGCATCTGGCCGGCAATACCCaa                <  2:621837/228‑1 (MQ=255)                                                                                                                                                                                                            caTTTGACGGAATATCCCCGGCCAC‑TGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCCCCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCTTGCGCACAGCGCCCGTCGTTCCACAAAATGGCAGGGCATAACACCCGTTGCTGAGCATCCAGCAAGGTTGCTCCGTGCATCTGGCCGGCAATACCCaa                <  1:13841/227‑1 (MQ=255)                                                                                                                                                                                                             aTTTGACGGAATATCTCCGGCTCA‑TGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCTTGCGCACAGCGCCCGTCGTTCCACAACATGGCAGGGCGTACCACCCGTTGCTGACCATCCAGCAAGGTTGCTCCGTGCATCTGGCCGGCAATACCCaa                >  1:601350/1‑227 (MQ=255)                                                                                                                                                                                                               ttGACGGAATATCTCCGGCTCA‑TGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCTTGCGCACAGCGCCCGTCGTTCCACAAAATGGCAGGGCGTAACACCCGTTGCTGAGCATCCAGCAAGGATGCTCCGTGCATCTGGCCGGCAATACCaaa                >  1:467780/1‑225 (MQ=255)                                                                                                                                                                                                               ttGACGGAATATCTCCGGCTCA‑TGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCTTGCGCACAGCGCCCGTCGTTCCACAAAATGGCAGGGCGTAACACCCGTTGCTGAGCATCCACCAAGGTTGCTCCGTCCATCTGGCCGGCATTACCCAATg              >  1:225181/1‑227 (MQ=255)                                                                                                                                                                                                                         ttctcCGGCTCA‑TGCCGCTAAACCCATAGCAATTTAGGCGCAGGAAATCCGGTCATCATCAGGTTGCCGGTAATCCCCCGCGATTGCGGAACTCGCTCTTCCAGCAAAGTGCACTCTTGCGCACAGCGCCCGTCGTTCCACAAAATGGCAGGGCGTAACACCCGTTGCTGAGCATCCAGCAAGGTTGCTCCGTGCATCTGGCCGGCAATACCCAATGCTTTAACGTc    <  1:404453/226‑1 (MQ=255)                                                                                                                                                                                                                          cctcCGGCTCA‑TGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCTTGCGCACAGCGCCCGTCGTTCCACAAAATGGCAGGGCGTAACACCCGTTGCTGAGCATCCAGCAAGGTTGCTCCGTGCATCTGGCCGGCAATACCCAATGCTTTAACGTCCt  <  1:613599/227‑1 (MQ=255)                                                                                                                                                                                                                            |                                                                                                                                                                                                                                   TAACAAAGCACCAGTAATTTCGCTGCCTT‑CGTATAATGCGGGCATCTGGTCACGAGATAAGTCGCAAGCCTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGTGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGTAATAATA‑CTTTGTCGATTTGACGGAATATCTCCGGCTCA‑TGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCTTGCGCACAGCGCCCGTCGTTCCACAAAATGGCAGGGCGTAACACCCGTTGCTGAGCATCCAGCAAGGTTGCTCCGTGCATCTGGCCGGCAATACCCAATGCTTTAACGTCCT  >  NC_000913/3728733‑3729175

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