Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	13,932	A→G	E590E (GAA→GAG)	dnaK →	chaperone protein DnaK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	13,932	0	A	G	96.8%	186.2 / ‑2.7	64	E590E (GAA→GAG)	dnaK	chaperone protein DnaK
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (22/39);  total (23/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.08e-02

ACGAAGATGAAATCCAGAAAATGGTACGCGACGCAGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCT  >  NC_000913/13685‑14173                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                  aCGAAGAAGACATCCAGAAAAGTGTACGCGACGCACAAGCTAACGCCGAAGCTGCCCGTAAGGTTGAATAGCTGGTACAGCCGCGCAACCGGTGCGACCATCTGCTGCACAGCACACGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCGCTGACTGCACTGGAAACTGCGCGGAAAGGTTAAGACAAATCCGCTATCGAAGCGAAAATGCAGGATCTg                                                                                                                                                                                                                                                <  2:353354/251‑1 (MQ=255)   gAAGATGAAATCCAGAAAATGGTACGCGACGCAGAAGCTAACGCCGAAGCTGACCGTAAGATTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGATGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGc                                                                                                                                                                                                                                              >  1:659647/1‑251 (MQ=255)     aGATGAAATCCAGAAAATGGAACGCGCCGCAGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGCCGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGcac                                                                                                                                                                                                                                            >  1:39778/1‑251 (MQ=255)      gATGAAATCCAGAAAATGGTACGCGACGCAGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGcaca                                                                                                                                                                                                                                           <  1:513955/251‑1 (MQ=255)       aTGAAATCCAGAAAATGGTACGCGACGCAGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGcac                                                                                                                                                                                                                                            >  1:419009/1‑249 (MQ=255)              ccAGAAAAGGGTACGCGACGCAGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTccc                                                                                                                                                                                                                                   <  1:370012/251‑1 (MQ=255)                   aaaTGGTACGCGACGCAGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGaa                                                                                                                                                                                                                               >  1:635245/1‑250 (MQ=255)                               aCGCAGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATc                                                                                                                                                                                                                  <  1:478475/251‑1 (MQ=255)                                   aGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGccc                                                                                                                                                                                                              >  1:261512/1‑251 (MQ=255)                                              ccGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACAt                                                                                                                                                                                                   >  1:398963/1‑251 (MQ=255)                                               cGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACAtc                                                                                                                                                                                                  >  1:582836/1‑250 (MQ=255)                                                 aaggtgacctgaaGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGcc                                                                                                                                                                                                <  2:189840/240‑1 (MQ=255)                                                  aGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCAGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGccc                                                                                                                                                                                               <  2:84141/251‑1 (MQ=255)                                                     tGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCcagc                                                                                                                                                                                            >  1:324845/1‑251 (MQ=255)                                                      gACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTcc                                                                                                                                                                                                                                    <  1:589569/210‑1 (MQ=255)                                                       aCCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGtt                                                                                                                                                                                                                                       >  2:202782/1‑206 (MQ=255)                                                       aCCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGtt                                                                                                                                                                                                                                       <  1:202782/206‑1 (MQ=255)                                                             aggttGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCCGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACt                                                                                                                                                                                    <  2:194659/248‑1 (MQ=255)                                                                ttGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGc                                                                                                                                                                                  <  1:566718/250‑1 (MQ=255)                                                                ggtaaGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGCAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGc                                                                                                                                                                                  <  2:658765/247‑1 (MQ=255)                                                                             tcAGACTCGCATCCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCt                                                                                                                                                                     <  2:582836/249‑1 (MQ=255)                                                                                   ccGCAACCAGGGCGACCATCTGCTGCACTGCACCCGTAAGCAGGTTGAAGAAGCAGGCGCCAAACGGCCGTGTGACGACAAAACTGGTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCTAAAATGCAGGAGCTGGCGCAGGTGTCCCAGAAGCTGATGGAAACCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGGGCTGATGTTTCTGCaa                                                                                                                                                              <  2:39778/250‑1 (MQ=255)                                                                                        aCCAGGGCGACCATCTGCTTCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACTACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAaacaac                                                                                                                                                         <  1:3075/251‑1 (MQ=255)                                                                                             ggCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGa                                                                                                                                                     >  1:449286/1‑250 (MQ=255)                                                                                                gACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGa                                                                                                                                                 <  1:254025/251‑1 (MQ=255)                                                                                                 aCCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGaaa                                                                                                                                                                                                                              >  1:373150/1‑173 (MQ=255)                                                                                                   cATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATg                                                                                                                                               >  1:50874/1‑250 (MQ=255)                                                                                                           gcACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGAtgt                                                                                                                                        <  1:258306/249‑1 (MQ=255)                                                                                                              caGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCg                                                                                                                                   <  1:173139/251‑1 (MQ=255)                                                                                                              caGCACCCGTAAGCAGGGTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAAATGCTATCGAGTCTGCACTTACTGCACGGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCg                                                                                                                                   <  2:419009/251‑1 (MQ=255)                                                                                                                  aCCCGTAAGCAGGTTAAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGAAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGc                                                                                                                               <  2:171408/251‑1 (MQ=255)                                                                                                                   cccGAAAGCAGGTTGAAGTAGCATGCGACAAACAGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAATCGAAAATGCAGGAGCTGTCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCATCAACATGCCCAGCAGCAGGCTGCCGGTGCTGATGCTTCTGCAAACAACGCTAAAGATGGCGATGTTTTCTACGCt                                                                                                                              <  2:68820/251‑1 (MQ=255)                                                                                                                     cGTAAGCAGGTTGAAGAATCATGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTTCACGGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGa                                                                                                                            <  2:449286/251‑1 (MQ=255)                                                                                                                           cAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTg                                                                                                                       >  1:322042/1‑250 (MQ=255)                                                                                                                            aGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTgaa                                                                                                                     <  1:190698/251‑1 (MQ=255)                                                                                                                             ggTTGAAGAAGCAGGCGACAAACTGCCGGCGGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCGTCTGCAAACAACGCGAAAGATGACGAGGTTGTCGACGCTGAATTTgaag                                                                                                                    <  2:289884/251‑1 (MQ=255)                                                                                                                                      agCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGAATGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCAGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCACCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGATGATGCCTATGCAAACAACGCGAACGATGACGCTGGTGTCGTCCCTGAATTTGAAGAAGTcca                                                                                                             >  2:319034/1‑247 (MQ=255)                                                                                                                                           gCGAGAAACTGCCGGCTGACGACAAAACGGCTATCGAGTCTGCCCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTGCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACaaa                                                                                                      <  2:13123/251‑1 (MQ=255)                                                                                                                                           gCGACAAACTGCCGGCTGACGCCAAAACTGCTATCGAGTCTGCACTTACTGCACTGGAAACTGCGCTGAAAGGAGAAGACAAAGCCGGTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCTAAGACa                                                                                                        <  2:599766/249‑1 (MQ=255)                                                                                                                                                  aCTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCa                                                                                                                                                               <  2:292296/187‑1 (MQ=255)                                                                                                                                                  aCTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCa                                                                                                                                                               >  1:292296/1‑187 (MQ=255)                                                                                                                                                  aCTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAaataat                                                                                               >  1:281858/1‑251 (MQ=255)                                                                                                                                                         gCTGACGACAAAACTGCTATCGAGTCTGCACTGGCTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCta                                                                                        <  1:356104/251‑1 (MQ=255)                                                                                                                                                            gacgacAAAACTGGAATCGAGTCTTCACCGACTGCACTGGAAACTGCTCTGAAAGTTTAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAAGTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATaa                                                                                     <  2:580144/251‑1 (MQ=255)                                                                                                                                                            gacgacAAAACTGCTATCGAGCCTGCACTGACTGCACTGGAAACTGCACTGAAAGGTCAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCGGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGGTTCTGCAAACAACGCGAAAGATGACGATGTGGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAAGCGCCCTATaa                                                                                     <  2:456660/251‑1 (MQ=255)                                                                                                                                                             acgacAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGAAACCACCGCGAAAGAGGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAATTAATAGCCCTATaa                                                                                     >  2:411133/1‑250 (MQ=255)                                                                                                                                                             acgaaaaaccttcTATCTAGTCTGCACTTACTGCACTGGAAACTGCTTTGAAAGGTGAAGACAAAGCCGCTATCGTAGCGAAAATGCAGTAGCTTGCACAGGTTTCTCAGAAACTGATGGAAATCTCCCAGCAGCAACATGCCCAGCAGCAGACTTCCGGTGCTGATTCTTCTTCAAACAACTCGAAAGATGACGATGGTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATaaa                                                                                    <  2:659647/239‑1 (MQ=255)                                                                                                                                                                   aaaCTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCg                                                                                                                                   <  1:105895/198‑1 (MQ=255)                                                                                                                                                                   aaaCTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCg                                                                                                                                   >  2:105895/1‑198 (MQ=255)                                                                                                                                                                        gCTATCGAGTCTGCACTGATTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAAtt                                                                           <  1:470773/249‑1 (MQ=255)                                                                                                                                                                          tATCGAGTCTGCGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATAc                                                                       >  1:432586/1‑251 (MQ=255)                                                                                                                                                                              gAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGac                                                                   >  1:436204/1‑251 (MQ=255)                                                                                                                                                                                    gcacggACTGCACTGGAAACTGCTCGTAAAGGTTAAGACAAAGCCGCTAACGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACg                                                                <  2:497235/243‑1 (MQ=255)                                                                                                                                                                                    gcacTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACgg                                                               <  1:648142/246‑1 (MQ=255)                                                                                                                                                                                                              gAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGccc                                      <  1:469350/248‑1 (MQ=255)                                                                                                                                                                                                                     gAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCAt                            <  1:201051/251‑1 (MQ=255)                                                                                                                                                                                                                        gACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCt                          >  1:107992/1‑250 (MQ=255)                                                                                                                                                                                                                           aaaGCGGCTATCGAAGCGAAAATGCAGGAGCTGGCGCAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGTGTAATTATACTGACACGTGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAggg                      <  2:324845/251‑1 (MQ=255)                                                                                                                                                                                                                              gCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTGCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCGATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCa                   <  1:444915/251‑1 (MQ=255)                                                                                                                                                                                                                                   tATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATGGCCCAGCAGCACCATGCCCAGCAGCAGCCTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCACGTGCATTCATCTAggggacattt               >  1:642489/1‑244 (MQ=255)                                                                                                                                                                                                                                    aTCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTa              <  1:206264/250‑1 (MQ=255)                                                                                                                                                                                                                                              aaaTGCAGTAGCTGGCACAGCGTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCAGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGTCGCTGAATTTgaagaa                                                                                                                  <  2:328517/140‑1 (MQ=255)                                                                                                                                                                                                                                              aaaTGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTgaagaa                                                                                                                  >  1:328517/1‑140 (MQ=255)                                                                                                                                                                                                                                               aaTGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCt  <  1:578830/251‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                  ACGAAGATGAAATCCAGAAAATGGTACGCGACGCAGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCT  >  NC_000913/13685‑14173

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