Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	67,759	G→A	G193G (GGC→GGT)	araA ←	L‑arabinose isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	67,759	0	G	A	88.5%	133.6 / ‑1.5	52	G193G (GGC→GGT)	araA	L‑arabinose isomerase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (42/4);  minor base C (1/5);  total (43/9)
Fisher's exact test for biased strand distribution p-value = 2.70e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTTGCCGTGGATTTGTGTGGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCC‑AGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCGCCATCGGTG‑ACCGCCACTTCACGCATGTTATCGCCAAATCGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACCACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGG  >  NC_000913/67509‑68004                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                       tggtgAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTTGCCGTGGATTTGCGTGGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATAGACCAGCGCGTTAACATCGCCGTCGCTGATGGGGTTCACCACCTGCACCAGATCGCCAACCGCCC‑AGGTATTGACGGGGAAACCGAACTTGAGCTGTGCGGCAACTTTATCg                                                                                                                                                                                                                                                        >  2:277001/1‑251 (MQ=255)      aaCGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCC‑AGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATc                                                                                                                                                                                                                                                   <  1:498733/251‑1 (MQ=255)       aCGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCC‑AGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCg                                                                                                                                                                                                                                                  <  1:164334/251‑1 (MQ=255)            tGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCC‑AGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTG‑Ac                                                                                                                                                                                                                                            >  1:55350/1‑251 (MQ=255)             ggAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCC‑AGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTG‑Acc                                                                                                                                                                                                                                           <  1:520251/251‑1 (MQ=255)                    cccTTTTTCCAGGAAACGCTTCACCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCACCCGCCC‑AGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCACCGGTG‑ACCGCCACtt                                                                                                                                                                                                                                    <  1:600691/250‑1 (MQ=255)                        tGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCC‑AGGTATTGACGGAGAAACCGAACTTGATCTGTGCTGCAACTTTATCAACATCGGGG‑ACCGCCACTTCACg                                                                                                                                                                                                                                >  2:197128/1‑251 (MQ=255)                        tGTTACAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCCGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCACGTTAACATCGCCGTCGCTGCTGGAGTTCACCACCTGCACCAGATCGCCACCCGCCC‑AGGTATTGGCGGAGAAACCGAACTTTATCTGTGCGGCACCTTTATAACAAACGATG‑ACCGCCACTTCACg                                                                                                                                                                                                                                >  2:106805/1‑251 (MQ=255)                         gTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCTGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGCTCGCCAACCGCCC‑AGGTATTGACGGAGAAACCGACCTTGATCTGTGCGGCAACTTTATCCACATCGATG‑ACCGCCACttaccgc                                                                                                                                                                                                                               >  2:261346/1‑246 (MQ=255)                             cAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCC‑AGGAGTTGACGGAGAAACCGCACTTGATCTGTGCGGCAACTTTCTCACCATCGGTG‑ACCGCCGCTTCACGCATGt                                                                                                                                                                                                                           >  2:265365/1‑251 (MQ=255)                                 aaccgcttcaccccaagtcaaatccgcgCCGTTTCCAGCATTTTCTGTCTTTTTTCGCCGTGGATTTTTGTTGCAGCTGTCAGGGTGTAGCAGTTTTCGTCCTCACCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTCCCCCACCTGCCCCAGACCCCCACCCGCCA‑AGGTATTGACGGAGAAACCTAACTTGATCTGGGCGGCAACTTTATCCCCATCGGTG‑ACCGCCACTTCACGCATGTTATc                                                                                                                                                                                                                       <  1:421967/227‑1 (MQ=255)                                 aaaCGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCC‑AGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTG‑ACCGCCACTTCACGCATGTTATc                                                                                                                                                                                                                       >  1:416179/1‑251 (MQ=255)                                         cATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCC‑AGGTATTGACGGAGAAACCGCACTTGATCTGTGCGGCAACTTTATCACCATCGGTG‑ACCGCCACTTCACGCATGTTATCGCCaaaac                                                                                                                                                                                                               >  2:270389/1‑249 (MQ=255)                                              ccAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCC‑AGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTG‑ACCGCCACTTCACGCATGTTATCGCCAAAACGGc                                                                                                                                                                                                            >  1:109411/1‑249 (MQ=255)                                                     aaTACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGAAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCC‑AGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTG‑CCCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTc                                                                                                                                                                                                    >  2:701615/1‑250 (MQ=255)                                                         cgcgCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCC‑AGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTG‑ACCGCCACTTCACGCATGTTATCGCCaaaac                                                                                                                                                                                                               >  1:84604/1‑233 (MQ=255)                                                                         cGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTAACCACCTGCACCAGATAGCCAACCGCCC‑AGGTATTGCCGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTG‑ACAGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGttt                                                                                                                                                                                >  1:468504/1‑250 (MQ=255)                                                                             ttgttgttttcccccggggATTTGGTTGGCAGGCGTCTTGGTGTAGCGCCTTCCGTCCCAACCGACCAGCGCGTTAACATCGCGGTCGCGGAGGGATTTCCCCCCCTGCCCAAGACCCCCACCCCCCC‑AGGTATGGACGGAAAAACCGATCTTTATCTGTGGGGAAATTTTACCCCCATCGGTG‑ACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTagaa                                                                                                                                                                            <  1:197128/233‑2 (MQ=255)                                                                                             tGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCC‑AGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTG‑ACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATcc                                                                                                                                                            >  1:155385/1‑250 (MQ=255)                                                                                                aTTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCC‑AGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTG‑ACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACg                                                                                                                                                                 >  1:399680/1‑242 (MQ=255)                                                                                                        tgCAGGCGTCATGTTGGAGCAGCTTTCGTCCTCTGCGCCCAGCGCTTTAACACCGCCTTCGCTGGTGGAGTCCCCCACCTGCACCAGATCGCCACCCGCCC‑AGTTATTGAGGGAAAACCCGATTTTGATCTGGGCGACAATTTTACCCCCACCGTGG‑ACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATa                                                                                                                                                 <  1:701615/249‑1 (MQ=255)                                                                                                         gCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGAGTTAACATCGCCGTCGCTGATGGCGTTCACCACCTGCACCAGATCGCCAACCGCCC‑AGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTG‑ACCGCCACTTCCCGCATGTTATCGCCAAAACGGCCGACTTTCAGATGACGGGTATCCTATTTAGAAACCGCCTGACGCATCCAGGAGCCGATa                                                                                                                                                 >  2:82696/1‑249 (MQ=255)                                                                                                          cAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCC‑AGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTG‑ACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGc                                                                                                                                              >  1:686438/1‑251 (MQ=255)                                                                                                              cGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCC‑AGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTG‑ACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGCATCCTGTTTCGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCAt                                                                                                                                          >  2:559144/1‑251 (MQ=255)                                                                                                               gTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCC‑AGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTG‑ACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATg                                                                                                                                         >  1:523037/1‑251 (MQ=255)                                                                                                                cctgggTGTAGCAGTTTTCTTACTCACCGCCCAGCGCGTTACCACCGCCGTCGTGGAGGGTGTCCCCCCCCGGCCCAAGA‑CCCCAACCCCCCCAGTTATTGGCGGAAAAACCAAATTTGATTGGTGCGGCATTTTTACCCCCACCGTTG‑ACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATgg                                                                                                                                        <  1:488486/247‑1 (MQ=255)                                                                                                                  aTGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCC‑AGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTG‑ACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACAGCCTGACGCATCCAGGAGCCGATACGCTCATGGGc                                                                                                                                      >  2:10377/1‑251 (MQ=255)                                                                                                                         agcagcTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCC‑AGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTG‑ACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAACCCGCCTGACGCCTCCAGGAGCCGATACGCTCATGGGCTTGttat                                                                                                                               >  2:534571/1‑249 (MQ=255)                                                                                                                            agcTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCC‑AGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTG‑ACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATcc                                                                                                                            >  1:625651/1‑251 (MQ=255)                                                                                                                                  cGTACTCATAAACCAACGCGTTAACATCGCCTTCGCTGCTGGAGTTCACCACCTGCACAAGAACGCCAACCGCCC‑AGGTATTGACGGAGAAACCGCACTTGATCTGAGCGGCAACTTTCTCACCATCGGCG‑ACCGCCACTTCACGCAGGTGATCGCCAGAACGGCAGCCTTTCAGATGACGGGTATCCTGTTTAGAACCCGCCTGACGCATCCAGGAGCCGATACGCACAGGGGATTGTTTagccgccc                                                                                                                        >  2:105493/1‑242 (MQ=255)                                                                                                                                     aCTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCC‑AGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTG‑ACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTg                                                                                                                    >  1:136147/1‑250 (MQ=255)                                                                                                                                       tCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCC‑AGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTG‑ACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGAcc                                                                                                                 >  1:532561/1‑251 (MQ=255)                                                                                                                                              ccAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCC‑AGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTG‑ACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAAcg                                                                                                          >  1:257020/1‑250 (MQ=255)                                                                                                                                                  cgcgTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCC‑AGGTATTGACGGAGCAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTG‑ACCGCCACTTCACGCATGTTATCGCCAAACCGGCAGACTTTCAGATGCCGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCAGATCCGCTCATGGGCTTGTTTATCCTGCCAGTGGCCGGTaccgacgg                                                                                                      >  2:342914/1‑244 (MQ=255)                                                                                                                                                   gcgTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCC‑AGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTG‑ACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGc                                                                                                     >  1:698340/1‑251 (MQ=255)                                                                                                                                                     gTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCC‑AGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTG‑ACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCAt                                                                                                   >  1:686016/1‑251 (MQ=255)                                                                                                                                                      ttAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCC‑AGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTG‑ACCGCCACTTCACGCATGTTATCGCCAACACGGCAGCCTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGAAACGACGGCATg                                                                                                  >  2:531437/1‑251 (MQ=255)                                                                                                                                                        aaCATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCC‑AGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTG‑ACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCg                                                                                                                >  1:530722/1‑235 (MQ=255)                                                                                                                                                             cGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCC‑AGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTG‑ACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTg                                                                                            >  1:300909/1‑250 (MQ=255)                                                                                                                                                                     cTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCC‑AGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTG‑ACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATAcg                                                                                   >  1:391452/1‑251 (MQ=255)                                                                                                                                                                      tGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCC‑AGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTG‑ACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATAcgc                                                                                  >  1:529608/1‑251 (MQ=255)                                                                                                                                                                              ttCACCACCTGCACCAGATCGCCAACCGCCC‑AGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTG‑ACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGCTGACGGGTATCCTGTTTAGAAACCGCCTGACGCa                                                                                                                                                               >  2:353683/1‑166 (MQ=255)                                                                                                                                                                                      cTGCACCAGATCGCCAACCGCCC‑AGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCACCTTTATCACCATCGGTG‑ACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAACCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGGCCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCg                                                                   >  1:478874/1‑250 (MQ=255)                                                                                                                                                                                      cTGCACCAGATCGCCAACCGCCC‑AGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTG‑ACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCg                                                                   >  1:581523/1‑250 (MQ=255)                                                                                                                                                                                            cAGATCGCCAACCGCCC‑AGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTG‑ACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGCCCGGTAACGACGGCATGTTGCTGACGCATAccgatcgggagagcgtcgtggcggg                                                            >  2:490554/1‑227 (MQ=255)                                                                                                                                                                                                      aCCGCCC‑AGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTG‑ACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCAtg                                                  >  1:504977/1‑251 (MQ=255)                                                                                                                                                                                                        cGCCC‑AGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTG‑ACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGACACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCAtgtg                                                >  2:337602/1‑251 (MQ=255)                                                                                                                                                                                                                                    cTTGATCTGAGCGGCAACTTTATCACCATCGGTG‑ACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCCTCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGTCGCATACGCGCGCCAATGAAGCCGAAATCGCGACCGCCATGTGCAGTCTGGGTCCGGTTCATAAAgacc                      >  2:215187/1‑247 (MQ=255)                                                                                                                                                                                                                                               gcgccaCTTTAACCCCAACGGTGGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACt           <  1:337602/246‑1 (MQ=255)                                                                                                                                                                                                                                               cGGCAACTTTATCACCATCGGTG‑ACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCAt                                                                                                                                          >  2:364753/1‑123 (MQ=255)                                                                                                                                                                                                                                                  cAACTTTATCACCATCGGTG‑ACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGCCCATATCGATACTGTcc       >  1:382516/1‑251 (MQ=255)                                                                                                                                                                                                                                                       ttATCACCATCGGTG‑ACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACgg  >  1:280154/1‑251 (MQ=255)                                                                                                                                                                                                                                                       ttATCACCATCGGAG‑ACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAACGTCCATATCGATACTGTCCCACgg  >  2:709493/1‑251 (MQ=255)                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                       TGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTTGCCGTGGATTTGTGTGGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCC‑AGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCGCCATCGGTG‑ACCGCCACTTCACGCATGTTATCGCCAAATCGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACCACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGG  >  NC_000913/67509‑68004

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