Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	13,539	A→G	L459L (CTA→CTG)	dnaK →	chaperone protein DnaK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	13,539	0	A	G	98.7%	244.8 / ‑4.7	76	L459L (CTA→CTG)	dnaK	chaperone protein DnaK
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (34/40);  total (35/41)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.43e-01

GTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGT‑TCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTAGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAATGGTACGCGACGCAGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCG  >  NC_000913/13303‑13780                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                  gtgcggACTGGTGCCGAGAAAGACTTACTGTTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAATCCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGT‑TCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAAcc                                                                                                                                                                                                                                      <  2:323691/245‑1 (MQ=255)  ttCTGACTGGTGACGTGACAGACGTACTGCTGCTGGCCGTTACCCCACTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATAACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGAGT‑TCTCTACCGCTGAAGACAACCAGACTGCGGTAAACATCCATGTGCTGCAGGGTGAACTTAACCATGCGACTGATAACAAATCGCTGGGTCAGTTCAACCTGGATGGTATCAACCCg                                                                                                                                                                                                                                    >  2:372016/1‑251 (MQ=255)      gACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTAGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGT‑TCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCAGGCAc                                                                                                                                                                                                                                >  1:423935/1‑251 (MQ=255)       aCTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGT‑TCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCAcc                                                                                                                                                                                                                               <  1:7427/251‑1 (MQ=255)          ggggacttgaaaGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGTGATCGAAACCATGGGCGGTGTGATGACGACGCTGAATGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGT‑TCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACcg                                                                                                                                                                                                                              <  2:588647/240‑1 (MQ=255)            tGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGCCGCTGATCGCGAAAAACCCCACTATCCCGACCAAGCACAGCCAGGTGT‑TCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTCAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACcacgcg                                                                                                                                                                                                                          >  2:179848/1‑251 (MQ=255)            ggACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGT‑TCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACcgcgcg                                                                                                                                                                                                                          <  1:288669/250‑1 (MQ=255)               cgtgaacgACGTACAGCTGCTGAACGTTACCCCGCTGTCTCTGGGTCTCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATACCAAACCAGCACAGCCAGGTGT‑TCTCTCCCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGATGCAGGGTGAACGTAAACGTGCGGCTGATAACAACTATCTGTGTCAGTTCAACCAGGATGGTATCaaaccagacacccacaga                                                                                                                                                                                                                        >  2:101772/8‑234 (MQ=255)                 cgaaagacgttcggcgggtggccGTTACCCCGCTGTCGCTGGGTAGCGAAACCAGGGGCGGTGTGATGCCGGCGCTGATCGCGAAAAACCCCACGAACCCGACCAAGCACAGCCAGG‑GTGTCTCTAGCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCGGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCTCGGCAt                                                                                                                                                                                                                      <  2:405668/228‑1 (MQ=255)                  gaaaGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGT‑TCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATg                                                                                                                                                                                                                     <  1:25315/249‑1 (MQ=255)                  gaaaGACGGACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCAGGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTAACCCGACCAAGCACAGCCAGGTGT‑TCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGc                                                                                                                                                                                                                    <  1:352015/250‑1 (MQ=255)                     aGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCCGTGTGATGACGACGCTGATCGCGAAAAACACCCATATCCCAACCAAGCACAGCCAGGGGT‑TCTCTACCGATGAAGACAACCAGGCTGCGGTAACCATCCATGTGCTGCAGGGGGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCTACCCGGCACCGCGCGGCATGACGc                                                                                                                                                                                                                 >  1:317077/1‑251 (MQ=255)                                  ctgGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGT‑TCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAATTAc                                                                                                                                                                                                    >  2:368832/1‑251 (MQ=255)                                     gacggttccccGCTGTCTCTGGGTATCGAAACCATGGGCGGTGGGATGACGGAGCTGATCGCGAAAAACACCACTAGCCCGATCAAGCACAGCCAAGTGT‑TCTCTACCGCTGAAGACAACCAGTCTGCGGGAACCAGCCATGTTCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAGTCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCTCAGATCGAATTTACCtt                                                                                                                                                                                                 <  2:234517/244‑1 (MQ=255)                                       cGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGT‑TCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCg                                                                                                                                                                                               <  1:490870/251‑1 (MQ=255)                                              ccGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGT‑TCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCg                                                                                                                                                                                         >  1:189840/1‑250 (MQ=255)                                                   gtcgctGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCTAAAAACACCACCATCCCGACCAAGCACAGCCAGGTGT‑TCTCTACCGCTGAAGACAACCAGTCTGCGGTATCCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCTCCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTg                                                                                                                                                                                   <  2:411207/247‑1 (MQ=255)                                                   gTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGT‑TCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATgttg                                                                                                                                                                                   <  1:354885/251‑4 (MQ=255)                                                         gggTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGT‑TCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACg                                                                                                                                                                                <  1:358975/248‑1 (MQ=255)                                                         gggTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGT‑TCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTa                                                                                                                                                                             <  1:161960/251‑1 (MQ=255)                                                         gggTAGCGAAACCATGGGCGGTTTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAATCACAGCCAGGTGT‑TCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACTTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTGCAACCTGGATGTTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACg                                                                                                                                                                                <  1:83759/248‑1 (MQ=255)                                                          ggTATCGAAACCATGGGAGGTGTGATGACGACGCTGATCGCGAATAACACCACTATCCCGACCAAGCACAGACAGGAGT‑TCTCTACCGCTGAAGACAACCAGTCTGCGGAAACCATCCCTGTGCTGCAGGGGGAACGTAAACGTGCGGCTGATCACAAATCACTGGGTCAGTTCAACCTGGATGGTATCAAACCGACACCGAGCGGCATGCCGCAGATCGATGTTCCCTTCGATATCAGTGCTGACGGTAt                                                                                                                                                                            >  2:29959/1‑251 (MQ=255)                                                            tATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGT‑TCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATcc                                                                                                                                                                          >  1:492993/1‑251 (MQ=255)                                                             agcGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACCCCACTATCCCGACCAAGCACAGCCAGGTGT‑TCTCGACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCTGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCt                                                                                                                                                                         <  2:99807/249‑1 (MQ=255)                                                             aTCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGT‑TCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGCTGCTGACGGTATCCt                                                                                                                                                                         >  1:428420/1‑251 (MQ=255)                                                                 aaaCCATGTGCGGTGTGATGAGGACGCTGATCGCGAAAAACACCACTAGCCCGACCAAGCACAGCCAGGTGT‑TCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCAc                                                                                                                                                                     <  2:182797/251‑1 (MQ=255)                                                                 aaaCCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGT‑TCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCATGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCAc                                                                                                                                                                     <  1:280492/251‑1 (MQ=255)                                                                     cATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGT‑TCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGttt                                                                                                                                                                 <  1:462337/251‑1 (MQ=255)                                                                        gggCGGTGGGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGG‑TCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGGATCTTGCAGGTTTCcg                                                                                                                                                              <  2:261042/251‑1 (MQ=255)                                                                              tgtgATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGT‑TCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGaa                                                                                                                                                          >  1:379636/1‑249 (MQ=255)                                                                                tgatgaCGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGT‑TCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGaa                                                                                                                                                          >  1:580144/1‑247 (MQ=255)                                                                                  atgaCGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGT‑TCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATaa                                                                                                                                                    <  1:181717/251‑1 (MQ=255)                                                                                          gCTGATCGCGAAAAACACTACTATCCCGACCAAGCACAGCCAGGTGT‑TCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTa                                                                                                                                                                                                     >  2:118031/1‑194 (MQ=255)                                                                                          gCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGT‑TCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTa                                                                                                                                                                                                     <  1:118031/194‑1 (MQ=255)                                                                                           ctggtCGCGAAAACCACCACTATCCCGACCAAGCACAGCCAGGTGT‑TCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGATGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACATTCGATATCGATGCACACGAGATCCTGCGCGTTTCCGCGAAAAATAAAAGCAGCgg                                                                                                                                           >  2:469350/5‑251 (MQ=255)                                                                                              aTCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGT‑TCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTa                                                                                                                                         >  1:194659/1‑250 (MQ=255)                                                                                                  cgAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGT‑TCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGATACCTTCGATAACGATACTGACGGGAACCTGCACGTTTCCGCGAAAGATAAAAACAGCGGGAAag                                                                                                                                      >  2:370012/1‑249 (MQ=255)                                                                                                   gAAAAACACACCTATCCCGACCAAGCACAGCCAGGTGT‑TCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCAACCATGTTCTGCAGGGTGAACGTAAACGGGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCTCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGc                                                                                                                                   <  2:567137/251‑1 (MQ=255)                                                                                                   aaaaacCACCCCTAACCAGACCAAGCTCAGCCAGGTGT‑TCTCTACCGCTGAAGACCACCAGTCTGCTGGAACCATCCATGTGCTGCAGGGTGAACGTAAACGGGCGGCTGATAACAAATCTCTGGGGCAGTTCAACCTGGATGGTTTCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGTTGACGGTAGCCGGCACGTTTCCGAGAAAGATAAAATCAGTGGTAAAGAGc                                                                                                                                   <  2:85554/250‑1 (MQ=255)                                                                                                             acTATCCCGACCAAGCACAGCCAGGTGT‑TCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCAc                                                                                                                         <  1:572391/251‑1 (MQ=255)                                                                                                                   ccGACCAAGCACAGCCAGGTGT‑TCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCaa                                                                                                                   <  1:14237/251‑1 (MQ=255)                                                                                                                       ccAATCACAGCCAGGTTT‑TCGCTACCGCGGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGct                                                                                                               <  2:556482/251‑1 (MQ=255)                                                                                                                               aGCCAGGTGT‑TCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGGGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGAc                                                                                                                                                                                 <  1:441277/177‑1 (MQ=255)                                                                                                                                           ctctACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGa                                                                                              >  1:658765/1‑249 (MQ=255)                                                                                                                                             ctACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCa                                                                                                                          >  1:418981/1‑219 (MQ=255)                                                                                                                                             ctACCGCGAAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCTCCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCa                                                                                                                          <  2:418981/219‑1 (MQ=255)                                                                                                                                             ccaCCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAg                                                                                            <  2:313116/247‑1 (MQ=255)                                                                                                                                               aCCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATg                                                                                         <  1:160359/250‑1 (MQ=255)                                                                                                                                                ccGCTGACGTCAACCAGTCTGCTGTAACCATCCAGGTGCTGCAGGGTGAACGTAAACGCGCGGCCGATAACAAATCTCTGGGTCAGTTCAACCTTGATGGTATCAACCCGGCCCCGCGCGGCATGCCGCAGATCGAAGTGACCTTCGATATCGATGCTGACGGTATCCTGCTCGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGaa                                                                                       <  2:492993/251‑1 (MQ=255)                                                                                                                                                    tGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGAGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCCCCGAACGGCATGCAGCAGATCGAAGGTACCTTCGACATCGATGCTGCAGGTATCCTGCACGTTTCCGCGACAGATAAAAACAGAGATAAAGAGCAGAAGATCACCATCAAGGCTTATTCTGGGCTGAACGAAGATAAAATcc                                                                                   >  2:648142/1‑251 (MQ=255)                                                                                                                                                         aCAACCAGTCTGCGGTAACCATACATGTGCTGCAGGGTGAACGTAAACGTGCTACTGCTAACAACTCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCTCAGATCGACGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTACAGAGCAAAAGATCACCATCAAGGCTTCTTCTGTTCTGAACGCAGATGacctacagat                                                                               >  2:30160/1‑241 (MQ=255)                                                                                                                                                          cAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGaaa                                                                              >  1:553198/1‑250 (MQ=255)                                                                                                                                                              cAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAATGGt                                                                         >  2:142086/1‑251 (MQ=255)                                                                                                                                                               aGTCTGCGGTACCCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCgg                                                                                                                                                                                                                         >  2:30552/1‑106 (MQ=255)                                                                                                                                                               aGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCgg                                                                                                                                                                                                                         <  1:30552/106‑1 (MQ=255)                                                                                                                                                                         aaCCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAATGGTACGCGACGc                                                                >  1:84141/1‑249 (MQ=255)                                                                                                                                                                               ccatGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAATGGTACGCGACGCAGAAGCTa                                                        <  1:192175/251‑1 (MQ=255)                                                                                                                                                                                    tgctgcAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCt                                                                                                    >  1:296300/1‑202 (MQ=255)                                                                                                                                                                                    tgctgcAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTGCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTGCTGGTCt                                                                                                    <  2:296300/202‑1 (MQ=255)                                                                                                                                                                                     gctgcAGGGGGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTAAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAATGGTACGCGACGCAGAAGCTAACGcc                                                   <  2:477010/250‑1 (MQ=255)                                                                                                                                                                                       tgcAGGGAGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGACGATCACCATCAAGGCTTCTTCTGGACTGAACTAAGATGAAATCCAGAAAATGGTACGCGACGCCGCAGCTAACGccaa                                                 >  2:206264/1‑248 (MQ=255)                                                                                                                                                                                          aGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGcttc                                                                                                             >  1:471722/1‑187 (MQ=255)                                                                                                                                                                                          aGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCAACAAGGcttc                                                                                                             <  2:471722/187‑1 (MQ=255)                                                                                                                                                                                           gggTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGCAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGCTAAAAACAGCGGTAAAGAGCAGAAGATCAACATCAAGGCTTCTTATGGTCTGAACGAAGATGAAATCAAGAAAATGGTACGCGACGCAGAAGCTAACGCCGTAGCTg                                            >  2:513955/1‑251 (MQ=255)                                                                                                                                                                                             gTGACCGTGAACGTGCGTCTGATAACAAATCTCTGGGACAGTTCAACCTGGATGGTATCAACCCGGCACCGAGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCTCGTTTCCGCGAAAGATAAAAACAGCGGTAACGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGACCGAAGATGAGATCCAGAAAATGGGACGCGACGCAGACGCTAACGCCGACGCttac                                          >  2:190698/1‑248 (MQ=255)                                                                                                                                                                                             gTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGATGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTACAGAGCAGCAGCTCACCATCAAGGCTTCTTCTGGTCTGAACGAAGCTGCAATCCAGAAACTGGGACGCGCCGCAGAAGCTAACGCCGAAACTGAc                                          >  2:478475/1‑251 (MQ=255)                                                                                                                                                                                                aaCGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAATGGTACGCGACGCAGAAGCTAACGCCGAAGCTGACCGt                                       <  1:256667/251‑1 (MQ=255)                                                                                                                                                                                                       aCGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAATGGTACGCGACGCAGAAGCTAACGCCGAAGCTGACCGTAAGTTTg                                <  1:101772/251‑1 (MQ=255)                                                                                                                                                                                                       aCGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAATGGTACGCGACGCAGAAGCTAACGCCGAAGCTGACCGTAAGTTTg                                <  1:142086/251‑1 (MQ=255)                                                                                                                                                                                                        cGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAATTGTACGCGACGAAGAAGCTAACGCCGAAGCTGAcc                                         >  1:138400/1‑241 (MQ=255)                                                                                                                                                                                                        cGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGGATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAATGGTACGCGACGCAGAAGCTAACGCCGAAGATGACCGTACGTTTg                                >  2:254025/1‑250 (MQ=255)                                                                                                                                                                                                        cGTGCGGCTGATAACAAATCTCTGGGTCAGGTCAACCTGGATGGTATCCACCCGGCACCGCGCGGCATGCCGCATATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCGTCTGGTCTGAACGAAGATGAAATCCAGAAAATTGTACGCGACGCAGAAGCTAACGCCGAAGCTGAcc                                         <  2:138400/241‑1 (MQ=255)                                                                                                                                                                                                           gCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAATGGTACGCGACGCAGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAaga                            <  1:311732/251‑1 (MQ=255)                                                                                                                                                                                                              gCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGaa                                                                                             >  1:101328/1‑183 (MQ=255)                                                                                                                                                                                                                                 gTCAGTTCAACCTAGATGGTATCAAACCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCAACATAAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAACTGGTACTCGCCGCAGCAGCTAACGCCTAAGCTGACCGTAAGTTTGGACAGCTGGTAAGGACTCGCaccc        >  2:173139/1‑246 (MQ=255)                                                                                                                                                                                                                                      ttCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAATGGTACGCGACGCAGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCg  >  1:289884/1‑250 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                  GTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGT‑TCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTAGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAATGGTACGCGACGCAGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCG  >  NC_000913/13303‑13780

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