Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,639,635	Δ1 bp	6.9%	coding (941/1275 nt)	hisS ←	histidyl tRNA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,639,635	0	A	.	6.9%	204.5 / 6.1	58	coding (941/1275 nt)	hisS	histidyl tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base A (26/28);  new base . (2/2);  total (28/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

CTGCCGTTTGCTCACCAGAGCGCAAATCCTTCACTACTGCTGTGCCGTTAGCCACTTCAGACTCACCCAGCACCACAGCAACGCGGGCACCCCATTTATCAGCACGGGCAAACTGTTTCTTAAAGTTGCCGCCGCCGTGGTTGGTCATCAATTTCACGCCCGGTAATTCATCACGCAGACGCTCAGCTAATGCCATAGCCGCAGATTGTGTATCAGCACCTGAAGCCACCAGGTATATATCGACAACAGGATCGGCTTTAAATTCCGGATTAACGGCCTGTACTAACAATACAAGACGTTCGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTCCACAAGACCGTCATAACGACCGCCTGCACACACGGTGCCCTGGGAGCCGAGACTGTTAGTCACCCACTCGAAAACGGTACGGTTGTAGTAATCCAGACCACGCACCAGACGCTGGTTTACGGTGTAAGCGATCCCCGCGCTCTCCAGCAGTTTGCACAGACCGGCAAAATGCTCACGAGATTCCTCGTCCAGATAGTCACCTAATGCCGGAGCGTCG  >  NC_000913/2639348‑2639916                                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                                          cTGCCGTTTGCTCACCAGAGCGCAAATCCTTCACTACTGCTGTGCCGTTAGCCACTTCAGACTCACCCAGCACCACAGCAACGCGGGCACCCCATTTATCAGCACGGGCAAACTGTTTCTTAAAGTTGCCGCCGCCGTGGTTGGTCATCAATTTCACGCCCGGTAATTCATCACGCAGACGCTCAGCTAATGCCATAGCCGCAGATTGTGTATCAGCACCTGAAGCCACCAGGTATATATCGACAACAGGATCGGCTTTAAATTCCGGATTAACGGCCTGTACTAACaa                                                                                                                                                                                                                                                                                          <  1:471998/289‑1 (MQ=255) cTGCCGTTTGCTCACCAGAGCGCAAATCCTTCACTACTGCTGTGCCGTTAGCCACTTCAGACTCACCCAGCACCACAGCAACGCGGGCACCCCATTTATCAGCACGGGCAAACTGTTTCTTAAAGTTGCCGCCGCCGTGGTTGGTCATCAATTTCACGCCCGGTAATTCATCACGCAGACGCTCAGCTAATGCCATAGCCGCAGATTGTGTATCAGCACCTGAAGCCACCAGGTATATATCGACAACAGGATCGGCTTTAAATTCCGGATTAACGGCCTGTACTAACaa                                                                                                                                                                                                                                                                                          >  1:336553/1‑289 (MQ=255)                   agCGCAAATCCTTCACTACTGCTGTGCCGTTAGCCACTTCAGACTCACCCAGCACCACAGCAACGCGGGCACCCCATTTATCAGCACTGGCAAACTGTTTCTTAAAGTTGCCGCCGCCGTGGTTGGTCATCAATTTCACGCCCGGTAATTCATCACGCAGACGCTCAGCTAATGCCATAGCCGCAGATTGTGTAACAGCACCTGAAGCCACCAGGTATATATCGAAAACAGGATCAGCTTTAAATTCCGGATTAACGGCCTGTACTAACAATACAAGACCTTCGAGGc                                                                                                                                                                                                                                                                         >  1:689285/1‑288 (MQ=255)                                          gtgCCGTTAGCCACTTCAGACTCACCCAGCACCACAGCAACGCGGGCACCCCATTTATCAGCACGGGCAAACTGTTTCTTAAAGTTGCCGCCGCCGTGGTTGGTCATCAATTTCACGCCCGGTAATTCATCACGCAGACGCTCAGCTAATGCCATAGCCGCAGATTGTGTATCAGCACCTGAAGCCACCAGGTATATATCGACAACAGGATCGGCTTTAAATTCCGGATTAACGACCTGTACTAACAATACAAGACGTTCGAGGCCCATAGCAAAACCGACAGCCGgtg                                                                                                                                                                                                                                                 >  1:442279/1‑289 (MQ=255)                                                ttAGCCACTTCAGACTCACCCAGCACCACAGCAACGCGGGCACCCCATTTATCAGCACGGGCAAACTGTTTCTTAAAGTTGCCGCCGCCGTGGTTGGTCATCAATTTCACGCCCGGTAATTCATCACGCAGACGCTCAGCTAATGCCATAGCCGCAGATTGTGTATCAGCACCTGAAGCCACCAGGTATATATCGACAACAGGATCGGCTTTAAATTCCGGATTAACGGCCTGTACTAACAATACAAGACGTTCGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCAc                                                                                                                                                                                                                                           >  1:462266/1‑289 (MQ=255)                                                                  cccAGCACCACAGCAACGCGGGCACCCCATTTATCAGCCCGGGCAAACTGTTTCTTAAAGTTGCCGCCGCCGTGGTTGGTCATCAATTTCACGCCCGGTAATTCATCACGCAGACGCTCAGCTAATGCCATAGCCGCAGATTGTGTATCAGCACCTGAAGCCACCAGGTCTATATCCACAACAGGATCGGCTTTAAATTCCGGATTAACGGCCTGTCCTACCAATACAAGACGTTCGAGGCCCATAGCCAAACCGACAGCCGGCGTTGCCCGACCGCCCAGTCGTTcc                                                                                                                                                                                                                          >  2:153560/1‑288 (MQ=255)                                                                  cccAGCACCACAGCAACGCGGGCACCCCATTTATCAGCACGGGCAAACTGTTTCTTAAAGTTGCCGCCGCCGTGGTTGGTCATCAATTTCACGCCCGGTAATTCATCACGCAGACGCTCAGCTAATGCCATAGCCGCAGATTGTGTATCAGCACCTGAAGCCACCAGGTATATATCGACAACAGGATCGGCTTTAAATTCCGGATTAACGGCCTGTACTAACAATACAAGACGTTCGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTcc                                                                                                                                                                                                                          >  1:277026/1‑288 (MQ=255)                                                                        accacAGCAACGCGGGCACCCCATTTATCAGCACGGGCAAACTGTTTCTTAAAGTTGCCGCCGCCGTGGTTGGTCATCAATTTCACGCCCGGTAATTCATCACGCAGACGCTCAGCTAATGCCATAGCCGCAGATTGTGTATCAGCACCTGAAGCCACCAGGTATATATCGACAACAGGATCGGCTTTAAATTCCGGATTAACGGCCTGTACTAACAATACAAGACGTTCGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTCCACaa                                                                                                                                                                                                                      <  1:351720/286‑1 (MQ=255)                                                                        accacAGCAACGCGGGCACCCCATTTATCAGCACGGGCAAACTGTTTCTTAAAGTTGCCGCCGCCGTGGTTGGTCATCAATTTCACGCCCGGTAATTCATCACGCAGACGCTCAGCTAATGCCATAGCCGCAGATTGTGTATCAGCACCTGAAGCCACCAGGTATATATCGACAACAGGATCGGCTTTAAATTCCGGATTAACGGCCTGTACTAACAATACAAGACGTTCGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTCCACAAGAc                                                                                                                                                                                                                   <  2:277026/289‑1 (MQ=255)                                                                          cacaGCAACGCGGGCACCCCATTTATCAGCACGGGCAAACTGTTTCTTAAAGTTGCCGCCGCCGTGGTTGGTCATCAATTTCACGCCCGGTAATTCATCACGCAGACGCTCAGCTAATGCCATAGCCGCAGATTGTGTATCAGCACCTGAAGCCACCAGGTATATATCGACAACAGGATCGGCTTTAAATTCCGGATTAACGGCCTGTACTAAC‑ATACAAGACGTTCGAGGCCCATAGCAAAACCGACAGc                                                                                                                                                                                                                                                      <  1:151919/251‑1 (MQ=255)                                                                          cacaGCAACGCGGGCACCCCATTTATCAGCACGGGCAAACTGTTTCTTAAAGTTGCCGCCGCCGTGGTTGGTCATCAATTTCACGCCCGGTAATTCATCACGCAGACGCTCAGCTAATGCCATAGCCGCAGATTGTGTATCAGCACCTGAAGCCACCAGGTATATATCGACAACAGGATCGGCTTTAAATTCCGGATTAACGGCCTGTACTAAC‑ATACAAGACGTTCGAGGCCCATAGCAAAACCGACAGc                                                                                                                                                                                                                                                      >  2:151919/1‑251 (MQ=255)                                                                                           cccATTTATCAGCACGGGCAAACTGTTTCTTAAAGTTGCCGCCGCCGTGGTTGGTCATCAATTTCACGCCCGGTAATTCATCACGCAGACGCTCAGCTAATGCCATAGCCGCAGATTGTGTATCAGCACCTGAAGCCACCAGGTATATATCGACAACAGGATCGGCTTTAAATTCCGGATTAACGGCCTGTACTAACAATACAAGACGTTCGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTCCACAAGACCGTCATAACGACCGCCTg                                                                                                                                                                                                 >  2:541771/1‑288 (MQ=255)                                                                                           cccATTTATCAGCACGGGCAAACTGTTTCTTAAAGTTGCCGCCGCCGTGGTTGGTCATCAATTTCACGCCCGGTAATTCATCACGCAGACGCTCAGCTAATGCCATAGCCGCAGATTGTGTATCAGCACCTGAAGCCACCAGGTATATATCGACAACAGGATCGGCTTTAAATTCCGGATTAACGGCCTGTACTAACAATACAAGACGTTCGAGGCCCATAGCAAAACCGACAGCCGGTGGTGCACGACCGCCCAGGTGTTCCACAAGACAGTTATACCGAACGCCTg                                                                                                                                                                                                 >  2:36631/1‑288 (MQ=255)                                                                                                       cACGGGCAAACTGTTTCTTAAAGTTGCCGCCGCCGTGGTTGGTCATCAATTTCACGCCCGGTAATTCATCACGCAGACGCTCAGCTAATGCCATAGCCGCAGATTGTGTATCAGCACCTGAAGCCACCAGGTATATATCGACAACAGGATCGGCTTTAAATTCCGGATTAACGGCCTGTACTAACAATACAAGACGTTCGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCCCGACCGCCCAGTTGTTCCACAAGACCGTCATAACGACCGCCTGCCCACACGGTGc                                                                                                                                                                                     >  2:203111/1‑288 (MQ=255)                                                                                                       cACGGGCAAACTGTTTCTTAAAGTTGCCGCCGCCGTGGTTGGTCATCAATTTCACGCCCGGTAATTCATCACGCAGACGCTCAGCTAATGCCATAGCCGCAGATTGTGTATCAGCACCTGAAGCCACCAGGTATATATCGACAACAGGATCGGCTTTAAATTCCGGATTAACGGCCTGTACTAACAATACAAGACGTTCGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTCCACAAGACCGTCATAACGACCTCCTGCACACCCGGTGc                                                                                                                                                                                     >  1:595148/1‑288 (MQ=255)                                                                                                                  tGTTTCTTAAAGTTGCCGCCGCCGTGGTTGGTCATCAATTTCACGCCCGGTAATTCATCACGCAGACGCTCAGCTAATGCCATAGCCGCAGATTGTGTATCAGCACCTGAAGCCACCAGGTATATATCGACAACAGGATCGGCTTTAAATTCCGGATTAACGGCCTGTACTAACAATACAAGACGTTCGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTCCACAAGACCGTCATAACGACCGCCTGCACACACGGTGCCCTGGGAGCCga                                                                                                                                                                         <  1:553290/289‑1 (MQ=255)                                                                                                                    tttCTTAAAGTTGCCGCCGCCGTGGGTGGTCATCAATTTCACGCCCGGTAATTCATCACGCAGACGCTCAGCTAATGCCATAGCCGCAGATTGTGTATCAGCACCTGAAGCCACCAGGTATATATCGACAACAGGATCGGCTTTAAATTCCGGATTAACGGCCTGTACTAACAATACAAGACGTTCGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTCCACAAGACCGTCATAACGACCGCCTGCACACACGGTGCCCTGGGAGCCgaga                                                                                                                                                                       <  1:350397/289‑1 (MQ=255)                                                                                                                     ttCTTAAAGTTGCCGCCGCCGTGGTTGGTCATCAATTTCACGCCCGGTAATTCATCACGCAGACGCTCAGCTAATGCCATAGCCGCAGATTGTGTATCAGCACCTGAAGCCACCAGGTATATATCGACAACAGGATCGGCTTTAAATTCCGGATTAACGGCCTGTACTAACAATACAAGACGTTCGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTCCACAAGACCGTCATAACGACCGCCTGCACACACGGTGCCCTGGGAGCCgaga                                                                                                                                                                       <  1:622098/288‑1 (MQ=255)                                                                                                                      tCTTAAAGTTGCCGCCGCCGTGGTTGGTCATCAATTTCACGCCCGGTAATTCATCACGCAGACGCTCAGCTAATGCCATAGCCGCAGATTGTGTATCAGCACCTGAAGCCACCAGGTATATATCGACAACAGGATCGGCTTTAAATTCCGGATTAACGGCCTGTACTAACAATACAAGACGTTCGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTCCACAAGACCGTCATAACGACCGCCTGCACACACGGTGCCCTGGGAGCCgag                                                                                                                                                                        >  1:101984/1‑286 (MQ=255)                                                                                                                                gccgccgccGTGGTTGGTCATCAATTTCACGCCCGGTAATTCATCACGCAGACGCTCAGCTAATGCCATAGCCGCAGATTGTGTATCAGCACCTGAAGCCACCAGGTATATATCGACAACAGGATCGGCTTTAAATTCCGGATTAACGGCCTGTACTAACAATACAAGACGTTCGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTCCACAAGACCGTCATAACGACCGCCTGCACACACGGTGCCCTGGGAGCCGAGACTGTTAGTCAc                                                                                                                                                            >  1:109869/1‑288 (MQ=255)                                                                                                                                 ccgccgccgTGGTTGGTCATCAATTTCACGCCCGGTAATTCATCACGCAGACGCTCAGCTAATGCCATAGCCGCAGATTGTGTATCAGCACCTGAAGCCACCAGGTATATATCGACAACAGGATCGGCTTTAAATTCCGGATTAACGGCCTGTACTAACAATACAAGACGTTCGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTCCACAAGACCGTCATAACGACCGCCTGCACACACGGTGCCCTGGGAGCCGAGACTGTTAGTCAcc                                                                                                                                                           <  2:336553/288‑1 (MQ=255)                                                                                                                                        cgTGGTTGGTCATCAATTTCACGCCCGGTAATTCATCACGCAGACGCTCAGCTAATGCCATAGCCGCAGATTGTGTATCAGCACCTGAAGCCACCAGGTATATATCGACAACAGGATCGGCTTTAAATTCCGGATTAACGGCCTGTACTAACAATACAAGACGTTCGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTCCACAAGACCGTCATAACGACCGCCTGCACACACGGTGCCCTGGGAGCCGAGACTGTTAGTCACCCACTCGaa                                                                                                                                                   <  1:357338/289‑1 (MQ=255)                                                                                                                                                gtCATCAATTTCACGCCCGGTAATTCATCACGCAGACGCGCAGCTAATGCCATAGCCGCAGATTGTGTATCAGCACCTGAAGCCACCAGGTATATATCGACAACAGGATCGGCTTTAAATTCCGGATTAACGGCCTGTACTAACAATACAAGACGTTCGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTCCACAAGACCGTCATAACGACCGCCTGCACACACGGTGCCCTGGGAGCCGAGACTGTTAGTCACCCACTCGAAAacggta                                                                                                                                            <  1:153560/288‑1 (MQ=255)                                                                                                                                                     caATTTCACGCCCGGTAATTCATCACGCAGACGCTCAGCTAATGCCATAGCCGCAGATTGTGTATCAGCACCTGAAGCCACCAGGTATATATCGACAACAGGATCGGCTTTAAATTCCGGATTAACGGCCTGTACTAACAATACAAGACGTTCGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTCCACAAGACCGTCATAACGACCGCCTGCACACACGGTGCCCTGGGAGCCGAGACTGTTAGTCACCCACTCGAAAACGGTACGGTTg                                                                                                                                      <  2:476924/289‑1 (MQ=255)                                                                                                                                                      aaTTTCACACCCGGTAATTCACCCCGCAGACGCTCAGCTAATGCAATAGCCGCAGATTGTGTATCAGCACCGGAAGCCACCAGGTATATATCGACAACAGGATCGGCTTTAAATTCCGGATTAACGGCCTGTACTAACAATACAAGACGTTCGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTCCGCAAGACCGTCATAACGACCGCCTGCACACACGGTGCCCTGGGAGCCGAGACTGTTAGTCACCCACTCGAAAacggtacggt                                                                                                                                        <  1:169376/286‑1 (MQ=255)                                                                                                                                                       aTTTCACGCCCGGTAATTCATCACGCAGACGCTCAGCTAATGCCATAGCCGCAGATTGTGTATCAGCACCTGAAGCCACCAGGTATATATCGACAACAGGATCGGCTTTAAATTCCGGATTAACGGCCTGTACTAACAATACAAGACGTTCGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTCCACAAGACCGTCATAACGACCGCCTGCACACACGGTGCCCTGGGAGCCGAGACTGTTAGTCACCCACTCGAAAACGGTACGGTTgta                                                                                                                                    <  1:36631/289‑1 (MQ=255)                                                                                                                                                       aTTTCACGCCCGGTAATTCATCACGCAGACGCTCAGCTAATGCCATAGCCGCAGATTGTGTATCAGCACCTGAAGCCACCAGGTATATATCGACAACAGGATCGGCTTTAAATTCCGGATTAACGGCCTGTACTAACAATACAAGACGTTCGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTCCACAAGACCGTCATAACGACCGCCTGCACACACGGTGCCCTGGGAGCCGAGACTGTTAGTCACCCACTCGAAAACGGTACGGTTgta                                                                                                                                    <  1:541771/289‑1 (MQ=255)                                                                                                                                                                  ggTAATTCATCACGCAGACGCTCAGCTAATGCCATAGCCGCAGATTGTGTATCAGCACCTGAAGCCACCAGGTATATATCGACAACAGGATCGGCTTTAAATTCCGGATTAACGGCCTGTACTAACAATACAAGACGTTCGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTCCACAAGACCGTCATAACGACCGCCTGCACACACGGTGCCCTGGGAGCCGAGACTGTTAGTCACCCACTCGAAAACGGTACGGTTGTAGTAATCCAGa                                                                                                                          >  1:137127/1‑288 (MQ=255)                                                                                                                                                                              cGCAGACGCTCAGCTAATGCCATAGCCGCAGATTGTGTATCAGCACCTGAAGCCACCAGGTATATATCGACAACAGGATCGGCTTTAAATTCCGGATTAACGGCCTGTACTAAC‑ATACAAGACGTTCGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTCCACAAGACCGTCATAACGACCGCCTGCACACACGGTGCCCTGGGAGCCGAGACTGTTAGTCACCCACTCGAAAacggtacggt                                                                                                                                        >  1:39475/1‑261 (MQ=255)                                                                                                                                                                              cGCAGACGCTCAGCTAATGCCATAGCCGCAGATTGTGTATCAGCACCTGAAGCCACCAGGTATATATCGACAACAGGATCGGCTTTAAATTCCGGATTAACGGCCTGTACTAAC‑ATACAAGACGTTCGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTCCACAAGACCGTCATAACGACCGCCTGCACACACGGTGCCCTGGGAGCCGAGACTGTTAGTCACCCACTCGAAAacggtacggt                                                                                                                                        <  2:39475/261‑1 (MQ=255)                                                                                                                                                                               gCAGACGCTCAGCTAATGCCATAGCCGCAGATTGTGTATCAGCACCTGAAGCCACCAGGTATATATCGACAACAGGATCGGCTTTAAATTCCGGATTAACGGCCTGTACTAACAATACAAGACGTTCGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTCCACAAGACCGTCATAACGACCGCCTGCACACACGGTGCCCTGGGAGCCGAGACTGTTAGTCACCCACTCGAAAACGGTACGGTTGTAGTAATCCAGACCACGCACCAGACg                                                                                                            <  1:479361/289‑1 (MQ=255)                                                                                                                                                                                       tCAGCTAATGCCATAGCCGCAGATTGTGTATCAGCACCTGAAGCCACCAGGTATATATCGACAACAGGATCGGCTTTAAATTCCGGATTAACGGCCTGTACTAACAATACAAGACGTTCGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTCCACAAGACCGTCATAACGACCGCCTGCACACACGGTGCCCTGGGAGCCGAGACTGTTAGTCACCCACTCGAAAACGGTACGGTTGTAGTAATCCAGACCACGCACCAGACGCTGGTTTa                                                                                                    <  1:353954/289‑1 (MQ=255)                                                                                                                                                                                        cAGCTAATGCCATAGCCGCAGATTGTGTATCAGCACCTGAAGCCACCAGGTATATATCGACAACAGGATCGGCTTTAAATTCCGGATTAACGGCCTGTACTAACAATACAAGACGTTCGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTCCACAAGACCGTCATAACGACCGCCTGCACACACGGTGCCCTGGGAGCCGAGACTGTTAGTCACCCACTCGAAAACGGTACGGTTGTAGTAATCCAGACCACGCACCAGACGCTGGtt                                                                                                      <  2:653354/286‑1 (MQ=255)                                                                                                                                                                                        cAGCTAATGCCATAGCCGCAGATTGTGTATCAGCACCTGAAGCCACCAGGTATATATCGACAACAGGATCGGCTTTAAATTCCGGATTAACGGCCTGTACTAACAATACAAGACGTTCGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTCCACAAGACCGTCATAACGACCGCCTGCACACACGGTGCCCTGGGAGCCGAGACTGTTAGTCACCCACTCGAAAACGGTACGGTTGTAGTAATCCAGACCACGCACCAGACGCTGGtt                                                                                                      >  1:653354/1‑286 (MQ=255)                                                                                                                                                                                                   aTAGCCGCAGATTGTGTATCAGCACCTGAAGCCACCAGGTATATATCGACAACAGGATCGGCTTTAAATTCCGGATTAACGGCCTGTACTAACAATACAAGACGTTCGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTCCACAAGACCGTCATAACGACCGCCTGCACACACGGTGCCCTGGGAGCCGAGACTGTTAGTCACCCACTCGAAAACGGTACGGTTGTAGTACTCCAGACCACGCACCAGACGCTGGTTTACGGTGTAAgaga                                                                                        >  1:659710/1‑289 (MQ=255)                                                                                                                                                                                                     aGCCGCAGATTGTGTATCAGCACCTGAAGCCACCAGGTATATATCGACAACAGGATCGGCTTTAAATTCCGGATTAACGGCCTGTACTAACAATACAAGACGTTCGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTCCACAAGACCGTCATAACGACCGCCTGCACACACGGTGCCCTGGGAGCCGAGACTGTTAGTCACCCACTCGAAAACGGTACGGTTGTAGTAATCCAGACCACGCACCAGACGCTGGTTTACGGTGTAAGCGAt                                                                                       >  1:145847/1‑288 (MQ=255)                                                                                                                                                                                                                 tgtATCAGCACCTGAAGCCACCAGGTATATATCGACAACAGGATCGGCTTTAAATTCCGGATTAACGGCCTGTACTAACAATACAAGACGTTCGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTCCACAAGACCGTCATAACGACCGCCTGCACACACGGTGCCCTGGGAGCCGAGACTGTTAGTCACCCACTCGAAAACGGTACGGTTGTAGTAATCCCGACCACGCACCAGACGCTGGTTTACGGTGTAAGCGATCCCCGCCCcctc                                                                           >  2:616286/1‑288 (MQ=255)                                                                                                                                                                                                                           ccTGAAGCCACCAGGTATATATCGACAACAGGATCGGCTTTAAATTCCGGATTAACGGCCTGTACTAACAATACAAGACGTTCGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTCCACAAGACCGTCATAACGACCGCCTGCACACACGGTGCCCTGGGAGCCGAGACTGTTAGTCACCCACTCGAAAACGGTACGGTTGTAGTAATCCAGACCACGCACCAGACGCTGGTTTACGGTGTAAGCGATCCCcgcg                                                                                >  2:518716/1‑273 (MQ=255)                                                                                                                                                                                                                           ccTGAAGCCACCAGGTATATATCGACAACAGGATCGGCTTTAAATTCCGGATTAACGGCCTGTACTAACAATACAAGACGTTCGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTCCACAAGACCGTCATAACGACCGCCTGCACACACGGTGCCCTGGGAGCCGAGACTGTTAGTCACCCACTCGAAAACGGTACGGTTGTAGTAATCCAGACCACGCACCAGACGCTGGTTTACGGTGTAAGCGATCCCcgcg                                                                                <  1:518716/273‑1 (MQ=255)                                                                                                                                                                                                                            cTGAAGCCACCAGGTATATATCGACAACAGGATCGGCTTTAAATTCCGGATTAACGGCCTGTACTAACAATACAAGACGTTCGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTCCACAAGACCGTCATAACGACCGCCTGCACACACGGTGCCCTGGGAGCCGAGACTGTTAGTCACCCACTCGAAAACGGTACGGTTGTAGTAATCCAGACCACGCACCAGACGCTGGTTTACGGTGTAAGCGATCCCCGCGCTCTCCAGCAGTTTGca                                                               <  1:257011/289‑1 (MQ=255)                                                                                                                                                                                                                                                  gACAACAGGATCGGCTTTAAATTCCGGATTAACGGCCTGTACTAACAATACAAGACGTTCGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTCCACAAGACCGTCATAACGACCGCCTGCACACACGGTGCCCTGGGAGCCGAGACTGTTAGTCACCCACTCGAACACGGTACGGTTGTAGTAATCCAGACCACGCACCAGACGCTGGTTTACGGTGTAAGCGATCCCCGCGCCCTCCAGCAGTTTGCACAGACCGGCAAAATGCTCACga                                         >  1:611942/1‑289 (MQ=255)                                                                                                                                                                                                                                                   acaacaGGATCGGCTTTAAATTCCGGATTAACGGCCTGTACTAACAATACAAGACGTTCGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTCCACAAGACCGTCATAACGACCGCCTGCACACACGGTGCCCTGGGAGCCGAGACTGTTAGTCACCCACTCGAAAacggta                                                                                                                                            >  2:399846/1‑189 (MQ=255)                                                                                                                                                                                                                                                   acaacaGGATCGGCTTTAAATTCCGGATTAACGGCCTGTACTAACAATACAAGACGTTCGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTCCACAAGACCGTCATAACGACCGCCTGCACACACGGTGCCCTGGGAGCCGAGACTGTTAGTCACCCACTCGAAAacggta                                                                                                                                            <  1:399846/189‑1 (MQ=255)                                                                                                                                                                                                                                                   acaacaGGATCGGCTTTAAATTCCGGATTAACGGCCTGTACTAACAATACAAGACGTTCGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTCCACAAGACCGTCATAACGACCGCCTGCACACACGGTGCCCTGGGAGCCGAGACTGTTAGTCACCCACTCGAAAACGGTACGGTTGTAGTAATCCAGACCACGCACCAGACGCTGGTTTACgg                                                                                                 <  1:504167/232‑1 (MQ=255)                                                                                                                                                                                                                                                   acaacaGGATCGGCTTTAAATTCCGGATTAACGGCCTGTACTAACAATACAAGACGTTCGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTCCACAAGACCGTCATAACGACCGCCTGCACACACGGTGCCCTGGGAGCCGAGACTGTTAGTCACCCACTCGAAAACGGTACGGTTGTAGTAATCCAGACCACGCACCAGACGCTGGTTTACgg                                                                                                 >  2:504167/1‑232 (MQ=255)                                                                                                                                                                                                                                                   acaacaGGATCGGCTTTAAATTCCGGATTAACGGCCTGTACTAACAATACAAGACGTTCGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTCCACAAGACCGTCATAACGACCGCCTGCACACACGGTGCCCTGGGAGCCGAGACTGTTAGTCACCCACTCGAAAACGGTACGGTTGTAGTAATCCAGACCACGCACCAGACGCTGGTTTACGGTGTAAGCGATCCCcg                                                                                  >  1:186102/1‑247 (MQ=255)                                                                                                                                                                                                                                                   acaacaGGATCGGCTTTAAATTCCGGATTAACGGCCTGTACTAACAATACAAGACGTTCGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTCCACAAGACCGTCATAACGACCGCCTGCACACACGGTGCCCTGGGAGCCGAGACTGTTAGTCACCCACTCGAAAACGGTACGGTTGTAGTAATCCAGACCACGCACCAGACGCTGGTTTACGGTGTAAGCGATCCCcg                                                                                  <  2:186102/247‑1 (MQ=255)                                                                                                                                                                                                                                                   acaacaGGATCGGCTTTAAATTCCGGATTAACGGCCTGTACTAACAATACAAGACGTTCGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTCCACAAGACCGTCATAACGACCGCCTGCACACACGGTGCCCTGGGAGCCGAGACTGTTAGTCACCCACTCGAAAACGGTACGGTTGTAGTAATCCAGACCACGCACCAGACGCTGGTTTACGGTGTAAGCGATCCCCGCGCTCTCCAGCAGTTTGCACAGACCGGCAAAATGCTCACgag                                        >  1:306771/1‑289 (MQ=255)                                                                                                                                                                                                                                                   acaacaGGATCGGCTTTAAATTCCGGATTAACGGCCTGTACTAACAATACAAGACGTTCGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTCCACAAGACCGTCATAACGACCGCCTGCACACACGGTGCCCTGGGAGCCGAGACTGTTAGTCACCCACTCGAAAACGGTACGGTTGTAGTAATCCAGACCACGCACCAGACGCTGGTTTACGGTGTAAGCGATCCCCGCGCTCTCCAGCAGTTTGCACAGACCGGCAAAATGCTCACgag                                        >  1:456469/1‑289 (MQ=255)                                                                                                                                                                                                                                                    caacaGGGTCGGCTTTAAATTCCGGATTAACGGCCTGTACTAACAATGCAAGACGTTCGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTCCACAAGACCGTCATAACGACCGCCTGCACACACGGTGCCCTGGGAGCCGAGACTGTTAGTCACCCACTCGAAAACGGTACGGTTGTAGTAATCCAGACCACGCACCAGACGCTGGTTTACGGTGTAAGCGATCCCCGCGCTCTCCAGCAGTTTGCACAGACCGGCAAAATGCTCACgag                                        <  2:76330/288‑1 (MQ=255)                                                                                                                                                                                                                                                               gCTTTAAATTCCGGATTAACGGCCTGTACTAACAATACAAGACGTTCGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTCCACAAGACCGTCATAACGACCGCCTGCACACACGGTGCCCTGGGAGCCGAGACTGTTAGTCACCCACTCGAAAACGGTACGGTTGTAGTAATCCAGACCACGCACCAGACGCTGGTTTACGGTGTAAGCGATCCCCGCGCTCTCCAGCAGTTTGCACAGACCGGCAAAATGCTCACGAGATTCCTCGTCCa                            >  1:486321/1‑289 (MQ=255)                                                                                                                                                                                                                                                                    aaaTTCCGGATTAACGGCCTGTACTAACAATACAAGACGTTCGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTCCACAAGACCGTCATAACGACCGCCTGCACACACGGTGCCCTGGGAGCCGAGACTGTTAGTCACCCACTCGAAAACGGTACGGTTGTAGTAATCCAGACCACGCACCAGACGCTGGTTTACGGTGTAAGCGATCCCCGCGCTCTCCAGCAGTTTGCACAGACCGGCAAAATGCTCACGAGATTCCTCGTCCAGATAg                       <  1:616286/289‑1 (MQ=255)                                                                                                                                                                                                                                                                     aaTTCCGGATTAACGGCCTGTACTAACAATACAAGACGTTCGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTCCACAAGACCGTCATAACGACCGCCTGCACACACGGTGCCCTGGGAGCCGAGACTGTTAGTCACCCACTCGAAAACGGTACGGTTGTAGTAATCCAGACCACGCACCAGACGCTGGTTTACGGTGTAAGCGATCCCCGCGCTCTCCAGCAGTTTGCACAGACCGGCAAAATGCTCACGAGATTCCTCGTCCAGATAGt                      <  1:104217/289‑1 (MQ=255)                                                                                                                                                                                                                                                                     aaTTCCGGATTAACGGCCTGTACTAACAAGACAAGACGTTCGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTCCACAAGACCGTCATAACGACCGCCTGCACACACGGTGCCCTGGGAGCCGAGACTGTTAGTCACCCACTCGAAAACGGTACGGTTGTAGTAATCCAGACCACGCACCAGACGCTGGTTTACGGTGTAAGCGATCCCCGCGCTCTCCAGCAGTTTGCACAGACCGGCAAAATGCTCACGAGATTCCTCGTCCAGATAGt                      <  1:384818/289‑1 (MQ=255)                                                                                                                                                                                                                                                                         ccGGATTAACGGCCTGTACTAACAATACAAGACGTTCGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTCCACAAGACCGTCATAACGACCGCCTGCACACACGGTGCCCTGGGAGCCGAGACTGTTAGTCACCCACTCGAAAACGGTACGGTTGTAGTAATCCAGACCACGCACCAGACGCTGGTTTACGGTGTAAGCGATCCCCGCGCTCTCCAGCAGTTTGCACAGACCGGCAAAATGCTCACGAGATTCCTCGTCCAGATAGTCAcc                  <  2:145847/289‑1 (MQ=255)                                                                                                                                                                                                                                                                                 aCGGCCTGTACTAACAATACAAGACGTTCGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTCCACAAGACCGTCATAACGACCGCCTGCACACACGGTGCCCTGGGAGCCGAGACTGTTAGTCACCCACTCGAAAACGGTACGGTTGTAGTAATCCAGACCACGCACCAGACGCTGGTTTACGGTGTAAGCGATCCCCGCGCTCTCCAGCAGTTTGCACAGACCGGCAAAATGCTCACGAGATTCCTCGTCCAGATAGTCACCTAATGCCg          <  1:268359/289‑1 (MQ=255)                                                                                                                                                                                                                                                                                     ccTGTACTAACAATACAAGACGTTCGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTCCACAAGACCGTCATAACGACCGCCTGCACACACGGTGCCCTGGGAGCCGAGACTGTTAGTCACCCACTCGAAAACGGTACGGTTGTAGTAATCCAGACCACGCACCAGACGCTGGTTTACGGTGTAAGCGATCCCCGCGCTCTCCAGCAGTTTGCACAGACCGGCAAAATGCTCACGAGATTCCTCGTCCAGATAGTCACCTAATGCCGGAg       <  2:109869/288‑1 (MQ=255)                                                                                                                                                                                                                                                                                      cTGTACTAACAATACAAGACGTTCGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTCCACAAGACCGTCATAACGACCGCCTGCACACACGGTGCCCTGGGAGCCGAGACTGTTAGTCACCCACTCGAAAACGGTACGGTTGTAGTAATCCAGACCACGCACCAGACGCTGGTTTACGGTGTAAGCGATCCCCGCGCTCTCCAGCAGTTTGCACAGACCGGCAAAATGCTCACGAGATTCCTCGTCCAGATAGTCACCTAATGCCGGAGcg     <  2:101984/289‑1 (MQ=255)                                                                                                                                                                                                                                                                                       tGTACTAACAATACAAGACGTTCGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTCCACAAGACCGTCATAACGACCGCCTGCACACACGGTGCCCTGGGAGCCGAGACTGTTAGTCACCCACTCGAAAACGGTACGGTTGTAGTAATCCAGACCACGCACCAGACGCTGGTTTACGGTGTAAGCGATCCCCGCGCTCTCCAGCAGTTTGCACAGACCGGCAAAATGCTCACGAGATTCCTCGTCCAGATAGTCACCTAATGCCGGAGcg     >  2:298444/1‑288 (MQ=255)                                                                                                                                                                                                                                                                                         tACTAACAATACAAGACGTTCGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTCCACAAGACCGTCATAACGACCGCCTGCACACACGGTGCCCTGGGAGCCGAGACTGTTAGTCACCCACTCGAAAACGGTACGGTTGTAGTAATCCAGACCACGCACCAGACGCTGGTTTACGGTGTAAGCGATCCCCGCGCTCTCCAGCAGTTTGCACAGACCGGCAAAATGCTCACGAGATTCCTCGTCCAGATAGTCACCTAATGCCGGAGcgtcg  >  2:223740/1‑289 (MQ=255)                                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                                          CTGCCGTTTGCTCACCAGAGCGCAAATCCTTCACTACTGCTGTGCCGTTAGCCACTTCAGACTCACCCAGCACCACAGCAACGCGGGCACCCCATTTATCAGCACGGGCAAACTGTTTCTTAAAGTTGCCGCCGCCGTGGTTGGTCATCAATTTCACGCCCGGTAATTCATCACGCAGACGCTCAGCTAATGCCATAGCCGCAGATTGTGTATCAGCACCTGAAGCCACCAGGTATATATCGACAACAGGATCGGCTTTAAATTCCGGATTAACGGCCTGTACTAACAATACAAGACGTTCGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTCCACAAGACCGTCATAACGACCGCCTGCACACACGGTGCCCTGGGAGCCGAGACTGTTAGTCACCCACTCGAAAACGGTACGGTTGTAGTAATCCAGACCACGCACCAGACGCTGGTTTACGGTGTAAGCGATCCCCGCGCTCTCCAGCAGTTTGCACAGACCGGCAAAATGCTCACGAGATTCCTCGTCCAGATAGTCACCTAATGCCGGAGCGTCG  >  NC_000913/2639348‑2639916

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑