Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	22,945	T→A	V185V (GTT→GTA)	ileS →	isoleucine‑‑tRNA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	22,945	0	T	A	95.5%	187.4 / 1.6	67	V185V (GTT→GTA)	ileS	isoleucine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base A (34/29);  total (36/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.84e-02

AGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCAGTTCACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTC‑CAGGCAGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTAATCTGGACCACCACGCCGTGGACTCTGCCTGCCAACCGCGCAATCTCTATTGCACCAGATTTCGACTATGCGCTG  >  NC_000913/22698‑23179                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                            aGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCAGTTCAc                                                                                                                                                                                                                                          >  1:237385/1‑251 (MQ=255)  gCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACt                                                                                                                                                                                                                                         >  1:32924/1‑251 (MQ=255)      aaaggcgcgCAAGAATACGGTAAGCCGTGGGTGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGCCCCAGGTTGACGGGCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCAGCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACTAAGGCGCGAAGCCGGTACACTGGTg                                                                                                                                                                                                                                     <  2:350929/243‑1 (MQ=255)       aaGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGc                                                                                                                                                                                                                                    >  1:500961/1‑251 (MQ=255)           cGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTg                                                                                                                                                                                                                                <  1:385030/251‑1 (MQ=255)               cAAGAATACGGTAAGCCGGGTGAGAAATTAACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGAGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCAGCGCGCTGGGCAAAATAATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGAGCGTGGACTg                                                                                                                                                                                                                            >  1:359829/1‑251 (MQ=255)                        gttaagccgtgtgcgaattcccccgccgccgAGTTCCGCGCCCAGTGCCTCGAATACGCGGCGACCCAGGTTACCGGTCATCGCAAAGACTTTATCCGTTTGGTCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATTGACGTCAAATCTGAAGCCAACCCCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCTCCTGCACAAAGGCTCGAAGCCGGTACTCTTGTGCGTTTACTGCCCTTCTgc                                                                                                                                                                                                                   <  2:299647/229‑1 (MQ=255)                         gTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACACAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTgcg                                                                                                                                                                                                                  >  1:220551/1‑251 (MQ=255)                            aGCCGGGTGAGAAATTAACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGCCGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACACGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTg                                                                                                                                                                                                               >  1:503381/1‑251 (MQ=255)                              ccGGGTGAGAAATTCACCGCCGCCGAGTTCCGCTCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCAGCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGgc                                                                                                                                                                                                             <  2:484699/251‑1 (MQ=255)                              ccGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGgc                                                                                                                                                                                                             >  1:45915/1‑251 (MQ=255)                              ccGGGGGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTAACGGTCAACGCAAAGACTTTATCCGTCTGGAAGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCAGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACAATGGTGCGTTGGCTGCCGATCTGCGCTggg                                                                                                                                                                                                             >  2:11877/1‑250 (MQ=255)                                gggTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGACAACGCAGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGCCTGCCGTTCTGCGCTggggg                                                                                                                                                                                                           >  2:328995/1‑248 (MQ=255)                                  gTGAGAACGTCACCGCCGTCGAGTTCCGCGCCAAGTGCCGCGAATACGCCGCGACCCAGGTTTACGGGCAACGCAAAGACTGTATCCGTCTGGGTGTGCTGGGCGACTGGTCGCACCAGGACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGAAAAATCATCGGCAACGCTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGgcggaa                                                                                                                                                                                                         <  2:32110/251‑1 (MQ=255)                                    aataaaTTCACCGCCGCCGAGTTCCGCGCCAATTGCCGCGAATACGCGGCTACCCAGGTTGTCGTTCAACGCAATGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCAGCGGCAACGGTCACCTGCACAAATTCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGgcggaagc                                                                                                                                                                                                       <  1:484975/248‑1 (MQ=255)                                               cgccgccgAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGt                                                                                                                                                                                            >  1:412669/1‑251 (MQ=255)                                                   gccgAGTTCCGCGCCACGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCTCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTa                                                                                                                                                                                        <  1:451241/251‑1 (MQ=255)                                                   gacgAGTTCCCCGCCAAGTGCCGCGAATACGGGGCGCCCCAGGTTGACGGTCAACGCAAAGACTTTCTCCGTCGGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCAAGGACTTCAAAACTGAAGCCAACATCACCCGCGCGCTGGGCAAAACCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTa                                                                                                                                                                                        <  1:156259/249‑1 (MQ=255)                                                     cgAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACg                                                                                                                                                                                      >  1:397639/1‑251 (MQ=255)                                                            cgcgCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGGTCTGCGCTGGCGGAAGCGGAAGTTGAGGATTACGACAAAAc                                                                                                                                                                               <  2:72077/251‑1 (MQ=255)                                                                   aGTGCCGCGAAGACGCGGCGGCCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGTCTTCAAAACTGAAGCCAACAGCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAACCTTCTCCg                                                                                                                                                                        <  2:160819/251‑1 (MQ=255)                                                                     tGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTc                                                                                                                                                                      >  1:351610/1‑251 (MQ=255)                                                                           gAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAACCTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCCCAAAGGCGCGCAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGGATTACGCCAAAACTTCTCCATcccacgc                                                                                                                                                                >  2:278557/1‑246 (MQ=255)                                                                             aTACGCGTCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGGCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACg                                                                                                                                                              <  1:18707/251‑1 (MQ=255)                                                                               ccgcgGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGTAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGtt                                                                                                                                                            <  1:427118/250‑1 (MQ=255)                                                                                  cggcgACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGGGTATTACGACAAAACTTCTCCGTCCATCGCCGGTGCt                                                                                                                                                         >  2:362273/1‑251 (MQ=255)                                                                                   ggcgACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCtt                                                                                                                                                        <  1:484801/251‑1 (MQ=255)                                                                                        cccgggTTGACGGTCAACCCAAAGCCTTTTGCCGTGTGGGCGTGCTGGCCGCATGGTCGCACCCGTACCTGACCATGGGCGTCAACACTGAGTCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCATCGGTCACCTGCACAAAGGCGCGAACCGAGTGCACTGGTGCGTTGACGGCCGCTCTGCGCTGGCGGAAGCGGAAGTGGAGTATTACGTCAAAACGTCACCGCCCATCGACGtttctcgccag                                                                                                                                                   <  2:234874/247‑10 (MQ=255)                                                                                         ccAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAgg                                                                                                                                                 >  2:265486/1‑251 (MQ=255)                                                                                              ttGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGGGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAggtggtc                                                                                                                                            >  2:356569/1‑246 (MQ=255)                                                                                                taaggagatggCAAAGACTTTATCCGTCTGGGCGTGCTGTGCGACTGGTCGCACCCGTACCTGACCATGGTCTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGATCACCTGCACAAAGGCGCGAAGCCGGTACCCTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTCCGCCAAAACTTCTCCGTCCATCGACATTGCttttcaagctgtcgg                                                                                                                                           >  2:267733/11‑239 (MQ=255)                                                                                                    gTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGAGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATc                                                                                                                                        >  2:448187/1‑249 (MQ=255)                                                                                                                   ttATCCGGCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGc                                                                                                                        <  2:356107/250‑1 (MQ=255)                                                                                                                     aTCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCCTGGACTTCAAAACTGAAGCCACCATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACCCTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAACGGACGTTGAGTCTTCCGACAAAACTTCTCCGTCCATCTACGTTGCTTTT‑CAGTCGGTCGATCAGGATGCACTGAAAGCaa                                                                                                                      >  2:243926/1‑250 (MQ=255)                                                                                                                     aTCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCaa                                                                                                                      >  1:425430/1‑250 (MQ=255)                                                                                                                         gTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTg                                                                                                                                                                                                                            >  1:357607/1‑145 (MQ=255)                                                                                                                         gTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTg                                                                                                                                                                                                                            <  2:357607/145‑1 (MQ=255)                                                                                                                             gggCGTGCTGGGCGACTGGTTGCACTCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTg                                                                                                                <  2:292874/248‑1 (MQ=255)                                                                                                                                 gTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGACGCCAACATCATCCGCGCGCTGGGCAAACTCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGAGCGTTGACTGCCGTTCTGCGCTGGCGGCAACGGACGTTGAGTCTTACGACAACACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACCGAAAGCAACAtttacctta                                                                                                           >  2:320997/1‑243 (MQ=255)                                                                                                                                 gTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTa                                                                                                           >  1:229796/1‑249 (MQ=255)                                                                                                                                     tGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCGAAACTGAAGCCACCATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCAGTTCACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGCAGTTGAGTATTACGACAACACTTCTCCGTCCATCGACGTTGCTTTC‑CAGGCCGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAAc                                                                                                     >  2:255105/1‑251 (MQ=255)                                                                                                                                     tGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAAc                                                                                                     >  1:120092/1‑251 (MQ=255)                                                                                                                                      gggcgtatGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGGGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCaa                                                                                                      <  2:348562/242‑1 (MQ=255)                                                                                                                                      gggCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAAc                                                                                                     <  1:303790/250‑1 (MQ=255)                                                                                                                                                    aCCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGATGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTTGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAATCAAAATTTGCCGTAAGCAACGTTAACGGCCCAAtc                                                                                      <  1:369689/251‑1 (MQ=255)                                                                                                                                                      ccGTACCTGACCATGGACTTCAAAACGGAAGCCAACTTCAGCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCGGCACAAAGGCGCGAAGCCGGTACACGGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTGATCGGCCCAAtcgc                                                                                    <  2:231977/251‑3 (MQ=255)                                                                                                                                                        gTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGt                                                                                                   >  1:55836/1‑234 (MQ=255)                                                                                                                                                        gTACCTGACCATGGACTTCAAAACGGAAGCCAACATCATCCGCGCGCTGGGCAAAATCAGCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGt                                                                                                   <  2:55836/234‑1 (MQ=255)                                                                                                                                                                 cATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGAAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCt                                                                         >  1:281383/1‑251 (MQ=255)                                                                                                                                                                    ggaattCAAAACTGAAGCCAACATCAACCGCGCGCTGGCCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGACGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGAGCTGGa                                                                      <  2:269656/247‑1 (MQ=255)                                                                                                                                                                     gACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGTCAAAAGCATCGGCAACGGTCACCTGCCCAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATGTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGac                                                                     <  2:500961/251‑1 (MQ=255)                                                                                                                                                                        ttCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGaccac                                                                  >  2:217186/1‑251 (MQ=255)                                                                                                                                                                         tCAAAACGGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACGTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGaccac                                                                  <  2:393554/250‑1 (MQ=255)                                                                                                                                                                                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)                                                                                                                                                                                                                 catcGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTggggat                                                                           >  2:29508/1‑197 (MQ=255)                                                                                                                                                                                                                 catcGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGtgat                                                                           <  1:29508/201‑4 (MQ=255)                                                                                                                                                                                                                                 cTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACt          >  1:482778/1‑250 (MQ=255)                                                                                                                                                                                                                                   tcacaAAGGCGCGAAGCCGGTACACTGGTGCGTTGTCTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTATCCGTCCATCGACGGTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGGTATCGGCCCAAGCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCGGCCTGCGAACCGCGCAATCTCGTTTGCACTGGATTGCGACTAt        <  2:503381/249‑1 (MQ=255)                                                                                                                                                                                                                                       aaaGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCt   >  2:115229/1‑251 (MQ=255)                                                                                                                                                                                                                                       aaaGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCt   >  1:379488/1‑251 (MQ=255)                                                                                                                                                                                                                                        aaGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCCCCTGATTTCGCCTATGCGCTg  >  2:125876/1‑251 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                            AGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCAGTTCACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTC‑CAGGCAGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTAATCTGGACCACCACGCCGTGGACTCTGCCTGCCAACCGCGCAATCTCTATTGCACCAGATTTCGACTATGCGCTG  >  NC_000913/22698‑23179

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