Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,909,227	A→C	15.9%	G146G (GGT→GGG)	pyrG ←	CTP synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,909,227	0	A	C	15.9%	120.6 / 3.5	59	G146G (GGT→GGG)	pyrG	CTP synthetase
Reads supporting (aligned to +/- strand):  ref base A (31/16);  new base C (2/9);  total (33/26)
Fisher's exact test for biased strand distribution p-value = 6.08e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCA‑GGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAA‑CGGCAAGGATTCGATATCACCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCAC‑GCGCTCTTTGATTGCGTTAGTGATGTGCGGAATAACCTGCACGGTTGCGCCGAGGTAGTCACCGCGGCGTTCTTTACGCAGAACGTCAG‑AGTAGATACGACCCGTGGTGAAGTTGTTGCGGCGGCTCATTTTGGTACGAATGAAACGCTCGTAGTGCCCCAGGTCCAGGTCGGTTTCA  >  NC_000913/2908990‑2909464                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                tGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCA‑GGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCACCAGCCATCTGCCGAATGGCTTCCAGGAA‑CGGCAAGGATTCGATATCa                                                                                                                                                                                                                                                 >  1:487661/1‑238 (MQ=255) tGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCA‑GGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAA‑CGGCAAGGATTCGATATCa                                                                                                                                                                                                                                                 >  1:214619/1‑238 (MQ=255) tGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCA‑GGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAA‑CGGCAAGGATTCGATATCa                                                                                                                                                                                                                                                 >  1:286656/1‑238 (MQ=255) tGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCA‑GGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAA‑CGGCAAGGATTCGAAATCa                                                                                                                                                                                                                                                 >  1:321661/1‑238 (MQ=255) tGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCA‑GGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGGGCATAAACAGAGTGTGCTACCGGCCAATTTCAAAAGCCATCTGGCGAATCCCTTAGGGGAA‑CGGCAAGGATTCCATAacac                                                                                                                                                                                                                                                >  2:340669/1‑239 (MQ=255) tGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCA‑GGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGGCAGGTGCATAAACAGAGTGTGCTCACGGCCACTTTCAACAGCCATCTGGCGAATCGCTTCGAGGGA‑CGGCAAGGATTCGCTATCa                                                                                                                                                                                                                                                 >  1:305869/1‑238 (MQ=255)    attcttcgccttgtcgtgcgCCGGAAAGGGGCGATGTGAACGAAAAACCA‑GGTGGTCAGGGTGGATCCCGAGGGCGTGCAGTTCTTTTCCGGGGTGCTGGGCCGTTTGGGTTTTGCTTTCCCCAGCCGCCGCCATGTACGGCACCATCGTCAGGTGCAAAAACGGAGTGTGCTCACGGCCAATTTCAACAGCCATCGGGCGAATCGCTGCGAGGAA‑CGGCAAGGATTCGATATCACCTa                                                                                                                                                                                                                                             <  2:487661/226‑1 (MQ=255)           gCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCA‑GGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAA‑CGGCAAGGATTCGATATCACCTACTGTAc                                                                                                                                                                                                                                       >  1:556239/1‑238 (MQ=255)            cctgttcgttcGCTGGAACAGCGCGATCTGAACGACAAATCA‑GTATGTCAGGCTGTATACCGATGGAGAGCAGCTCTTTTTCAGAGGGCTGAGCCGGTTTGGTTTTGACTTCGCCAGACGCGGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCTAATCGCTTCGAGGAA‑CGGCAAGGATTCGATATCACCTACGGTAccg                                                                                                                                                                                                                                     <  1:412885/236‑1 (MQ=255)                ttcgttcgGCGGAACAGCGCGATCTGAACGACAAATCA‑GGATGTCAGGCTGGATACCGATGGAGAGCAGTTCTTTTACAGAGTGCTGGGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAA‑CGGCAAGGATTCGATATCACCTACTGTAccgc                                                                                                                                                                                                                                    <  1:159863/236‑1 (MQ=255)                ttcgttcgCCGGAACAGCGCGATCTGAACGACAAATCA‑GGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCCCCAGCGTCAGGTGCATAACCAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGCCGAATCGCTTCGAGGAA‑CGGCAAGGATTCGATATCACCTACTGTACCGCCGa                                                                                                                                                                                                                                 >  1:357185/1‑239 (MQ=255)                ttcgttcgCCGGAACAGCGCGATCTGAACGACAAATCA‑GGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGAGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCCCTTCGAGTAA‑CGGCAAGGATTCCATTTCACCTACTGTAcctcc                                                                                                                                                                                                                                   >  1:187204/1‑237 (MQ=255)                      ccgcGGAACAGGGCGATCTGGACG‑TAAATCAGGGGTGCTAGGGTGGGTACCGGTGGAGAGCAGCTCTTTTAAAGAGTGCGGAGCCGGTTTGGTTTTGACTTCCCCAGACGCTGCAATGTCCGGACCCAGCGGCAGGTGCATAAACAGAGTGTGCTCACGGCAAATTTCAACAGCCATCTGGCGAATCGCTTGGAGGAA‑CGGCAAGGATTCGATATCCCCTACTGTACCGCCGATTTCTa                                                                                                                                                                                                                           <  2:214619/236‑1 (MQ=255)                       ccgggaaacggggtctCTGAACGAAAAATCA‑GGGTGTCGGGCTGGATCCGGATGGAGAGCACCTCTTTTAAAGAGTGCTGGGGCGTTTTGGTTTTGACTCCACCAGACGGTGGCAGGTACGGCCCCAGGGTCAGGGGCATAAACAGAGTGTGCCCACGGACAATTTCAACAGCCATCTGGCGAATCGTTTCGAGGAA‑CGGCAAGGATTCGATATCTCCTACTGTACCGCCGATTTCTAc                                                                                                                                                                                                                          <  2:892035/224‑1 (MQ=255)                                 cgcgATCTGAACGACAAATCA‑GGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGCATCGCTTCGGGGAA‑CGGCAAGGATTCGCAATTCCCCACTTTACCCCCGATTTCTACCAGaacaca                                                                                                                                                                                                                 >  2:130803/1‑235 (MQ=255)                                             gaaaaaTCC‑GGATGTCAGGCTGGATACCGAGGGAGGGCACCTTTTTTACAGAGTGCTGATTCGGTTTGGTTTTGACTTCCCCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAA‑CGGCAAGGATTCGATATCACCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCg                                                                                                                                                                                                    <  1:659632/236‑1 (MQ=255)                                                             ggctcggtattccGGTGGAGGCCCGCTTTTTTACAGGGTGCTGGGTCGGTTTGGTTTTGACTTCCCCAGACGCTGCCATGTCGGGCACCAGCGTCAGGGGCATAAACAGGGTGTGCCCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAA‑CGGCAAGGATTCGATATCACCTCCTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCAc‑g                                                                                                                                                                                    <  2:321661/228‑1 (MQ=255)                                                             cAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAA‑CGGCAAGGATTCGATATCACCTACTGTACCGCCGATTTCTACCAGGACTACGTCCTGACCTTCGCCACCTTCCAGCAc                                                                                                                                                                                      >  1:475838/1‑238 (MQ=255)                                                                        ccGAGGGAGAGCCGCTTGTTTACAGAGTTCGTAGTCGGTTTGGTTGTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAA‑CGGCAAGGATTCGATATCACCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCAC‑GCGCTCTTTGAt                                                                                                                                                                         <  1:557489/239‑1 (MQ=255)                                                                           aTGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGGGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAA‑CGGCACGGATTCGCAATCACCTCCTTTAACGCCGATTTCTTCCAGACACACGCCATGACCTCCGCCCCCTccccgacccgcgcctttagttgg                                                                                                                                                                       >  2:479020/1‑219 (MQ=255)                                                                                   cagcTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAA‑CGGCAAGGATTCGATATCACCTACTGTACCGCCGCTTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCCC‑GCCCTCTTTGATTGGGttgggga                                                                                                                                                              >  2:812605/1‑234 (MQ=255)                                                                                   cagcTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAA‑CGGCAAGGATTCGATATCACCTACTGTACCCCCGATTTCTACCCGAACTACTGCCTTACCTTCCCCACCTTCCAGCCC‑CCGCCCTTTTATTGCGGTAGTgg                                                                                                                                                              >  2:486642/1‑238 (MQ=255)                                                                                      ctctTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAA‑CGGCAAGGATTCGAAATCACCTACTGTACCGCCGATTTCTACCGGAACTACGTCATGACCTTCCCCACCTTCCTGCAC‑GCCCCCTTTGGTTGCGTTATTGAtgt                                                                                                                                                           >  2:586293/1‑239 (MQ=255)                                                                                       tctTTGACTGAGTGCTGGTTCGGTTTGGTTTTGACTTCCCCAGACGCTGCCATGTACGCCTCCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCAAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAA‑CGGCAAGGATTCGATATCACCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCGC‑GCGCTCTTTGATTGCGTTAGTGAt                                                                                                                                                             <  1:160649/236‑1 (MQ=255)                                                                                       ctttttAAAGAGTGCTGGGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAA‑CGGCAAGGATTCGATATCACCTCCTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCAC‑GCGCTCTTTGATTGCGTTAGTGAtgt                                                                                                                                                           <  1:161274/236‑1 (MQ=255)                                                                                           ttACGGGGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGCCGCGGCCATGTACGGCCCCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAA‑CGGCAAGGATTCGATATCACCTACGGGACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCAC‑GCGCTCTTTGATTGCGTTAGTGATGTGCGGa                                                                                                                                                      <  1:641595/239‑1 (MQ=255)                                                                                           ttACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAA‑CGGCAAGGATTCGATATCACCTACTGTACCGCCGATTTCTACCAGTACTACGTCATTACCTTCGCCACCTTCCAGCAC‑GCGCCCTTTGATTGCGTTAGTGATGTGCGGa                                                                                                                                                      >  1:422358/1‑239 (MQ=255)                                                                                                cgggggtGAGCCGGTTTGGTTTTGCCCTCACCAGACGCTGCCATGTACGGCCCCAGCGTCGGGTGCATAAACAGGGTGTGCTCACGGCAAATTTCAACAGCCATCTGGCAAATCGCTTCGAGGAA‑CGGCAAGGATTCGATATCCCCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCAC‑GCGCTCTTTGATTGCGTTAGTGATGTGCGGAATAAc                                                                                                                                                 <  1:843425/233‑1 (MQ=255)                                                                                                   tGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAA‑CGGCAAGGATTCGATATCACCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCCGCCC‑GCGCCCTTTGATTACGTTAGTGATGTGGGGAATAACCt                                                                                                                                               >  2:245473/1‑238 (MQ=255)                                                                                                   tGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAA‑AGGCAAGGATTCGATATCCCCCACTGTACCGCCGATTTTTACCCGTACTACGCCATGACCTTCGCCACCTTCCCGCCC‑GCGCCCTTTGTTTGCGTTAGTGAtgggggcaaaaacag                                                                                                                                               >  2:872823/1‑225 (MQ=255)                                                                                                           cGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAA‑CGGCAAGGATTCGATATCACCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCTCCACCTTCCAGCAC‑GCGCTCTTTTATTGCGTTAGTGAGGTGCGGAGAAACCTGCACgggt                                                                                                                                       >  2:659991/1‑236 (MQ=255)                                                                                                             gtttggtttTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCCTAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCAACTGGCGAATCGCTTCGAGGAA‑CGGCAAGGATTCGGAATCACCTACTGTCCCGCCGTTTTCTACCACAACTACGACACGACCTTCGCACCCATCCAGCCC‑CCCCTCTTTGATTGAGTTGGGCATAGGCGGAAAAActcagccgggtggt                                                                                                                                    >  2:843425/1‑226 (MQ=255)                                                                                                             gtttggtttTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGGA‑AGGCAAGGATTCGATATCACCTACGGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCCCCTTCCA‑CACAGCCCTCTTTGTTTGCGTTATGTATGTGCGGGAAAACCTGCCAGGTTgg                                                                                                                                     >  2:160159/1‑237 (MQ=255)                                                                                                             gtttggtttTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAA‑CGGCAAGGATTCGATATCACCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCAC‑GCGCTCTTTGATTGCGTTAGTGATGTGCGGAATAACCTGCACGGTTgcg                                                                                                                                    >  2:14604/1‑239 (MQ=255)                                                                                                             gtttggtttTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAA‑CGGCAAGGATTCGATATCACCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCAC‑GCGCCCTTTGATTGCGTTAGTGATGTTCGGAATAACCTGCCCGGTTgc                                                                                                                                     >  1:774017/1‑238 (MQ=255)                                                                                                               ttggtttTGACTCCACCAAGCGCTGCCATGTCCGGCACCAGCGTCAGGTGCATAAACAGAGTTGGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAA‑CGGCAAGGATTCGATATCACCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCAC‑GCGCTCTTTGATTGCGTTAGTGATGTGCGGAATAACCTGCACGGTTGCGcc                                                                                                                                  <  1:340669/239‑1 (MQ=255)                                                                                                                      tgccttccccaggccgttccctGTCCGGCACCAGCGTCAGGTGCATAAACATATGTTGCTCACGGCAAATTACAACGTCCATCGGGCGAATCGCTTGGAGGAA‑CGGAAAGGATTCGATACCCCCTCCTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCAC‑GCGCTCTTTGATTGCGTTAGTGATGTGCGGAATAACCTGCACGGTTGCGCCGAGGTAg                                                                                                                           <  2:286656/218‑1 (MQ=255)                                                                                                                                   ccgtGCCATGTACGGCCCCAGCGTCAGGTGCATAAACAGGGGGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAA‑CGGCAAGGATTCGATATCACCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCAC‑GCGCTCTTTGATTGCGTTAGTGATGTGCGGAATAACCTGCACGGTTGCGCCGAGGTAGTCACCgcggc                                                                                                                 <  1:130803/233‑1 (MQ=255)                                                                                                                                   cGCTGCCATGTACGGCCCCAGCGTCAGGTGCATAAACGGGGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGTTTCGAGGAA‑CGGCAAGGATTCGATATCCCCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCAC‑GCGCTCTTTGATTGCGTTAGTGATGTGCGGAATAACCTGCACGGTTGCGCCGAGGTAGTCACCgcggc                                                                                                                 <  1:116895/236‑1 (MQ=255)                                                                                                                                      tGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAA‑CGGCAAGGATTCGATATCACCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCAC‑GCGCCCTTTGGTTGCGTTAGTGGTGTGCGGAATAACCTGCACGGTTGCGCCGAGGACGTCACCGCGGGGTCCtt                                                                                                           >  2:59555/1‑239 (MQ=255)                                                                                                                                                  cACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAA‑CGGCAAGGATTCGATATCACCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCAC‑GCGCTCTTTGATTGCGTTAGTGATGTGCGGAATAACCTGCACGGTTGCGCCGAGGCAGTCACCGCGGCGTTCTTTACGCAGaagt                                                                                                >  1:114939/1‑236 (MQ=255)                                                                                                                                                    cccgcgcCAGGTTCTAAAACAGGGTGTGCTCACGGCCAATTCCAACAGCCATCTGGCGAATCGCTTCGAGGAA‑CGGAAAGGATTCGATACCCCCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCAC‑GCGCTCTTTGATTGCGTTAGTGATGTGCGGAATAACCTGCACGGTTGCGCCGAGGTAGTCACCGCGGCGTTCTTTACGCAGAACGTCa                                                                                             <  2:422358/236‑1 (MQ=255)                                                                                                                                                      atgttcgtgTGCAAAAACCGGGTGGGCTCACGGCCACTTGCAACACCCTGCTGGCGTATCGTTGCGAGGAA‑CGGCAAGGTGTGGCTATCACCTCCTGTCCCGCCGATTTCTACCAGTACTACGTCATGACCTGCGCCCCCTTCCAGCAC‑GCGCTCTTTGATTGCGTTAGTGATGTGCGGAATAACCTGCACGGTTGCGCCGAGGTAGTCACCGCGGCGTTCTTTACGCAGACCGGCag‑a                                                                                          <  2:604514/231‑1 (MQ=255)                                                                                                                                                       ggtccgggggcataaaagggtggtgCCCACGGCAATTTTAAACCGCCTTCTGGC‑CATCGCTTTGGGGAACCGGCAGGATTTCGAGACCCCCTGCTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCAC‑GCGCTCTTTGATTGCGTTAGTGATGTGCGGAATAACCTGCACGGTTGCGCCGAGGTAGTCACCGCGGCGTTCTTTACGCAGAACGTCag‑ag                                                                                         <  2:114939/217‑1 (MQ=255)                                                                                                                                                       gCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAA‑CGGCAAGGATTCGATATCACCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCAC‑GCGCTCTTTGATTGCGTTAGTGATGTGCGGAATAACCTGCACGGTTGCGCCGAGGTAGTCACCGCGGCGTTCTTTACGCAGAACGTCcg‑ag                                                                                         >  1:753963/1‑239 (MQ=255)                                                                                                                                                       gCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAA‑CGGCAAGGATTCGATATCACCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGAAC‑GCGCTCTTTGATTGCGTTAGTGATGTGCGGAATAACCTGCACGGTTGCGCCGAGGTAGTCACCGCGGCGTTCTTTACGCCGAACGTGag‑ag                                                                                         >  1:542762/1‑239 (MQ=255)                                                                                                                                                                 cATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAA‑CGGCAAGGATTCGATATCACCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCAC‑GCGCTCTTTGATTGCGTTAGTGATGTGCGGAATAACCTGCACGGTTGCGCCGGGGTAGTCACCGCGGCGTTCTTTCCGCAGAACG‑CAGAAGTAGATACGa                                                                                >  1:857188/1‑238 (MQ=255)                                                                                                                                                                         gagTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAA‑CGGCAAGGATTCGATATCACCTACTGTACCGCCGATTTCTAACAGTACTACGTCATGACCTTCGCCACCTTCCAGAAC‑GCGCTCTTTGATTGCGTTAGTGATGTGCGGAATAACCTGCACGGTTGCGCCGAGGTAGTCGCCGCGGCGATCTTTCCGCAGAAGGTCAG‑GTTCGGTACGACCCgtgggg                                                                       >  1:885355/1‑237 (MQ=255)                                                                                                                                                                         gagTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAACCGCTTCGAGGAA‑CGGCAAGGATTCGATATCACCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCAC‑GCGCTCTTTGATTGCGTTAGTGATGTGCGGAATAACCTGCACGGTTGCGCCGAGGGAGTCACCGCGGCGTTCTTTACGCAGAACGTCAG‑AGAAGATACGAccccggggg                                                                       >  2:754504/1‑233 (MQ=255)                                                                                                                                                                              tgCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAA‑CGGCAAGGATTCGATATCACCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCAC‑GCGCTCTTTGATTGCGTTAGTGATGTGCGGAATAACCTGCACGGTTGCGCCGAGGTAGTCACCGCGGCGTTCTTTACGCAGAACGTCAG‑AGTAGATACGACCCGTGGTGAAgtt                                                                  >  1:147448/1‑239 (MQ=255)                                                                                                                                                                                cTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAA‑CGGCAAGGATTCGATATCACCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCAC‑GCGCTCTTTGATTGCGTTTGTGATGTGCGGAATAACCTGCACGGTTGCGCCGGGGTAGTCACCGCGGCGTTCTTTACGCAGAACGTCAG‑AGTAGATACGACCCGTGGTGAAgttgt                                                                >  1:660665/1‑239 (MQ=255)                                                                                                                                                                                cTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAA‑CGGCAAGGATTCGATATCACCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCAC‑GCGCTCTTTGATTGCGTTAGTGATGTGCGGAATAACCTGCACGGTTGCGCCGAGGTAGTCACCGCGGCGTTCTTTACGCAGACCGTCAG‑AGTAGATACGACCCGGGGTGAAgttgt                                                                >  1:490282/1‑239 (MQ=255)                                                                                                                                                                                      gCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAA‑CGGCAAGGATTCGATATCACCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCAC‑GCGCTCTTTGATTCCGTTAGTGATGTGCGGACTAACCTGCACGGTTGCGCCGAGGCAGTCACCGCGGCGTTCTTTACGCAGAACGTCCG‑GGTAGATACGACCCGTGGTGAAGTTGTTccggc                                                          >  1:127823/1‑239 (MQ=255)                                                                                                                                                                                        cAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAA‑CGGCAAGGATTCGATATCACCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCAC‑GCGCTCTTTGATTGCGTTAGTGATGTGCGGAATAACCTGCACGGTTGCGCCGAGGTAGTCACCGCGGCGTTCTTTACGCAGAACGTCAG‑AGTAGATACGACCCGTGGTGAAGTTGTTgcggcg                                                         <  1:14604/238‑1 (MQ=255)                                                                                                                                                                                         aaTTTCAACCGCCATCTGGCGAATCGCTTCGAGGAA‑CGGCAAGGATGCGATATCACCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCAC‑GCGCTCTTTGATTGCGTTAGTGATGTGCGGAATAACCTGCACGGTTGCGCCGAGGTAGTCACCGCGGCGTTCTTTACGCAGAACGTCAG‑AGTAGATACGACCCGTGGTGAAGTTGTTGCggcggc                                                       <  1:678437/239‑1 (MQ=255)                                                                                                                                                                                          aTTTCAACAGCCATTTGGCGAATCGCTTCGAGGAA‑CGGCAAGGATTCGATATCACCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCAC‑GCGCTCTTTGATTGCGTTAGTGATGTGCGGAATAACCTGCACGGTTGCGCCGAGGTAGTCACCGCGGCGTTCTTTACGCAGAACGTCAG‑AGTAGATACGACCCGTGGTGAAGTTGTTGCGGCGGCt                                                      <  1:245473/239‑1 (MQ=255)                                                                                                                                                                                          aTTTCAACAGCCATCTGGCGAATCGCTTCGAGGAA‑CGGCAAGGATTCGATATCACCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCAC‑GCGCTCTTTGATTGCGTTAGTGATGTGCGGAATAACCTGCACGGTTGCGCCGAGGTAGTCACCGCGGCGTTCTTTACGAAGAACGTCAG‑GGAAAATACGACCCGTGGGGGAGGTGTTgggggtggc                                                      >  2:425433/1‑234 (MQ=255)                                                                                                                                                                                          aTTTCAACAGCCATCTGGCGAATCGCTTCGAGGAA‑CGGCAAGGATTCGATATCACCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCAC‑GCGCTCTTTGATTGCGTTAGTGATGTGCGGAATAACCTGCACGGTTGCGCCGAGGTAGTCACCGCGGCGGTCTTTACGCAGAACGTCAG‑AGTAGATCCGACCCGTGGTGAAGTTGTTGCggcggc                                                       >  1:353797/1‑238 (MQ=255)                                                                                                                                                                                                    cctccttgcgactcgttggggggAA‑CGGAAGGGTTTCGATTCCCCCCACTGTCCCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCAC‑GCGCTCTTTGATTGCGTTAGTGATGTGCGGAATAACCTGCACGGTTGCGCCGAGGTAGTCACCGCGGCGTTCTTTACGCAGAACGTCAG‑AGTAGATACGACCCGTGGTGAAGTTGTTGCGGCGGCTCATTTTGGTa                                            <  2:774017/218‑1 (MQ=255)                                                                                                                                                                                                    ccATCGTGCGAATCGCTTCGAGGAA‑CGGAAAGGTTTCGATACCACCTACTGTCCCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCAC‑GCGCTCTTTGATTGCGTTAGTGATGTGCGGAATAACCTGCACGGTTGCGCCGAGGTAGTCACCGCGGCGTTCTTTACGCAGAACGTCAG‑AGTAGATACGACCCGTGGTGAAGTTGTTGCGGCGGCTCATTTTGGTa                                            <  1:739401/239‑1 (MQ=255)                                                                                                                                                                                                         tGGCGAATCGCTTCGAGGAA‑CGGCAAGGATTCGATATCACCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCAC‑GCGCTCTTTGATTGCGTTAGTGATGTGCGGAATAACCTGCACGGTTGCGCCGAGGTAGTCACCGCGGCGTTCTTTACGCAGAACGTCAG‑AGTAGATACGCACCGTGGGGAAGGTGTTGGGGCGGCTCATTTTggggaggaa                                       >  2:337603/1‑232 (MQ=255)                                                                                                                                                                                                                  aaagtacggtaccGGAAAGGTTTCGATTTCCCCTTCTGTCCCGCCGATTTCTCCCAGTACGACGTCATGGCCTTCGCCACCTTCCAGCAC‑GCGCTCTTTGATTGCGTTAGTGATGTGCGGAATAACCTGCACGGTTGCGCCGAGGTAGTCACCGCGGCGTTCTTTACGCAGAACGTCAG‑AGTAGATACGACCCGTGGTGAAGTTGTTGCGGCGGCTCATTTTGGTACGAATGAAACGCt                               <  2:406184/227‑1 (MQ=255)                                                                                                                                                                                                                                   aGGATTCGAATTCCCCTACGGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCAC‑GCGCTCTTTGATTGCGTTAGTGATGTGCGGAATAACCTGCACGGTTGCGCCGAGGTAGTCACCGCGGCGTTCTTTACGCAGAACGTCAG‑AGTAGATACGACCCGTGGTGAAGTTGTTGCGGCGGCTCATTTTGGTACGAATGAAACGCTCGTAGTGCCccaggtcc              <  1:754504/239‑1 (MQ=255)                                                                                                                                                                                                                                      aTTCGATATCACCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCAC‑GCGCTCTTTGATTGCGTTAGTGATGTGCGGAATAACCTGCACGGTTGCGCCGAGGTAGTCACCGCGGCGTTCTTTACGCAGAACGTCAG‑AGTAGATACGACCCGTGGTGAAGTTGTTGCGGCGGCTCATTTTTGTACGAATGAAACGCTCGTATGGccccccgccttct           >  1:406184/1‑230 (MQ=255)                                                                                                                                                                                                                                               cACCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCAC‑GCGCTCTTTGATTGCGTTAGTGATGTGCGGAATAACCTGCACGGTTGCGCCGAGGTAGTCACCGCGGCGTTCTTTACGCAGAACGTCAG‑AGTAGATACGACCCGTGGTGAAGTTGTTGCGGCGGCTCCTTTTGGGACGGATAAAACGCCCGTAGTGCCCCCGGTCCCGGGCGGTTTCa  >  1:150357/1‑239 (MQ=255)                                                                                                                                                                                                                                                aCCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCAC‑GCGCTCTTTGATTGCGTTAGTGATGTGCGGAATAACCTGCACGGTTGCGCCGAGGTAGTCACCGCGGCGTTCTTTACGCAGACGGTCAG‑AGTAGATACGACCCGGGGTGAAGTTGTTGCGGGGGCCCATTTTGGTCCGAATTAAACGCTCGTAGTGCCCCAGGGCACGGGCGGTTTCa  >  1:238787/1‑238 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                TGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCA‑GGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAA‑CGGCAAGGATTCGATATCACCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCAC‑GCGCTCTTTGATTGCGTTAGTGATGTGCGGAATAACCTGCACGGTTGCGCCGAGGTAGTCACCGCGGCGTTCTTTACGCAGAACGTCAG‑AGTAGATACGACCCGTGGTGAAGTTGTTGCGGCGGCTCATTTTGGTACGAATGAAACGCTCGTAGTGCCCCAGGTCCAGGTCGGTTTCA  >  NC_000913/2908990‑2909464

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