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	96.2%	142.8 / ‑2.3	52	V386V (GTA→GTG)	dnaK	chaperone protein DnaK
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (35/15);  total (36/16)
Fisher's exact test for biased strand distribution p-value = 5.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.40e-01

AAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTAGATGGTATCAACCCGGCACCGCGC  >  NC_000913/13079‑13563                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    aaaGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCCGGGTGGTGTTCTGACTGGTGACGTGAAAGACGt                                                                                                                                                                                                                                             >  1:232103/1‑250 (MQ=255)        ggTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAACACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTACGGGTGGTGTTCTGACTGGTGACGTGAATGACGTActgatgc                                                                                                                                                                                                                                     >  2:24967/1‑251 (MQ=255)            gAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGACAGTCGTACTGCTGCTggg                                                                                                                                                                                                                                 >  2:23124/1‑250 (MQ=255)                             ttCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGAAGCAATCGGTACTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCAGGACGTTACCCCGCTGactc                                                                                                                                                                                                                >  1:4742/1‑251 (MQ=255)                                   tGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGATGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGCCGTTACCCCGCTGTCTCTGGGt                                                                                                                                                                                                           >  2:419265/1‑250 (MQ=255)                                                   gCACGGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCgg                                                                                                                                                                                          <  1:470615/251‑1 (MQ=255)                                                     aCTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGgtg                                                                                                                                                                                        >  1:465092/1‑251 (MQ=255)                                                      cTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGgtgt                                                                                                                                                                                       >  1:65131/1‑251 (MQ=255)                                                         cAGGACGCTGGCCTGTCCGTATCTGATATCGACGTCGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTCAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGAGGTGTTCTGACGGGTGACGTGAAAGACGTACTGCTGCTGAACGTTACCCCGCTGCCTCTGGGTATCGAAACCATGGGGGGTGtgat                                                                                                                                                                                    >  2:315278/1‑251 (MQ=255)                                                          aGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAACGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGAGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGCCGTTACCCCGCTGTCTCTGGGTCTCGAAACCATGGGCGGTGtgatg                                                                                                                                                                                   >  1:170276/1‑251 (MQ=255)                                                               gCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATgacgac                                                                                                                                                                              <  1:259451/251‑1 (MQ=255)                                                                  ggCCTGTCCGTTACTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGTCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCt                                                                                                                                                                           <  2:181632/251‑1 (MQ=255)                                                                        tCCGTATCTGATATCGACGACGTTATAATCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATcgc                                                                                                                                                                     >  1:437706/1‑251 (MQ=255)                                                                         ccggatcggATATCGACGACGATATCCTCGTTGGTGGTCAGACACGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGGTGTAGCAATCGGTGCTGCTGTTCAGGGTGGCGTTCTGACTGGTGACGCGAAAGACGTACGGCTGCTGGACGTTACCCAGGTGTCTCTGGGTATCGCAAACATGGTCGGTGTGAGTACGACACTGacaacg                                                                                                                                                                    >  2:417147/9‑246 (MQ=255)                                                                            tATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGCCTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGTTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGGCGACGCTGATCGCGaa                                                                                                                                                                  >  1:280888/1‑250 (MQ=255)                                                                                       gacgacGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTTGTAAAGAGCCGCGTAAAGACGTGAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTAt                                                                                                                                                      <  2:217912/251‑1 (MQ=255)                                                                                       gacgacGTTATCCGCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACCCCACTAt                                                                                                                                                      <  1:410941/251‑1 (MQ=255)                                                                                        acgacgTTAACCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCGTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATc                                                                                                                                                     <  1:250957/251‑1 (MQ=255)                                                                                          gacgTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATccc                                                                                                                                                   <  1:99433/251‑1 (MQ=255)                                                                                                 tCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGATGTTCAGGGTGGTGGTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACTCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTAACCCGGCCaa                                                                                                                                             >  1:138005/1‑250 (MQ=255)                                                                                                  ccTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGc                                                                                                                                           >  1:446386/1‑251 (MQ=255)                                                                                                   cTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGCACGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAATCATGGGCGGTGTGATGACGCCGCTGGTCGCGAAAGACACCACTATCCCGACCacgc                                                                                                                                           >  2:176207/1‑247 (MQ=255)                                                                                                         ggtggtCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTTCCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCa                                                                                                                                    <  1:295432/251‑1 (MQ=255)                                                                                                                   cTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGGCGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTctct                                                                                                                          >  1:170854/1‑251 (MQ=255)                                                                                                                    tCGTATGCCAATGGTTCAGCAGAAAGTTGCTGAGTTCTTTGGTACAGAGCCGCGTAAAGACGTTAACCCGGCCGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGCCGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACACTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTctct                                                                                                                          >  2:432176/1‑250 (MQ=255)                                                                                                                     cGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTctct                                                                                                                          <  2:112853/249‑1 (MQ=255)                                                                                                                      gtaagCCAATGGTTGATACGAAAGATGATGAGTTCTTTTATAAACAGACACGTAATGACGTTAAACCGGACGAAGCTGTCGCAGTCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGGAAGACGTACTGCTGCCGGCCGATACCCAGATGTCTCTGGGTATTGAAACCATGGGCGGTGTGCTAACGAAGCTGATCGCGAAAAACACCACTACCCCGCCCAAGCTCAGCCAGGTGTTCTGTAcc                                                                                                                       >  2:497685/5‑251 (MQ=255)                                                                                                                            cAATGGTTCAGAAGAAAGCTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGa                                                                                                                  >  1:197380/1‑250 (MQ=255)                                                                                                                                 gTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCg                                                                                                                                                                                                       >  2:30830/1‑160 (MQ=255)                                                                                                                                 gTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCg                                                                                                                                                                                                       <  1:30830/160‑1 (MQ=255)                                                                                                                                            aGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAGAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGCGGTGTTCTGACTGGTGACGTGCAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCgg                                                                                                 >  1:65800/1‑251 (MQ=255)                                                                                                                                              ttGCTGAGTTCTTTGGAAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTAAGGGAGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGAACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACCCCACTATCCCGACCACGCACAGCCAGGGGTTCTCACCCGCAGAAGACAACCaagcggcggt                                                                                                >  2:427152/1‑241 (MQ=255)                                                                                                                                                  tGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCAGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGGGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAAcc                                                                                            >  1:373511/1‑250 (MQ=255)                                                                                                                                                  ggAGTTCTATGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTTCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAATACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGTTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCTCTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTACcca                                                                                           <  2:327497/250‑1 (MQ=255)                                                                                                                                                     gTTCTTGGGTAAAGAGCCGCGTAAAGTCGTTAACCCGGTCGAAGCTGTAGCAATCGGTTCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCACTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCATGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAAccatcc                                                                                        <  1:495900/251‑1 (MQ=255)                                                                                                                                                        cTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCAtg                                                                                     >  1:184368/1‑251 (MQ=255)                                                                                                                                                        cTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGAAGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACACCCAGGTGTTCTCTACCGCTGAAGACAACCCGTCTGCGGCAAccacccagg                                                                                     >  2:328123/1‑249 (MQ=255)                                                                                                                                                            ggTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGtgc                                                                                  <  1:215729/250‑1 (MQ=255)                                                                                                                                                                    gCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTg                                                                         <  1:432176/251‑1 (MQ=255)                                                                                                                                                                       gcgTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAAc                                                                      >  1:135813/1‑251 (MQ=255)                                                                                                                                                                            aaGACGTTTACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCGGCTGGACGTTACCCCGTTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTGGACCGCGGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTaaa                                                                 <  1:329335/251‑1 (MQ=255)                                                                                                                                                                                           aCGAAGCTGTAGCAATAGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAAc                                                  >  1:445558/1‑251 (MQ=255)                                                                                                                                                                                             gAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACa                                                 >  1:156271/1‑250 (MQ=255)                                                                                                                                                                                             gAAGCTGTAGCAATCGGAGCTGCTGTTCAGGGTGGTGTTCTGACTGGGGACGTGAAAGAAGTCCTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTATACCGCCTAAGACAACCAGTCTGCGGGAACCATCCAGGTGCTGCAGGGGGAACGTAAACGTGAGGATGATAACa                                                 >  2:392273/1‑250 (MQ=255)                                                                                                                                                                                               aGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGATCTGACTGGTGACGTAAAAGACGTACTGCGGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACAAAGCACAGCCAGGAGTTCTCTACCGCTGAAGACAACCAGTCTGCGGCAACCATCCATGTGCAGCATGGTGCACGTAAACGTCCGGCTGATAACAAAt                                              >  1:357407/1‑251 (MQ=255)                                                                                                                                                                                                gCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAAtc                                             >  2:458503/1‑251 (MQ=255)                                                                                                                                                                                                 cTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAAtct                                            >  1:505249/1‑251 (MQ=255)                                                                                                                                                                                                   gTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAAtctc                                           <  1:23124/250‑1 (MQ=255)                                                                                                                                                                                                            ggTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCCTGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGtt                                 >  2:58139/1‑251 (MQ=255)                                                                                                                                                                                                              tgctgcTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTc                                >  1:369629/1‑250 (MQ=255)                                                                                                                                                                                                                                          ggTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCCAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACcgcg   >  1:217284/1‑251 (MQ=255)                                                                                                                                                                                                                                           gTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACACGGCACcgcgc  >  1:218724/1‑251 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    AAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTAGATGGTATCAACCCGGCACCGCGC  >  NC_000913/13079‑13563

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