Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,640,045	+C	coding (13/534 nt)	rrrQ ←	putative lysozyme

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,640,045	1	.	C	100.0%	167.2 / NA	44	I5V (ATC→GTC)	rrrQ	putative lysozyme
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (14/30);  total (14/30)

CATGCCAGCGCCTTATCACGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGAT‑TTTTGTGTTCATAGATATTTATGATGAGGACGCTCGTGCTTATTGGCAGGATTTTCAATCTTAAAGGAGTACTGATGCTGCAGATAAGACTCAACTTTTTCTGACAATTTTTCTGCTACTTCCAGGAAGACTTGCCGGACGCTCCTTCTGGCTGCTGCCTCATAAAACTCCAGCGCAGCTCCTTCAACACGGTCCATGGCGACATCCAGGTCAAAAATTTCACCGTCAAAGCGTTC  >  NC_000913/1639795‑1640281                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                             cATGCCAGCGCCTTATCACGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGAt                                                                                                                                                                                                                                               <  1:270353/251‑1 (MQ=255)   tGCCAGCGCATTATCACGCTCACTGGCGTTGACCTGGTCGCATTTTTCCTTCGACCGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGCTGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGAtct                                                                                                                                                                                                                                             >  2:427496/1‑249 (MQ=255)     ccAGCGCCTTATCACGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGAt‑cttt                                                                                                                                                                                                                                          >  2:24804/1‑251 (MQ=255)                      tCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgc                                                                                                                                                                                                                          >  1:181060‑M1/1‑242 (MQ=255)                       cAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGCAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCTGTATCTGATCTTTTGTGTTCAcccgcccg                                                                                                                                                                                                                           <  2:523989‑M1/249‑9 (MQ=255)                          tGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcg                                                                                                                                                                                                                      >  1:255068‑M1/1‑238 (MQ=255)                            gCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtc                                                                                                                                                                                                                    <  1:480163‑M1/251‑16 (MQ=255)                             cGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtct                                                                                                                                                                                                                   <  1:427496‑M1/251‑17 (MQ=255)                             cGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtct                                                                                                                                                                                                                   >  2:487035‑M1/1‑235 (MQ=255)                                tGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtctta                                                                                                                                                                                                                 <  2:539168‑M1/250‑19 (MQ=255)                                   ccTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGCTCTTTTGTGTTCAcccgcccgcttcgtcttacgat                                                                                                                                                                                                             >  2:490297‑M1/1‑229 (MQ=255)                                                 tCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGTCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaac                                                                                                                                                                                                <  1:365442‑M1/250‑36 (MQ=255)                                                   cgtcgACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACGGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaacttt                                                                                                                                                                                             <  2:35485‑M1/249‑39 (MQ=255)                                                              tCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaacttttggatcatac                                                                                                                                                                                   <  1:554428‑M1/250‑49 (MQ=255)                                                              gcaTATTGGGAAAAACGTTTTTCCCATCCACCACTTGGGCCCCCCGTCAGATGGTCCATATGCCAGAACCATCCCGGTATGCCATTGTGTGGTTCCCTTCTTTTTCGTCCAGAAACTGGTTAAGTATCTGAGGAGCAGATTCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCGTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatac                                                                                                                                                                                   <  2:537098‑M1/249‑49 (MQ=255)                                                               cATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactg                                                                                                                                                                                 <  1:13141‑M1/251‑51 (MQ=255)                                                                  atTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaag                                                                                                                                                                              <  1:24804‑M1/251‑54 (MQ=255)                                                                        aaaaaCGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctcccc                                                                                                                                                                        <  2:224897‑M1/251‑60 (MQ=255)                                                                         aaaaCGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccg                                                                                                                                                                       <  1:384221‑M1/251‑61 (MQ=255)                                                                          aaaCGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaacttttggatcatactgaagctccccgc                                                                                                                                                                      >  2:101020‑M1/1‑190 (MQ=255)                                                                          aaaCGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgc                                                                                                                                                                      <  1:432124‑M1/251‑62 (MQ=255)                                                                           aaCGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgcc                                                                                                                                                                     >  2:106297‑M1/1‑189 (MQ=255)                                                                                       atccaCCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctg                                                                                                                                                         >  1:259105‑M1/1‑177 (MQ=255)                                                                                                   ggCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcaacgc                                                                                                                                             <  2:73519‑M1/251‑87 (MQ=255)                                                                                                    cccccccggctgatggtccatatGCCTGCGCCATCGCGGTATGCCATTGTGGGTTTACCTTCTTTTCCGTCAAGAAACGGGTCTAGTTTCTGAGGAGCAGATGCGCCGGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTTATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatgctgaagctccccgccagatgcctcctgtagacgcatcgcg                                                                                                                                            <  2:392677‑M1/240‑88 (MQ=255)                                                                                                            aCAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgctt                                                                                                                                                                                                                        >  2:52498‑M1/1‑156 (MQ=255)                                                                                                            aCAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgctt                                                                                                                                                                                                                        <  1:52498‑M1/167‑12 (MQ=255)                                                                                                                         atGCCAGACCCATCGCGTTATGCAATTGTGGGGTTCCCAGCTTTTTCGTCAAGAATCTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCCCCAATCAGCGCCAGAACGGCCGCCGTCAGGCCGTATCTGATCTTTTGTGTTCCcccgcccgcttcgtcttacgatagtcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcc                                                                                                                       <  2:220364‑M1/251‑109 (MQ=255)                                                                                                                                      cgcgGTATGCCATTGTGTGGTTACCTTCTTTTTCGGCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCTCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTAGCTGATCTTTTGGGTTGAcccgcgcgcttcgtcttaagatattcatcataagctttgggatcatactgaagctccccgccagatgcctcctgtcgacgcatcgcgcgaccttcgggaactaaagcccctttccagctat                                                                                                          >  2:233151‑M1/1‑130 (MQ=255)                                                                                                                                          gTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaag                                                                                                      <  1:365884‑M1/251‑126 (MQ=255)                                                                                                                                            aTGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcg                                                                                                    >  2:117812‑M1/1‑124 (MQ=255)                                                                                                                                                    gtgtgGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaa                                                                                             >  2:369835‑M1/1‑116 (MQ=255)                                                                                                                                                                   ttttCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgct                                                                                <  1:406629‑M1/248‑148 (MQ=255)                                                                                                                                                                    tttCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACGGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgct                                                                             <  1:411286‑M1/250‑151 (MQ=255)                                                                                                                                                                              aaaCGGGTAAAGTTTCGGAGGAGCCGATGGGCCAGCCCCAATCATCGCCTGAACGGCAGCCGACAGGCCATATCTGTTCTTCTGTGTTCAcccgcccgcttcggcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttt                                                                  <  2:513625‑M1/251‑162 (MQ=255)                                                                                                                                                                               aaCTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagtttt                                                                 >  2:195579‑M1/1‑89 (MQ=255)                                                                                                                                                                               aaCTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagct                                                                                                            <  1:213562‑M1/208‑120 (MQ=255)                                                                                                                                                                                   ggTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgttt                                                             <  1:469478‑M1/251‑167 (MQ=255)                                                                                                                                                                                       aaGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCGGCCGACAGGCCGTATCTGATCTTTTGTTTGCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttga                                                         <  2:538700‑M1/251‑171 (MQ=255)                                                                                                                                                                                         gTATCTGATGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaac                                                       <  2:32249‑M1/251‑173 (MQ=255)                                                                                                                                                                                         gTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaac                                                       <  1:564‑M1/251‑173 (MQ=255)                                                                                                                                                                                               gaggagCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaa                                                   <  1:64581‑M1/249‑177 (MQ=255)                                                                                                                                                                                                        aTGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaaaagag                                        <  1:344202‑M1/251‑188 (MQ=255)                                                                                                                                                                                                                            gcCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaaaagagcgtcagttttaagcattta                     >  1:529705‑M1/1‑44 (MQ=255)                                                                                                                                                                                                                              cAGAACGGCAGCCGACAGGCCGTATCTGATCTCTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataacctttgggatcatcctgaagctctccgccagatgcctcctgtagacgcatcgcgcgacctccgggaactaaatcccctttccagctataaagcgaagccaaaccaattcctgctgcttgtgcaagttttgtttttgaaccgaaatacaaaagagcgtcagttttaagcatttaaaa                  <  1:233151‑M1/251‑210 (MQ=255)                                                                                                                                                                                                                                            aCAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaaaagagcgtcagttttaagcatttaaaacacctttattgtta    <  1:7630‑M1/251‑224 (MQ=255)                                                                                                                                                                                                                                              tggcCGTATCTGATCTTTTGTGTTCAcccgcctgcttcgtgttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcacctttccagctataacacgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaaaagagcgtcaattttaagcatttaaaacacctttattgttagt  >  2:179729‑M1/2‑26 (MQ=255)                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                             CATGCCAGCGCCTTATCACGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGAT‑TTTTGTGTTCATAGATATTTATGATGAGGACGCTCGTGCTTATTGGCAGGATTTTCAATCTTAAAGGAGTACTGATGCTGCAGATAAGACTCAACTTTTTCTGACAATTTTTCTGCTACTTCCAGGAAGACTTGCCGGACGCTCCTTCTGGCTGCTGCCTCATAAAACTCCAGCGCAGCTCCTTCAACACGGTCCATGGCGACATCCAGGTCAAAAATTTCACCGTCAAAGCGTTC  >  NC_000913/1639795‑1640281

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