Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	23,137	C→G	A249A (GCC→GCG)	ileS →	isoleucine‑‑tRNA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	23,137	0	C	G	93.9%	192.9 / 5.5	67	A249A (GCC→GCG)	ileS	isoleucine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base C (2/2);  new base G (28/34);  total (31/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 = 1.53e-01

CGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCAGTTCACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTCCAGGCAGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTAATCTGGACCACCACGCCGTGGACTCTGCCTGCCAACCGCGCAATCTCTATTGCACCAGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAAAAGGTGCGGAGCTTGAGCTGCTGCGCTTT‑ACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTACCCTGGATGCCGGTACC  >  NC_000913/22893‑23377                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                  cgcTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTAAGCCGTGGACTCTGCCTGCGAACcgc                                                                                                                                                                                                                                             >  1:496499/1‑251 (MQ=255)     gggCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCAGTTCACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTCCAGGCAGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCC‑AATCTCGCTGGTAATCTGGACCACCACGCCGTGGACTCTGCCTGCCAACCGCGCAAt                                                                                                                                                                                                                                        <  2:237385/251‑1 (MQ=255)      ggCAAAATCATCGGCATCTGGCCCCTGCACAAAGGCGCGATGCCGGTACACTGGTGCGGTGACTGCCGTTCTGCGCTGGCGGAAGCTGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGGCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACTTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAATCGCGCAAt                                                                                                                                                                                                                                        <  2:400878/251‑1 (MQ=255)                              cTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACt                                                                                                                                                                                                                 >  1:482778/1‑250 (MQ=255)                                tcacaAAGGCGCGAAGCCGGTACACTGGTGCGTTGTCTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTATCCGTCCATCGACGGTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGGTATCGGCCCAAGCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCGGCCTGCGAACCGCGCAATCTCGTTTGCACTGGATTGCGACTAt                                                                                                                                                                                                               <  2:503381/249‑1 (MQ=255)                                    aaaGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCt                                                                                                                                                                                                          >  1:379488/1‑251 (MQ=255)                                    aaaGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCt                                                                                                                                                                                                          >  2:115229/1‑251 (MQ=255)                                     aaGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCCCCTGATTTCGCCTATGCGCTg                                                                                                                                                                                                         >  2:125876/1‑251 (MQ=255)                                                 cggtacACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACAGAAAGCAAAATTTGCCGTCAGAAAAGTTAACGGCCCAATCTCGCTGGTGATCTGTACCACTACGCCGTGGCATCTGCCTGCAAACCGCGCAATCTCTATTCCACCTGATTTCGACAATGCTCTGGAGCAGTTCGAc                                                                                                                                                                                             >  2:274235/6‑251 (MQ=255)                                                  ggtacACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACg                                                                                                                                                                                            <  1:265486/247‑1 (MQ=255)                                                  ggtacACTGGTGCGTTGACTGCCGTTCGGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACg                                                                                                                                                                                            <  2:149249/247‑1 (MQ=255)                                                       aCTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGCCTATGCGATGGTGCAGATCGACGGTc                                                                                                                                                                                         >  1:368787/1‑249 (MQ=255)                                                                   aCTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTg                                                                                                                                                                           <  2:86781/251‑1 (MQ=255)                                                                   aCTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGAGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTg                                                                                                                                                                           >  2:306621/1‑251 (MQ=255)                                                                           tCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAg                                                                                                                                                                    <  1:251561/250‑1 (MQ=255)                                                                              gcgcGGGCTGAAGCTGAAGTTTAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTTAAAGCAAAATTTGCCTTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTTTATTGCATCTGATTTCGACTATGCGCTGGTGCAGATCGTCGTTCAGGTCGTGATTCTGGCGAAAGATCt                                                                                                                                                                <  1:360000/251‑1 (MQ=255)                                                                                        aagcggaagTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAACATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGc                                                                                                                                                      >  1:230879/1‑251 (MQ=255)                                                                                             gaagTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGCTTTCGACTATGCGCTGGTGCAGATCGCCGGTCAGGCCGTGATTCTGGCAGAAGATCTGTTTGaagacgcaat                                                                                                                                                 >  2:450379/1‑243 (MQ=255)                                                                                               agTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGc                                                                                                                                               >  1:117650/1‑251 (MQ=255)                                                                                                    attactACGACAAATCTTTTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGTATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTTGTGATCTGGACCACTACGCCGTGGACTCTGCCTTCGAACCGCGCAATCGCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGt                                                                                                                                          <  1:54215/246‑1 (MQ=255)                                                                                                    aggaTTACGACAAAACTTCTCCGTCCAACGTCGTTGCTGTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCCCTGGTGAGCTGGACCACTACGCCGTGGACTCTGCCGGCGAACCGCGCAATCTCTATTGCACCTGATTTCGAATATGCGCTGGTGCAGATCGACGGTCAGGTCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTACTGCAGCGt                                                                                                                                          <  2:351610/248‑1 (MQ=255)                                                                                                    aGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGt                                                                                                                                          <  2:425430/251‑1 (MQ=255)                                                                                                                aaaCTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTggcc                                                                                                                              >  1:228682/1‑248 (MQ=255)                                                                                                                        ccggccATCGACGTTTCTTTTCAGGCGGACGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGACCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCGGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGGATCTGCGTGACCGATTac                                                                                                                        <  2:225698/245‑1 (MQ=255)                                                                                                                         cGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACAcc                                                                                                                     <  1:78738/251‑1 (MQ=255)                                                                                                                          gtccctctaCGTTGCTTTTCAGGCGGTCGACCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTCTTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCTAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCa                                                                                                                    <  2:281383/243‑1 (MQ=255)                                                                                                                            ccattggcGTTGCTTTTCAGGCGGTCGATTAGGATGCACAGAAAGCAAAATTTTCCGTAAGCAACGTTAACGGCCCAATCTATCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCTCCGGAGTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCAt                                                                                                                   <  2:278228/243‑1 (MQ=255)                                                                                                                             cATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTTACCGATTACAcccttt                                                                                                                 >  2:6420/1‑247 (MQ=255)                                                                                                                              agcGACGGTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATtct                                                                                                                <  1:267733/249‑1 (MQ=255)                                                                                                                                cGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATtct                                                                                                                <  1:356569/249‑1 (MQ=255)                                                                                                                                cGACGTTGCGTTTCAGTCGGTCGATCAGGATGCACTGAAAGCAAAAGTTTCCGTAAGAAACGTTAACGGCCCAATCTCGCGGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCAACTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACCCCATtctca                                                                                                              <  2:460447/251‑2 (MQ=255)                                                                                                                                  aCGTTGCTTGTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCGGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATtctcagc                                                                                                            <  1:74963/251‑4 (MQ=255)                                                                                                                                       gCTTTTCAGGCGGTCGATCAGGATGCCCTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACgg                                                                                                        <  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)                                                                                                                                                                                                             cccgccgctcttccgatcttacGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTAAAAGGTGCGGAGCTTGAGTTGCTGCGCTTT‑ACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTc                                <  2:33242/231‑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)                                                                                                                                                                                                                                tacGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTAAAAGGTGCGGAGCTTGAGTTGCTGCGCTTT‑ACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTc                                >  1:33242/2‑232 (MQ=255)                                                                                                                                                                                                                                      gTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCCCCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAACGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCCGCACGGTAAAAGGTGCGGAGCTTGAGGTGCTGCGCTTT‑ACCCATCCGTTTATGGGCTTCGACGTTCCGGAAATCCCGGGCGATCACGTTAACCTGGATGCCg       >  2:24608/1‑251 (MQ=255)                                                                                                                                                                                                                                       tGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTAAAAAGTGCGGAGCTTGAGTTGCTGCGCTTT‑ACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTACCCTGGATGCCgg      >  1:359047/1‑251 (MQ=255)                                                                                                                                                                                                                                        ggACTCTGCCTGCGAACCGCGCAATCTCTAGTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTAAAAGGTGCGGAGCTTGAGTTGCTGCGCTTT‑ACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTACCCTGGATGCCGGt     <  2:379488/251‑1 (MQ=255)                                                                                                                                                                                                                                           ccctgcctgcGAACCGCGCAATCGCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGAGCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTAAAAGGTGCGGAGCTTGAGTTGCTGCGCTTT‑ACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTACCCTGGATGCCGGTAcc  <  2:412669/249‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                  CGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCAGTTCACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTCCAGGCAGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTAATCTGGACCACCACGCCGTGGACTCTGCCTGCCAACCGCGCAATCTCTATTGCACCAGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAAAAGGTGCGGAGCTTGAGCTGCTGCGCTTT‑ACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTACCCTGGATGCCGGTACC  >  NC_000913/22893‑23377

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: ‑