Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	23,161	A→T	P257P (CCA→CCT)	ileS →	isoleucine‑‑tRNA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	23,161	0	A	T	95.6%	203.5 / NA	71	P257P (CCA→CCT)	ileS	isoleucine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base A (2/1);  major base T (34/31);  minor base G (0/3);  total (36/35)
Fisher's exact test for biased strand distribution p-value = 2.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.17e-01

CTGCACAAAGGCGCGAAGCCAGTTCACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTCCAGGCAGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTAATCTGGACCACCACGCCGTGGACTCTGCCTGCCAACCGCGCAATCTCTATTGCACCAGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAAAAGGTGCGGAGCTTGAGCTGCTGCGCTTT‑ACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTACCCTGGATGCCGGTACCGGTGCCGTTCACACCGCGCCTGGCCACG  >  NC_000913/22922‑23405                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                      cTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACt                                                                                                                                                                                                                                             >  1:482778/1‑250 (MQ=255)   tcacaAAGGCGCGAAGCCGGTACACTGGTGCGTTGTCTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTATCCGTCCATCGACGGTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGGTATCGGCCCAAGCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCGGCCTGCGAACCGCGCAATCTCGTTTGCACTGGATTGCGACTAt                                                                                                                                                                                                                                           <  2:503381/249‑1 (MQ=255)       aaaGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCt                                                                                                                                                                                                                                      >  2:115229/1‑251 (MQ=255)       aaaGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCt                                                                                                                                                                                                                                      >  1:379488/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)                                                                                                                                                                                                                          aCCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTAAAAGGTGCGGAGCTTGAGTTGCTGCGCTTT‑ACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTACCCTGGATGCCGGTACCGGTGCCGTTcac                  <  1:217186/251‑1 (MQ=255)                                                                                                                                                                                                                                cAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACAATTCTCAGCACGGTAAAAGGGGCGGGGCTTGAGTTGATGCGCTTT‑ACCCATCCGGTTATGGGCTTCGACGTTCCGGCAATCCTAGGCGATCACGTTTCCCTGGATGCCGGTACCGGTGCCGTTCACACcgcc            >  2:80285/1‑250 (MQ=255)                                                                                                                                                                                                                                    ctctATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGGAAAAGGTGCGGGGCTTGAGTTGCTGCGCTTT‑ACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTACCCTGGATGCCGGTACCGGTGCCGTTCACACCGCGCCTg        >  1:338702/1‑251 (MQ=255)                                                                                                                                                                                                                                     tctATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGGTTACACCATTCTCAGCACGGTAAAAGGTGCGGAGCTTGAGTTGCTGCGCTTT‑ACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTACCCTGGATGCCGGTACCGGTGCCGTTCACACCGCGCCTgg       >  1:419814/1‑251 (MQ=255)                                                                                                                                                                                                                                     tctATTGCACCTGATTTCGACTATGCGCTGGCGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGCCCGATTCCCCCATTATCAGCACGGTAAAAGGTGCGGGGCATGAGTTGCTGCGCTTT‑ACCCACCCGTTTATGGGCTTCGACCTTCACGCCATCCTCGGCGCTCACGCTACCCTGGATGCCGGCAACGGTGCCGTTCCCACCGCGCCt         >  2:423945/1‑249 (MQ=255)                                                                                                                                                                                                                                      ctATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTAAAAGGTGCGGAGCTTGAGTTGCTGCGCTTT‑ACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTACCCTGGATGCCGGTACCGGTGCCGTTCACACCGCGCCTGGc      >  1:392036/1‑251 (MQ=255)                                                                                                                                                                                                                                          tGCACCTGATTTCGACTATGCGCTGGTGCAGATCTCCGGTGAGGCCGTGTTTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGAACGATTACAACATTCTCAGCACGGTAAAAGGTGCGGCGCTTTCGTTGCTGCGCTTT‑ACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCCCGCCGATCACGTTACCCTGGTTGCCGATACCGGTGCCGGTCACACAGCGCCCGGccccg  >  2:475728/1‑248 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                      CTGCACAAAGGCGCGAAGCCAGTTCACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTCCAGGCAGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTAATCTGGACCACCACGCCGTGGACTCTGCCTGCCAACCGCGCAATCTCTATTGCACCAGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAAAAGGTGCGGAGCTTGAGCTGCTGCGCTTT‑ACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTACCCTGGATGCCGGTACCGGTGCCGTTCACACCGCGCCTGGCCACG  >  NC_000913/22922‑23405

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