Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,796,648	A→C	19.2%	V46G (GTG→GGG)	ygaU ←	uncharacterized protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,796,648	0	A	C	19.2%	126.7 / 2.1	61	V46G (GTG→GGG)	ygaU	uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (21/26);  new base C (7/6);  total (28/33)
Fisher's exact test for biased strand distribution p-value = 7.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGACTTAACGGTATAAAACTGGCTGGCAGTGGCTGGTGTCGCCGTTTTCACCTGATCATCGACACTGGCAATACCGGAAATATTCCCCACCGCAACAAGGATTTTCTCCTTCGCCTCCT‑GACTCAGGCCGTCACCAGTGACCGTCGCTTTGCCGTCGGCAATTTGAATATTCACTTTATCGGCATCCGGTATACCGGTTTTGTTCAGATGCTCCTGCACCTTCTTCGCCTGATCGTCTTTATCGTGCTGACCTGTAACCGCGTCCCAGAGTTTTTCTCCGGCATCTTTCACAAAATTGAACAGACCCATAGCTACCTCATTCGTTAACGGAAACACTCAGAAAGTGTAGCAGAGGAACGCGGCGCGTTAACGCCGCGGGGATGAGCAC  >  NC_000913/2796428‑2796863                                                                                                                                                                                                                              |                                                                                                                                                                                                                        aTTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGACTTAACGGTATAAAACTGGCTGGCAGTGGCTGGTGTCGCCGTTTTCACCTGATCATCGACACTGGCAATACCGGAAATATTCCCCACCGCAACCACGATTTTCTCCTTCGCCCCCT‑GAATCAGGCCGGCACCAGTGAACGTCGCTTTTCCGTAGGGCATTTTGATATTACCTTTAta                                                                                                                                                                                                                  >  1:164987/1‑227 (MQ=255)  ttAGCGTTACCGTAGACCTGGTTGGAAATGGCCCTCAGAGTGTCGCCAGCCTTAACGGTATAAAACTGGCTGGCAGTGGCTGGTGTCGCCGTTTTCACCTTATCATCGACACTGGCAATACCGGAAATATTCCCCACCGCAACAAGGATTTTCTCCTTCGCCTCCT‑GACTCAGGCCGTCACCAGTGACCGTCGCTTTGCCGTCGGCAATTTGAATATTCACTTTATCg                                                                                                                                                                                                                 <  1:312093/228‑1 (MQ=255)      cGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGACTTAACGGTATAAAACTGGCTGGCAGTGGCTGGTGTCGCCGTTTTCACCTGATCATCGAAACTGGCAATACCGGAAATATTCCCCACCGCAACAAGGATTTTCTCCTTCGCCTCCT‑ACCTCAACCCGCCACCAGAGCCCGTCGCTTTCCCCTCGCCAATTTCAAAATTCCAATTATCGGCAt                                                                                                                                                                                                             >  1:76437/1‑228 (MQ=255)       gTTACCGTCGACCTGTTTGGAAATGGCACCCAGGGTGTCGCCCGACTTAACGGTATAAAACTGGCTGGCAGTGGCTGGGGTCGCCTTTTTCCCCTGATCATCGACACTGGCAATACCGGAAATATTCCCCACCGCAACAAGGATTTTCTCCTTCGCCTCCT‑GACTCAGGCCGTCACCAGTGACCGTCGCTTTGCCGTCGGCAATTTGAATATTCACTTTATCGGCATc                                                                                                                                                                                                            <  1:647696/228‑1 (MQ=255)        ttACCGTAGACCTGTTTGGAAATGGCCCTCAGAGTGTCGCCAGACTTAACGGTATAAAACTGGCTGGCAGTGGCTGGTGTCGCCGTTTTCACCTGATCATCGACACTGGCAATACCGGAAATATTCCCCACCGCAACAAGGATTTTCTCCTTCGCCTCCT‑GACTCAGGCCGTCACCAGTGACCGTCGCTTTGCCGTCGGCAATTTGAATATTCACTTTATCGGCATc                                                                                                                                                                                                            <  2:76437/227‑1 (MQ=255)         tACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGACTTAACGGTATAAAACTGGCTGGCAGTGGCTGGTGTCGCCGTTTTCACCTGATCATCAACACTGGCAATACCGGAAATATTCCCCACCGCAACAAGGATTTTCTCCTTCGCCTCCT‑GACTCCGGCCGTCACCCGTGACCGTCGCTTTGCCGTCGGCAATTTGAATATTACCTTAATCGCCATCCg                                                                                                                                                                                                          >  1:478818/1‑228 (MQ=255)           ccGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGACTTAACGGTATAAAACTGGCTGGCAGTGGCTGGTGTCGCCGTTTTCACCTGATCATCGCCACTGGCAATACCGGGAATATTCCCCACCGCCACAAGGCTTTTCTCCTTCGCCTCCT‑GACTCAGGCCGGCACCCGTGACCGTCGCTTTTCCCTCGGGCATTTGAAAATTACCTTTATCGGccacc                                                                                                                                                                                                           >  1:261710/1‑221 (MQ=255)           ccGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGACTTAACGGTATAAAACTGGCTGGCAGTGGCTGGTGTCGCCGTTTTCACCTGATCATCGACACTGGCAATACCGGAAATATTCCCCCCCCCAACAAGGATTTTCTCCTTCGCCCCCT‑GACTCAGGCCGTCACCAGAGACCGTCCCTTTGCCGCCGGCAATTTGAAAATTCACTTTTACGGCAACCGGt                                                                                                                                                                                                        >  1:37030/1‑228 (MQ=255)           ccGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGACTTAACGGTATAAAACTGGCTGGCAGTGGCTGGTGTCGCCGTTTTCACCTGATCATCGACACTGGCAATACCGGAAATATTCCCCACCGCAACAAGGATTTTCTCCTTACCCCCCT‑GACTCAGGCCGTCACCAGTGCCCGTCGCTTTTCCGCCCGCAATTTGAATATTCCCTTTTACGGCCTCCGGt                                                                                                                                                                                                        >  1:75442/1‑228 (MQ=255)               agccttGTTTGGAAATGGCCGTCAGGGTGTCGCCAGACTTAACGGTCTAAAACTGGCTGGCATTGGCTGGTGTCGCCGTTTTCACCTGATCATCGACACTGGCAATCCCGGAAATATTCCCCACCGCAACAAGGATTTTCTCCTTCGCCTCCT‑GACTCAGGCCGTCACCAGTGACCGTCGCTTTGCCGTCGGCAATTTGAATATTCCCTTTATCGGCATCCGGTATAc                                                                                                                                                                                                    <  1:488679/223‑1 (MQ=255)                gACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGACTTAACGGTATAAAACTGGCTGGCAGTGGCTGGTGTCGCCGTTTTCACCTGATAATCGACACTGGCAATACCGGAAATATTCCCCACCGCAACAAGGATTTTCTCCTTCGCCTCCT‑GACTCAGGCCGTCACCAGTGACCGTCGCTTTGCCGTCGGCAATTTGAATATTCACTTTATCGGCATCCGGttta                                                                                                                                                                                                     >  2:488679/1‑226 (MQ=255)                    tttttGGAAAGGGCACCCGGTGGTGCGCCCGACTTAAAGGTATAAAACTGGCTGGCAGTGGCTGGTGTCGCCGTTTTCACCTGATCATCGACCCTGGCAATACCGGAAATATTCCCCACCGCAACAAGGATTTTCTCCTTCGCCTCCT‑GACTCAGGCCGTCACCAGTGACCGTCGCTTTGCCGTCGGCAATTTGAATATTCACTTTATCGGCATCCGGTATACCGGtt                                                                                                                                                                                               <  1:481315/226‑1 (MQ=255)                        tGGAAATGGCACTCAGAGTGTCGCCAGACTTTACGGTATAAAACTGGCTGGCAGTGGCTGGTGTCGCCGTTTTCACCTGATCATCGACACTGGCAATACCGGAAATATTCCCCACCGCAACAAGGATTTTCTCCTTCGCCTCCT‑GACTCAGGCCGCCACCAGTGACCGTCGCTTTGCCGTCGGCCATTTGAATATTCACTTTATCGGCATCCGGCATACCGGctta                                                                                                                                                                                             >  2:234832/1‑225 (MQ=255)                             aTGGCCCTCAGGGTGTCCCCAAACTTAACGGTATAAAACTGGCTGGCAGTGGCTGGGGTCGCCGTTTTCACCTGATCAGCGACACTGGCAATCCCGGAAATATTCCCCACCGCAACAAGGATTTTCTCCTTCGCCTCCT‑GACTCAGGCCGTCACCAGTGACCGTCGCTTTGCCGTCGGCAATTTGAATATTCACTTTATCGGCATCCGGTATACCGGTTTTGTTCAGa                                                                                                                                                                                      <  1:329261/228‑1 (MQ=255)                              tGGCACCCAGGGTGTCGCCAGACTTAACGGTATAAAACTGGCTGGCAGTGGCTGGGGTCGCCGTTTTCACCTGATCATCGACACTGGCAATACCGGAAATATTCCCCACCGCAACAAGGATTTTCTCCTTCGCCTCCT‑GACTCAGGCCGTCACCAGTGACCGTCGCTTTGCCGTCGGCAATTTGAATATTCACTTTATCGGCATCCGGTATACCGGTTTTGTTCAGa                                                                                                                                                                                      <  2:75442/227‑1 (MQ=255)                                                  gACTTAACGGTATAAAACTGGCTGGCAGTGGCTGGTGTCGCCGTTTTCACCTGATCATCGACACTGGCAATACCGGAAATATTCCCCACCGCAACAAGGATTTTCTCCTTCGCCTCCT‑GACTCAGGCCGTCACCAGTGACCGTCGCTTTGCCGTCGGCAATTTGAATATTCCCTTCATCGGCCTCCGGTATACCGGTTTTGTTCCGATGCTCCTGCAACTTCTTCGcc                                                                                                                                                                 >  2:214130/1‑228 (MQ=255)                                                      tAACGGTATAAACGTGGGTTGCAGGGGCTGGGGTCCCCGTTTTCCCCTGATCAACGACACTGGCAATACCGGAAATATTCCCCACCGCAACAAGGATTTTCTCCTTCGCCTCCT‑GACTCAGGCCGTCACCAGTGACCGTCGCTTTGCCGTCGGCAATTTGAATATTCACTTTATCGGCATCCGGTATACCGGTTTTGTTCAGATGCTCCTGCACCTTCTTCGCCTGAt                                                                                                                                                             <  1:280002/228‑1 (MQ=255)                                                       aaCGGTATAAAACTGGCTGGCAGTGGCTGGTGTCGCCGTTTTCACCTGATCATCGACACTGGCAATACCGGAAATATTCCCCACCGCAACAAGGATTTTCTCCTTCGCCTCCT‑GACTCAGGCCGTCACCAGTGACCGTCGCTTTGCCGTCGGCAATTTGAATATTCACTTTATCGGc                                                                                                                                                                                                               <  2:490366/177‑1 (MQ=255)                                                         cGGTATAAAACTGGCTGGCAGTGGCTGGTGTCGCCGTTTTCACCTGATCATCGACACTGGCAATACCGGAAATATTCCCCACCGCAACAAGGATTTTCTCCTTCGCCTCCT‑GACTCAGGCCGTCACCAGTGACCGTCGCTTTGCCGTCGGCAATTTGAATATTCACTTTATCGGCATCCGGTATACCGGTTTTGTTCAGATGCTCCTGCACCTTCTTCGCCTGATCGt                                                                                                                                                          <  2:344261/228‑1 (MQ=255)                                                              ttaacctggtggggcggtgttGGTGTCCCCGTTTTCCCCGGACTCTCGACCTTGGCAATACCGGAAAGCTCCCCCACCGCAACAAGGATTTTCTCCTTCGCCTCCT‑GACTCAGGCCGCCACCAGTGACCGTCGTTTTGCCGTCGGCAATTTGAATATTCACTTTATCGGCATCCGGTATACCGGTTTTGTTCAGATGCTCCTGCACCTTCTTCGCCTGATCGTCTTTa                                                                                                                                                     <  1:20402/208‑1 (MQ=255)                                                              taAAACTGGCTGGCAGTGGCTGGTGTCGCCGTTTTCACCTGATCATCGACACTGGCAATACCGGAAATATTCCCCACCGCAACAAGGATTTTCTCCTTCGCCTCCT‑GACTCAGGCCGTCACCAGTGACCCTCGCTTTGCCGTCGGCAATTTGAATATTCACTTTATCGGCCTCCGGTTTACCGGTTTTTGTCAGATGCTCCTGCCCCTTCTTCGCCTGATCGCCtt                                                                                                                                                       >  2:213068/1‑226 (MQ=255)                                                              taAAACTGGCTGGCAGTGGCTGGTGCTGCCGTTTTCACCTAATCGTCGACACCGGCAATACCGGAAATATCCCCCACCGCAACAAGGATTTTCTCCTTCGCCTCCT‑TACTCAGGCCGTCACCAGTGACCGTCGCTTTGCCGTCGGCAATTTGAATATTCACTTTATCGGCATCCGGTATACCGGTTTTGTTCAGATGCTCCTGCACCTTCTTCGCCTGATCGTCTTTa                                                                                                                                                     <  1:476495/228‑1 (MQ=255)                                                                     ggctggcAGTGGCTGGTGTCGCCGTTTTCACCTGATCATCGACACTGGCAATACCGGAAATATTCCCCACCGCAACAAGGATTTTCTCCTTCGCCTCCT‑GACTCAGGCCGTCACCAGTGACCGTCGCTTTGCCGTCGGCAATTTGAATATTCACTTTATCGGCATCCGGTATACCGGTTTTGTTCAGATGCTCCTGCACCTTCTTCGCCTGATCGACTTTATCGTGc                                                                                                                                               >  2:97273/1‑227 (MQ=255)                                                                      gctggcAGTGGCTGGTTTCGCCGTTTTCACCTGATCATCGACACTGGCAATACCGGCAATATTCCCCACCGCAACACGGATTTTCTCCTTCGCCTCCT‑GACTCCGGCCGTCACCAGTGACCGTCGCTTTGCCGTCGGCAATTTGAATATTCCCTTACTCGGCATCCGGCTTACCGGTTTTGTTCAGATGCTCCTGCACCTTCTTCGCCTGCTCGTCTTTATCGTGCt                                                                                                                                              >  2:272839/1‑227 (MQ=255)                                                                      gctggcAGTGGCTGGTGTCGCCGTTTTCACCTGATCATCGACACTGGCAATACCGGAAATATTCCCCACCGCAACAAGGATTTTCTCCTTCGCCTCCT‑GACTCAGGCCGTCACCAGTGACCGTCGCTTTGCCGTCGGCAATTTGAATATTCACTTTATCGGCATCCGGTATACCGGTTTTGTTCAGATGCTCCTTCACCTTCTTTGCCGGTTCGCCTTTATCGTGCt                                                                                                                                              >  2:7727/1‑227 (MQ=255)                                                                       gtggcAGGGGTTGGTGGCCCCTTTTTCACCTGATCATCGACACGGGAAATCCCGGAAATTTTCCCCCCCTCAACAAGGATTTTCTCCTTCGCCTCCT‑GACTCAGGCCGTCACCAGTGACCGTCGCTTTGCCGTCGGCAATTTGAATATTCACTTTATCGGCATCCGGTATACCGGTTTTGTTCAGATGCTCCTGCACCTTCTTCGCCTGATCGTCTTTATCGTGCTg                                                                                                                                             <  1:214130/226‑1 (MQ=255)                                                                           cAGTGGCTGGTGTCGCCGTTTTCACCTGATCATCGACACTGGCAATACCGGAAATATTCCCCACCGCAACAAGGATTTTCTCCTTCGCCTCCT‑GACTCAGGCCGTCACCAGTGACCGTCGCTTTGCCGTCGGCAATTTGAATATTCACTTTATCGGCATCCGGTATACCGGTTTTGTTCAGATGTCCCTGCACCTTCTTCGCCTGCTCGTCTTTTTCGTGCttgcct                                                                                                                                         >  2:408063/1‑222 (MQ=255)                                                                                      tttgCCTTTTTCCCCTGCTCACCGACACTGTCAATACCGGAAATATTCCCCCCCCCAACAAGGATTTTCTCCTTCGCCTCCT‑GAATTAGGCCGTCACCAGGGACCGTCGCTTTGCCGTCGGCAATTTGAATATTCACTTTATCGGCATCCGGTATACCGGTTTTGTTCAGATGCTCCTGCACCTTCTTCGCCTGATCGTCTTTATCGTGCTGACCTGTAACCGCGTcc                                                                                                                             <  1:343760/225‑1 (MQ=255)                                                                                        ccccGTTTTCACCTGATCCTCGCACCTGGCAATACCGGAAATATTCCCCACCGCAACAAGGATTTTCTCCTTCGCCTCCT‑GACTCAGGCCGTCACCAGTGACCGTCGCTTTGCCGTCGGCAATTTGAATATTCACTTTATCGGCATCCGGTATACCGGTTTTGTTCAGATGCTCCTGCACCTTCTTCGCCTGATCGTCTTTATCGTGCTGACCTGTAACCGCGTCCCa                                                                                                                           <  2:240782/226‑1 (MQ=255)                                                                                          ccGTTTTCACCTGATCATCGACACTGGCAATACCGGAAATATTCCCCACCGCAACAAGGATTTTCTCCTTCGCCTCCT‑GACTCAGGCCGTCACCAGTGACCGTCGCTTTGCCGTCGGCAATTTGAATATTCACTTTATCGGCATCCGGTATACCGGTTTTGTTCAGATGCTCCTGCACCTTCTTCGCCTGATCGTCTTTATCGTGCTGACCTGTAACCGCGTCCCaga                                                                                                                         >  2:360491/1‑228 (MQ=255)                                                                                                            cGACACTGGCAATACCGGAAATATTCCCCACCGCAACAAGGATTTTCTCCTTCGCCTCCT‑GACTCAGGCCGTCACCAGTGACCGTCGCTTTGCCGTCGGCAATTTGAATATTCACTTTATCGGCATCCGGTATACCGGTTTTGTTCAGATGCTCCTGCACCTTCTTCGCCTGATCGTCTTTATCGTGCTGACCTGTACCCGCGTCCCAGAGTTTTTCTCCGGCATCtt                                                                                                       >  2:379067/1‑228 (MQ=255)                                                                                                                  tGGCAATACCGGAAATATTCCCCACCGCAACAAGGATTTTCTCCTTCGCCTCCT‑GACTCAGGCCGTCACCAGTGACCGTCGCTTTTCCGTCGGCAATTTGAATATTCACTTTATCGGCATCCGGTATACCGGTTTTGGTCAGATGCTCCTGCACCTTCTTCGCCTGATAGTCTTTATCTTGCTGACGTGTATCCGCGGCCCAGGGGTTTTTTCCGGCATCttttccc                                                                                                  >  1:643018/1‑223 (MQ=255)                                                                                                                  tGGCAATACCGGAAATATTCCCCACCGCAACAAGGATTTTCTCCTTCGCCTCCT‑GACTCAGGCCGTCACCAGTGACCGTCGCTTTGCCGTCGGCAATTTGAATATTCACTTTATCGGCATCCGGTATACCGGTTTTGTTCAGATGCTCCTGCACCTTCTTCGCCTGATCGTCTTTATCTGGCTGACCTGTAACCACAGCCCGAAGTTTTTCTCCGGAATCTTTcaca                                                                                                  >  1:197790/1‑227 (MQ=255)                                                                                                                    gCAATACCGGAAATATTCCCCACCGCAACAAGGATTTTCTCCTTCGCCTCCT‑GACTCAGGCCGTCACCAGTGACCGTCGCTTTGCCGTCGGCAATTTGAATATTCACTTTATCGGCATCCGGCATACCGGTTTTGTTCAGATGCTCCTGCACCTTCTTCGCCTGATCGTCTTTATCGTACTGACCTGTAAGCGACGCCCAGGGTCTTTATCCTGCAGCATTcacgcaa                                                                                               >  2:279101/1‑224 (MQ=255)                                                                                                                         cccGGAAATTTCCCCCACCGAAACAAGGATTTTCTCCTTCGCCTCCT‑GGCTCAGGCCGTCACCAGTGACCGTCGCTTTGCCGTTGGCAATTTGAATATTCACTTTATCGGCATCCGGTATACCGGTTTTGTTCAGATGCTCCTGCACCTTCTTCGCCTGATCGTCTTTATCGTGCTGACCTGTAACCGCGTCCCAGAGTTTTTCTCCGGCATCTTTCACAAAATTGaa                                                                                          <  1:213068/227‑1 (MQ=255)                                                                                                                               aaTTTTCCCCCCGGCAAAAGGGGTTTTATTCTTCTCCTCCTGGCCTACGGCAGTC‑CCAGTGGCCGTAGGTTTGCCGGCGGGAATTTGAATATTCACTTTATCGGCATCCGGTAGACCGGTTTTGTTCAGATGCTCCGGCACCTTCTTCGCCTGATCGTCTTTATCGTGCTGACCTGTATCCGCGTCCCAGAGTTTTTCTCCGGCATCTTTCACAAAATTGAACAGAcc                                                                                    <  1:234832/228‑1 (MQ=255)                                                                                                                                 tatTCCCCGCCGGCACAAGGATTTTCTCCTTCGCCTTCT‑GACTCAGGCCGTCACCAGTGACCGTCGCTTTGCCGTCGGCAATTTGAATATTCACTTTATCGGCATCCGGTATACCGGTTTTGTTCAGATGCTCCTGCACCTTCTTCGCCTGATCGTCTTTATCGTGCTGACCTGTAGCCGCGTCCCAGAGTTTTTCTCCGGCATCTTTCACAAAATTGAACAGACCCa                                                                                  <  2:197790/228‑1 (MQ=255)                                                                                                                                 tatCCCCCCCCGCAACAAGGATTTTCCCCTTCCCCTCCT‑GACTCAGGCCGTCCCCAGTGACCGTCGCTTTGCCGGCGGCAATTTGAATATTCCCTTTATCGGCATCCGGTATACCGGTTTTGTTCAGATGCTCCTGCACCTTCTTCGCCTGATCGTCTTTATCGTGCTGACCTGTAACCGCGTCCCAGAGTTTTTCTCCGGCATCTTTCACAAAATTGAACAGACCCa                                                                                  <  1:7727/228‑1 (MQ=255)                                                                                                                                   ttCCCCACCGCAACAAGGATTTTCTCCTTCGCCTCCT‑GACTCAGGCCGTCACCAGTGACCGTCGCTTTGCCGTCGGCAATTTGAATATTCACTTTATCGGCATCCGGTATACCGGTTTTGTTCAGATGCTCCTGCACCTTCTTCGCCTGATCGTCTTTATCTTGCTGACCTGTAACCGCGTCCCAGAGTTTTTTTCCGGCACATTTCCACAAATTTAACAAAccacta                                                                                >  1:208910/1‑224 (MQ=255)                                                                                                                                       ccACCGCAACAAGGATTTTCTCCTTCGCCTCCT‑GACTCAGGCCGTCACCAGTGACCGTCGCTTTGCCGTCGGCAATTTGAATATTCACTTTATCGGCATCCGGTATACCGGTTTTGTTCAGATGCTCCTGCACCTTCTTCGCCTGATCGTCTTTATCGTGCTGACCTGTAACCGCGTCCCAGAGTTTTTTTCCGGCATCTTTCCCAAAATTAAACAGAccccaaacaa                                                                            >  1:18277/1‑221 (MQ=255)                                                                                                                                           cGCAACAAGGATTTTCTCCTTCGCCTCCT‑GACTCAGGCCGTCACCAGTGACCGTCGCTTTGCCGTCGGCAATTTGAATATTCACTTTATCGGCATCCGGTATACCGGTTTTGTTCAGATGCTCCTGCACCTTCTTCGCCTGATCGTCTTTATCGTGCTGACCTGTAACCGCGTCCCCGAGTTTTTCTCCGGCCTCTTTCCCAAAATTTAACAAACCCAAACCCACCTc                                                                        >  1:607587/1‑228 (MQ=255)                                                                                                                                           cGCAACAAGGATTTTCTCCTTCGCCTCCT‑GACTCAGGCCGTCACCAGTGACCGTCGCTTTGCCGTCGGCAATTTGAATATTCACTTTATCGGCATCCGGTATACCGGTTTTGTTCAGATGCTCCTGCACCTTCTTCGCCTGATCGTCTTTATCGTGCTGACCTGTAACCGCGTCCCAGAGTTTTTCTCCGGCATCTTTCACAAAATTGAACAAACCCAAAGCAACCTc                                                                        >  1:602994/1‑228 (MQ=255)                                                                                                                                             caacaaGGATTTTCTCCTTCGCCTCCT‑GACTCAGGCCGTCACCAGTGACCGTCGCTTTGCCGTCGGCAATTTGAATATTCACTTTATCGGCATCCGGTATACCGGTTTTGTTCAGATGCTCCTGCACCTTCTTCGCCTGATCGTCTTTATCGTGCTGACCTGTAACCGCGTCCCAGAGTTTTTCTCCGGCATCTTTCACAAAATTGAACAGACCCATAGCTACCTCAt                                                                      >  2:569574/1‑228 (MQ=255)                                                                                                                                                        ttCTCCTTCGCCTCCT‑GACTCAGGCCGTCACCAGTGACCGTCGCTTTGCCGTCGGCAATTTGAATATTCACTTTATCGGCATCCGGTATACCGGTTTTGTTCAg                                                                                                                                                                                       >  1:48460/1‑104 (MQ=255)                                                                                                                                                         tctcCTTCGCCTCCT‑GACGCAGGCCGTCACCAGTGACCGTCGCTTTGCCGTCGGCAATTTGAATATTCACTTTATCGGCATCCGGTATACCGGTTTTGTTCAGa                                                                                                                                                                                      <  2:48460/104‑1 (MQ=255)                                                                                                                                                          ctcCTTCCCCTCCT‑GACTCAGGCCGTCACCAGTGCCCGGCGCTTTGCCGTCGGCAATTTGAATATTCACTTTATCGGCATCCGGTATACCGGTTTTGTTCAGATGCTCCTGCACCTTCTTCGCCTGATCGTCTTTATCGTGCTGACCTGTAACCGCGTCCCAGAGTTTTTCTCCGGCATCTTTCACAAAATTGAACAGACCCATAGCTACCTCATTCGTTAACGGaaa                                                         <  1:360491/228‑1 (MQ=255)                                                                                                                                                          cgccTTCGCCTCCT‑GACCCATGCCGTCACCAGTGACCGTCGCTTTGCCGTCGGCAATTTGAATATTCACTTTATCGGCATTCGGTATACCGGTTTTGTTCAGATGCTCCTGCGCCTTCTTCGCCTGATCGTCTTTATCGTGCTGACCTGTAACCGCGTCCCAGAGTTTTTCTCCGGCATCTTTCACAAAATTGAACAGACCCATAGCTTCCTCATTCGTTAACGGaaa                                                         <  2:164987/226‑1 (MQ=255)                                                                                                                                                                         gACTCAGGCCGTCACCAGTGACCGTCGCTTTGCCGTCGGCAATTTGAATATTCACTTTATCGGCATCCGGTATACCGGTTTTGTTCAGATGCTCCTGCACCTTCTTCGCCTGATCGTCTTTATCGTGCTGACCTGTAACCGCGTCCCAGAGTTTTTCTCCGGCATCTTTCACAAAATTGAACAGACCCATAGCTACCTCATTCGTTAACGGAAACACTCAGAAAgtgt                                           <  2:37030/228‑1 (MQ=255)                                                                                                                                                                           cTCAGGCCGTCACCAGTGACCGTCGCTTGGCCGTCGGCAATTTGAATATGCACTTTAGCGGCATCCGGTATACCGGTTTTGTTCAGATGCTCCTGCACCTTCTTCGCCTGATCGTCTTTATCGTGCTGACCTGTAACCGCGTCCCAGAGTTTTTCTCCGGCATCTTTCACAAAATTGAACAGACCCATAGCTACCTCATTCGTTAACGGAAACACTCAGAAAGTGTAg                                         <  2:627868/228‑1 (MQ=255)                                                                                                                                                                             cAGGCCGTCTCCAGTGACCGTCGCTTTGCCGTCGGCAATTTGAATATTCACTTTATCGGCATCCGGTATACCGGTTTTGTTCAGATGCTCCTGCACCTTCTTCGCCTGATCGTCTTTATCGTGCTGACCTGTAACCGCGTCCCAGAGTTTTTCTCCGGCATCTTTCACAAAATTGAACAGACCCATAGCTACCTCATTCGTTAACGGAAACACTCAGAAAGTGTAGCa                                       <  2:208910/228‑1 (MQ=255)                                                                                                                                                                                gCCGTCACCAGTGACCGTCGCTTTGCCGTCGGCAATTTGAATATTCACTTTATCGGCATCCGGTATACCGGTTTTGTTCAGATGCTCCTGCACCTTCTTCGCCTGATCGTCTTTATCGTGCTGACCTGTAACCGCGTCCCAGAGTTTTTCTCCGGCATCTTTCACAAAATTGAACAGACCCATAGCTACCTCATTCGTTAACGGAAACACTCAGAAAGTGTAGCag                                      >  2:190442/1‑226 (MQ=255)                                                                                                                                                                                   gTCACCAGTGACCGTCGCTTTGCCGTCGGCAATTTGAATATTCACTTTATCGGCATCCGGTATACCGGTTTTGTTCAGATGCTCCTGCACCTTCTTCGCCTGATCGTCTTTATCGTGCTGACCTGTAACCGCGTCCCAGAGTTTTTCTCCGGCATCTTTCTCAAAATTGAACAGACCCATAGCTACCTCCTTCGTTCACGGAAAACCTCAGAAAGGGACGCAGAGGaa                                 >  2:110441/1‑228 (MQ=255)                                                                                                                                                                                        cagttgccttccgttttCCGTCGGCAATTTGAATATTCCCTTTATCGGCATCCGGTATACCGGTTTTGTTCAGATGCTCCTGCACCTTCTTCGCCTGATCGTCTTTATCGTGCTGACCTGTAACCGCGTCCCAGAGTTTTTCTCCGGCATCTTTCACAAAATTGAACAGACCCATAGCTACCTCATTCGTTAACGGAAACACTCAGAAAGTGTAGCAGAGGAACgcgg                            <  1:97273/215‑1 (MQ=255)                                                                                                                                                                                          gTGACCGTCGCTTTGCCGTCGGCAATTTGAATATTCACTTTATCGGCATCCGGTATACCGGTTTTGGTCAGATGCTCCTGCACCTTCTTCGCCTGATCGTCTTTATCGTGCTGACCTTTAACCGCGTCCCAAAGTTTTTCTCCGGCATCTTTCACAAAATTGCACAGACCCATAGCTACCTCATTCGTTAACGGAAACACCCAGAAAGAGTAGCAGAGGAACccgggc                          >  1:158652/1‑226 (MQ=255)                                                                                                                                                                                                gttgCTTTGCCGTCGACAATTTGAATATTCGCTTTATCGGCATCCGGTATACCGGTTTTGTTCAGATGCTCCTGCCCCTTCTTCGCCTGATCGTCTTTATCGTGCTGACCTGTAACCGCGTCCCAGAGTTTTTCTCCGGCATCTTTCACAAAATTGAACAGACCCATAGCTACCTCATTCGTTAACGGAAACACTCAGAAAGTGTAGCAGAGGAACGCGGCGCGTTaa                    <  2:478818/225‑1 (MQ=255)                                                                                                                                                                                                          ctgcGGCATTTTGAATTTTCCTTTTATCGGTCTCCGGTATACCGTTTTTGTTCAGCTGCTCCTGCACCTTCTTCGCCTGATCGTCTTTATCGTGCTGACCTGTAACCGCGTCCCAGAGTTTTTCTCCGGCATCTTTCACAAAATTGAACAGACCCATAGCTACCTCATTCGTTAACGGAAACACTCAGAAAGTGTAGCAGAGGAACGCGGCGCGTTAACGCCGCgggg          <  1:525504/225‑1 (MQ=255)                                                                                                                                                                                                                cAATTTGAATATTCACTTTATCGGCATCCGTTATACCGGTTTTGTTCAGATGCTCCTGCACCTTCTTCGCCTGATCGTCTTTATCGTGCTGACCTGTAACCGCGTCCCAGAGTTTTTCTCCGGCATCTTTCACAAAATTGAACAGACCCATAGCTACCTCATTCGTTAACGGAAACACTCAGAAAGTGTAGCAGAGGAACGCGGCGCGTTAACGCCGCGGGGATGAGc    <  2:496999/228‑1 (MQ=255)                                                                                                                                                                                                                 atttttAATATTCCTTTTCTCGGCATCCGTGTTACCGTTTTTGTTCAGATGCTCCTGCACCTTCTTCGCCTGACCGGCTTTATCGTGCTGACCTGTAACCGCGTCCCAGAGTTTTTCTCCGGCATCTTTCACAAAATTGAACAGACCCATAGCTACCTCATTCGTTAACGGAAACACTCAGAAAGTGTAGCAGAGGAACGCGGCGCGTTAACGCCGCGGGGATGAGca   <  1:408063/226‑1 (MQ=255)                                                                                                                                                                                                                 aaTTTGAATATTCCCTTTAACCGCCTTCGGTAATCCGGTTTTGTCCATGTGCTCCTGCACCTTCGTCGCCTGTTCGTCCTTATCGTGCTGACCTGTAACCGCGTCCCAGAGTTTTTCTCCGGCATCTTTCACAAAATTGAACAGACCCATAGCTACCTCATTCGTTAACGGAAACACTCAGAAAGTGTAGCAGAGGAACGCGGCGCGTTAACGCCGCGGGGATGAGca   <  2:582707/228‑1 (MQ=255)                                                                                                                                                                                                                 aaTTTGAATATTCACTTTATCGGAATCCGGTATACCGGTTTTGTTCAGATGCTCCTGCACCTTCTTCGCCTGATCGTCTTTATCGTGCTGACCTGTAACCGCGTCCCAGAGTTTTTCTCCGGCATCTTTCACAAAATTGAACAGACCCATAGCTACCTCATTCGTTAACGGAAACACTCAGAAAGTGTAGCAGAGGAACGCGGCGCGTTAACGCCGCGGGGATGAGca   <  2:261710/228‑1 (MQ=255)                                                                                                                                                                                                                  ttttgatttaccACTTTATCGGCTACCGGATTACCGGTTTTAGTCTTATGCTCCTGCCCTGTCTTCTCCTGTTTGTCTTTATCGTTCTGACCTGTAATCCCGGCCCAGAGTTTTTCGCCGGCATCTTTCACAAAATTGAACAGACCCATAGCTACCTCATTCGTTAACGGAAACACTCAGAAAGTGTAGCAGAGGATCTCGGCGCGTTAACGCCGCGGGGATGAGcac  <  1:279101/217‑1 (MQ=255)                                                                                                                                                                                                                              |                                                                                                                                                                                                                        ATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGACTTAACGGTATAAAACTGGCTGGCAGTGGCTGGTGTCGCCGTTTTCACCTGATCATCGACACTGGCAATACCGGAAATATTCCCCACCGCAACAAGGATTTTCTCCTTCGCCTCCT‑GACTCAGGCCGTCACCAGTGACCGTCGCTTTGCCGTCGGCAATTTGAATATTCACTTTATCGGCATCCGGTATACCGGTTTTGTTCAGATGCTCCTGCACCTTCTTCGCCTGATCGTCTTTATCGTGCTGACCTGTAACCGCGTCCCAGAGTTTTTCTCCGGCATCTTTCACAAAATTGAACAGACCCATAGCTACCTCATTCGTTAACGGAAACACTCAGAAAGTGTAGCAGAGGAACGCGGCGCGTTAACGCCGCGGGGATGAGCAC  >  NC_000913/2796428‑2796863

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