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%	171.3 / NA	49	I5V (ATC→GTC)	rrrQ	putative lysozyme
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (21/28);  total (21/28)

GCCAGCGCCTTATCACGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGAT‑TTTTGTGTTCATAGATATTTATGATGAGGACGCTCGTGCTTATTGGCAGGATTTTCAATCTTAAAGGAGTACTGATGCTGCAGATAAGACTCAACTTTTTCTGACAATTTTTCTGCTACTTCCAGGAAGACTTGCCGGACGCTCCTTCTGGCTGCTGCCTCATAAAACTCCAGCGCAGCTCCTTCAACACGGTCCATGGCGACATCCAGGTCAAAAATTTCACCGTCAAAGCGTTCTTTGTC  >  NC_000913/1639798‑1640287                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                   gCCAGCGCCTTATCACGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTGCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGAtct                                                                                                                                                                                                                                                   <  1:139543/250‑3 (MQ=255)   cAGCGCCTTATCACGCTCAATGGCGTTGACCTGGTCGCATTCTTCCTTCGACAGTTTCATATTGGGCAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGCGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCGGAt‑ctttt                                                                                                                                                                                                                                               >  2:496220/1‑250 (MQ=255)           tATCACGCTCAATGTCGTTGACCTTGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGAGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCATCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTATTGTGTTCa                                                                                                                                                                                                                                         <  1:312447/250‑1 (MQ=255)           tATCACGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAc                                                                                                                                                                                                                                        >  1:164529‑M1/1‑250 (MQ=255)                       tGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcg                                                                                                                                                                                                                            >  2:397157‑M1/1‑238 (MQ=255)                          cGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCACGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCNGATCTTTTGTGTTCAcccgcccgcttagtct                                                                                                                                                                                                                         >  2:236113‑M1/1‑235 (MQ=255)                                           ttttCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCACGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCGTTGGGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatctta                                                                                                                                                                                                         >  1:216902‑M1/1‑218 (MQ=255)                                            tttCCTTCGACAGTTTCAGATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCGCCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATTCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataa                                                                                                                                                                                                        <  1:236113‑M1/250‑34 (MQ=255)                                                  tCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaacttt                                                                                                                                                                                                   <  1:252605‑M1/249‑39 (MQ=255)                                                     ccaGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCCCCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggat                                                                                                                                                                                              <  2:88292‑M1/250‑44 (MQ=255)                                                      cAGTTTCATATGGGGAAAAACGGTTTTTCCATCCACCACTGGGGCCCCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatc                                                                                                                                                                                             <  2:149353‑M1/251‑45 (MQ=255)                                                           tCATATTGGGAAAAACGGTTGTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCAATGGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGGATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaacgttgggatcatac                                                                                                                                                                                         <  1:257055‑M1/250‑49 (MQ=255)                                                            cATATTGGGAAAAACGGTTTTTCCATCCACCCCTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactg                                                                                                                                                                                       <  1:269183‑M1/251‑51 (MQ=255)                                                                ttGGGAAAAACGGTTTTTCCATCCGCCACTGTGGCACCCCGACGGATGGTCCATATGCCAGAACCATCGCGGGATGCAATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGTCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttactatattcatcataaattgtgggatcatactgaag                                                                                                                                                                                    <  2:410989‑M1/250‑54 (MQ=255)                                                                   ggAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctcc                                                                                                                                                                                <  2:164529‑M1/251‑58 (MQ=255)                                                                      aaaaCGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcata                                                                                                                                                                                          <  2:141640‑M1/238‑48 (MQ=255)                                                                      aaaaCGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcata                                                                                                                                                                                          >  1:141640‑M1/1‑191 (MQ=255)                                                                       aaaCGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgc                                                                                                                                                                            <  1:302661‑M1/251‑62 (MQ=255)                                                                        aaCGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGCTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgcc                                                                                                                                                                           >  2:284823‑M1/1‑189 (MQ=255)                                                                         aCGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGACATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgcc                                                                                                                                                                           >  1:307109‑M1/1‑188 (MQ=255)                                                                          cGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgcc                                                                                                                                                                           >  1:258475‑M1/1‑187 (MQ=255)                                                                                        accacTGTGGCACCCCGACAGATGGTCCATATGCCAGAACAATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGCATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgccacctataa                                                                                                                                                            >  2:462775‑M1/1‑173 (MQ=255)                                                                                          cacTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgaagccg                                                                                                                                                         >  2:262454‑M1/1‑171 (MQ=255)                                                                                              gtgGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGACCGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggctcatactgaagctccccgccagatgcctcctgtagacgcatc                                                                                                                                                     >  2:317713‑M1/1‑167 (MQ=255)                                                                                                ggCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTTTGCCATTGTGTGGTTACCTTCTTTTCCGTCCAGAAACTGGTCAAGTATCTGAGGAGCTGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGTCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgc                                                                                                                                                   <  2:443383‑M1/251‑87 (MQ=255)                                                                                                         aaagaTGGTCCATATGCCAGAACCATCGCGGTAGGCCAGGGTGTGGTTACCTTCTGTTTCGTCCAGAAACTGTTCAAGTAGCTGAGGAGCAGACGCGCCAGCTCCAATCAGCGCCAGAACGGCCGCCGACAGGCCGTATCTGATCTTTAGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatc                                                                                                                                                                                             <  2:194890‑M1/198‑45 (MQ=255)                                                                                                         aCAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAGCGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatc                                                                                                                                                                                             >  1:194890‑M1/1‑156 (MQ=255)                                                                                                         aCAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttc                                                                                                                                          <  1:284823‑M1/251‑96 (MQ=255)                                                                                                               ggTCCATATGCCAGAACCATCGCGTTATGCCATGTTGTGGTTCCCTTCTTTTTCGTCCAGAACCTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaacttttggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaac                                                                                                                                    <  2:429212‑M1/251‑102 (MQ=255)                                                                                                                  ccATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGAtct                                                                                                                                                                                                                                                   <  1:213120/137‑3 (MQ=255)                                                                                                                  ccATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGAtct                                                                                                                                                                                                                                                   >  2:213120/1‑135 (MQ=255)                                                                                                                            gAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgccccctgtagacgcatcccgcgaccttcggggactaaatcccctttc                                                                                                                       >  2:365214‑M1/1‑137 (MQ=255)                                                                                                                                    gcgttatgccttttggtggttaccttctttttcGCCCAGAAACTGGTCAGGTATCTGAGGGGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgttattcatcataaactttgggatcatactgaagccccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctata                                                                                                               <  2:316190‑M1/236‑123 (MQ=255)                                                                                                                                                          ccTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacct                                                                                         <  1:397157‑M1/251‑145 (MQ=255)                                                                                                                                                                tttacgtccAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCCCCAGAACGGCACCCTCCAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaataccttct                                                                                      <  2:55103‑M1/244‑148 (MQ=255)                                                                                                                                                                     gTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgc                                                                              >  1:399395‑M1/1‑96 (MQ=255)                                                                                                                                                                        cAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgca                                                                             >  1:318335‑M1/1‑93 (MQ=255)                                                                                                                                                                                ggTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgttt                                                                   >  2:21510‑M1/1‑85 (MQ=255)                                                                                                                                                                                        aTCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtg                                                                               >  2:154706‑M1/1‑77 (MQ=255)                                                                                                                                                                                        aTCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtg                                                                               <  1:154706‑M1/231‑155 (MQ=255)                                                                                                                                                                                                     aggcgccaGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaaaagag                                              <  1:151093‑M1/249‑188 (MQ=255)                                                                                                                                                                                                     aggcgccaGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggat                                                                                                                                                                                              <  2:171674‑M1/105‑44 (MQ=255)                                                                                                                                                                                                     aTGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggat                                                                                                                                                                                              >  1:171674‑M1/1‑64 (MQ=255)                                                                                                                                                                                                       gcgcCAGCACCACTCAGCGCCAGACCGGCGGCCGACAGGCCGTATCGGATCTTTTGGGTTCAcccgcccgctccgtcttacgatattcatcataaactttgggatcatactgaagccccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaaaagagcg                                            <  2:170948‑M1/251‑190 (MQ=255)                                                                                                                                                                                                            cgcACCAATCAGCGCCGGCACGGCAGCCGGCAGGCCGTATCTGATCTTTTGTGTTCCcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgccgcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgttggtgaaccgaaatacaaaagagcgtcagt                                       <  2:304189‑M1/250‑195 (MQ=255)                                                                                                                                                                                                             cccccaaaCAGCTCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaaaagagcgtcagtt                                      <  1:290144‑M1/248‑196 (MQ=255)                                                                                                                                                                                                                 cAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttcc                                                                                                                      <  2:335404‑M1/167‑116 (MQ=255)                                                                                                                                                                                                                 cAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttcc                                                                                                                      >  1:335404‑M1/1‑52 (MQ=255)                                                                                                                                                                                                                       gcgcCAGAACGGCGGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaaaagagcgtcagttttaagcattt                            <  1:42321‑M1/251‑206 (MQ=255)                                                                                                                                                                                                                        cgcCAGACCGGCAGCCGACAGGCCTTATCTGATCTTTTGTGTTCAcccgcccgctacgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatgcctgctgcttgtgcaagttttgtttttgaaccgaaatacaaaagagcgtcagttttaagcattta                           <  1:365214‑M1/251‑207 (MQ=255)                                                                                                                                                                                                                                    aGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcacgttttgtttttgaaccgaaatacaaaagagcgtcagttttaagcatttaaaacacctttat               >  2:24910‑M1/1‑33 (MQ=255)                                                                                                                                                                                                                                       cGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaaaagagcgtcagttttaagcatttaaaacacctttattgt            >  1:88758‑M1/1‑30 (MQ=255)                                                                                                                                                                                                                                                  aTCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaaaagagcgtcagttttaagcatttaaaacacctttattgttagtcataac  <  1:148234‑M1/250‑232 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                   GCCAGCGCCTTATCACGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGAT‑TTTTGTGTTCATAGATATTTATGATGAGGACGCTCGTGCTTATTGGCAGGATTTTCAATCTTAAAGGAGTACTGATGCTGCAGATAAGACTCAACTTTTTCTGACAATTTTTCTGCTACTTCCAGGAAGACTTGCCGGACGCTCCTTCTGGCTGCTGCCTCATAAAACTCCAGCGCAGCTCCTTCAACACGGTCCATGGCGACATCCAGGTCAAAAATTTCACCGTCAAAGCGTTCTTTGTC  >  NC_000913/1639798‑1640287

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