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	94.0%	200.3 / 1.8	68	L459L (CTA→CTG)	dnaK	chaperone protein DnaK
Reads supporting (aligned to +/- strand):  ref base A (3/1);  new base G (36/27);  total (39/29)
Fisher's exact test for biased strand distribution p-value = 6.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.00e-01

GGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTAGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAATGGTACGCGACGCAGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCAT  >  NC_000913/13312‑13785                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                       ggTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCCAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACcgcg                                                                                                                                                                                                                                 >  1:217284/1‑251 (MQ=255)  gTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACACGGCACcgcgc                                                                                                                                                                                                                                >  1:218724/1‑251 (MQ=255)        tgaaaGACGTACTGCTGCTGGACGTTACCCAGCTGTCTCTGGGTATAGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAGAAACACCACTATACCGCCCAAGAACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGACTGCGGTGACCATCCATGTGCTGCAGGGTGACCGTGAACGAGCGGCTGATAACAAATCTCTGGGTCAGTGCAACCTGGATGGGATCACACCGGCAcagcacgacatt                                                                                                                                                                                                                          >  2:313170/3‑240 (MQ=255)             gACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATAAACCCGGCACCGCGCGGCATGCCGc                                                                                                                                                                                                                      >  2:418271/1‑250 (MQ=255)              aCGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGc                                                                                                                                                                                                                      >  1:310608/1‑249 (MQ=255)                        gctgGACGTTACCCCGCTGTTTCTGGGTACCGAAACCATGTGCGGTGGGAGGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGATGTTCTCTACCGCTGAAGACATCCAGTCTCCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCTGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAGCCCGGCCCCGCGCGGCATGCCGCAGATCGAAGTTa                                                                                                                                                                                                          <  1:379739/251‑1 (MQ=255)                        gctgGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCCGATCGAAGtt                                                                                                                                                                                                           >  2:417084/1‑250 (MQ=255)                          tgGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGCAGTTAcc                                                                                                                                                                                                        >  2:386819/1‑251 (MQ=255)                          tgGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAAACTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTAcc                                                                                                                                                                                                        >  1:180831/1‑251 (MQ=255)                           gggcgtgcccccGCGGTCTCGGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTAGCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCAGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACTTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGCAGTTAc                                                                                                                                                                                                         <  2:170276/241‑1 (MQ=255)                            gACGTTACCACGCTGACGCTGGGTATCGAAACCATGGGCGGTGTGATGAGGACGCTGATCGCGAAAAACACCACTAGCCCGACCAAGCCCAGCCAGGTGGTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCAGCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTAccct                                                                                                                                                                                                      <  2:138005/251‑3 (MQ=255)                             ccggtACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCCCTACCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATTGTATCAACCCGGCACCGCGCGGCAGGCCGCAGATCGAAGTTACCTTc                                                                                                                                                                                                     <  2:217284/247‑1 (MQ=255)                                   ccccGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTTGGTCAGTTCAACCTAGATGGTATCAACCCGGCACCGCGCGGCATGTCGCAGATCGAAGTTACCTTCGATATc                                                                                                                                                                                               >  1:123827/1‑251 (MQ=255)                                           tctcTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGGGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCACCCCGGCACCGCACGGCATGCCGCAGATCGGAGTTAACTTCGTTATCGATGCTg                                                                                                                                                                                        >  2:452338/1‑250 (MQ=255)                                                gggTATCGAAACCATTGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACg                                                                                                                                                                                     <  1:147082/248‑1 (MQ=255)                                                  gTATCGAACCCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCg                                                                                                                                                                                                                       >  2:44715/1‑212 (MQ=255)                                                  gTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCg                                                                                                                                                                                                                       <  1:44715/212‑1 (MQ=255)                                                   tATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGt                                                                                                                                                                                                                                                  >  1:148475/1‑184 (MQ=255)                                                    aTCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATcc                                                                                                                                                                               <  1:58139/250‑1 (MQ=255)                                                        aaaCCATGGGGGGTGTGATGACGACGCTGATCGCGAAAAACACCACTAACCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTGCCTTCGATATCGATGCTGACGGTATCCTGCAc                                                                                                                                                                          <  2:197380/251‑1 (MQ=255)                                                         aaCCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACg                                                                                                                                                                         <  1:458503/251‑1 (MQ=255)                                                            cATGGGCGGTGTGATAACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCAAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGttt                                                                                                                                                                      >  1:73855/1‑251 (MQ=255)                                                               gggCGGTTTGATGACGACTCTTATCGCGAAAAACTCCACCACCCCGCACAATCACATCCATGAGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTATCCATCCATGTGCTGCAGTGTGAACGTAAACGTGCGGCTGATAACAAATCGCTGGGTCAGGTCAACCTGGATGGTATCAGCCCGGCTCCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTAGCCGGCACGGTTc                                                                                                                                                                     <  2:26650/249‑1 (MQ=255)                                                                              gacgCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAAc                                                                                                                                                     >  1:259183/1‑250 (MQ=255)                                                                                cgCTGATCGCGAAAAACACCACGACCCCGCCCAAGCACAGCCAGGGGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAGCGTGCGGCGGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCGCCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATAGCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGc                                                                                                                                                  <  2:259183/251‑1 (MQ=255)                                                                                 gCTGATCGCGAAAAACACCACTACCCCGACCAAGCACAGCCAGGGGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCg                                                                                                                                                 <  2:135813/251‑1 (MQ=255)                                                                                       cgcgAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTCCCGCGAAAGCTAAAAAAAGCGGTaa                                                                                                                                             >  2:337130/1‑249 (MQ=255)                                                                                        gcgAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAaga                                                                                                                                          <  1:289113/251‑1 (MQ=255)                                                                                                   cacTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCagaagctc                                                                                                                                >  1:450353/1‑247 (MQ=255)                                                                                                           cGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCa                                                                                                                         >  1:265260/1‑249 (MQ=255)                                                                                                            gACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAgg                                                                                                                      >  1:323087/1‑251 (MQ=255)                                                                                                             aCCAAGCCCAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGAGCACCATCAAGGc                                                                                                                     <  2:310608/251‑1 (MQ=255)                                                                                                             aCCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGAAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGCTGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGc                                                                                                                     >  2:91026/1‑251 (MQ=255)                                                                                                                   cacaGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCttcttc                                                                                                               >  2:500689/1‑251 (MQ=255)                                                                                                                   cacaGCCAGGTGGTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCttcttc                                                                                                               <  1:190602/251‑1 (MQ=255)                                                                                                                      aGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTgg                                                                                                            <  1:419265/251‑1 (MQ=255)                                                                                                                          aGGTGTTCTCTACCGCTGAAGGCAACCAGTCTGCGGTAACCATCCATTGGCTGCAGGGTGAACGTAAACTTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGGATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCATGGCTTCTTCTGGTCTg                                                                                                        <  2:232103/251‑1 (MQ=255)                                                                                                                           ggTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTg                                                                                                        >  1:261855/1‑250 (MQ=255)                                                                                                                            gtgtTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGa                                                                                                       >  1:206228/1‑250 (MQ=255)                                                                                                                            gtgtTCGCGACCGCTGAAGACACCCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGaa                                                                                                      <  2:261855/251‑1 (MQ=255)                                                                                                                               gtctcTTCCGCTGAAGACAACCAGTCTGCGGTAACCATTCATGTGCCGCAGGTTGAACGTAAACGTGCGGCTGATAACAAATCTCTTGGTCAGTTCAACCTGGATGGTATCAACCCGGCGCCGCGCGGCATGCCGCAGATCGAAGTTTCCTTCGATATCGATTCTGACGGTATCCTGCACGTCTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCTTCAAGGCTTCTTCTGGTCTGAACg                                                                                                    <  2:67704/249‑1 (MQ=255)                                                                                                                                      ccGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGa                                                                                             >  1:348242/1‑250 (MQ=255)                                                                                                                                       cGCTGCAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGCAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGaa                                                                                            >  2:184450/1‑250 (MQ=255)                                                                                                                                       cGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGaa                                                                                            >  2:421576/1‑250 (MQ=255)                                                                                                                                          tGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATcc                                                                                        >  1:178066/1‑251 (MQ=255)                                                                                                                                           gAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGAAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGCACGAAGATGAAATca                                                                                        >  2:399043/1‑249 (MQ=255)                                                                                                                                            aaGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATcc                                                                                        >  1:438471/1‑249 (MQ=255)                                                                                                                                                 aaCCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAAt                                                                                 >  1:382906/1‑251 (MQ=255)                                                                                                                                                    cAGTCTGCGGGAACCAGCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAATGGt                                                                              <  2:505249/251‑1 (MQ=255)                                                                                                                                                     aGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCGCTTGGTCAGTTCAACCTAGATGGTATCAACCCGGCACCGCGCGGCATGTCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATTAAATCCAGAAAATGGTa                                                                             <  2:123827/251‑1 (MQ=255)                                                                                                                                                        atGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAATGGTAcg                                                                           <  1:315278/249‑1 (MQ=255)                                                                                                                                                                  catccatGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAATGGTACGCGACGCCGAAg                                                                <  2:264975/251‑1 (MQ=255)                                                                                                                                                                   atccatGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGaaaa                                                                                  >  1:439511/1‑232 (MQ=255)                                                                                                                                                                   atccagGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGaaaa                                                                                  <  2:439511/232‑1 (MQ=255)                                                                                                                                                                         gtgCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATAAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAATGGTACGCGACGCAGAAGCTAACGc                                                         >  1:26015/1‑251 (MQ=255)                                                                                                                                                                                aGGGTTAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCATCCTGGATGGTATCATCCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAATGGTACGCGACGCAGAAGCTAACGCCGAAGCt                                                  <  1:358614/251‑1 (MQ=255)                                                                                                                                                                                aGGGTGAACGTAAACGGGCGGCTGATAACAAATCTCTGGGTCAGGTCAACCTTGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATAACGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAATAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAAGGGTACGCGACGCAGAAGCTAACGCCGAAGCt                                                  <  2:218724/251‑1 (MQ=255)                                                                                                                                                                                  ggTGACCGTAAGCGTGCGGCTGCTAACAAATCTCTGGATCCGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGATGTTACCTTCGATATCGATGCTGACGGTATCCTGCATGTTTACGCGAAATATAAAAACAGAGGTACAGAGCAGACGCGCACCATAAAGGACGATTCTGGTCAGAACGTAGATGACATCCAGAAACTGGTACTCGGCGAAGAAGATACCGCCGAAGc                                                   >  2:3331/1‑248 (MQ=255)                                                                                                                                                                                       aCGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAATGGTACGCGACGCAGAAGCTAACGCCGAAGCTGACCgc                                            >  1:314554/1‑249 (MQ=255)                                                                                                                                                                                         gTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAATGGTACGCGACGCAGAAGCTAACGCCGAAGCTGACCGTAAg                                         <  1:500689/251‑1 (MQ=255)                                                                                                                                                                                           aaaCGTGCGGCTGATAACAAATCTCTGGGTCAGTCTAACCTAGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAATGGTACGCGACGCAGAAGCTAACGCCGAAGCTGACCGTAAGtt                                       >  1:186366/1‑251 (MQ=255)                                                                                                                                                                                              cGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCGGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTGCTTCTGGTCTGAACGAAGATGAAATCCAGAAAATGGTACGCGACGCAGACGCTAACGCCGAAGCTGACCGTAAGTTTGa                                    <  2:180831/251‑1 (MQ=255)                                                                                                                                                                                                               aTCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAATGGTACGCGACGCAGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAAGGCTGGTACAGACTc                   >  2:242887/1‑251 (MQ=255)                                                                                                                                                                                                                           ggtaaaCCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCGGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAATGGTACGCGACGCAGAAGCTAACGCCGAAGCTGACCGTAAGTGTGAAGAGCTGGTACAGACTCGCAACCAGGGCg       <  2:450353/247‑1 (MQ=255)                                                                                                                                                                                                                             tCAACCTAGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAATGGTACGCGACGCAGAAGCTAACGCCAAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGAc     >  1:4409/1‑251 (MQ=255)                                                                                                                                                                                                                              cAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGCAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAATGGTACGCGACGCAGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGAcc    >  2:63067/1‑251 (MQ=255)                                                                                                                                                                                                                               aaCCTGGATGGTATCAACCCGACACCGCGCGGCATGCCGCATATCGAAGTTACCATCGATATCGATGCTGACTGTATCCTGCACGTTTCCGCGAAAGCTAAAACCAGCGGTACAGAGCAGAAAATCACCATCAAGGCTTCTTCTGGTCTGAACGACAATGCAATCCAGAAACTAGTACGAGACGCAGAAGCTCACGACGAAGCTGCCCGAAAGTTTATAGATCTGGTAAAGACTCGCaccccactcaacc    >  2:367741/1‑238 (MQ=255)                                                                                                                                                                                                                                aCCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAATGGTACGCGACGCAGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCAt  >  2:395941/1‑251 (MQ=255)                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                       GGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTAGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAATGGTACGCGACGCAGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCAT  >  NC_000913/13312‑13785

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