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	96.0%	180.9 / ‑0.6	50	coding (13/534 nt)	rrrQ	putative lysozyme
Reads supporting (aligned to +/- strand):  ref base . (0/2);  new base C (19/29);  total (19/31)
Fisher's exact test for biased strand distribution p-value = 5.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑T‑TTTTGTGTTCATAGATATTTATGATGAGGACGCTCGTGCTTATTGGCAGGATTTTCAATCTTAAAGGAGTACTGATGCTGCAGATAAGACTCAACTTTTTCTGACAATTTTTCTGCTACTTCCAGGAAGACTTGCCGGACGCTCCTTCTGGCTGCTGCCTCATAAAACTCCAGCGCAGCTCCTTCAACACGGTCCATGGCGACATCCAGGTCAAAAATTTCACCGTCAAAGCGTTCTTT  >  NC_000913/1639819‑1640284                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                aTGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttc                                                                                                                                                                                                                          <  1:389126‑M1/251‑13 (MQ=255)   ggtgttgACCTGTTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCTCCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGTTATGCCATTGTGTGGTTCCCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCGGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTTTGTTCAcccgcccgcttcgt                                                                                                                                                                                                                        <  2:403092‑M1/248‑15 (MQ=255)          aCCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacg                                                                                                                                                                                                                  >  1:459455‑M1/1‑230 (MQ=255)              ggTCGCCTTTTTCATTTGACAGTTTCAGATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGTCAGATTGTCCATATGCCAGAACCATCGCGGGATTCCATTGTGTGGTTACCTTCTTTTTCGTCCTGAAACGGGTCAAGTGTCTGAGGAGCAGATGCTCCAGCACCATTCAGCGCCAGAACGGCAGCCGACAGGCAGTATCTGATT‑TTTTGTGTTCAcccgcccgcttcttcttacgatat                                                                                                                                                                                                              <  2:36180‑M1/250‑25 (MQ=255)                ccgCATTTTTCCGTCGACAGTTTCATAGTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACATATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGTGCCACAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattca                                                                                                                                                                                                           <  1:423522‑M1/250‑28 (MQ=255)                       tttCCTTCGACAGTTTCAGATTAGGAAAAACGGTTTTTCCATCCACTACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGACACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAccctcccgcttcgtcttacgatattcatcataa                                                                                                                                                                                                     >  2:381640‑M1/1‑217 (MQ=255)                          ccTTCGCCAATTCCATATTGGGAAAACCGGTTTTCCCATCCTCCACTGTGGCTCCCCGACAGATGGTCCATATGCCGGAACCATCGCGGTATGCCATTGTGGTGTTGCCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGGGGAGCAGATGCGCCAGCACCAATCATCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTTGTCCcccgcccgcgtcgtcttacgatattcatcataaact                                                                                                                                                                                                  <  2:38583‑M1/250‑37 (MQ=255)                             tCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttg                                                                                                                                                                                               <  2:3412‑M1/250‑40 (MQ=255)                                aCAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggat                                                                                                                                                                                           >  2:327406‑M1/1‑208 (MQ=255)                                  aGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatca                                                                                                                                                                                         <  1:33892‑M1/251‑46 (MQ=255)                                   gTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcat                                                                                                                                                                                        >  1:139407‑M1/1‑205 (MQ=255)                                          atTGGGAAAAACGGTTTTTCCTTCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTTTCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaa                                                                                                                                                                                  <  2:81451‑M1/250‑53 (MQ=255)                                                  aaaCGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgc                                                                                                                                                                         <  1:313998‑M1/251‑62 (MQ=255)                                                    aCGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgca                                                                                                                                                                        >  1:100444‑M1/1‑188 (MQ=255)                                                          tttCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCCTCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggctcatactgaagctccccgccagatgcc                                                                                                                                                                 >  2:375597‑M1/1‑182 (MQ=255)                                                                 cccccactgtGGCACCCCGACAGATGGTCCATATGCCATAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAATTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaacttttggatcatactgaagctcccctccagatgcctcctgta                                                                                                                                                          <  1:160958‑M1/248‑77 (MQ=255)                                                                  caccacTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtag                                                                                                                                                         <  1:345939‑M1/251‑78 (MQ=255)                                                                    ccacTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagac                                                                                                                                                       <  1:370150‑M1/251‑80 (MQ=255)                                                                         gtgGCACCACGACAGCTGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACCGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatacttaagctccccgccagatgcctactgtagacgcatc                                                                                                                                                  >  2:163605‑M1/1‑167 (MQ=255)                                                                           ggCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgc                                                                                                                                                <  1:277553‑M1/251‑87 (MQ=255)                                                                           ggCACCCCGACAGATGGTCCATATGACAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgcagctccccgccagatacctcctgtagacgcctcgc                                                                                                                                                >  2:183892‑M1/1‑165 (MQ=255)                                                                                 ccGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataacctttgggatcatactgacgctccccgccagatgcctcctgaagacgcatatcgcgacc                                                                                                                                          >  2:93766‑M1/1‑159 (MQ=255)                                                                                                       gAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttc                                                                                                                    >  1:378962‑M1/1‑137 (MQ=255)                                                                                                         aCCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttcc                                                                                                                   >  1:151848‑M1/1‑135 (MQ=255)                                                                                                             gcgcggTATGCCATTGTGTGGTTGCCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagcta                                                                                                              <  2:100475‑M1/250‑121 (MQ=255)                                                                                                                    aTGCCCTTGTGTGGTTACCTTCTTTTTCGTCCAGCAAATGGTAAAGCATCTGAGGAGCAGATGCGCTAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTGTCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggctcatactgaagctcaccgccagatgcctcctgtagacgcatcgctcgaccttcggaaactaaatcccctttccagctataaagcg                                                                                                       >  2:379746‑M1/1‑124 (MQ=255)                                                                                                                           tgtttggttcCCTTCTTTTGCGTCCAAAACCGGGTCAAGTATCCGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGCCAGGCCGTATCTGA‑TCTTTTGTGTTCAccctcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagcca                                                                                                 <  2:111238‑M1/246‑134 (MQ=255)                                                                                                                              gtgGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCTCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaaccgaagccaaccg                                                                                             >  2:42435‑M1/1‑114 (MQ=255)                                                                                                                                ggTTACCTTCTTTTCCGTCCAGAAACTGGTCAAGTACCTGAGGAGCGGTTGCGCCACCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaa                                                                                           <  2:382258‑M1/251‑140 (MQ=255)                                                                                                                                           ttttCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgct                                                                                   <  1:217633‑M1/248‑148 (MQ=255)                                                                                                                                           ttttCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGAGGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgct                                                                                   <  1:77422‑M1/248‑148 (MQ=255)                                                                                                                                               cGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGACAGAACGGCAGACGAAAGGCCGTATCTGA‑TCTTTTGTGTTCAccagcccgcttcgtcttacgatataaatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgcggattgtg                                                                            >  1:294601‑M1/1‑97 (MQ=255)                                                                                                                                                        aCTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttg                                                                   <  1:98591‑M1/251‑164 (MQ=255)                                                                                                                                                           gttcaAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgttt                                                                <  2:273891‑M1/249‑167 (MQ=255)                                                                                                                                                              cAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttg                                                             <  1:46145‑M1/251‑170 (MQ=255)                                                                                                                                                                            gCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaaa                                               <  1:453754‑M1/251‑184 (MQ=255)                                                                                                                                                                             cAGATGCGCCAGCTCCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAccggcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctcc                                                                                                                                                                             <  2:190092‑M1/124‑58 (MQ=255)                                                                                                                                                                             cAGATGCGCCAGCACCAATCAGCGCCAGCAAGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgaccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctcc                                                                                                                                                                             >  1:190092‑M1/1‑67 (MQ=255)                                                                                                                                                                               gATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgcgatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagtgttgtttttgaaccgaaatacaaaaga                                            <  1:408307‑M1/251‑187 (MQ=255)                                                                                                                                                                               gATGCGCCAGCACCAATAAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaaacg                                             >  2:433501‑M1/1‑65 (MQ=255)                                                                                                                                                                                 tGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaaaagagc                                          >  1:279665‑M1/1‑63 (MQ=255)                                                                                                                                                                                         cccaaatcagCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgttcaagttttgtttttgaaccgaaatacaaaagagcgtcagttt                                  <  2:321565‑M1/247‑197 (MQ=255)                                                                                                                                                                                                 aGCGCCAGAACGGCAGCCGACAGGCCGTATCTGAGC‑TTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaaaagagcgtcagttttaagcatt                          <  1:266671‑M1/251‑205 (MQ=255)                                                                                                                                                                                                 aGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaaaagagcgtcagttttaagcatt                          <  1:327406‑M1/251‑205 (MQ=255)                                                                                                                                                                                                       aGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataacctttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaaaagagcgtcagttttaagcatttaaaac                    >  1:132212‑M1/1‑41 (MQ=255)                                                                                                                                                                                                          aCGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaaaagagcgtcagttgtaagcatttaaaacacc                 <  1:143484‑M1/251‑214 (MQ=255)                                                                                                                                                                                                           cGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaanagagcgtcagttttaagcatttaaaccacct                >  2:216273‑M1/1‑37 (MQ=255)                                                                                                                                                                                                            ggCAGCCGACAGGCCGTATCTTA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaaaagagcgtcagttttaagcatttaaaacacctt               <  1:433501‑M1/251‑216 (MQ=255)                                                                                                                                                                                                                 ccGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaaaagagcgtcagttttaagcatttaaaacacctttattg          <  1:264917‑M1/251‑221 (MQ=255)                                                                                                                                                                                                                         ccGTATCTGA‑TCTTTTGTGTTCCcccgcccgcatcggcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaaaagatcgtcagttttaagcatttaaaacacctttattgtttgtcat  <  2:425786‑M1/251‑229 (MQ=255)                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                ATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑T‑TTTTGTGTTCATAGATATTTATGATGAGGACGCTCGTGCTTATTGGCAGGATTTTCAATCTTAAAGGAGTACTGATGCTGCAGATAAGACTCAACTTTTTCTGACAATTTTTCTGCTACTTCCAGGAAGACTTGCCGGACGCTCCTTCTGGCTGCTGCCTCATAAAACTCCAGCGCAGCTCCTTCAACACGGTCCATGGCGACATCCAGGTCAAAAATTTCACCGTCAAAGCGTTCTTT  >  NC_000913/1639819‑1640284

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