Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	24,977	C→T	L863L (CTG→TTG)	ileS →	isoleucine‑‑tRNA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	24,977	0	C	T	98.2%	178.5 / ‑6.6	57	L863L (CTG→TTG)	ileS	isoleucine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (28/26);  total (29/28)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

CGAACGTGAAAAATACGTCTTCACCGGTGAGTGGTACGAAGGCCTGTTTGGCCTGGCAGACAGTGAAGCGATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAAC‑TGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGT‑TGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCTGCGGCCGCTGTGTCAGCAACGTCGCCGGTGACGGTGAAAAACGTAAGTTTGCCTGATGA  >  NC_000913/24730‑25210                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                           cGAACGTGAAAAATACGTCTTCACCGGTGAGTGGTACGAAGGCCTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTGGTCTTGt                                                                                                                                                                                                                                         <  2:157517/251‑1 (MQ=255)    aCGTGAAAAATACGTCTTCACCGGTGAGTGGTACGAAGGCCTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGT‑TGa                                                                                                                                                                                                                                     <  1:77281/251‑1 (MQ=255)      gTGAAAAATACGTCTTCACCGGTGAGTGGTACGAAGGCCTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGT‑TGAcc                                                                                                                                                                                                                                   <  1:82548/251‑1 (MQ=255)       tgcaaaaTACGTCTTCACCGGTGAGTGGTACGAAGGCCTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGCTTTGTCTTGT‑TGGCCt                                                                                                                                                                                                                                  >  2:3093/4‑251 (MQ=255)       tGAAAAATACGTCTTCACCGGTGAGTGGTACGAAGGCCTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAACGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGT‑TGACCt                                                                                                                                                                                                                                  >  2:17340/1‑251 (MQ=255)         aaaaaTACGTCTTCACCGGTGAGTGGTACGAAGGCCTGTTTGGTCTGGCCGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTACCCTTGTATGCAGAACCGGAAC‑TGGCGGCGCAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGT‑TTACCTcc                                                                                                                                                                                                                                >  1:130622/1‑251 (MQ=255)           aaaTACGTCTTCACCGGTGAGTGGTACGAAGGCCTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGT‑TGACCTCCgg                                                                                                                                                                                                                              >  2:92386/1‑251 (MQ=255)            actACGTCTTCACCGGTGAGTGGGACGAAGGCCTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGAGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACAGCGCTGGGCGATGAATTACGATTTGTC‑TGTGTGTCCTCCGgc                                                                                                                                                                                                                             >  2:6859/3‑251 (MQ=255)            aaTACGTCTTCACCGGTGAGTGGTACGAAGGACTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGAGATGAATTACGATTTGTCTTGT‑TGACCTCCGgc                                                                                                                                                                                                                             >  1:74642/1‑251 (MQ=255)            aaTACGTCTGCACCGGTGAGTGGTACGAAGGCCTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGT‑TGACCTCCGgc                                                                                                                                                                                                                             <  2:98410/251‑1 (MQ=255)              tACGTCTTCACCGGTGAGTGGTACGAAGGCCTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCGGGGCGATGAATTACGATTTGTCTTGT‑TGCCCTCCGgcgc                                                                                                                                                                                                                           >  1:166630/1‑251 (MQ=255)                    ttCACCGGTGAGTGGTACGAAGGCCTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGT‑TGACCTCCGGCGCTACCGt                                                                                                                                                                                                                     >  1:78463/1‑251 (MQ=255)                          ggTGAGTGGTACGAAGGCCTGTTTGGTCTGGCCGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGCAAGTGCGCGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGT‑TGACCTCCGGCGCTACCGTTGc                                                                                                                                                                                                                  >  2:194445/1‑248 (MQ=255)                                     cGAAGGCCTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGT‑TGGCCTCCGGCGCTACCGTTGGAGAGTATAacgacg                                                                                                                                                                                                    <  2:128181/251‑1 (MQ=255)                                        aGGCCTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGT‑TGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCAc                                                                                                                                                                                                 <  2:218829/251‑1 (MQ=255)                                           ccTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGT‑TGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACctg                                                                                                                                                                                              <  1:120392/251‑1 (MQ=255)                                             tGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGACGTGAACAAAGTCATTGAGCACGCGCGTGCCGACAAGACAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGCTGAATTACGATTTGTCTTGT‑TGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACctgct                                                                                                                                                                                            >  1:35913/1‑251 (MQ=255)                                                       ggCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCAGTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGT‑TGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTcag                                                                                                                                                                                   <  2:195948/250‑1 (MQ=255)                                                            aCAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCAAAGTGGACAAAGTAATTGAGCCAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGT‑TGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGc                                                                                                                                                                             >  1:39156/1‑251 (MQ=255)                                                                      aaTGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGT‑TGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTc                                                                                                                                                                    >  1:5508/2‑250 (MQ=255)                                                                          aaCGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGT‑TGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAgg                                                                                                                                                               <  2:68074/251‑1 (MQ=255)                                                                            cgatgctgtCTGGGACGAGCTGTGGAAAGTGCGTGGCGAAGTGAACAAAGTCAGTGAGCAAGGGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGT‑TGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGc                                                                                                                                                             <  2:173655/243‑1 (MQ=255)                                                                                  gTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGT‑TGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGaa                                                                                                                                                         >  1:215880/1‑249 (MQ=255)                                                                                       gggACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGT‑TGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTcgcg                                                                                                                                                  >  1:81710/1‑251 (MQ=255)                                                                                              gCTGTTGAAATTGCGGGGCGAAGTAAACAAAGTCATTGAGCAAGCGCGTGCCTGCAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGT‑TGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTa                                                                                                                                           <  2:80817/251‑1 (MQ=255)                                                                                                 gTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGT‑TGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTaa                                                                                                                                          >  2:56455/1‑249 (MQ=255)                                                                                                          ggggggggaaggtgaCAAAGTCAATGGGCAAGCGGGTGCCCGCAAGAAAGGGGGGGGGTCGGTGGGAGCGGCAGTAACCTTGTATGCAGGACCGGAACCTG‑CTGCGAAACTGGCCGCGCTGGGCCATGAATTACGGTTTGTCCTTT‑TGACCCCCGGCGCTTCCGTTGCAGACCATAACGACGCACCTTCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGCCGCGTTGAGTAAAGCCGAAGGTg                                                                                                                               <  2:168402/237‑1 (MQ=255)                                                                                                           cggggCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGT‑TGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTga                                                                                                                              <  2:173175/248‑1 (MQ=255)                                                                                                            gTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGT‑TGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACctgct                                                                                                                                                                                            >  1:59170/1‑188 (MQ=255)                                                                                                            gTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGT‑TGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACctgct                                                                                                                                                                                            <  2:59170/188‑1 (MQ=255)                                                                                                            gTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGT‑TGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTg                                                                                                                               >  1:243560/1‑249 (MQ=255)                                                                                                                 gAAGAGAACAAAGACATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAAACGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGCTGCATTACGATTTGTCTTGT‑TGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAAGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTg                                                                                                                        >  2:138854/1‑251 (MQ=255)                                                                                                                  aaGTGAACAAAGTCAGTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGT‑TGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGc                                                                                                                       <  2:118185/251‑1 (MQ=255)                                                                                                                   aGTGAACAAAGGCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGT‑TGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCGGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGcc                                                                                                                      <  2:123695/251‑1 (MQ=255)                                                                                                                    gTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGT‑TGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGccc                                                                                                                     >  1:48375/1‑251 (MQ=255)                                                                                                                         aaaaGGCATTGAGCAAGCGCGGGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTGGTATGCAGAACCGGAAC‑GGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCGTGT‑GGACCTCCGGCGAGAACGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACGCAAAGGGCTGAAAGGCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACgct                                                                                                                <  2:35913/250‑1 (MQ=255)                                                                                                                              tCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGT‑TGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTgg                                                                                                           >  2:70067/1‑251 (MQ=255)                                                                                                                                       caGCGCATGCCGAATAGCAAGTGGGTGGCTGGCTGGACGCGGCAGTCACCTTGTATGCAGAACCGGCAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATCACGATTTGTCTTGT‑TGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCCGAGCGACGTACTCCCACCGCTGAAAGTCGCGTTGAGTACCGCCGAAGTGGGGACGTGCCCAAGCTGCTGGCGCTacccc                                                                                                  >  2:44687/2‑248 (MQ=255)                                                                                                                                               gCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGT‑TGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATg                                                                                           <  1:92386/250‑1 (MQ=255)                                                                                                                                                                        tGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGT‑TGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGa                                                                  >  2:76020/1‑250 (MQ=255)                                                                                                                                                                         ggAAGCGGCAGTAACCTTGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGT‑TGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGaa                                                                 >  1:40707/1‑250 (MQ=255)                                                                                                                                                                            aGCGGCAGTAACCTTGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGT‑TGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCt                                                             <  1:235724/251‑1 (MQ=255)                                                                                                                                                                               ggCAGTAACCTTGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGT‑TGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCTGCg                                                          <  2:5508/251‑1 (MQ=255)                                                                                                                                                                                   gTAACCTTGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGT‑TGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCTGCGGCCg                                                      <  1:46458/251‑1 (MQ=255)                                                                                                                                                                                        cTTGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGT‑TGACCTCCGGCGCTACCGTTACAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGAGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACACAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCTGCGGCCGCtgtg                                                 >  1:163003/1‑251 (MQ=255)                                                                                                                                                                                          tGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGT‑TGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGCTGTCGGCAAGGTGGCGGAACACGCAGAAATCTGCGGCCGCTGTGTc                                               >  2:205037/1‑251 (MQ=255)                                                                                                                                                                                                cAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGT‑TGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCTGCGGCCGCTGTGTCAGCAAc                                         <  1:102205/251‑1 (MQ=255)                                                                                                                                                                                                 aGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGT‑TGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCTGCGGCCGCTGTGTCAGCAACg                                        <  1:56455/251‑1 (MQ=255)                                                                                                                                                                                                      cGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGT‑TGACCTCCGGCGCTACCGTTGAAGACTATAACGACGCACCTGCTGCTGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGACGCGTTGAGTAAAGCCGGAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGGAGGCAAGGAGGCGGAACACGAAGAAATCTGCGGCCGCTGTGGCAGCAACGTCgac                                   >  1:267743/1‑249 (MQ=255)                                                                                                                                                                                                              ggcGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGT‑TGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCTGCGTCCGCTGTGTCAGCAACGTCGCcggtgacg                            <  1:90616/249‑1 (MQ=255)                                                                                                                                                                                                                    gAAACTGACCGCGCTGGGCGATGACTTACGATTTGTGTTGT‑TGACCTCCGGCGCTACCGGTGCAGACTATAACGACGCACCTGCTGCTGCTCAGCAGCGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGCAGGTGAGCAGTGCCCAAGCTGCTGGCACTACACCCAGGATGTCGGCACGGTGGCGGAACACGCAGAAATCTGCGGCCGCTGTGTCCGCCACGTCGCCGGTGACGGTGaa                       >  2:128517/1‑250 (MQ=255)                                                                                                                                                                                                                     aaaCTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGT‑TGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCTGCGGCCGCTGTGTCAGCAACGTCGCCGGTGACGGTGaaaa                     <  2:48375/251‑1 (MQ=255)                                                                                                                                                                                                                          gACCGCGCTGGGCGATGAATTACGATTTGTCTTGT‑TGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCTGCGGCCGCTGTGTCAGCAACGTCGCCGGTGACGGTGAAAAACGt                 >  1:99010/1‑250 (MQ=255)                                                                                                                                                                                                                          gACCGCGCTGGGCGATGAATTACGATTTGTCTTGT‑TGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCTGCGGCCGCTGTGTCAGCAACGTCGCCGGTGACGGTGAAAAACGTa                <  1:71667/251‑1 (MQ=255)                                                                                                                                                                                                                          gACCGCGCTGGGCGATGAATGACGATTGGTCTTGT‑TGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGTCTGAAAGTCGCGTTGAGTAAAGCCGAATGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCTGCGGCCGCTGTGTCAGCAACGTCGCCGGTGACGGTGAAAAACGt                 <  2:130622/250‑1 (MQ=255)                                                                                                                                                                                                                                    ggCGATGAATTACGATTTGTCTTGT‑TGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGACGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCTGCGGCCGCTGTGTCAGCAACGTCGCCGGTGACGGTGAAAAACGTAAGTTTGCCtg      >  2:126333/1‑251 (MQ=255)                                                                                                                                                                                                                                        aTGAATTACGATTTGTCTTGT‑TTCCCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAGGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCTGCGGCCGCTTTGTCAGCAACGTCGCCGGTGACGGTGAAAAACGTAAGTTTGCCtgatga  <  2:12022/251‑1 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                           CGAACGTGAAAAATACGTCTTCACCGGTGAGTGGTACGAAGGCCTGTTTGGCCTGGCAGACAGTGAAGCGATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAAC‑TGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGT‑TGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCTGCGGCCGCTGTGTCAGCAACGTCGCCGGTGACGGTGAAAAACGTAAGTTTGCCTGATGA  >  NC_000913/24730‑25210

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