Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	23,116	C→T	T242T (ACC→ACT)	ileS →	isoleucine‑‑tRNA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	23,116	0	C	T	93.7%	176.0 / 0.2	66	T242T (ACC→ACT)	ileS	isoleucine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base C (2/2);  new base T (28/31);  total (30/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.71e-01

TGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCAGTTCACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTC‑CAGGCAGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTAATCTGGACCACCACGCCGTGGACTCTGCCTGCCAACCGCGCAATCTCTATTGCACCAGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAAAAGGTGCGGAGCTTGAGCTGCTGCGCTTT‑ACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTAC  >  NC_000913/22873‑23361                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       tGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGt                                                                                                                                                                                                                                                 >  1:8002/1‑251 (MQ=255)  gAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTg                                                                                                                                                                                                                                                <  1:65444/251‑1 (MQ=255)    aGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTgg                                                                                                                                                                                                                                               >  2:251561/1‑250 (MQ=255)                 cgcgcgCTGGGCAAACTCACCGGTAACGGTCACCTGTGCAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAACCTGCTCCGTCCATCGACGTTG‑TTTCTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCTCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTctgcctgcgaa                                                                                                                                                                                                                                 <  1:20772/251‑4 (MQ=255)                     cgcTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTAAGCCGTGGACTCTGCCTGCGAACcgc                                                                                                                                                                                                                             >  1:496499/1‑251 (MQ=255)                         gggCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCAGTTCACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTC‑CAGGCAGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCC‑AATCTCGCTGGTAATCTGGACCACCACGCCGTGGACTCTGCCTGCCAACCGCGCAAt                                                                                                                                                                                                                        <  2:237385/251‑1 (MQ=255)                          ggCAAAATCATCGGCATCTGGCCCCTGCACAAAGGCGCGATGCCGGTACACTGGTGCGGTGACTGCCGTTCTGCGCTGGCGGAAGCTGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGGCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACTTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAATCGCGCAAt                                                                                                                                                                                                                        <  2:400878/251‑1 (MQ=255)                                                  cTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACt                                                                                                                                                                                                 >  1:482778/1‑250 (MQ=255)                                                    tcacaAAGGCGCGAAGCCGGTACACTGGTGCGTTGTCTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTATCCGTCCATCGACGGTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGGTATCGGCCCAAGCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCGGCCTGCGAACCGCGCAATCTCGTTTGCACTGGATTGCGACTAt                                                                                                                                                                                               <  2:503381/249‑1 (MQ=255)                                                        aaaGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCt                                                                                                                                                                                          >  1:379488/1‑251 (MQ=255)                                                        aaaGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCt                                                                                                                                                                                          >  2:115229/1‑251 (MQ=255)                                                         aaGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCCCCTGATTTCGCCTATGCGCTg                                                                                                                                                                                         >  2:125876/1‑251 (MQ=255)                                                                     cggtacACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACAGAAAGCAAAATTTGCCGTCAGAAAAGTTAACGGCCCAATCTCGCTGGTGATCTGTACCACTACGCCGTGGCATCTGCCTGCAAACCGCGCAATCTCTATTCCACCTGATTTCGACAATGCTCTGGAGCAGTTCGAc                                                                                                                                                                             >  2:274235/6‑251 (MQ=255)                                                                      ggtacACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACg                                                                                                                                                                            <  1:265486/247‑1 (MQ=255)                                                                      ggtacACTGGTGCGTTGACTGCCGTTCGGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACg                                                                                                                                                                            <  2:149249/247‑1 (MQ=255)                                                                           aCTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGCCTATGCGATGGTGCAGATCGACGGTc                                                                                                                                                                         >  1:368787/1‑249 (MQ=255)                                                                                       aCTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTg                                                                                                                                                           <  2:86781/251‑1 (MQ=255)                                                                                       aCTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGAGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTg                                                                                                                                                           >  2:306621/1‑251 (MQ=255)                                                                                               tCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAg                                                                                                                                                    <  1:251561/250‑1 (MQ=255)                                                                                                  gcgcGGGCTGAAGCTGAAGTTTAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTTAAAGCAAAATTTGCCTTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTTTATTGCATCTGATTTCGACTATGCGCTGGTGCAGATCGTCGTTCAGGTCGTGATTCTGGCGAAAGATCt                                                                                                                                                <  1:360000/251‑1 (MQ=255)                                                                                                            aagcggaagTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAACATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGc                                                                                                                                      >  1:230879/1‑251 (MQ=255)                                                                                                                 gaagTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGCTTTCGACTATGCGCTGGTGCAGATCGCCGGTCAGGCCGTGATTCTGGCAGAAGATCTGTTTGaagacgcaat                                                                                                                                 >  2:450379/1‑243 (MQ=255)                                                                                                                   agTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGc                                                                                                                               >  1:117650/1‑251 (MQ=255)                                                                                                                        attactACGACAAATCTTTTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGTATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTTGTGATCTGGACCACTACGCCGTGGACTCTGCCTTCGAACCGCGCAATCGCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGt                                                                                                                          <  1:54215/246‑1 (MQ=255)                                                                                                                        aggaTTACGACAAAACTTCTCCGTCCAACGTCGTTGCTGTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCCCTGGTGAGCTGGACCACTACGCCGTGGACTCTGCCGGCGAACCGCGCAATCTCTATTGCACCTGATTTCGAATATGCGCTGGTGCAGATCGACGGTCAGGTCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTACTGCAGCGt                                                                                                                          <  2:351610/248‑1 (MQ=255)                                                                                                                        aGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGt                                                                                                                          <  2:425430/251‑1 (MQ=255)                                                                                                                                    aaaCTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTggcc                                                                                                              >  1:228682/1‑248 (MQ=255)                                                                                                                                            ccggccATCGACGTTTCTTTT‑CAGGCGGACGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGACCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCGGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGGATCTGCGTGACCGATTac                                                                                                        <  2:225698/245‑1 (MQ=255)                                                                                                                                             cGTCCATCGACGTTGCTTTT‑CAGGCGGTCGGTCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTctg                                                                                                                                                                                                                                         >  1:277680/1‑119 (MQ=255)                                                                                                                                             cGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACAcc                                                                                                     <  1:78738/251‑1 (MQ=255)                                                                                                                                             cGTCCAGCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTctg                                                                                                                                                                                                                                         <  2:277680/119‑1 (MQ=255)                                                                                                                                              gtccctctaCGTTGCTTTT‑CAGGCGGTCGACCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTCTTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCTAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCa                                                                                                    <  2:281383/243‑1 (MQ=255)                                                                                                                                                ccattggcGTTGCTTTT‑CAGGCGGTCGATTAGGATGCACAGAAAGCAAAATTTTCCGTAAGCAACGTTAACGGCCCAATCTATCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCTCCGGAGTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCAt                                                                                                   <  2:278228/243‑1 (MQ=255)                                                                                                                                                 cATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTTACCGATTACAcccttt                                                                                                 >  2:6420/1‑247 (MQ=255)                                                                                                                                                  agcGACGGTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATtct                                                                                                <  1:267733/249‑1 (MQ=255)                                                                                                                                                    cGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATtct                                                                                                <  1:356569/249‑1 (MQ=255)                                                                                                                                                    cGACGTTGCGTTT‑CAGTCGGTCGATCAGGATGCACTGAAAGCAAAAGTTTCCGTAAGAAACGTTAACGGCCCAATCTCGCGGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCAACTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACCCCATtctca                                                                                              <  2:460447/251‑2 (MQ=255)                                                                                                                                                      aCGTTGCTTGT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCGGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATtctcagc                                                                                            <  1:74963/251‑4 (MQ=255)                                                                                                                                                           gCTTTT‑CAGGCGGTCGATCAGGATGCCCTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACgg                                                                                        <  2:32924/250‑1 (MQ=255)                                                                                                                                                                ttcAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTaaa                                                                                    >  1:445739/3‑250 (MQ=255)                                                                                                                                                                  cAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTAAAAg                                                                                  <  1:243926/250‑1 (MQ=255)                                                                                                                                                                     gcggTCGATCAGGAGGCACTAAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACGACGTCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTGCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGTCGAAATATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTAAAAGGTGc                                                                              <  2:45915/248‑1 (MQ=255)                                                                                                                                                                           gATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTAAAAGGTGCGGAGCt                                                                        >  1:338687/1‑251 (MQ=255)                                                                                                                                                                              cAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTAAAAGGTGCGGAGCTTg                                                                      <  1:362273/250‑1 (MQ=255)                                                                                                                                                                                ggATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCCTTCTCAGCACGGGAAAAGGTGCGGAGCGTGAgt                                                                   >  1:265810/1‑250 (MQ=255)                                                                                                                                                                                    gCACTGAAAGCAAAATTTTCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCGGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGTCGGGCAGGCCGTGATTCTGTCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTAGCGGCGTGACCGATTACACCATTCTCAGCACGGTAAAAGGTTCGGAGCTTGAgttgct                                                               <  2:399939/251‑1 (MQ=255)                                                                                                                                                                                        tGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTAGCTGGTGATCTGGACCACTACGCAGTGGACTATGCCTGCGAACCGCGCACTCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTTGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATAGGCGTGCCCGACTACACCATTCTCAGCACGGCAAAAggggaggggcttggactgtcgcgc                                                           >  2:215072/1‑229 (MQ=255)                                                                                                                                                                                          aaaGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGCTTACACCATTCTCAGCACGGTAAAAGGTg                                                                               >  2:423572/1‑229 (MQ=255)                                                                                                                                                                                          aaaGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTAAAAGGTg                                                                               <  1:423572/229‑1 (MQ=255)                                                                                                                                                                                          aaaGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACGACGCCGTCGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTAAAATGTGCGGAGCTTGAGTTGCTGCGCtt                                                         <  2:123954/251‑1 (MQ=255)                                                                                                                                                                                             gCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTAATCTGGACCACCACGCCGTGGACTCTGCCTGCCAACCGCGCAGTCTCTATTGCACCAGATTTCGACTATGCGTTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAAAAGGTGCGGAGCTTGAGCTGCTGCGCTTT‑Ac                                                     >  1:47707/1‑251 (MQ=255)                                                                                                                                                                                              cAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTCCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGAGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGAACGGTAAAAGGTGCGGAGCTTGAGTTGCTGCGCTTT‑Acc                                                    >  1:225499/1‑251 (MQ=255)                                                                                                                                                                                              cAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTAATCTGGACCACCACGCCGTGGACTCTGCCTGTCAACCGCGCAATCTCTATTGCACCAGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAAAAGGTGCGGAGCTTGAGCTGCTGCGCTTT‑Acc                                                    <  1:255105/251‑1 (MQ=255)                                                                                                                                                                                               aaaaTTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGGTTACCCAATTTTCAGCACGGTAAAAGGGGCGGAGCTTGAGTTGCTGCGCTTT‑Accc                                                   >  2:101602/1‑251 (MQ=255)                                                                                                                                                                                                 aaTTTGCCGTAAGCAACGTTAACGGCCCAATCGCGCTGGTTATCTGTACCACTACGCTGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGTATGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACTGGAAAAGGTGCGGAGCTTGAGTTGCTGCGCTTT‑ACCCAt                                                 <  1:53061/251‑1 (MQ=255)                                                                                                                                                                                                        cGTAAGCAACGTTAACGGCCCAATCTCGCTGGCGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGCCTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTAAACCATTCTCAGCACGGTAAAAGGTGCGGGGCTTGAGTTGCTGCGCTTT‑ACCCATCCGttt                                           >  2:25972/1‑250 (MQ=255)                                                                                                                                                                                                                 cGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTAAAAGGTGCGGGGCTTGAGTTGCTGCGCTTT‑ACCCATCCGTTTATGGGCTTCg                                 >  1:423073/1‑251 (MQ=255)                                                                                                                                                                                                                    tAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACAATTCTCAGCACGGTAAAAGGTGCGGAGCtt                                                                       >  1:346162/1‑211 (MQ=255)                                                                                                                                                                                                                       cGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTAAAAGGTGCGGAGCTTGAGTTGCTGCGCTTT‑ACCCATCCGTTTATGGGCTTCGACGTTc                           <  1:278557/251‑1 (MQ=255)                                                                                                                                                                                                                       cGGCCCAATCTCGCTGGGGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTAAAAGGTGCGGAGCTTGAGTTGCTGCGCTTT‑ACCCATCCGTTTATGGGCTTCGACGTTc                           <  2:220551/251‑1 (MQ=255)                                                                                                                                                                                                                         gCCCAATCTCGCTGGTGAGCTGGACCCCTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTAAAAGGTGCGGAGCTTGAGTTGCTGCGCTTT‑ACCCATCCGTTTATGGGCTTCGACGTTCCg                         <  2:397928/251‑1 (MQ=255)                                                                                                                                                                                                                                  cGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGTCGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTAAAAGGTGCGGAGCTTGAGTTGCTGCGCTTT‑ACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTc                <  1:11877/251‑1 (MQ=255)                                                                                                                                                                                                                                   gCTGGTGATCTGGACCATGACGCCGTGTACGCTGCCTGCGAACCGCGCAATCTCTAGTGCTCCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTTTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTAAAAGGTGCGGAGCTTGAGTTGCTGCGCTTT‑ACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCg               <  1:328995/251‑1 (MQ=255)                                                                                                                                                                                                                                     tGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGCTTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGAATACACCATTCTCAGCACGGTACAAGGTGCGGAGCTTGAGTTGATGCGCTTTCACCC‑TCCGTTTATGGGCTTCGGCGATCCGGCAATCCTCGGc             >  2:395486/1‑251 (MQ=255)                                                                                                                                                                                                                                          aTCTGGACCACCACGCCGTGGACTCTGCCTGCCAACCGCGCAATCTCTATTGCACCAGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGAACGATTACCCCATTCTCGGCACGGTAAAAGGGGCGGAGCTTGAGCTGCTGCGCTTT‑ACCCATCCGTTTATGGGCTTCGGCGTTCGGGCAATCCCCGGCgctcc        >  1:111058/1‑247 (MQ=255)                                                                                                                                                                                                                                             tGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAACGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTACGCACGGTTAAAGGGGCGGAGGCTGGGTTGCTGCGCTTT‑ACCCATCCGTTTATGGGATTCGGCGTTCCGGCAATCCTCGGCGATCACGt     >  2:367104/1‑251 (MQ=255)                                                                                                                                                                                                                                                accactACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACCCCATTCTCAGCACGGAAAAAGGTGCGGAGCTTGAGTTGCTGCGCTTT‑ACCCATACGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCAATCACGTTAc  >  1:474640/1‑251 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       TGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCAGTTCACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTC‑CAGGCAGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTAATCTGGACCACCACGCCGTGGACTCTGCCTGCCAACCGCGCAATCTCTATTGCACCAGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAAAAGGTGCGGAGCTTGAGCTGCTGCGCTTT‑ACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTAC  >  NC_000913/22873‑23361

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