Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,245,696	C→A	A242S (GCC→TCC)	malE ←	maltose ABC transporter periplasmic binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,245,696	0	C	A	98.6%	235.8 / NA	71	A242S (GCC→TCC)	malE	maltose ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base A (31/39);  minor base T (0/1);  total (31/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

CTTTATTAACCGCTTCCAGACCTTCATCAGTCAGCAGATAGTTTTCGAGGAACTCTTTCGCCAGCTCTTTGTTCGGACTGGCGGCGTTAATACCTGCGCTCAGCACGCCAACGAACGGTTTGGATGGTTGACCCTTGAAGGTCGGCAGTACCGTTACACCATAATTCACTTTGCTGGTGTCGATGTTGGACCATGCCCACGGGCCGTTGATGGTCATCGCTGTTTCGCCTTTATTAAAGGCAGCTTCTGCGATGGAGTAATCGGTGTCTGCATTCATGTGTTTGTTTTTAATCAGGTCAACCAGGAAGGTCAGACCCGCTTTCGCGCCAGCGTTATCCACGCCCACGTCTTTAATGTCGTACTTGCCGTTTTCATACTTGAACGCATAACCCCCGTCAGCAGCAATCAGCGGCCAGGTGAAGTACGGTTCTTGCAGGTTGAACATCAGCGCGCTCTTACCTTTCGCTTTCAGTTCTTTATCCAGCGC  >  NC_000913/4245456‑4245942                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                       cTTTATTAACCGCTTCCAGACCTTCATCAGTCAGCAGATAGTTTTCGAGGAACTCTTTCGCCAGCTCTTTGTTCGGACTGGCGGCGTTAATACCTGCGCTCAGCACGCCAACGAACGGTTTGGCTGGTTGACCCTTGAAGGGCGGCAGTACCGTTACACCATAATTCACTTTGCTGGTGTCGATGTTGGACCATGCCCACGGGCCGTTGATGGTCATCGCTGTTTCGCCTTTATTAAAGGAAGCTTCTGCg                                                                                                                                                                                                                                              >  2:23848/1‑251 (MQ=255)           cGCTTCCAGACCTTCATCAGTCAGCAGATAGTTTTCGAGGAACTCTTTCGCCAGCTCTTTGTTCGGACTGGCGGCGTTAATACCTGCGCTCAGCACGCCAACGAACGGTTTGGATGGTTGACCCTTGAAGGTCGGCAGTACCGTTACACCATAATTCACTTTGCTGGTGTCGATGTTGGACCATGCCCACGGGCCGTTGATGGTCATCGCTGTTTCGCCTTTATTAAAGGAAGCTTCTGCGATGGAGTAAt                                                                                                                                                                                                                                    >  1:21702/1‑251 (MQ=255)                  aGACCTTCATCAGTCAGCAGATAGTTTTCGAGGAACTCTTTCGCCAGCTCTTTGTTCGGACTGGCGGCGTTAATACCTGCGCTCAGCACGCCAACGAACGGTTTGGATGGTTGACCCTTGAAGGTCGGCAGTACCGTTACACCATAATTCACTTTGCTGGTGTCGATGTTGGGCCATGCCCACGGGCCGTTGATGGTCATCGCTGTTTCGACTTTATTAAAGGAAGCTTCTGCGATGGAGTAATCGGGGTc                                                                                                                                                                                                                             >  2:145396/1‑251 (MQ=255)                   gACCTTCATCAGTCAGCAGATAGTTTTCGAGGAACTCTTTCGCCAGCTCTTTGTTCGGACTGGCGGCGTTAATTCCTGCGCTCAGCACGCCAACGAACGGTTTGGATGGTTGACCCTTGAAGGTCGGCAGTACCGTTACACCATAATTCACTTTGCTGGTGGCGATGTTGGACCATGCCCACGGGCCGTTGATGGTCATCGCTGTTTCGCCTTTATTAAAGGAAGCTTCTGCGATGGAGTAATCGGTGTCt                                                                                                                                                                                                                            <  1:429641/251‑1 (MQ=255)                         catcaGTCAGCAGATAGTTTTCGAGGAACTCTTTCGCCAGCTCTTTGTTCGGACTGGCGGCGTTAATACCTGCGCTCAGCACGCCAACGAACGGTTTGGATGGTTGACCCTTGAAGGTCGGCAGTACCGTTACACCATAATTCACTTTGCTGGTGTCGATGTTGGACCATGCCCACGGGCCGTTGATGGTCATCGCTGTTTCGCCTTTATTAAAGGAAGCTTCTGCGATGGAGTAATCGGTGTCTGCAtta                                                                                                                                                                                                                      <  1:486931/251‑2 (MQ=255)                              gtcagCAGATAGTTTTCGAGGAACTCTTTCGCCAGCTCTTTGTTCGGACTGGCGGCGTTAATACCTGCGCTCAGCACGCCAACGAACGGTTTGGATGGTTGACCCTTGAAGGTCGGCAGTACCGTTACACCATAATTCACTTTGCTGGTGTCGATGTTGGACCATGCCCACGGGCCGTTGATGGTCATCGCTGTTTCGCCTTTATTAAAGGAAGCTTCTGCGATGGAGTAATCGGTGTCTGCATTCAtgtg                                                                                                                                                                                                                 <  1:472243/251‑1 (MQ=255)                              gtcagCAGATAGTTTTCGAGGAACTCTTTCGCCAGCTCTTTGTTCGGACTGGCGGCGTTAATACCTGCGCTCAGCACGCCAACGAACGGTTTGGATGGTTGACCCTTGAAGGTCGGCAGTACCGTTACACCATAATTCACTTTGCTGGTGTCGATGTTGGACCATGCCCACGGGCCGTTGATGGTCATCGCTGTTTCGCCTTTATTAAAGGAAGCTTCTGCGATGGAGTAATCGGTGTCTGCATTCAtgtg                                                                                                                                                                                                                 <  1:604453/251‑1 (MQ=255)                               tcagCCGATAGTTTTCGAGGAACTCTTTCGCCAGCTCTTCGTTCGGACTGGCGGCGTTAATACCTGCGCTCAGCACGCCAACGAACGGTTTGGCTGGTTGAACCTTGAAGGTCGGCCGTACCGTTACACCATAATTCACTTTGCTGGTGTCGCTGTTGGACCATGCCCACGGGCCGTTGATGGTCATCGCTGTTTCGCCTTTATTAAAGGAAGCTTCTGCGATGGAGTAATCGGTGTCTGCATTCAttggt                                                                                                                                                                                                                >  2:702480/1‑247 (MQ=255)                                    agatcgatttCGAGGAACTCTTTCGCCAGCTCTTTGGTCGGACTGGCGGCGTTAATACCTGCGCTCAGCACGCCCACGAACGGTTTGGATGGTTGACCCTTGAAGGTCGGCAGTACCGTTACACCATAATTCACTTTGCTGGTGTCGTTTTTGTACCATGCCCACGGGCCGTTGATGGTCATAGCTGTTTAGCCTTTATTAAAGGAAGCTTCTGCGAATGAGTAATCGGGGTCTGCATTCATGtggttgtt                                                                                                                                                                                                           >  2:555446/8‑251 (MQ=255)                                      aTAGTTTTCGAGGAACTCTTTCGCCAGCTCTTTGTTCGGACTGGCGGCGTTAATACCTGCGCTCAGCACGCCAACGAACGGTTTGGATGGTTGACCCTTGAAGGTCGGCAGTACCGTTACACCATAATTCACTTTGCTGGTGTCGATGTTGGACCATGCCCACGGGCCGTTGATGGTCATCGCTGTTTCGCCTTTATTAAAGGAAGCTTCTGCTATGGAGTAATCGGTGTCTGCATTCATGTGTTTGtttt                                                                                                                                                                                                         <  2:520331/251‑1 (MQ=255)                                       aaGTTTTCGAGGAACTCTTTCGCCAGCTCTTTGTTCGGACTGGCGGCGTTAATACCTGCGCTCAGCACGCCAACGAACGGTTTGGATGGTTGACCCTTGAAGGTCGGCAGTACCGTTACACCATAATTCACTTTGCTGGTGTCGATGTTGGACCATGCCCACGTGCCGTTGATGGTCATCGCTGTTTCGCCTTTATTAAAGGAAGCTTCTGCGATGGAGTAATCGGTGTCTGCATTCATGTGTTTGttttt                                                                                                                                                                                                        <  1:95449/250‑1 (MQ=255)                                           gttCGAGGAACTCTGTCGCCAGCTCTTTGTTCGGTCTGGCGGCGTTAATACCTGCGCTCAGCACGCCCACGAACGGTTTGGATGGTTGACCCTTGAAGGTCGGCAGTACCGTTACCCCATAATTCCCGGTGCTGGTGTCGATGTTGGACCATGCCCACGGGCCGTTGATGGGCATCGCTGTTTCGCCTTTATTAAAGGTATCTTCTGCGATGGAGTAATCGGTGTCTGCATTCATGTGTTTGTTTTTAATc                                                                                                                                                                                                    <  2:260240/250‑1 (MQ=255)                                            ttCGAGGAACTCTTTCGCCAGCTCTTTGTTCGGACTGGCGGCGTTAATACCTGCGCTCAGCACGCCAACGAACGGTTTGGATGGTTGACCCTTGAAGGTCGGCAGTACCGTTACACCATAATTCACTTTGCTGGTGTCGATGTTGGACCATGCCCACGGGCCGTTGATGGTCATCGCTGTTTCGCCTTTATTAAAGGAAGCTTCTGCGATGGAGTAATCGGTGTCTGCATTCATGTGTTTGTTTTTAATc                                                                                                                                                                                                    >  1:8786/1‑250 (MQ=255)                                             tCGAGGAACTCTTTCGCCAGCTCTTTGTTCGGACTGGCGGCGTTAATACCTGCGCTCAGCACGCCAACGAACGGTTTGGATGGTTGACCCTTGAAGGTCGGCAGTACCGTTACACCATAATTCACTTTGCTGGTGTCGATGTTGGACCATGCCCACGGGCCGTTGATGGTCATCGCTGTTTCGCCTTTATTAAAGGAAGCTTCTGCGATGGAGTAATCGGTGTCTGCATTCATGTGTTTGTTTTTAATc                                                                                                                                                                                                    >  1:488717/1‑249 (MQ=255)                                              cGAGGAACTCTTTCGCCAGCTCTTTGTTCGGACTGGCGGCGTTAATTCCTGCGCTCAGCACGCCAACGAACGGTTTGGATTGTTGACCCTTGAAGGTCGGCAGTACCGTTACACCATAATTCACTTTGCTGGTGGCGATGTTGGACCATGCCCACGGGCCGTTGATGGTCATCGCTGTTTCGCCTTTATTAAAGGAAGCTTCTGCGATGGAGTAATCGGTGTCTGCATTCATGTGTGTGTTTTTAATCAgg                                                                                                                                                                                                 <  2:438514/251‑1 (MQ=255)                                                aGGAACTCTTTCGCCAGCTCTTTGTTCGGTCTGGCGGCGTTAATACCAGCGCTCAGCACGCCAACGACCGGTTTGGGTGGTTGACCCTTGAAGGTCGGCAGTACCGTTACACCATAATTCACTTTGCTGGTGTCGATGTTGGACCATGCCCACGGGCCGTTGATGGCCATCGCTGTTTCGCCTTTATTACAGGAAGCTTCTGCGATGGAGTCCTCGTTGTCTGCTTTTATGTGTTTGTTTTTAATCAGGtg                                                                                                                                                                                               >  2:236675/1‑250 (MQ=255)                                                 ggAACTCTTTCGCCAGCTCTTTGTTCGGACTGGCGGCGTTAATACCTGCGCTCAGCACGCCAACGAACGGTTTGGATGGTTGACCCTTGAAGGTCGGCAGTACCGTTACACCATAATTCACTTTGCTGGTGTCGATGTTGGACCATGCCCACGGGCCGTTGATGGTCATCGCTGTTTCGCCTTTATTAAAGGAAGCTTCTGCGATGGAGTAATCGGTGTCTGCATTCATGTGTTTGTTTTTAATCAGGTCa                                                                                                                                                                                              <  1:574845/251‑1 (MQ=255)                                                 ggAACTCTTTCGCCAGCTCTTTGTTCGGACTGGCGGCGTTAATACCTGCGCTCAGCACGCCAACGAACGGTTTGGATGGTTGACCCTGGAAGGTCGGCAGTACCGTTACACCATAATTCACTTTGCTGGTGTCGATGTTGGACCATGCCCACGGGCCGTTGATGGTCATCGCTGTTTCGCCTTTATTAAAGGAAGCTTCTGCGATGGAGTAATCGGTGTCTGCATTCATGTGTTTGTTTTTAATCAGGTCa                                                                                                                                                                                              <  1:102009/251‑1 (MQ=255)                                                  gaagtctTTCGCCAGCTCTTTGTTCGGACTGGCGGCGTTAATACCTGCGCTCAGCACGCCAACGAACGGTTTGGCTGGTTGACCCTTGAAGGTCGGCAGTACCGTTACACCATAATTCACTTTGCTGGTGTCGAGGTTGGACCATGCCCACGGGCCGTTGATGGTCATCGCTGTTTCGCCTTTATTAAAGGAAGCTTCTGCGATGGAGTAATCGCTGTCTGCCTTTATGTGTTTGTTTTTAATCAGGTCa                                                                                                                                                                                              >  2:396086/5‑250 (MQ=255)                                                   aaCTCTTTCGCCAGCTCTTTGTTCGGACTGGCGGCGTTAATACCTGCGCTCAGCACGCCAACGAACGGTTTGGATGGTTGACCCTTGAAGGTCGGCAGTACCGTTACACCATAATTCACTTTGCTGGTTTCGATGTTGGACCATGCCCACGGGCCGTTGATGGTCATCGCTGTTTCGCCTTTATTAAAGGAAGCTTCTGCGATGGAGTAATCGGTGTCTGCATTCATGTGTTTGTTTTTAATCAGGTCAAc                                                                                                                                                                                            <  1:655818/251‑1 (MQ=255)                                                   aaCTCTTTCGCCAGCTCTTTGTTCGGACTGGCGGCGTTAATACCTGCGCTCAGCACGCCAACGAACGGTTTGGATGGTTGACCCTTGAAGGTCGGCAGTACCGTTACACCATAATGCACTTTGCTGGTGTCGATGTTGGACCATGCCCACGGGCCGTTGATGGTCATCGCTTTTTCGCCTTTATTAAAGGAAGCTTCTGCGATGGAGTAATCGGTGTCTGCATTCATGTGTTTGTTTTTAATCAGGTCAAc                                                                                                                                                                                            <  2:273486/251‑1 (MQ=255)                                                           cGCCAGCTCTTTGTTCGGACTGGCGGCGTTAATACCTGCGCTCAGCACGCCAACGAACGGTTTGGATGGTTGACCCTTGAAGGTCGGCAGTACCGTTACACCATAATTCACTTTGCTGGTGTCGATGTTGGACCATGCCCACGGGCCGTTGATGGTCATCGCTGTTTCGCCTTTATTAAAGGAAGCTTCTGCGATGGAGTAATCGGTGTCTGCATTCATGTGTTTGTTTTTAATCAGGTCAACCAGGAAg                                                                                                                                                                                     <  1:213897/250‑1 (MQ=255)                                                            gCCAGCTCTTTGTTCGGACTGGCGGCGTTAATACCTGCGCTCAGCACGCCAACGAACGGTTTGGATGGTTGACCCTTGAAGGTCGGCAGTACCGTTACACCATAATTCACTTTGCTGGTGTCGATGTTGGACCATGCCCACGGGCCGTTGATGGTCATCGCTGTTTCGCCTTTATTAAAGGAAGCTTCTGCGATGGAGTAATCGGTGTCTGCATTCATGTGTTTGTTTTTAATCAGGTCAACCAGGAAGGt                                                                                                                                                                                   >  1:327798/1‑251 (MQ=255)                                                             ccAGCTCTTTGTTCGGACTGGCGGCGTTAATCCCTGCGCTCAGCACGCCAACGAACGGTTTGGATGGTTGACCCTTGAAGGTCGGCAGTACCGTTACACCATAATTCACTTTGCTGGTGTCGATGTTGGACCATGCCCACGGGCCGTTGATGGTCATCGCTGTTTCGCCTTTATTAAAGGAAGCTTCTGCGATGGAGTAATCGGTGTCTTCATTCATGTGTTTGTTTTTAATCAGGTCAACCATGAAGGTc                                                                                                                                                                                  <  2:297354/251‑1 (MQ=255)                                                             ccAGCTCTTTGTTCGGACTGGCGGCGTTAATACCTGCGCTCAGCACGCCAACGAACGGTTTGGATGGTTGACCCTTGAAGGTCGGCAGTACCGTTACACCATAATTCACTTTGCTGGTGTCGATGTTGGACCATGCCCACGGGCCGTTGATGGTCATCGCTGTTTCGCCTTTATTAAAGGAAGCTTCTGCGATGGAGTAATCGGTGTCTGCATTCATGTGTTTGTTTTTAATCAGGTCAACCAGGAAGGTc                                                                                                                                                                                  >  2:191987/1‑251 (MQ=255)                                                                   ctTTGTTCGTATTGGCGTCGGTAAACCCTGCGCTCATCCCTCCCACGAACGGTTTGTATTGTTTACCCTTGAATGTCGGCAGTACCGTTACACCATAATTCACTTTGCTTGTGTCGATGTTGGATCATGCCCCCGGGCCGTTGATGGTCATCGCTGTTTCGCCTTTATTAAAGGAAGCTTCTGCGATGGATTAATCGGTGTCTGCATTCAGGTGTTTGTTTTTAATCAGGTCAACCAGGAAgttccgac                                                                                                                                                                              <  2:223147/249‑8 (MQ=255)                                                                    tttGTTCGGACTGGCGGCGTTAATACCTGCGCTCAGCACGCCAACGAACGGGTTGGATGGTTGACCCTTGAAGGTCGGCAGGACCGTTACACCCTAATTCACTTTGCTGGTGTCGATGTTGGACCATGCCCACGGGCCGTTGATGGTCATCGCTGTTTCGCCTTTATTAAAGGAAGCTTCTGCGATGGAg                                                                                                                                                                                                                                        >  2:193313/1‑190 (MQ=255)                                                                    tttGTTCGGACTGGCGGCGTTAATACCCGCGCTCAGCACGCCAACGAACGGTTTGGATGGTTGACCCTTGAAGGTCGGCAGTACCGTTACACCATAATTCACTTTGCTGGTGTCGATGTTGGACCATGCCCGCGGGCCGTTGATGGTCATCGCTGTTTCGCCTTTATTAAAGGAAGCTTCTGCGATGGAg                                                                                                                                                                                                                                        <  1:193313/190‑1 (MQ=255)                                                                       gTTCGGACTGGCGGCGTTAATACCTGCGCTCAGCACGCCAACGAACGGTTTGGATGGTTGACCCTTGAAGGTCGGCAGTACCGTTACACCATAATTCACTTTGCTGGTGTCGATGTTGGACCATGCCCACGGGCCGTTGATGGTCATCGCTGTTTCGCCTTTATTAAAGGAAGCTTCTGCGATGGAGTAATCGGTGTCTGCATTCATGTGTTTGTTTTTAATCAGGTCAACCAGGAAGGTCAGACCCGCt                                                                                                                                                                         <  1:52428/250‑1 (MQ=255)                                                                                     cgTTAATACCTGCGCTCAGCACGCCCACGAACGGTTTGGATGGTTGACCCTTGAAGGTTGGCAGTACCGTTACACCATAATTCTCTTTGCTGGTGTCGATGTTGGACCATGCCCACGGGCCGTTGATGGTCATCGCTGTTTCGCCTTTATTAAAGGAAGCTTCTGCGATGGAGTAATCGGTGTCTGCATTCATGTGTTTGTTTTTAATCAGGTCAACCAGGAAGGTCAGACCCGCTTTCGCGCCAGCGTTa                                                                                                                                                          <  2:536140/251‑1 (MQ=255)                                                                                      gTTAATACCTGCGCTCAGCACGCCAACGAACGGTTTGGATGGTTGACCCTTGAAGGTCGGCAGTACCGTTACACCATAATTCACTTTGCTGGTGTCGATGTTGGACCATGCCCACGGGCCGTTGATGGTCATCGCTGTTTCGCCTTTATTAAAGGAAGCTTCTGCGATGGAGTAATCGGTGTCTGCATTCATGTGTTTGTTTTTAATCAGGTCAACCAGGAAGGTCAGACCCGCTTTCGCGCCAGCGTTAt                                                                                                                                                         >  2:156717/1‑251 (MQ=255)                                                                                           tACCTGCGCTCAGCACGCCAACGAACGGTTTGGATGGTTGACCCTTGAAGGTCGGCAGTACCGTTACACCATAATTCACTTTGCTGGTGTCGATGTTGGACCATGCCCACGGGCCGTTGATGGTCATCGCTGTTTCGCCTTTATTAAAGGAAGCTTCTGCGATGGAGTAATCGGTGTCTGCATTCCTGTGTTTGTTTTTAATCAGGTCAACCAGGAAGGTCAGACCCGCTTTCGCGCCAGCGTTATCCAcc                                                                                                                                                    >  2:225842/1‑250 (MQ=255)                                                                                           tACCTGCGCTCAGCACGCCAACGAACGGTTTGGATGGTTGACCCTTGAAGGTCGGCAGTACCGTTACACCATAATTCACTTTGCTGGTGTCGATGTTGGACCATGCCCACGGGCCGTTGATGGTCATCGCTGTTTCGCCTTTATTAAAGGAAGCTTCTGCGATGGAGTAATCGGTGTCTGCATTCATGtgtt                                                                                                                                                                                                               >  2:197730/1‑192 (MQ=255)                                                                                                gcgcTCAGCACGCCAACGAACGGTTTGGATTGTTGACCCTTGAAGGTCGGCAGTACCGTTACACCATAATTCACTTTGCTGGTGTCGATGTTGGACCATGCCCACGGGCCGTTGATGGTCATCGCTGTTTCGCCTTTATTAAAGGAAGCTTCTGCGATGGAGTAATCGGTGTCTGCATTCATGTGTTTGTTTTTAATCAGGTCAACCAGGAAGGTCAGACCCGCTTTCGCGCCAGCGTTATCCACGCCCAc                                                                                                                                               <  1:145396/251‑1 (MQ=255)                                                                                                 cgcTCAGCACGCCAACGAACGGTTTGGATGGTTGACCCTTGAAGGTCGGCAGTACCGTTACACCATAATTCACTTTGCTGGTGTCGATGTTGGACCATGCCCACGGGCCGTTGATGGTCATCGCTGTTTCGCCTTTATTAAAGGAAGCTTCTGCGATGGAGTAATCGGTGTCTGCATTCATGTGTTTGTTTTTAATCAGGTCAACCAGGAAGGTCAGACCCGCTTTCGCGCCAGCGTTATCCACGCCCACg                                                                                                                                              <  1:500973/251‑1 (MQ=255)                                                                                                   cTCAGCACGCCAACGAACGGTTTGGATGGTTGACCCTTGAAGGTCGGCAGTACCGTTACACCATAATTCACTTTGCTGGTGTCGATGTTGGACCATGCCCACGGGCCGTTGATGGTCATCGCTGTTTCGCCTTTATTAAAGGAAGCTTCTGCGATGGAGTAATCGGTGTCTGCATTCATGTGTTTGTTTTTAATCAGGTCAACCAGGAAGGTCAGACCCGCTTTCGCGCCAGCGTTATCCACGCCCACGTc                                                                                                                                            >  1:438710/1‑251 (MQ=255)                                                                                                    tCAGCACGCCAACGAACGGTTTGGATGGTTGACCCTTGAAGGTCGGCAGTACCGTTACACCATAATTCACTTTGCTGGTGTCGATGTTGGACCATGCCCACGGGCCGTTGATGGTCATCGCTGTTTCGCCTTTATTAAAGGAAGCTTCTGCGATGGAGTAATCGGTGTCTGCATTCATGTGTTTGTTTTTAATCAGGTCAACCAGGAAGGTCAGACCCGCTTTCGCGCCAGCGTTATCCACGCCCACGTCt                                                                                                                                           >  1:163621/1‑251 (MQ=255)                                                                                                     cAGCACGCCAACGAACGGTTTGGATGGTTGACCCTTGAAGGTCGGCAGTACCGTTACACCATAATTCACTTTGCTGGTGTCGATGTTGGACCATGCCCACGGGCCGTTGATGGTCATCGCTGTTTCGCCTTTATTAAAGGAAGCTTCTGCGATGGAGTAATCGGTGTCTGCATTCATGTGTTTGTTTTTAATCAGGTCAACCAGGAAGGTCAGACCCGCTTTCGCGCCAGCGTTATCCACGCCCACGTCt                                                                                                                                           <  1:595342/250‑1 (MQ=255)                                                                                                       gCACGCCAACGAACGGTTTTGATGGTTGACCCTTGAAGGTCGGCAGTACCGTTACACCATAATTCACTTTGCTGGTGTCGATGTTGGACCATGCCCACGGGCCGTTGATGGTCATCGCTGTTTCGCCTTTATTAAAGGAAGCTTCTGCGATGGAGTAATCGGTGTCTGCATTCATGTGTTTGTTTTTAATCAGGTCAACCAGGAAGGTCAGACCCGCTTTCGCGCCAGCGTTATCCACGCCCACGTCTTTa                                                                                                                                        <  2:454531/251‑1 (MQ=255)                                                                                                          cGCCAACGAACGGTTTGGATGGTTGACCCTTGAAGGTCGGCCGTACCGCTACACCATAATTCACTTTGCTGCACTCGGTGTTGGACCCTGCACACGGGCCGTTGATGGTCATCGCTGTTTCGCCTTTATTAGAGGAAGCTTCTGCGATGGAGTAATCGGCGTCTGCATTCATGTGTTTGTTTTTAATCAGGTCAACCAGGCAGGTCAGACCCGCTTTCGCGCCAGCGTTAt                                                                                                                                                         >  1:685116/1‑231 (MQ=255)                                                                                                             cAACGAACGGTTTGGATGGTTGACCCTTGAAGGTCGGCAGTACCGTTACACCCTAATTCACTTTGCTGGTGTCGATGTTGGACCATGCCCACGGGCCGTTGATGGTCATCGCTGTTTCGCCTTTATTAAAGGAAGCTTCTGCGATGGAGTAATCGGTGTCTGCATTCATGTGTTTGTTTTTAATCAGGTCAACCAGGAAGGTCAGACCCGCTTTCGCGCCAGCGTTATCCACGCCCACGTCTTTAATGTCg                                                                                                                                  >  1:688904/1‑251 (MQ=255)                                                                                                                cgaacgGTTTGGATGGTTGACCCGTGAAGGTCGGCAGGACCGTTACTCCATGATTCACTTTGCTGTTGTCGATGTTGGACCATGCCCTCGGGCCGTTTATGGTCATCGCTGTTTCGCCTTTATTAAAGGAAGCTTCTGCGATGGAGTAATCGGTGTCTGCATTCATGTGTTTGTTTTTGATCAGGTCAACCAGGAAGGTCAGACCCGCTTTCGCGCCAGCGTTCTCCACGCCCACGGCTTTAATGTCGTAc                                                                                                                               <  1:46049/251‑1 (MQ=255)                                                                                                                     ggTTTGGATGGTTGACCCTTGAAGGTCGGCAGTCCCGTTACACCATAATTCACTTTGCTGGTGTCGATGTTGGACCATGCCCACGGGCCGTTGATGGTCATCGCTGTTTCGCCTTTATTAAAGGAAGCTTCTGCGATGGAGTAATCGGTGTCTGCATTCATGTGTTTGTTTTTAATCAGGTCAACCAGGAAGGTCAGACCCGCTTTc                                                                                                                                                                      >  2:441003/1‑207 (MQ=255)                                                                                                                     ggTTTGGATGGTTGACCCTTGAAGGTCGGCAGTACCGTTACACCATAATTCACTTTGCTGGTGTCGATGTTGGACCATGCCCACGGGCCGTTGATGGTCATCGCTGTTTCGCCTTTATTAAAGGAAGCTTCTGCGATGGAGTAATCGGTGTCTGCATTCATGTGTTTGTTTTTAATCAGGTCAACCAGGAAGGTCAGACCCGCTTTc                                                                                                                                                                      <  1:441003/207‑1 (MQ=255)                                                                                                                     ggTTTGGATGGTTGACCCTTGAAGGTCGGCAGTACCGTTACACCATAATTCACTTTGCTGGTGTCGATGTTGGACCATGCCCACGGGCCGTTGATGGTCATCGCTGTTTCGCCTTTATTAAAGGAAGCTTCTGCGATGGAGTAATCGGTGTCTGCATTCATGTGTTTGTTTTTAATCAGGTCAACCAGGAAGGTCAGACCCGCTTTCGCGCCAGCGTTATCCACGCCCACGTCTTTAATGTCGTACTTGcc                                                                                                                          <  1:81678/251‑1 (MQ=255)                                                                                                                        ttGGATGGTTGACCCTTGAAGGTCGGCAGTACCGTTACACCATAATTCACTTTGCTGGTGTCGATGTTGGACCATGCCCACGGGCCGTTGATGGTCATCGCTGTTTCGCCTTTATTAAAGGAAGCTTCTGCGATGGAGTAATCGGTGTCTGCATTCATGTGTTTGTTTTTAATCAGGTCAACCAGGAAGGTCAGACCCGCTTTCGCGCCAGCGTTATCCACGCCCACGTCTTTAATGTCGTACTTGCCGt                                                                                                                        <  2:667312/250‑1 (MQ=255)                                                                                                                           gATGGTTGACCCTTGAAGGTCGGCAGTACCGTTACACCATAATTCACTTTGCTGGTGTCGATGTTGGACCATGCCCACGGGCCGTTGATGGTCATCGCTGTTTCGCCTTTATTAAAGGAAGCTTCTGCGATGGAGTAATCGGTGTCTGCATTCATGTGTTTGTTTTTAATCAGGTCAACCAGGAAGGTCAGACCCGCTTTCGCGCCAGCGTTATCCACGCCCACGTCTTTAATGTCGTACTTGCCGTTTTc                                                                                                                    >  2:513762/1‑251 (MQ=255)                                                                                                                                        tGAAGGTCGGCAGTACCGTTACACCATAATTCACTTTGCTGGTGTCGAGGTTGGACCATGCCCACGGGCCGTTGATGGTCATCGCTGTTTCGCCTTTATTAAAGGAAGCTTCTGCGATGGAGTAATCGGTGTCTGCATTCATGTGTTTGTTTTTAATCAGGTCAACCAGGAAGGTCAGACCCGCTTTCGCGCCAGCGTTATCCACGCCCACGTCTTTAATGTCGTACTTGCCGTTTTCATACTTGAACg                                                                                                         <  1:110165/249‑1 (MQ=255)                                                                                                                                             gTCGGCAGTACCGTTACACCATAATTCACTTTGCTGGTGTCGATGTTGGACCATGCCCACGGGCCGTTGATGGTCATCGCTGTTTCGCCTTTATTAAAGGAAGCTTCTGCGATGGAGTAATCGGTGTCTGCATTCATGTGTTTGTTTTTAATCAGGTCAACCAGGAAGGTCAGACCCGCTTTCGCGCCAGCGTTATCCACGCCCACGTCTTTAATGTCGTACTTGCCGTTTTCATACTTGAACGCATAAcc                                                                                                  >  1:355417/1‑251 (MQ=255)                                                                                                                                                ggCAGTACCGTTACACCATAATTCACTTTGCTGGTGTCGATGTTGGACCATGCCCACGGGCCGTTGATGGTCATCGCTGTTTCGCCTTTATTAAAGGAAGCTTCTGCGATGGAGTAATCGGTGTCTGCATTCATGTGTTTGTTTTTAATCAGGTCAACCAGGAAGGTCAGACCCGCTTTCGCGCCAGCGTTATCCACGCCCACGTCTTTAATGTCGTACTTGCCGTTTTCATACTTGAACGCATAAccccc                                                                                               <  1:228225/251‑1 (MQ=255)                                                                                                                                                       ccGTTACACCATAATTCACTTTGCTGGTGTCGATGTTGGAACATGCCCACGGGCCGTTGATGGTCATCGCTGTTTCGCCTTTATTAAAGGAAGCTTCTGCGATGGAGTAATCGGTGTCTGCATTCATGTGTTTGTTTTTAATCAGGTCAACCAGGAAGGTCAGACCCGCTTTCGCGCCAGCGGTATCCACGCCCACGTCTTTAATGTCGTACGTGCCGTTTTCATACTTGAACGCATAACCCACGgcagca                                                                                        <  2:374917/251‑1 (MQ=255)                                                                                                                                                             cacCATAATTCACTTTGCTGGTGTCGATGTTGGACCATGCCCACGGGCCGTTGATGGTCATCGCTGTTTCGCCTTTATTAAAGGAAGCTTCTGCGATGGAGTAATCGGTGTCTGCATTCATGTGTTTGTTTTTAATCAGGTCAACCAGGAAGGTCAGACCCGCTTTCGCGCCAGCGTTATCCACGCCCACGTCTTTAATGTCGTACTTGCCGTTTTCATACTTGAACGCATAACCCCCGTCAGCAGCAATc                                                                                  >  1:572542/1‑251 (MQ=255)                                                                                                                                                             cacCATAATTCACTTTGCTGGTGTCGATGTTGGACCATGCCCACGGGCCGTTGATGGTCATCGCTGTTTCGCCTTTATTAAAGGAAGCTTCTGCGATGGAGTAATCGGTGTCTGCATTCATGTGTTTGTTTTTAATCAGGTCAACCAGGAAGGTCAGACCCGCTTTCGCGCCAGCGTTATCCACGCCCACGTCTTTAATGTCGTACTTGCCGTTTTCATACTTGAACGCATAACCCCCGTCAGCAGCAATc                                                                                  >  1:487913/1‑251 (MQ=255)                                                                                                                                                              acCATAATTCACTTTGCTGGTGTCGATGTTGGACCATGCCCACGGGCCGTTGATGGTCATCGCTGTTTCGCCTTTATTAAAGGAAGCTTCTGCGATGGAGTAATCGGTGTCTGCATTCATGTGTTTGTTTTTAATCAGGTCAACCAGGAAGGTCAGACCCGCTTTCGCGCCAGCGTTATCCACGCCCACGTCTTTAATGTCGTACTTGCCGTTTTCATACTTGAACGCATAACCCCCGTCAGCAGCAATc                                                                                  >  2:326983/1‑250 (MQ=255)                                                                                                                                                                               tGGTGTCGATGTTGGACCATGCCCACGGGCCGTTGATGGTCATCGCTGTTTCGCCTTTATTAAAGGAAGCTTCTGCGATGGAGTAATCGGTGTCTGCATTCATGTGTTTGTTTTTAATCAGGTCAACCAGGAAGGTCAGACCCGCTTTCGCGCCAGCGTTATCCACGCCCACGTCTTTAATGTCGTACTTGCCGTTTTCATACTTGAACGCATAACCCCCGTCAGCAGCAATCAGCGGCCAGGTGAAGTAc                                                                >  1:444761/1‑251 (MQ=255)                                                                                                                                                                                ggTGTCGATGTTGGACCATGCCCACGGGCCGTTGATGGTCATCGCTGTTTCGCCTTTATTAAAGGAATCTTCTGCGATGGAGTAATCGGTGTCTGCATTCATGTGTTTGTTTTTAATCAGGTCAACCAGGAAGGTCAGACCCGCTTTCGCGCCAGCGTTATCCACGCCCACGTCTTTAATGTCTTACTTGCCGTTTTCATACTTGAACGCATAACCCCCGTCAGCAGCAATCAGCGGCCAGGTGAAGTAc                                                                <  1:399895/250‑1 (MQ=255)                                                                                                                                                                                     cgaggTTGGACCATGCCCACGGGTCGTTGATGGTCATCGCTGGTTCGCCTTTATTAAAGGAAGCTTCTGCGATGGAGTAATCGGTGTCTGCATTCATGTGTTTGTTTTGAATCAGGTCAACCAGGAAGGGCAGACCCGCTTGCGCGCCAGCGTTATCCACGCCCACTTCTTTAATGTCGTACTTGCCGTGTTCATACTTGAACGCATAACCCCCGTCAGCAGCAATCAGCGGCCAGGTGAAGTACGGTTc                                                           <  2:427138/246‑1 (MQ=255)                                                                                                                                                                                             gACCATGACCACGGGCCGTTGCTGGTCATCGCTGTTTCGCCTTTATTAAAGGAAGCTTCTGCGATGGAGTAATCGGTGTCTGCATTCATGTGTTTGTTTTTAATCAGGTCAACCAGGAAGGTCAGACCCGCTTTCGCGCCAGCGTTATCCACGCCCACGTCTTTAATGTCGTACTTGCCGTTTTCATACTTGAACGCATAACCCCCGTCAGCAGCAATCAGCGGCCAGGTGAAGTACGGTTCTTGCAGGtt                                                  >  2:643565/1‑251 (MQ=255)                                                                                                                                                                                              aCCATGCCCACGGGCCGTTGATGGTCATCGCTGTTTCGCCTTTATTAAAGGAAGCTTCTGCGATGGAGTAATCGGTGTCTGCATTCATGTGTTTGTTTTTAATCAGGTCAACCAGGAAGGTCAGACCCGCTTTCGCGCCAGCGTTATCCACGCCCACGTCTTTAATGTCGTACTTGCCGTTTTCATACTTGAACGCATAACCCCCGTCAGCAGCAATCAGCGGCCAGGTGAAGTACGGTTCTTGCAGGTTg                                                 <  1:156717/251‑1 (MQ=255)                                                                                                                                                                                                 aTGCCCACGGGCCGTTGATGGTCATCGCTGTTTCGCCTTTATTAAAGGAAGCTTCTGCGATGGAGTAATCGGTGTCTGCATTCATGTGTTTGTTTTTAATCAGGTCAACCAGGAAGGTCAGACCCGCTTTCGCGCCAGCGTTATCCACGCCCACGTCTTTAATGTCGTACTTGCCGTTTTCATACTTGAACGCATAACCCCCGTCAGCAGCAATCAGCGGCCAGGTGAAGTACGGTTCTTGCAGGTTGAAc                                              <  1:513762/251‑1 (MQ=255)                                                                                                                                                                                                     ccACGGGCCGTTGATGGTCATCGCTTTTTCGCCTTTATTAAAGTAAGCTTCTGCGATGGAGTAATCGGTGTCTGCATTCATGTGTTTGTTTTTAATCAGGTCAACCAGGAAGGTCAGACCCGCTTTCGCGCCAGCGTTATCCACGCCCACGTCTTTAATGTCGTACTTGCCGTTTTCATACTTGAACGCATATCCCCCTTCAGCAGCAATCAGCGGCCAGGTGAAGTACGGTTCTTGCAGGTTGAtcatca                                          <  2:619926/251‑1 (MQ=255)                                                                                                                                                                                                     ccACGGGCCGTTGATGGTCATCGCTGTTTCGCCTTTATTAAAGGAAGCTTCTGCGATGGAGTAATCGGTGTCTGCATTCATGTGTTTGTTTTTAATCAGGTCAACCAGGAAGGTCAGACCCGCTTTCGCGCCAGCGTTATCCACGCCCACGTCTTTAATGTCGTACTTGCCGTTTTCATACTTGAACGCATAACCCCCGTCAGCAGCAATCAGCGGCCAGGTGAAGAACGGTTCTTGCAGGTTGAACATc                                           >  1:279496/1‑250 (MQ=255)                                                                                                                                                                                                          ggacggtgagggTCATCGCTGTTTCGCCTTTATTAAAGGAAGCTTCTGCGATGGAGTAATCGGTGTCTGCATTCATGTGTTTGTTTTTAATCAGGTCAACCAGGAAGGTCACACCCGCTTTCGCGCCAGCGTTATCCACGCCCACGTCTTTAATGTCGTACTTGCCGTTTTCATACTTGAACGCATAACCCCCGTCAGCAGCAATCAGCGGCCAGGTGAAGTCCGGTTCTTGCAGGTTGAACATCAgcgcg                                     >  2:1150/11‑251 (MQ=255)                                                                                                                                                                                                          ggCCGTTGATGGTCATCGCTGTTTCGCCTTTATTAAAGGAAGCTTCTGCGATGGAGTAATCGGTGTCTGCATTCATGTGTTTGTTTTTAATCAGGTCAACCAGGAAGGTCAGACCCGCTTTCGCGCCAGCGTTATCCACGCCCACGTCTTTAATGTCGTACTTGCCGTTTTCATACTTGAACGCATAACCCCCGTCAGCAGCAATCAGCGGCCAGGTGAAGTACGGTTCTTGCAGGTTGAACATCAgcgcg                                     <  1:191987/251‑1 (MQ=255)                                                                                                                                                                                                              gTTGATGGTCATCGCTGTTTCGCCTTTATTAAAGGAAGCTTCTGCGATGGAGTAATCGGTGTCTGCATTCATGTGTTTGTTTTTAATCAGGTCAACCAGGAAGGTCAGACCCGCTTTCGCGCCAGCGTTATCCACGCCCACGACTTTAATGTCGTACTTGCCGTTTTCATACTTGAACGCATAACCCCCGTCAGCAGCAATCAGCGGCCAGGTGAAGTACGGTTCTTGCAGGTTGAACATCAGCGCGctct                                 <  1:684671/251‑1 (MQ=255)                                                                                                                                                                                                                     gTCATCGCTGTTTCGCCTTTATTAAAGGAAGCTTCTGCGATGGAGTAATCGGGGTCTGCATTCATGTGTTTGTTTTTAATCAGGTCAACCAGGAAGGTCAGACCCGCTTTCGCGCCAGCGTTATCCACGCCCACGTCTTTAATGTCGTACTTGCCGTTTTCATACTTGAACGCATAACCCCCGTCAGCAGCAATCAGCGGCCAGGTGAAGTACGGTTCTTGCAGGTGGAACATCAGCGCGCTCTTACCttt                          <  2:488717/251‑1 (MQ=255)                                                                                                                                                                                                                       cagcgctgtgtggccggtattaAAGGAAGATTCTGCTATGGAGTAAGCGGTGTCTGCATTCATGTGTTTGTTTGTAATCAGGTCAACCAGGAAGGTCAGACCCGCTGTCGCGGCAGCGTGAGCCACGCGTACGTCGTTAAGGTCGTACTTGGCGTTTTCAGTCTGGAACTCAGAGGCCCCGTCAGCAGCTATCAGCGGCCAGGGTAAGTACGGTTATGGGAGGTTGAGTATGAGCGCTCTCGTAGCTTTCg                        <  2:630226/234‑1 (MQ=255)                                                                                                                                                                                                                       cATCGCTGTTTCGCCTTTATTAAAGGAAGCTTCTGCGATGGAGTAATCGGTGTCTGCATTCATGTGTTTGTTTTTAATCAGGTCAACCAGGAAGGTCAGACCCGCTTTCGCGCCAGCGTTATCCACGCCCACGTCTTTAATGTCGTACTTGCCGTTTTCATACTTGAACGCATAACCCCCGTCAGCAGCAATCAGCGGCCAGGTGAAGTACGGTTCTTGCAGGTTGAACATCAGCGCGCTCTTACCTTTc                         <  2:550309/250‑1 (MQ=255)                                                                                                                                                                                                                                   gCCTTTATTAAAGGAAGCTTCTGCGATGGAGTAATCGGTGTCTGCAGTCATGTGTTTGTTTTTAATCAGGTCAACCAGGAAGGTCAGACCCGCTTTCGCGCCAGCGTTATCCACGCCCACGTCTTTAATGTCGTACTTGCCGTTTTCATACTTGAACGCATAACCCCCGTCAGCAGCAATCAGCGGCCAGGTGAAGTACGGTTCTTGCAGGTTGAACATCAGCGCGCTCTTACCTTTCGCTTGCAGTTCtt            <  2:568372/251‑1 (MQ=255)                                                                                                                                                                                                                                    ccTTTATTAAAGGAAGCTTCTGCGATGGAGTAATCGGTGTCTGCATTCATGTGTTTGTTTTTAATCAGGTCAACCAGGAAGGTCAGCCCCGCTTTCGCGCCAGCGTTATCCACGCCCACGTCTTTAATGTCGTACTTGCCGTTTTCATACTTGAACGCATAACCCCCGTCAGCAGCAATCAGCGGCCAGGTGAAGTACGGTTCTTGCAGGTTGAACATCAGCGCGCTCTTACCTTTCGCTTTCAGTTCttt           >  1:109438/1‑251 (MQ=255)                                                                                                                                                                                                                                             aaGGAAGCTTCTGCGATGGAGTAATCGGTGTCTGCATTCATGTGTTTGTTTTTAATCAGGTCAACCAGGAAGGTCAGACCCGCTTTCGCGCCAGCGTTATCCACGCCCACGTCTTTAATGTCGTACTTGCCGTTTTCATACTTGAACGCATAACCCCCGTCAGCAGCAATCAGCGGCCAGGTGAAGTACGGTTCTTGCAGGTTGAACATCAGCGCGCTCTTACCTTTCGCTTTCAGTTCTTTATCCAgcgc  >  1:493197/1‑251 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                       CTTTATTAACCGCTTCCAGACCTTCATCAGTCAGCAGATAGTTTTCGAGGAACTCTTTCGCCAGCTCTTTGTTCGGACTGGCGGCGTTAATACCTGCGCTCAGCACGCCAACGAACGGTTTGGATGGTTGACCCTTGAAGGTCGGCAGTACCGTTACACCATAATTCACTTTGCTGGTGTCGATGTTGGACCATGCCCACGGGCCGTTGATGGTCATCGCTGTTTCGCCTTTATTAAAGGCAGCTTCTGCGATGGAGTAATCGGTGTCTGCATTCATGTGTTTGTTTTTAATCAGGTCAACCAGGAAGGTCAGACCCGCTTTCGCGCCAGCGTTATCCACGCCCACGTCTTTAATGTCGTACTTGCCGTTTTCATACTTGAACGCATAACCCCCGTCAGCAGCAATCAGCGGCCAGGTGAAGTACGGTTCTTGCAGGTTGAACATCAGCGCGCTCTTACCTTTCGCTTTCAGTTCTTTATCCAGCGC  >  NC_000913/4245456‑4245942

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