Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,773,173	Δ1 bp	coding (893/1914 nt)	mtlA →	mannitol‑specific PTS enzyme IICBA component

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,773,173	0	G	.	96.8%	270.0 / 0.7	64	coding (893/1914 nt)	mtlA	mannitol‑specific PTS enzyme IICBA component
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base . (33/28);  total (34/30)
Fisher's exact test for biased strand distribution p-value = 2.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

GAAGCTAACCCAGGTCCAGGTATGGGCGTGCTGCTGGCGTACATGTTCTTTGGTCGTGGTAGCGCTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGGTTCTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCATCTCCGCTGTCTGCTGGCGATGTGACTAACGACCTGAGCCACGTACGTAAAATC  >  NC_000913/3772932‑3773420                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                        gAAGCTACCCCAGGTCCAGGTATGGGCGTGCTGCTGGCGTACATGTTCTTTGGTCGTGGTAGCGCTAAACAGTCTGCGGGCGGTGCGGCAATCATCCCCTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCtt                                                                                                                                                                                                                                               <  2:378756/251‑1 (MQ=255) gAAGCTAACCCAGGTCCAGGTATGGGCGTGCTGCTGGCGTACATGTTCTTTGGTCGTGGTAGCGCTAACCAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCt                                                                                                                                                                                                                                                <  2:168578/250‑1 (MQ=255)        aCCCAGGTCCAGGTATGGGCGTGCTGCTGGCGTACATGTTCTTTGGTCGTGGTAGCGCTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGGCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTAc                                                                                                                                                                                                                                        >  2:192544/1‑251 (MQ=255)         cccAGTTCCAGGGATGGGCGGGCGGCTGGCGTACATGTTCTTTGGTCGTGGTAGCGCTAAACAGGCTGCGGGCGGTGTGGCAATCATCCACGGCCTTGGGGGTATCCACGAAGTCTACTTCCCGTATGTGCTGATGAATGCTCGTCTGATCCTCGCAGTCAGCCGTGGCGGTAGGACTGGCGTGTTCACGCTGACTATCTTGGGGGGTGGTCTGGTTTCTCAGGCATCGCCGGGGT‑TATCCTTGCTGTACt                                                                                                                                                                                                                                       <  1:356396/251‑1 (MQ=255)             ggtccaggtATGGGCGTGCTGCTGGCGTACATGTTCTTTGGTCGTGGTAGCGCTAACCAGTCTGCGGGCGGTGCGGCAATCATCCACTTCATGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGGGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCGCCGG‑TTTTATCCTTGCTGTCCTGGCg                                                                                                                                                                                                                                   >  2:28824/1‑251 (MQ=255)                   ggtATGGGCGTGCTGCTGGCGTACATGTTCTTTGGTCGTGGTAGCGCTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGaca                                                                                                                                                                                                                             <  1:112457/251‑1 (MQ=255)                      aTGGGGGTGCTGCTGGCGTACATGTTCATTGGCCGTGGTAACGCTAAACAGCCTGCGGGCGCTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCACGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGAATGGCGTGTTCACGCTAACTATACTGGGCGTTGGCCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACAcc                                                                                                                                                                                                                           >  2:326352/1‑250 (MQ=255)                        gggCGTGCTGCTGGCGGACATGTTCTTTGGTCGTGGTAGCGCTAAACAGTCTGCGGGTGGTGCGGCAATCATCCACTTACTGGGGGGTATCCACGAAATCTACTTCCCGGATGTGCTGATGAATCCGCGTCTGAGCCTCGCAGGCATCCTCGGCGGTATGACGGTCGTGTGCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTGGCTGTACTGGCGATGACCCCaaa                                                                                                                                                                                                                        <  2:467050/251‑1 (MQ=255)                                         aCATGTTCTTTGGTCGTGGTAGCGCTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCt                                                                                                                                                                                                        <  2:465284/250‑1 (MQ=255)                                                  ttGGTCGTGGTAGCGCTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCCCGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGGACTGGCGATGACACCAAAAGGGGCTTACTTCGCTAACATCTCgg                                                                                                                                                                                              >  1:374708/1‑251 (MQ=255)                                                    ggttggtgTAGCGCTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTGCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCTCAGTCATCCTCGGCGTTATGACTGGCGTGTTCACGCTGACTATTCTGGGCGGTGGTCTGGTTTCGCCGGCATCTTCGG‑TTGTATCCTTGCTGTACTGGCGATGACATCAAAAGGTGCTTACTTCGCTAACATTGCggt                                                                                                                                                                                             <  2:360224/243‑2 (MQ=255)                                                       cGTGGTAGCGCTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGgtgt                                                                                                                                                                                          <  1:443489/250‑1 (MQ=255)                                                            tAGCGCTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGACATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCACAAGGTGCTTACTTCGCTAACATCGCGGTTGTGTTTgc                                                                                                                                                                                    >  1:296848/1‑251 (MQ=255)                                                            tAGCGCTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATACTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAACAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTgc                                                                                                                                                                                    >  1:168512/1‑251 (MQ=255)                                                             aGCGCTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATACGCGTCTGATACTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGCTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTgcg                                                                                                                                                                                   >  1:336726/1‑251 (MQ=255)                                                               cgcTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGTCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCCTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTgcggc                                                                                                                                                                                 <  1:40526/251‑1 (MQ=255)                                                                gcttaaCAGTCTGCGGGCGGGGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGTGTGTGTGTgcggcg                                                                                                                                                                                <  2:494685/247‑1 (MQ=255)                                                                   aaaCAGTCTGCGGGCGGTGCGGCATTCATCCACTTCCTGTGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTAGCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGGTT‑TATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTCCTTCGCTAACATCGCGGGTGTGTGTGCGGCGAt                                                                                                                                                                              <  1:527506/250‑1 (MQ=255)                                                                        ggtggggggcgggtcGGCAATCCTTCACTTCCTGGGGGGGTGCCACGAAATCCACCTTCCGGATGGGCTGATGAATCCGCGTCTGGTCCTCGGAGTTATCCTCGGCGGTATGACCGGCGGGTTCCCGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTTCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGGGCGGCGATGGCTGt                                                                                                                                                                        <  2:16831/237‑1 (MQ=255)                                                                         tCTGCGAGCGGTGCGGCGATCATCCACTTCCTGGGGGGTATACACGAAATCTACTTCCCGTATGTGCTGCTGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGCTGAAACCCAAAGGTGCTTACTTCGCTAACATTGCGGGTGTGTGTGCGGCGATGGCTGtt                                                                                                                                                                       >  2:38612/1‑250 (MQ=255)                                                                          cTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGCATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTCTCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCACAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGtcc                                                                                                                                                                      >  2:227876/1‑250 (MQ=255)                                                                              gggcggTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGGTGGCTGTCTCCtt                                                                                                                                                                  >  1:349898/1‑251 (MQ=255)                                                                               ggcggTGCGGCAATCAGCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTc                                                                                                                                                                 <  1:387953/251‑1 (MQ=255)                                                                                   gTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTg                                                                                                                                                             <  1:192544/251‑1 (MQ=255)                                                                                        gCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCCGTCTCCTTCGTTGTctct                                                                                                                                                        >  1:106106/1‑251 (MQ=255)                                                                                        gCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCGCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTGCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTctct                                                                                                                                                        <  2:523344/251‑1 (MQ=255)                                                                                          gatcatcCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGAATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTAGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACAGCACAAGGGGCTTACTTCGTTTACATCGCGGGTGGGTGTGCGGCGATGGCGGTCTCCCTCGGTTTCTCTGc                                                                                                                                                      >  2:137881/2‑251 (MQ=255)                                                                                           atcatcCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGCCTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTTCTTCGCTCACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCt                                                                                                                                                     >  1:30441/1‑251 (MQ=255)                                                                                                     tCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTg                                                                                                                                            <  1:473766/250‑1 (MQ=255)                                                                                                        tGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCAGCGTCTGATCCTCGCAGTCATCCTCGGCGGTCTGCCTGGCGTGTTCCCGCTGACTATCCTGGGAGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGAAACCAAAAGGTGCTTACTTCGCTCACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGaaa                                                                                                                                         >  1:255902/1‑250 (MQ=255)                                                                                                          gggggTATCCACGACATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAAcc                                                                                                                                      >  1:423536/1‑251 (MQ=255)                                                                                                                   cACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTg                                                                                                                             >  1:271058/1‑251 (MQ=255)                                                                                                                     cGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACGGGCGTGTTCACGCTGACTATCCTGGGCTGTGGTCTTGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACTCCAAAAGGTGCTTACTGCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTTTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGaa                                                                                                                           <  1:489563/251‑1 (MQ=255)                                                                                                                        aaTCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGaa                                                                                                                           >  1:532836/1‑248 (MQ=255)                                                                                                                             aCTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACGGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAaagaagaa                                                                                                                    <  2:405215/250‑1 (MQ=255)                                                                                                                                tCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAgat                                                                                                                 <  1:458895/250‑1 (MQ=255)                                                                                                                                 cccGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAgatga                                                                                                               <  1:482955/251‑1 (MQ=255)                                                                                                                                     tATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATAtt                                                                                                           >  1:303981/1‑251 (MQ=255)                                                                                                                                          gCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAag                                                                                                       <  1:479203/250‑1 (MQ=255)                                                                                                                                           cTGATGAATCCGCGACTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAagc                                                                                                      >  1:90421/1‑250 (MQ=255)                                                                                                                                                gaATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAAc                                                                                                <  1:149899/251‑1 (MQ=255)                                                                                                                                                  aTCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCCAAAGGTGCTTACTTCGCTAACATCGCGGGTGGGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACtc                                                                                              >  1:512663/1‑251 (MQ=255)                                                                                                                                                       cgTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTcgtcgt                                                                                         >  2:416874/1‑251 (MQ=255)                                                                                                                                                             aTCCTCGCAGTCCTCCTAGGCGGTAAGACTGGCGTGTTCACGCTGACTATCCTGGTAGGGGGTCTGGTTTCTTCGGCATCACCGGCTT‑TATCCTTGCTGTACTGGCGATGACCCCAAAAGGTGCTTACTTCGCTAACCTCGCGGGCGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGc                                                                                     >  2:543423/1‑249 (MQ=255)                                                                                                                                                             aTCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCCCGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCac                                                                                   <  2:303981/251‑2 (MQ=255)                                                                                                                                                             aTCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATACTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGGGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGc                                                                                     >  1:345450/1‑249 (MQ=255)                                                                                                                                                                cTCGCAGTCATCCTCGGCGGTACGACTGGCGTGTTCAAGCGGACTTTCCTGGGCGGTGGTCTGGTTTCTCCGGGATCTCCGG‑TTTTATCCTTGCCGTACTGGCGATGACACCAACAGGTGATTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGGCGTATGCAGGAc                                                                                >  2:389110/1‑251 (MQ=255)                                                                                                                                                                   gCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATg                                                                             <  1:416874/251‑1 (MQ=255)                                                                                                                                                                      gTCATCCTCGGCGGTATGACTGGCGTGTTCTCGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTATCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGaaa                                                                          <  1:336442/251‑1 (MQ=255)                                                                                                                                                                         aTCCTCGGCGGTATGACTGGCGTGTTCACGCTGCCTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCt                                                                       >  1:100864/1‑251 (MQ=255)                                                                                                                                                                              cggcggTATGACTGGCGTGTTAACGGTGGCTATCATGGGCGGTGGGCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGGTGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCCGCCAAGTGAAAGAAGAAGATTATATTGAAACAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTc                                                                  >  2:527900/1‑251 (MQ=255)                                                                                                                                                                                    tATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAgg                                                            >  2:163353/1‑251 (MQ=255)                                                                                                                                                                                      tGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCCAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGgcg                                                          >  2:250125/1‑251 (MQ=255)                                                                                                                                                                                       gACTGGAGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATACTTGCTGTACTGGCGATGACACACAAAGGTGCTTACTTCGCTAACATCGCGGGTGAGTGTGCGGCGATGGCTGTCTCCTTCGTAGTCTCTGCTTTTTTGATGAAAACCAGCAACGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGacgc                                                         >  1:549135/1‑251 (MQ=255)                                                                                                                                                                                          tGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCAtc                                                      >  1:397553/1‑251 (MQ=255)                                                                                                                                                                                           ggCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCAtct                                                     <  2:168176/251‑1 (MQ=255)                                                                                                                                                                                                      gCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCATCTCCGCTGTctgc                                          >  2:373302/1‑251 (MQ=255)                                                                                                                                                                                                         gACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCATCTCCGCTGTctgct                                         <  2:439741/249‑1 (MQ=255)                                                                                                                                                                                                         gACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCATCTCCGCTGTCTGCTgg                                       <  2:397553/251‑1 (MQ=255)                                                                                                                                                                                                             aTCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCATCTCCGCTGTCTGCTGGCGAt                                   >  2:88547/1‑251 (MQ=255)                                                                                                                                                                                                                     cGGTGGTCTGGTTTCTACGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCATCTCCGCTGTCTGCTGGCGATGTGACTa                            >  2:252936/1‑250 (MQ=255)                                                                                                                                                                                                                       gtggtCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCATCTCCGCTGTCTGCTGGCGATGTGACTAACg                         <  2:297654/251‑1 (MQ=255)                                                                                                                                                                                                                       gtggtCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCATATCCGCTGTCTGCTGGCGATGTGACTAACg                         <  1:544567/251‑1 (MQ=255)                                                                                                                                                                                                                                  tctcCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACAACAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCATTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAACGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCATCTCCGCTGTCTGCTGGCGATGTGACTAACGACCTGAGCCac              >  1:61369/1‑251 (MQ=255)                                                                                                                                                                                                                                               ccggTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGTACATGAAAGCTGAGTCTAAAGGCGCATCTCCGCTGTCTGCTGGCGATGTGACTAACGACCTGAGCCACGTACGTAAAatc  <  1:527900/249‑1 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                        GAAGCTAACCCAGGTCCAGGTATGGGCGTGCTGCTGGCGTACATGTTCTTTGGTCGTGGTAGCGCTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGGTTCTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCATCTCCGCTGTCTGCTGGCGATGTGACTAACGACCTGAGCCACGTACGTAAAATC  >  NC_000913/3772932‑3773420

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