Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,948,824	C→T	A567T (GCA→ACA)	aspS ←	aspartate‑‑tRNA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,948,824	0	C	T	98.2%	169.9 / ‑4.2	57	A567T (GCA→ACA)	aspS	aspartate‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (25/30);  total (26/31)
Fisher's exact test for biased strand distribution p-value = 4.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.80e-01

GACCACTAAGATCGAAACGGGACGCTTATACACTTTATCCTTCACGCTGCCTCTTCGTTGACTGCCTTCGCTCATCCCATTCACATAGTTATCTATGCTCATGGGAGTTCACTCAGTTGCCGCCTCGATGCAACGCGAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGCTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAG‑CGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACACCGTCTGCTGCATATCACCATTATGGATA  >  NC_000913/1948588‑1949056                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                          gACCACTAAGATCGAAACGGGACGCTTATACACTTTATCCTTCACGCTGCCTCTTCGTTGACTGCCTTCGCTCATCCCATTCACATAGTTATCTATGCTCATGGGAGTTCACTCAGGTGCCGACTCGAGGCAACGCGAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCTGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGac                                                                                                                                                                                                                             <  1:443594/251‑1 (MQ=255)              gAAACGGGACGCTTATACACTTTATCCTTCACGCTGCCTCTTCGTTGACTGCCTTCGCTCATCCCATTCACATAGTTATCTATGCTCATGGGAGTTCACTCAGTTGCCGCCTCGATGCAACGCGAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGt                                                                                                                                                                                                                >  1:566529/1‑251 (MQ=255)                 aCGGGACGCTTATACACTTTATCCTTCACGCTGCCTCTTCGTTGACTGCCTTCGCTCATCCCATTCACATAGTTATCTATGCTCATGGGAGTTCACTCAGTTGCCGCCTCGATGCAACGCGAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCCGACACGCTGCCGCCGTGGt                                                                                                                                                                                                             >  2:562927/1‑251 (MQ=255)                  cGGGACGCTTATACACCTTATCCTTCCCGCTGCCTCTTCGTTGACTGCCTTCGCTCATCCCATTCACATATTTATCTATGCTCATGGGAGTTCACTCAGTTGCCGCCTCGATGCAACGCGAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGATGCCGCCGTGGtt                                                                                                                                                                                                            <  2:566529/251‑1 (MQ=255)                                 cTTTATCCTTCACGCTGCCTCTTCGTTGACTGCCTTCGCTCATCCCATTCACATAGTTATCTATGCTCATGGGAGTTCACTCAGTTGCCGCCTCGATGCAACGCGAATGTTTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCAGTGTTTCAGTCATCAGAAAAGCTGCCGCCGTGGTTTTCGGGAAGGCGAt                                                                                                                                                                                              >  2:405330/1‑250 (MQ=255)                                    tATCCTTCACGCTGCCTCTTCGTTGACTGCCTTCGCTCATCCCATTCACATAGTTATCTATGCTCATGGGAGTTCACTCAGTTGCCGCCTCGATGCAACGCGAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAAc                                                                                                                                                                                           <  1:412578/250‑1 (MQ=255)                                      tccctCACGCTGCCTCTTCGTTGACTGCCTTCGCTCATCCCATTCACATAGTTATCTATGCTCATGGGAGTTCACTCAGTTGCCGCCTCGATGCAACGCGAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGt                                                                                                                                                                                         <  2:406314/246‑1 (MQ=255)                                       ccTTCACGCTGCCTCTTCGTTGACTGCCTTCGCTCATCCCATTCACATAGTTATCTATGCTCATGGGAGTTCACTCAGTTGCCGCCTCGATGCAACGCGAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGa                                                                                                                                                                                                                              >  1:166998/1‑212 (MQ=255)                                       ccTTCACGCTGCCTCTTCGTTGACTGCCTTCGCTCATCCCATTCACATAGTTATCTATGCTCATGGGAGTTCACTCAGTTGCCGCCTCGATGCAACGCGAATGATTTCGTGTATTTGAGTCATATCAGTTAGTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGa                                                                                                                                                                                                                              <  2:166998/212‑1 (MQ=255)                                           cACGCTGCCTCTTCGTTGACTGCCTTCGCTCATCCCATTCACATAGTTATCTATGCTCATGGGAGTTCACTCAGTTGCCGCCTCGATGCAACGCGAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGTTGCTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACgg                                                                                                                                                                                    >  1:168714/1‑250 (MQ=255)                                              gCTGCCTCTTCGTTGACTGCCTTCGCTCATCCCAATCACATAGTTATCTATGCTCATGGGAGTTCACTCAGTTGCCGCCTCGATGCATCGCGAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGatat                                                                                                                                                                                <  2:215998/251‑1 (MQ=255)                                               cTGCCTATTAGTTGACTGCATTCGCTAATCCCATTCAAATAGATATCTATGCTCATGGGAGTTCACTCAGTTGCAGCCTCGATGCCACGCGAATGATTTTGAGTATTTGAGTCATCTCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTAGGGTTAGGAACGCTCGGTGTTTCAGTCATCAGAAAAGCTGCCGCCGTGGTTTTCGGGAGGGCGGAAACGTCACCGataat                                                                                                                                                                               >  2:524504/1‑249 (MQ=255)                                                cgCCTCTTCGTTGACTGCCTTCGCTCATCCCATTCACATAGTTATCTATGCTCATGGGAGTTCACTCAGTTGCCGCCTCGATGCAACGCGAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACATCCTGAATGCTCAGCTCAGCCAGGGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTg                                                                                                                                                                              <  1:44383/250‑1 (MQ=255)                                                     ctTCGTTGACTGCCTTCGCTCATCCTATTCACATAGTTATCTATGCTCATGGGAGTTCACTCAGTTGACGCCTCGATGCAACGCGAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTGCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGGGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCg                                                                                                                                                                           <  2:392901/249‑1 (MQ=255)                                                         gTTGACTGCCTTCGCTCATCCCATTCACATAGTTATCTATGCTCATGGGAGTTCACTCAGTTGCCGCCTCGATGCAACGCGAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCg                                                                                                                                                                     >  1:518158/1‑251 (MQ=255)                                                                    ccccccATCCCATTCAAATCGTTATCTATGATCATGGGAGTGCACTCAGTTGCCGCCTCGATGCAACGCGAATGATTTCGTGTATTTGAGGCATATCAGTTATTCTCAGCCTTCTTCACACCCTGAATGCTCAGCTCAGCCAGTGCCGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGATACGCGGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTcagcaga                                                                                                                                                            <  2:379104/244‑2 (MQ=255)                                                                         cATCCCATTCACATAGTTATCTATGCTCATGGGAGTTCACTCAGTTGCCGCCTCGATGCAACGCGAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTc                                                                                                                                                      >  1:268743/1‑250 (MQ=255)                                                                         cATCCCATTCACATAGTTATCTATGCTCAGGGGAGTTCACCCATTTGCCGCCTCGATGCAACGCGAATGATTTCGGGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTACACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCGCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCa                                                                                                                                                     <  2:348534/251‑1 (MQ=255)                                                                          aTCCCATTCACATAGTTATCTATGCTCATGGGAGTTCACTCAGTTGCCGCCTCGATGCAACGCGAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTc                                                                                                                                                      >  1:533267/1‑249 (MQ=255)                                                                            cccATTCACAAAGTTATCTATGCTCATGGGAGTGCACACAGTTGCCGCCTCGATGCAACGCGAATGATTTCGTGGATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTGCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGAc                                                                                                                                                  <  2:533267/251‑1 (MQ=255)                                                                             ccATTCACATAGTTATCTATGCCCATGGGATTTCACTCAGTTGCCGCCTCGATGCAACGCGAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACg                                                                                                                                                 <  2:95703/251‑1 (MQ=255)                                                                              cATTCACATAGTTATCTATGCTCATGGGAGTTCACTCAGTTGCCGCCTCGATGCAACGCGAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACgg                                                                                                                                                >  1:480616/1‑251 (MQ=255)                                                                                     aTAGTTATCTATGCTCATGGGAGTTCACTCAGTTGCCGCCTCGATGCAACGCGAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGAc                                                                                                                                         >  1:275121/1‑251 (MQ=255)                                                                                        gTTATCTATGCTCATGGGAGTTCACTCAGTTGCCGCCTCGATGCAACGCGAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGAGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCg                                                                                                                                       >  1:77495/1‑250 (MQ=255)                                                                                           aTCTATGCTCATGGGAGTTCACTCAGTTGCCGCCTCGATGCAACGCGAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAAt                                                                                                                                    <  1:542132/250‑1 (MQ=255)                                                                                                 gCTCATGGGAGTTCACTCAGTTGCCGCCTCGATGCAACGCGAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGa                                                                                                                              <  1:259367/250‑1 (MQ=255)                                                                                                           gTTCACTCAGTTGCCGCCTCGATGCAACGCGAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTgcgg                                                                                                                   >  1:96205/1‑251 (MQ=255)                                                                                                              cccTCAGTTGCCTCCCCGATGCAACGCGAATGATTTCGTGTATTGGAGTCATAGCAGTTAATCTCAGCCTTCTTCACAACCGGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTgcggcg                                                                                                                 <  2:275121/248‑1 (MQ=255)                                                                                                                  cAGTTGCCGTCTCGGTGCAACTCGAATGATTTCTTGTATTTTAGTCATATCAGTGATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCGCAGCCAGTGCAGTCGGGTTAGCAAAGCGCGGTGTTTCAGACATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTa                                                                                                            <  2:286301/251‑1 (MQ=255)                                                                                                                   aGTTGCCGCCTCGATGCAACGCGAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGAGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCCGACCTGCGTGCGGGGGAGTAc                                                                                                           >  1:67475/1‑251 (MQ=255)                                                                                                                       gccgccTCGATGCAACGCGAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACACCCTGAATGCTCAGCCCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGt                                                                                                        <  1:104586/250‑1 (MQ=255)                                                                                                                       gccgccTCGATGCAACGCGAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCTGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGAAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGt                                                                                                        >  1:101740/1‑250 (MQ=255)                                                                                                                           ccTCGATGCAACGCGAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTCTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTc                                                                                                   >  1:362505/1‑251 (MQ=255)                                                                                                                              cGATGCAACGCGAATGATTTCTTGTATTTGAGTCATATCAGTTATTCTCAGCCTTTTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACCCGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATAGTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCa                                                                                                  <  2:354966/249‑1 (MQ=255)                                                                                                                                attcAACGCGAATGATTTCGTGTATTTGATTCATATCAGTTATTCTCAGCCTCCTCCACAACCTGAATGCGCAGCTCAGCCAGTGCAGGCGGTTTAGCAAATCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCagag                                                                                               <  2:401745/247‑1 (MQ=255)                                                                                                                                   cAACGCGAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAAt                                                                                                                                    <  1:329945/210‑1 (MQ=255)                                                                                                                                        ggAATGATTTCTAGTATTTGATCCATTCCAGTTATGCTCAGCCGTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGGTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTTCCGCCGTGGTGTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTCCAGAG‑CGTCGAGc                                                                                      <  2:386818/249‑1 (MQ=255)                                                                                                                                         gAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAG‑CGTCGAGc                                                                                      >  1:5485/1‑249 (MQ=255)                                                                                                                                                   gtgtATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAG‑CGTCGAGCAGGAAGCCGAAt                                                                          <  1:489101/251‑1 (MQ=255)                                                                                                                                                       aTTTGAGTCATATCAGTTATTATCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGCTAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGGCAAGGCCTAAATCCAGACCTGCGTGCGGCGGTGTACCGTATTT‑ATAGACGTCGAGCCGGAAACCGAAttttt                                                                      >  2:298166/1‑249 (MQ=255)                                                                                                                                                            aGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAG‑CGTCGAGCAGGAAGCCGAATTTCTCgcgc                                                                 >  1:106273/1‑251 (MQ=255)                                                                                                                                                              tCATATAAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTAAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTAAGGCATCAGACACGCTGCCGCCGTGGATTTCGGGAAGGCGTAAAAGTCACGGATATTGACGGTGACGGGCAGCAGCATGGTCAGACGGCCAGGACCCGATGCCAGACCAGCGTGAGGAGGAGAACCGGATTTCAGAG‑CGTCGGGCAGGAACACGAGTTTCTCGCGCTg                                                               >  2:433427/1‑251 (MQ=255)                                                                                                                                                                    cAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAG‑CGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCct                                                          <  1:399414/250‑1 (MQ=255)                                                                                                                                                                      gTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAG‑CGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCtt                                                       >  1:271655/1‑251 (MQ=255)                                                                                                                                                                              cAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAG‑CGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGa                                                 <  1:291161/249‑1 (MQ=255)                                                                                                                                                                              cAGCCTTCTTCACAACCGGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCCGTCATCAGACTCGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAG‑CGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGAt                                                <  2:5113/250‑1 (MQ=255)                                                                                                                                                                               aGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAG‑CGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATAc                                              <  1:455583/251‑1 (MQ=255)                                                                                                                                                                                gCCTTCTACACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCGGCCGCCGTGGTTTTCGGGAAGGCGATAACGGCTCGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAG‑CGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATAcc                                             <  1:524504/251‑1 (MQ=255)                                                                                                                                                                                   ttcttcACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCCGACCTGCGGGCGGCGGAGTAACGTATTTCAGAG‑CGTCCAGCAGGAAGCCGAACTTCTCGCGCTCGTCCTCTTCCTTGATAccc                                            >  2:133113/1‑249 (MQ=255)                                                                                                                                                                                       tactaacctgaaggcgcgaccgagcCCGTGCAGTCGGGGTAGCAAATCTCGGGGAGTCATTCATCAGACACGCGGGCGGCGTGGTTTTCTGGAATGCGATAACGTCACGGATAGTGTCGGGGCCGGTCAGCAGCATGGTGAGACGGTCAAGACCGAACGCCAGTCGTGCGTGCGGCGGAGTACCGTATTTCAGAG‑CGTCGAGCAGGAAGCCGAATTTCTCGCTCTGTTCCTCTGCGTTGATACCCATAAt                                       <  2:315714/228‑1 (MQ=255)                                                                                                                                                                                         acaACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCCGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAG‑CGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATAcc                                    >  2:142189/1‑251 (MQ=255)                                                                                                                                                                                           aaCCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAG‑CGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCaa                                  <  1:133113/251‑1 (MQ=255)                                                                                                                                                                                                 aaTGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCACGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAG‑CGTCGAGCAGGAAGCCGAATTTATCGCGCTGTTCCTCTTCGTTGATAACCAGAATACCACACAccc                            >  1:18471/1‑250 (MQ=255)                                                                                                                                                                                                 aaTGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCACGGTTTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGGCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAG‑CGTCGAGCAGTAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACACCg                            <  1:68370/251‑1 (MQ=255)                                                                                                                                                                                                                           ggTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAG‑CGTCGCGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGCATACCAAACACCGTCTGCTGCATATCACCATTAttagt   >  2:57604/1‑246 (MQ=255)                                                                                                                                                                                                                           ggTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAG‑CGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACACCGTCTGCTGCATATCTCCATTATGGATa  <  1:60042/251‑1 (MQ=255)                                                                                                                                                                                                                                    aGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAG‑CGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTg                                                  >  1:347933/1‑194 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                          GACCACTAAGATCGAAACGGGACGCTTATACACTTTATCCTTCACGCTGCCTCTTCGTTGACTGCCTTCGCTCATCCCATTCACATAGTTATCTATGCTCATGGGAGTTCACTCAGTTGCCGCCTCGATGCAACGCGAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGCTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAG‑CGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACACCGTCTGCTGCATATCACCATTATGGATA  >  NC_000913/1948588‑1949056

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