Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	13,320	A→G	V386V (GTA→GTG)	dnaK →	chaperone protein DnaK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	13,320	0	A	G	98.0%	150.0 / ‑4.4	50	V386V (GTA→GTG)	dnaK	chaperone protein DnaK
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (24/25);  total (24/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.52e-01

AACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCAT‑CCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTAGATGGTATCAACCCGGCACCGCG  >  NC_000913/13073‑13562                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    aaCTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTgaaa                                                                                                                                                                                                                                                  <  1:293155/251‑1 (MQ=255)         aGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATCTCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCTGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCGGGGTGGTGTTCTGCCTGGTGACGTGAAAGACGTAct                                                                                                                                                                                                                                          >  1:465807/1‑251 (MQ=255)         aGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAg                                                                                                                                                                                                                                                 <  1:7913/244‑1 (MQ=255)         aGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAg                                                                                                                                                                                                                                                 >  2:7913/1‑244 (MQ=255)                 tGAAGATCTGGTAAAACGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTgg                                                                                                                                                                                                                                  >  1:116303/1‑251 (MQ=255)                          ggTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTAcc                                                                                                                                                                                                                          <  1:1661/250‑1 (MQ=255)                                ccGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGt                                                                                                                                                                                                                   <  2:175538/251‑1 (MQ=255)                                 cGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGtc                                                                                                                                                                                                                  >  1:30164/1‑251 (MQ=255)                                   ttCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTAct                                                                                                                                                                                                                                          <  1:57797/225‑1 (MQ=255)                                   ttCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGGAct                                                                                                                                                                                                                                          >  2:57797/1‑225 (MQ=255)                                               gCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGaaa                                                                                                                                                                                                    <  1:201929/251‑1 (MQ=255)                                                             tGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGtgtg                                                                                                                                                                                      >  1:268782/1‑251 (MQ=255)                                                                             ggccGTAGCTGATATCGACGACGTTAGCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGAt                                                                                                                                                                        <  2:96435/247‑1 (MQ=255)                                                                             gTCCGTATCTGATATCGACGACGTTATCCTCGTTGGAGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGCTGCCGACGCTTATcg                                                                                                                                                                      >  1:278345/1‑251 (MQ=255)                                                                                          atCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCCGAAGAAAGTTGCTGAGTTATTTGGTAAAGAGCCGCGTAAAGACGTTAACACGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTATGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCGGGGTATCGAAACCATGGGCGGAGTGATGACGACGCTGATCGCGAAAAAcaccac                                                                                                                                                         >  1:157982/1‑251 (MQ=255)                                                                                          atCGACGACGTTATCCGCGTGGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTTGTGTCCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAAcaccac                                                                                                                                                         <  1:21313/251‑1 (MQ=255)                                                                                            cgacgaCGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAAGCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCAGGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACt                                                                                                                                                        <  2:54778/250‑1 (MQ=255)                                                                                                         cTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGCAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGATCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACACAGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGAACAAGc                                                                                                                                           >  1:191258/1‑250 (MQ=255)                                                                                                            gTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAg                                                                                                                                       <  1:222587/251‑1 (MQ=255)                                                                                                                        aCTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTtctc                                                                                                                           <  1:306880/251‑1 (MQ=255)                                                                                                                         cTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGAGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTctct                                                                                                                          >  1:460115/1‑251 (MQ=255)                                                                                                                                   aaTGGTTCAGAAGAAAGTTGCTGATTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTACCCCGACCAAGCACAGCCAGGTGTTctct                                                                                                                          <  2:484471/241‑1 (MQ=255)                                                                                                                                   aaTGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCAAGGTGTTctct                                                                                                                          >  1:484471/1‑241 (MQ=255)                                                                                                                                         tCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGCAAAACACCACTATCCCGACCAAGCACAGCACGGTGTTCTCTACCGCTGAAGACAAcc                                                                                                          >  1:139541/1‑251 (MQ=255)                                                                                                                                          cagtagaaAGTTGCTGAGTTCTTGGGTAAAGAGCCGCGTAAAGACGGTAACCCGGACGAAGCTGTAGCAATCGGGGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGGCTCTGGGTATCGAAACCATGGGCGGTGGGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAAcc                                                                                                          <  2:111446/246‑1 (MQ=255)                                                                                                                                          cAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAAcc                                                                                                          <  2:473316/250‑1 (MQ=255)                                                                                                                                             aagaaATTTGCTGAGGTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGGTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTc                                                                                                      <  2:147348/251‑1 (MQ=255)                                                                                                                                                    ttGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGAGGTGTTATGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGt                                                                                                >  1:55765/1‑250 (MQ=255)                                                                                                                                                                     aaaGAGCCGCGTAAAGAAGTTAACCCGGACGAAACAGTAGCAATCGGTGCTGCTGTTCAGGGAGGTGTTCTGACTGGTGAAGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGAACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGAGGTAACCAT‑CCATGtgctac                                                                              >  1:134474/1‑248 (MQ=255)                                                                                                                                                                        gagCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAACGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGCTCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGCCAACCAGTCTGCGGTAACCAT‑CAATGTGCTGCTggg                                                                          >  2:127253/1‑251 (MQ=255)                                                                                                                                                                            cgcgTAAAGACGTTAAACCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGACAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGCAACCCTGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCAGACCAAGCACCGCCAGGTGTTCTCTTCCGCGGCGTACAACCAGTCTGCGGTAAC‑ATACCATGTGGTGCAGGGTGa                                                                       >  2:108389/1‑250 (MQ=255)                                                                                                                                                                                 aaaGACGTTAGCCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGGCTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCGCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCAT‑CCATGTGCTGCAGGGTGAACGTaa                                                                 <  2:183981/251‑1 (MQ=255)                                                                                                                                                                                  aaGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGGGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACAAGTCTGCGGTAACCAT‑CCATGGGCTGCAGGGTGAACgaaa                                                                 >  2:355001/1‑247 (MQ=255)                                                                                                                                                                                  aaGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTATCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGCACGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCAA‑CAATGTGCTGCAGGGGGAACgaaa                                                                 >  2:434922/1‑247 (MQ=255)                                                                                                                                                                                         tAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGGTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGGCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCAT‑CCATGTGCTGCAGGGTGAACGTAAACGTGCg                                                          <  2:154582/250‑1 (MQ=255)                                                                                                                                                                                          aaCCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCAT‑CCATGTGCTGCAGGGTGAACGTAAACGTGCGGc                                                        >  1:85442/1‑251 (MQ=255)                                                                                                                                                                                          aaCCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCAT‑CCATGTGCTGCAGGGTGAACGTAAACGTGCGGc                                                        <  1:222941/251‑1 (MQ=255)                                                                                                                                                                                             ccGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCAT‑CCATGTGCTGCAGGGTGCACGTAAACGTGCGGCTg                                                      >  1:330459/1‑250 (MQ=255)                                                                                                                                                                                                   gAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCAT‑CCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACaa                                               <  1:421139/251‑1 (MQ=255)                                                                                                                                                                                                    aaGCTGTAGCAATCTGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTTCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATTGTCGGTGTGATTACGACGCTGATCGCGAAAAACACCACTATccc                                                                                                                                                   <  1:205359/151‑1 (MQ=255)                                                                                                                                                                                                    aaGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATccc                                                                                                                                                   >  2:205359/1‑151 (MQ=255)                                                                                                                                                                                                       gtGTAGCAATCCGGGCTGCTGTTGAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTTCTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGCCCAAGCAGAGCCAGTTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCAT‑CCATGTGCTGCCGGGTGAACGTAAACGGGCGGCTGATAGTAAAtct                                           <  2:200910/250‑1 (MQ=255)                                                                                                                                                                                                             aaagcGGTGCTGTGGTTCAGGTGGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCATGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCAT‑CCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGt                                     <  1:455467/247‑1 (MQ=255)                                                                                                                                                                                                                 cGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCAT‑CCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAg                                  <  1:154980/250‑1 (MQ=255)                                                                                                                                                                                                                     gctgctGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCAT‑CCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCa                              >  1:364523/1‑250 (MQ=255)                                                                                                                                                                                                                        gctgTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACCCCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCAT‑CCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAAc                            <  1:243782/249‑1 (MQ=255)                                                                                                                                                                                                                              cAGGGTGGTGGTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCAT‑CCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATgg                    <  1:277184/251‑1 (MQ=255)                                                                                                                                                                                                                                 ggttgtgtttggACTGGTGACGTGAAAGCCGTACTGCTGCGGGACGTTACCCCGCTGTCTCTGGGTATCGAAACTATTGGCGGGGTGATGAGGACGCTGATCGCGAAAAGCACCACTAACCCGACCAAGCACAGCTAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCAT‑CCATGTGCTGCAGGGTGAACGGAAACTTGCGGCTGATAACAAATCTCTGTGTCAGTTCAACCTGGAGGGTAt                 <  2:434962/240‑1 (MQ=255)                                                                                                                                                                                                                                  gtggtgTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCAT‑CCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATc                <  1:139841/251‑1 (MQ=255)                                                                                                                                                                                                                                             aCTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCAT‑CCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGAAcc     >  1:438885/1‑251 (MQ=255)                                                                                                                                                                                                                                                ggTGACGTGAAAGACGTACTGCTGCTGCACGTTACCCCGATGTCTCTGGGCATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCAT‑CCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACcgcg  >  1:126932/1‑251 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    AACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCAT‑CCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTAGATGGTATCAACCCGGCACCGCG  >  NC_000913/13073‑13562

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