Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	67,612	G→T	A242A (GCC→GCA)	araA ←	L‑arabinose isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	67,612	0	G	T	96.2%	176.1 / NA	53	A242A (GCC→GCA)	araA	L‑arabinose isomerase
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base T (25/26);  total (25/28)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCAGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCAGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTTGCCGTGGATTTGTGTG‑GCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCGCCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAATCGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAGACCGCCTGACGCATCC  >  NC_000913/67368‑67850                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                cGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGGCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTT‑GCAGGc                                                                                                                                                                                                                                          <  1:5518/251‑1 (MQ=255)       aCATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTT‑GCAGGCGTCATg                                                                                                                                                                                                                                    <  1:338822/251‑1 (MQ=255)        cATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTT‑GCAGGCGTCATgg                                                                                                                                                                                                                                   <  2:266455/251‑1 (MQ=255)        cATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGGGTT‑GCAGGCGTCATgg                                                                                                                                                                                                                                   <  1:137468/251‑1 (MQ=255)            aCCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTT‑GCAGGCGTCATGGTGTa                                                                                                                                                                                                                               <  1:116503/251‑1 (MQ=255)                     atgCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTT‑GCAGGCGTCATGGTGTAGCAGCTTTc                                                                                                                                                                                                                      >  2:38570/1‑251 (MQ=255)                                   ggcggTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTT‑GCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGAcc                                                                                                                                                                                                         >  2:233800/5‑250 (MQ=255)                                   ggcggTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTT‑GCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGAcc                                                                                                                                                                                                         >  2:3728/5‑250 (MQ=255)                                               gTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTT‑GCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACAt                                                                                                                                                                                            >  1:210762/1‑251 (MQ=255)                                                 cGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTT‑GCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATc                                                                                                                                                                                           <  1:224497/250‑1 (MQ=255)                                                   ccTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTT‑GCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGcc                                                                                                                                                                                        >  2:83978/1‑251 (MQ=255)                                                      tCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTCTTCGCCGTGGATTTGTGTT‑GCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCg                                                                                                                                                                                       <  2:186395/249‑1 (MQ=255)                                                        gCCCGCAAATCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAATGTTGTGGTGAACGCGTGGAAGCCACCTTGTGCCAGGAAACGCTTCATTCCCAGCTCAATACGCGCCGCTTCCAGCACGGTCTTTCGTTTTTCGCCGTGGATTTGTGTG‑GCATGCGTCATGGTGTATCAGCTTTCTTACTCATCGACCAGCGCGTTAACATCGCCGTCGc                                                                                                                                                                                   <  1:172605/251‑1 (MQ=255)                                                                 gCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTT‑GCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGt                                                                                                                                                                          <  1:290455/251‑1 (MQ=255)                                                                 gCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTT‑GCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGt                                                                                                                                                                          <  1:233800/251‑1 (MQ=255)                                                                               gcATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTT‑GCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCAc                                                                                                                                                             <  1:358089/250‑1 (MQ=255)                                                                                    aGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTT‑GCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATc                                                                                                                                                       <  1:301566/251‑1 (MQ=255)                                                                                        gCTGTACGCCCAAACCTGGAAGCTGTTTCAGACCGAGCAACTCTTCAGAGGAGCTGGTGAACGCGTGGCAGCCACCTTGTCCCAAAAAACGCTTCCTCCCCAGCTCAACACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTT‑GCAGGAGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACCTCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGcc                                                                                                                                                    >  2:366441/1‑250 (MQ=255)                                                                                           gTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTT‑GCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAAcc                                                                                                                                                >  2:228358/1‑251 (MQ=255)                                                                                              cGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTT‑GCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGcc                                                                                                                                             <  1:360827/251‑1 (MQ=255)                                                                                                  cAGACCTGGAAGCTGTTTCATACCATGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCTACCTTGGTCCAGGAATCCCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTTGTTTTGCGCCGTGGATTTGTTTT‑GCAGGCGTTATGGTGTGGCAGCTTTCGTACTCATCGACCAGCGCGTTAACTTCGCCGTCGCGGATGGAGATCACCACCTGCACCAGATCGCCAACCGCCCAgg                                                                                                                                         <  1:353381/251‑1 (MQ=255)                                                                                                    gACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTT‑GCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGt                                                                                                                                        >  1:228402/1‑250 (MQ=255)                                                                                                      cctggAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTT‑GCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTAtt                                                                                                                                     <  2:197980/248‑1 (MQ=255)                                                                                                         ggAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTT‑GCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGTAGTTCACCGCCGGCACCAGAGCGCCAACCGCCCAGGTATTg                                                                                                                                    <  2:21498/249‑1 (MQ=255)                                                                                                                  ttCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTT‑GCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAAc                                                                                                                          <  2:23816/250‑1 (MQ=255)                                                                                                                   tCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTT‑GCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTca                                                                                                                                                                       >  2:105245/1‑204 (MQ=255)                                                                                                                   tCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTT‑GCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTca                                                                                                                                                                       <  1:105245/204‑1 (MQ=255)                                                                                                                    cAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTT‑GCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCg                                                                                                                        <  2:161481/250‑1 (MQ=255)                                                                                                                     aGACCGTGCAAATCTTCAAAGGTGGTGGTGCACGCGTGGAAACCAACTTCTTCCAGGAAACGCTTCATCCCCAGCTCCATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTT‑GCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATAGCCGTCGATGATGGCGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGa                                                                                                                       >  1:148047/1‑250 (MQ=255)                                                                                                                         cGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTT‑GCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACAGCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTg                                                                                                                  <  2:54707/251‑1 (MQ=255)                                                                                                                              aaaTCTTCAAAGGTTGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCGTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTT‑GCAGGCCTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCtg                                                                                                             <  2:204129/251‑1 (MQ=255)                                                                                                                               aaTCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTT‑GCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCtgt                                                                                                            >  1:371077/1‑251 (MQ=255)                                                                                                                                  cgtCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTT‑GCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCg                                                                                                         <  2:210762/249‑1 (MQ=255)                                                                                                                                     cAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCTCGTTCTGTCGTTTTTCGCCGTGGATTTGTGTT‑GCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGTTGATGGAGTTCACCACCTGCTCCTGATCGCCAACCGCCCAGGTATTGACGGAGAAACCTAACTTGATCTGTGCGGc                                                                                                       <  2:179996/250‑1 (MQ=255)                                                                                                                                          gtggtggtgAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTT‑GCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGACACCGAACTTGATCTGTGCGGCAACttt                                                                                                 >  2:339094/1‑251 (MQ=255)                                                                                                                                          gtggtggtgAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTT‑GCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGAGCTGTGCGGCAACttt                                                                                                 <  2:146977/251‑1 (MQ=255)                                                                                                                                            ggtggtgAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTT‑GCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTAt                                                                                               <  2:204851/251‑1 (MQ=255)                                                                                                                                                                 aCCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTT‑GCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACtt                                                                          >  1:198436/1‑251 (MQ=255)                                                                                                                                                                  ccTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTT‑GCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACTTc                                                                         >  1:179127/1‑251 (MQ=255)                                                                                                                                                                             gAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTT‑GCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTAt                                                              >  1:18990/1‑251 (MQ=255)                                                                                                                                                                             gAAACGCTTCATCCACCGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTT‑GCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTAt                                                              >  1:112996/1‑251 (MQ=255)                                                                                                                                                                                  gCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTT‑GCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGcc                                                          >  2:81523/1‑250 (MQ=255)                                                                                                                                                                                      cATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTT‑GCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCaaaac                                                     >  2:132372/1‑249 (MQ=255)                                                                                                                                                                                         ccccAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTT‑GCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATcacc                                                                                           >  2:232930/1‑207 (MQ=255)                                                                                                                                                                                         ccccAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTT‑GCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATcacc                                                                                           <  1:232930/210‑4 (MQ=255)                                                                                                                                                                                         ccccAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTT‑GCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGc                                                  >  1:161347/1‑251 (MQ=255)                                                                                                                                                                                                      cgcgCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTT‑GCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATg                                     >  1:375494/1‑251 (MQ=255)                                                                                                                                                                                                                 cggccctttCTGTCGTTTTTCGCCGTGG‑TTTTTGTTTGCAGGCGTCATGGTGTAGCATCTTTCGTACTCATCGACCAGCCCGTAAACATCGCCGTCGCTGATGGAGTTCCCCCCCTGCCCAAGTTCCCCACCCGCCCGGGTATTGCCGGAAAACCCGAACTTGATCTGTGCGGCAACTTTTTCCCCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCt                          <  1:56088/244‑1 (MQ=255)                                                                                                                                                                                                                     aCGTTCTGTCGTTTTTCGCCGTGGATTTGTGTT‑GCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGttt                      >  1:123889/1‑251 (MQ=255)                                                                                                                                                                                                                         tCTGTGGTTTTTCGCCGGGGGCTTGTGTT‑GCAGACGTCATGGTGTAGCGGCTTTCGTACTCATCGACCAGCGCGTTAACACCGCCGTCGCTGGTGGAGTTCCCCACCTGCACCAGATCGCCAACCGCCCGGGTATTGACGGAGAACCCGCACTTGATCTGTCCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCACAACGGCATACTTTCAGATGACGGGTATCCTGTTTaga                   >  2:16595/1‑250 (MQ=255)                                                                                                                                                                                                                              cgattgtcgCCGCGGATTTGTGTT‑GCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCAAACCGCCCAGGCATTGCCGGAGAAACCGAACTTGATCTGCGCGGCAACTTTATCACCATCGGGGACCGCCACTTCACGCATGTTATCGCCAAAACGGCCGACTTTCCGATGACGGGTCTCCTGTTTCGAAACCGc             >  2:331528/9‑251 (MQ=255)                                                                                                                                                                                                                                         gTGGATTTGTGTG‑GCAGGCGTCTGGGTGTCGCAGTTTTCGTCCTCATCCACCAGCGCTTTAACACCGCCGGCGCTGATGGAGTCCCCCCCCTCCACCAGATCGCCAACCGCCCAGGTATTGCCGGAAAAACCGAACTTGCTCTGTGCGCCACCTTTACCACCATCGTTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATcc  <  2:259673/251‑1 (MQ=255)                                                                                                                                                                                                                                          tGGATTTGTGTT‑GCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATcc  >  1:325181/1‑250 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                CGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCAGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCAGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTTGCCGTGGATTTGTGTG‑GCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCGCCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAATCGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAGACCGCCTGACGCATCC  >  NC_000913/67368‑67850

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