Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	67,900	C→G	V146V (GTG→GTC)	araA ←	L‑arabinose isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	67,900	0	C	G	97.3%	221.8 / ‑0.6	75	V146V (GTG→GTC)	araA	L‑arabinose isomerase
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base G (43/30);  total (44/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.10e-01

GCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCGCCATCGGTG‑ACCGCCACTTCACGCATGTTATCGCCAAATCGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACCACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATAC‑TGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGCCACACCACCAGACCAGCGCAACGATCGTCGTAATTCGCGTCGCG  >  NC_000913/67655‑68145                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       gcgTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTG‑ACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGc                                                                                                                                                                                                                                                   >  1:698340/1‑251 (MQ=255)   gTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTG‑ACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCAt                                                                                                                                                                                                                                                 >  1:686016/1‑251 (MQ=255)    ttAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTG‑ACCGCCACTTCACGCATGTTATCGCCAACACGGCAGCCTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGAAACGACGGCATg                                                                                                                                                                                                                                                >  2:531437/1‑251 (MQ=255)           cGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTG‑ACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTg                                                                                                                                                                                                                                          >  1:300909/1‑250 (MQ=255)                   cTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTG‑ACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATAcg                                                                                                                                                                                                                                 >  1:391452/1‑251 (MQ=255)                    tGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTG‑ACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATAcgc                                                                                                                                                                                                                                >  1:529608/1‑251 (MQ=255)                                    cTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCACCTTTATCACCATCGGTG‑ACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAACCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGGCCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCg                                                                                                                                                                                                                 >  1:478874/1‑250 (MQ=255)                                    cTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTG‑ACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCg                                                                                                                                                                                                                 >  1:581523/1‑250 (MQ=255)                                          cAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTG‑ACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGCCCGGTAACGACGGCATGTTGCTGACGCATAccgatcgggagagcgtcgtggcggg                                                                                                                                                                                                          >  2:490554/1‑227 (MQ=255)                                                    aCCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTG‑ACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCAtg                                                                                                                                                                                                >  1:504977/1‑251 (MQ=255)                                                      cGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTG‑ACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGACACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCAtgtg                                                                                                                                                                                              >  2:337602/1‑251 (MQ=255)                                                                                 cTTGATCTGAGCGGCAACTTTATCACCATCGGTG‑ACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCCTCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGTCGCATACGCGCGCCAATGAAGCCGAAATCGCGACCGCCATGTGCAGTCTGGGTCCGGTTCATAAAgacc                                                                                                                                                                    >  2:215187/1‑247 (MQ=255)                                                                                     gtcggggtgcgacatctatacacgtCGGTGCACCGCCACTCCACGCATGTTACCGCCAAACCGGCAGACTTTTAGATGTCGGGTATCCTGTTTAGAACCCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTAGCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACGCGTGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCata                                                                                                                                                                 <  2:686016/226‑1 (MQ=255)                                                                                            gcgccaCTTTAACCCCAACGGTGGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATAC‑t                                                                                                                                                        <  1:337602/246‑1 (MQ=255)                                                                                               cAACTTTATCACCATCGGTG‑ACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGCCCATATCGATAC‑TGTcc                                                                                                                                                    >  1:382516/1‑251 (MQ=255)                                                                                                    ttATCACCATCGGTG‑ACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATAC‑TGTCCCACgg                                                                                                                                               >  1:280154/1‑251 (MQ=255)                                                                                                    ttATCACCATCGGAG‑ACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAACGTCCATATCGATAC‑TGTCCCACgg                                                                                                                                               >  2:709493/1‑251 (MQ=255)                                                                                                      atcaccATCGGTG‑ACCGCCACTTCACGCATGTTAACGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCATGGCGAATCCAGGAGACGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACTCGCGCCAATGAAGCCTAAATCGCGACCGGCATGGGCGGTGTGGGTCAGGTTCATAAAGTCCATAGTGATAC‑TGGCCCCCGGcg                                                                                                                                             >  2:76282/5‑250 (MQ=255)                                                                                                        gcacctggggtaaCCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATAC‑TGTCCCACGGCAg                                                                                                                                            <  1:314996/239‑1 (MQ=255)                                                                                                           cATCGGTG‑ACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATAC‑TGTCCCACGGCAGcgcc                                                                                                                                        >  1:10155/1‑251 (MQ=255)                                                                                                              ggtgaatccgcccttacacccATGTTAGCGCCAAAACGGCAGACTTTCAGATGACGGTTATCCTGTTTAGAACCCGCCTGACGCATCCAGGAGCCGATCCGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCCATAC‑TGTCCCACGGCAGcgccgc                                                                                                                                      <  2:8114/231‑1 (MQ=255)                                                                                                              cGGTG‑ACCGCCACTTCACGAATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCCTCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTACCGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGCCCGCCATGTGCAGTCTGGTTCGGGTTCATACAGTCCATATCGATAC‑TGTCCCACGGccgcgccgcg                                                                                                                                     >  2:76445/1‑251 (MQ=255)                                                                                                                  g‑ACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGAAGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCGGGTTCAGGTTCCTAACGTCCATAACGA‑ACATGTCCAACGGCAGCGCCGCGTTgc                                                                                                                                 >  2:520812/1‑250 (MQ=255)                                                                                                                     ccGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATAC‑TGTCCCACGGCAGCGCCGCGTTGAAc                                                                                                                               >  1:560672/1‑251 (MQ=255)                                                                                                                             tCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATAC‑TGTCCCACGGCAGCGCCGCGTTGAACTGGGggtg                                                                                                                       >  2:640599/1‑251 (MQ=255)                                                                                                                              cACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGt                                                                                                                                                                                                                                               >  1:697234/1‑131 (MQ=255)                                                                                                                              cACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCAGCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGt                                                                                                                                                                                                                                               <  2:697234/131‑1 (MQ=255)                                                                                                                                  ccgtttATCGACAAAACGGAAGACTTTCAGATGATGGGTAGCCTTTTAAGAAACCGCCGGACGCAGCCAGGAGCCGATACGCTCATGGGCTGGTTTATCCGGCCAGGGACCGGTAACGAGGGCATGTTGCTGACGCATACGCGCGCCATGGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATAC‑TGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAAtt                                                                                                                  <  2:686438/247‑1 (MQ=255)                                                                                                                                    ggTTATCGCAAAACCGGCAGACTTTCAAGTGCGGGGTTTCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATCCGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATAC‑TGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTg                                                                                                                 <  2:581523/249‑1 (MQ=255)                                                                                                                                       tATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATAC‑TGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTgcagc                                                                                                             >  1:113142/1‑251 (MQ=255)                                                                                                                                       tATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGGTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGAGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATAC‑TGTCCCACGGCGGCGCCGCGTTGAACTGGGTGTGGAATTgcagc                                                                                                             >  2:171739/1‑251 (MQ=255)                                                                                                                                        aTCGCCAACACGGCAGACCTTCAGATGACGGGAACCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGAAACGACGGCATGTTGCTGAAGCGTACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCCGTCTGGTTCAGGTTCATAAAGTCCCTATCGCTAC‑TGACCCCCGGCAGCGCCGCGTTGTACTGGGTGTGGAATTgcaga                                                                                                             >  2:54400/1‑249 (MQ=255)                                                                                                                                        aTCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATAC‑TGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTgcagc                                                                                                             <  1:640599/250‑1 (MQ=255)                                                                                                                                              aacacGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCCTGGGCTTGTTTATCCTGCCCGTGACCGGTAACGACGGCATGTTGCTGACACATACGCGCGCCAATGACGCCGAACTCGCGCCCGCCATGTGCAGTCTGGTTCAGGTTCAGAACGTCCAAATCGATAC‑TGTCCCAGGGCAGCGCAGCGCTGAACTGGGCGTGGacatacagctacggtt                                                                                                      >  2:273393/5‑236 (MQ=255)                                                                                                                                                 acgcgagacttacgaaTGACGGGTACCCTGTTTAAAACCCCCCTGCCGCATCCAGGAGCCGATACGCGCATGGTCTTGTTTATCCTGCCAGTTACCGGTAACGACGGCATGTTGCTGACGCTTACGCGCGCTAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATAC‑TGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGTAACGGTttgt                                                                                                   <  2:155385/236‑1 (MQ=255)                                                                                                                                                     cAGACTTTCATATGCCGGGTATCTTGTTTAGAAACCGCCTGACGCATCCAGGGGCCGATACGCTAATGGGCTTTTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATTC‑TGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAg                                                                                               <  2:504977/251‑1 (MQ=255)                                                                                                                                                       gctggtaaggtGACGTGTAGCCGGTTAAGAAACCGCCTGCCGCAGCAAGGACGCGGTATGCTCAGGGGTTGGTTTATCCTGCCAGGGACCGGTAACGACGGCATGTTGCTGATGCGTACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATAC‑TGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTTAGCa                                                                                             <  2:55350/241‑1 (MQ=255)                                                                                                                                                           ttCAGATGACGGGTACCCTTTTGAGAACCCGCCGGACGCACCCAGTAGCTGATACGTTCAGGGGCTTGTTTATCTTGCCAGTGACCGGTAACGCCGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATAC‑TGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATgg                                                                                          <  2:539311/250‑1 (MQ=255)                                                                                                                                                                 ggACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGTATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATAC‑TGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGcc                                                                                   <  1:625593/250‑1 (MQ=255)                                                                                                                                                                     gggtttcatGTTTAGAATCCGCCTGAAGCATCCCGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTATCGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATAC‑TGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTg                                                                               <  2:523037/243‑1 (MQ=255)                                                                                                                                                                     gggTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGaa                                                                                                                                                                                                                     >  2:245116/1‑118 (MQ=255)                                                                                                                                                                     gggTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGaa                                                                                                                                                                                                                     <  1:245116/118‑1 (MQ=255)                                                                                                                                                                        tATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCACATCGATAC‑TGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATc                                                                            >  1:66116/1‑251 (MQ=255)                                                                                                                                                                         aTCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATAC‑TGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAgg                                                                                     >  1:562843/1‑241 (MQ=255)                                                                                                                                                                         aTCCTGTTTAGAAACCGCCTGACGCACCCAGGGTCCGATTCGCTCAGGGGCTTGTTTGTCCGGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATAC‑TGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAgg                                                                                     <  2:562843/241‑1 (MQ=255)                                                                                                                                                                           ccTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATAC‑TGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCcac                                                                         <  1:20964/251‑1 (MQ=255)                                                                                                                                                                                ttAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACCACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATAC‑TGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGGATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACAtttt                                                                    >  1:493274/1‑251 (MQ=255)                                                                                                                                                                                ttAGAAACCGCCTGACGCACCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATAC‑TGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACAttt                                                                     <  1:76445/250‑1 (MQ=255)                                                                                                                                                                                   aaaaCCGCCTGACGCTTCCAGAGGCCGATACGCTCATGTGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATAC‑TGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGc                                                                 <  1:390110/248‑1 (MQ=255)                                                                                                                                                                                       ccGCCTGACGCATCCCGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCTACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATAC‑TGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGa                                                                                                                     >  1:200794/1‑195 (MQ=255)                                                                                                                                                                                       ccGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCGGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAAGGAAGCCGAACTCGCTACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATAC‑TGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGa                                                                                                                     <  2:200794/195‑1 (MQ=255)                                                                                                                                                                                        cGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATAC‑TGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCgggg                                                            <  1:606337/251‑1 (MQ=255)                                                                                                                                                                                              cccaaTCCGGGAGCCGATACGCCCATGGGTTTGTTTTTCCTGCCAGTGACCGGTACCGACGGCATGTTGCTGATGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATAC‑TGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAgg                                                      <  2:300909/247‑1 (MQ=255)                                                                                                                                                                                                 cATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATAC‑TGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGgtgt                                                   <  1:531437/251‑1 (MQ=255)                                                                                                                                                                                                    ccAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGACTGGTTCAGGTTCATAAAGTCCATATCGATAC‑TGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGGGTGcg                                                >  2:702573/1‑250 (MQ=255)                                                                                                                                                                                                       ggAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGAATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATAC‑TGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGcc                                             >  1:288117/1‑251 (MQ=255)                                                                                                                                                                                                       ggAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACCACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATAC‑TGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGcc                                             <  1:277001/251‑1 (MQ=255)                                                                                                                                                                                                            cGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTACATATCGATAC‑TGTCCCACGGCAGCGCCGCGTTGAACTGGGGGTGGAATTGCAGCAACGGTTTGTTGAGCATGGACAGGCCGTTGATCCACATTTTGGCCGGGGCGAAGGTGTGCAGCCAcacc                                        >  1:26620/1‑251 (MQ=255)                                                                                                                                                                                                            cGATACGCTCATGGGCTTGTT‑ATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATAC‑TGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGCCAcacca                                       <  1:534571/251‑1 (MQ=255)                                                                                                                                                                                                             gATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTAATAAAGTCCATATCGATAC‑TGTCCCACGGCAGCGCCGCGTTGAACTGGGAGTGGAATTGCAGCAACGGTTTGTTGGGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGGAGGTGTGCAGCCAcaccc                                       >  2:596014/1‑250 (MQ=255)                                                                                                                                                                                                                  gCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGACTGGTTCAGGTTCATAAAGTCCATATCGATAC‑TGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGCGAAGGTGTGCAGCCACaccacc                                     >  1:290240/1‑248 (MQ=255)                                                                                                                                                                                                                        gggCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATAC‑TGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGCCACACCACCAGACCAgc                             <  1:106805/250‑1 (MQ=255)                                                                                                                                                                                                                          gCTTGTTTATCCTTCCAGTGACCGGTAACGACGGCATGTGGCTGACGCATACGCGCGCTAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCAGAAAg                                                                                                                                                                       <  1:451831/111‑1 (MQ=255)                                                                                                                                                                                                                          gCTTGTCTATCCTGCCAGTGACCGGTAACGCCGACATGTAGCTGTCGAATCCGCGCGCCCATGAAGCCGCACTAGCGCCCGCCATGTGCAGTCTGGTTCAGGTTCATCAAg                                                                                                                                                                       >  2:451831/1‑111 (MQ=255)                                                                                                                                                                                                                             tGTTTATCCGGCCAGTGACCGGTAACGACGGCATTTTGCTGACACATACGCGCGCCAATGAAGCGGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATAC‑TGTCCCACGGCAGCGCCGCGGTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGCCACACCACCAGACCAGCGCatc                        <  2:257020/250‑3 (MQ=255)                                                                                                                                                                                                                                   tacTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATAC‑TGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGATTTGCAGCAACGGGTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGTGCGGAGGTGTGCAGCCACACCACCAGACCATCGCATCGAtcgtc                 >  2:224368/3‑251 (MQ=255)                                                                                                                                                                                                                                    ccTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCCTACGCGCGCCCATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAACGACCATATCGATAC‑TGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGCCACACCACCAGACCAGCGCAACGAtcgtcg                >  2:296808/1‑251 (MQ=255)                                                                                                                                                                                                                                        ccAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATAC‑TGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGCCACACCACCAGACCAGCGCAACGATCGTCGTAAt            >  1:435423/1‑251 (MQ=255)                                                                                                                                                                                                                                         cAGTGTCCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGTTAC‑TGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGATGGTTGAGCATGGTCAGGCCGTTGATCCACCTTTTGGACGAGGAGAGGGTGTGCAGCCACACCACCAGACCAGCGCAACGGTCGTCGTAAt            >  2:384456/1‑250 (MQ=255)                                                                                                                                                                                                                                         cAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATAC‑TGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGCCACACCACCAGACCAGCGCAACGATCGTCGTAAtt           <  1:10377/251‑1 (MQ=255)                                                                                                                                                                                                                                         cAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATAC‑TGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGCCACACCACCAGACCAGCGCAACGATCGTCGTAAt            <  1:384456/250‑1 (MQ=255)                                                                                                                                                                                                                                          aGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATAC‑TGTCACACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAAGTGTGCAGCCACACCACCAGACCAGCGCAACGATCGTCGTAATtc          >  1:71136/1‑251 (MQ=255)                                                                                                                                                                                                                                               ccGGTAACGACGGCTTGTTGCTGACGCATACGCGCGCCAATTAAGCCGAACTCGCGACCGCCATGTGCAGTCGTGTTCAGGTTCATAAAGTCCATATCGATAC‑TGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGCCACACCACCAGCCCAGCGCAACTATCGTCGTAATtcgcgtc     <  1:206387/251‑1 (MQ=255)                                                                                                                                                                                                                                                 ggTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGCGCAGTCTGGCGCAGGATCATACAGTCCATATCGATAC‑CGTCCCCCGGCAGCGCCGCGGTGAACTGGGTGAGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGCTGATCCACATTTTGGCAGGGGAGAAGGTGTGCAGCCATACCACCAGACCAGCGCAACGATCGTCGTAATtcgcgtcg    >  1:99724/1‑250 (MQ=255)                                                                                                                                                                                                                                                  gTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATAC‑TGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGCCACACCACCAGACCAGCGCAACGATCGTCGTAATtcgcgtcgcg  >  1:551315/1‑251 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       GCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCGCCATCGGTG‑ACCGCCACTTCACGCATGTTATCGCCAAATCGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACCACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATAC‑TGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGCCACACCACCAGACCAGCGCAACGATCGTCGTAATTCGCGTCGCG  >  NC_000913/67655‑68145

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