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	97.8%	156.2 / ‑4.3	45	coding (13/534 nt)	rrrQ	putative lysozyme
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (23/21);  total (23/22)
Fisher's exact test for biased strand distribution p-value = 4.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCGCCTTATCACGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑T‑TTTTGTGTTCATAGATATTTATGATGAGGACGCTCGTGCTTATTGGCAGGATTTTCAATCTTAAAGGAGTACTGATGCTGCAGATAAGACTCAACTTTTTCTGACAATTTTTCTGCTACTTCCAGGAAGACTTGCCGGACGCTCCTTCTGGCTGCTGCCTCATAAAACTCCAGCGCAGCTCCTTCAACACGGTCCATGGCGACATCCAGGTCAAAAATTTCACCGTCAAAGCGTTCTTT  >  NC_000913/1639801‑1640284                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                aGCGCCTTATCACGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTtg                                                                                                                                                                                                                                            >  2:71800/1‑251 (MQ=255)   cgcCTTATCACGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTtgtg                                                                                                                                                                                                                                          <  1:84652/251‑1 (MQ=255)           cACGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAccc                                                                                                                                                                                                                                   <  1:398507‑M1/250‑4 (MQ=255)             cGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCTATAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcc                                                                                                                                                                                                                                <  1:46714‑M1/251‑7 (MQ=255)               cTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGATTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccg                                                                                                                                                                                                                              >  1:527532‑M1/1‑243 (MQ=255)                      gCGTTTACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGACAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGGATCTGA‑TCTTTTGTGTTCCcccgcccccttcgtc                                                                                                                                                                                                                       >  2:303882‑M1/1‑236 (MQ=255)                       cGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtct                                                                                                                                                                                                                      >  2:398827‑M1/1‑235 (MQ=255)                               tagtcGCATTTTTCCTTCGACAGTTTCATATTGGGAAAACCGGTTTTTCCATCCACCACTGTGGCACCCCGACAGAAGGTCAATATGCCAGAACCATCGCGGTATGCCATTGTGAGGTTACCTTCTTTTTCGTCACGAAACTCGTCAAGTATCTGAGGAGCAGATGCGCCAACACCAATCAGCGCCAGAACGGCCGCCGACAGGCCGAATCTGA‑TCTTTTGTGTTCAcacgcccgcttcgtcttccgatat                                                                                                                                                                                                              >  2:345369‑M1/3‑227 (MQ=255)                                    gCATTTCTCCTTCTACAGTTTCATATTGGGAAAACAGGCTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCCGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatc                                                                                                                                                                                                         >  2:213066‑M1/1‑222 (MQ=255)                                              ttCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttg                                                                                                                                                                                               >  1:9082‑M1/1‑212 (MQ=255)                                                   cAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACACCGACAGATGGTCCATATGCCAGAACCATCGCGGCATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGCCAGGCCGAATCTGA‑TCTTTTGTATTCAcccgcccgcttcgtcttacgatataccacagatactttgtgatc                                                                                                                                                                                          >  2:295614‑M1/1‑207 (MQ=255)                                                        tCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatac                                                                                                                                                                                      <  1:398827‑M1/250‑49 (MQ=255)                                                            atTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCATCATCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaa                                                                                                                                                                                  <  1:71800‑M1/250‑53 (MQ=255)                                                               gggAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctc                                                                                                                                                                              >  1:32278‑M1/1‑195 (MQ=255)                                                                            tttCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcc                                                                                                                                                                 >  1:374697‑M1/1‑182 (MQ=255)                                                                              tccatccaCCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTAAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcct                                                                                                                                                                <  1:380503‑M1/250‑71 (MQ=255)                                                                                 atccaCCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctg                                                                                                                                                            >  1:346807‑M1/1‑177 (MQ=255)                                                                                  tcccccactgtgTCCCCCCGACAGATGGCCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGGCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcct                                                                                                                                                             <  2:287531‑M1/245‑74 (MQ=255)                                                                                   ccaccaCTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgt                                                                                                                                                           >  1:255992‑M1/1‑175 (MQ=255)                                                                                   ccaccaCTGTGGCACCCCGAAAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgt                                                                                                                                                           >  1:114855‑M1/1‑175 (MQ=255)                                                                                           gtgGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcat                                                                                                                                                   <  1:116349‑M1/250‑84 (MQ=255)                                                                                                cccccgACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTGTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCGcccgcccgcttcgtcttacgatattcatcataaactttgggatcatgctgaagctccccgccagatgcctcctgtagacgcatcgcgcg                                                                                                                                             <  1:466648‑M1/250‑90 (MQ=255)                                                                                                     gACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGAAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgacctt                                                                                                                                        >  2:409158‑M1/1‑157 (MQ=255)                                                                                                      ccaGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttc                                                                                                                                       <  1:89210‑M1/250‑96 (MQ=255)                                                                                                       cAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGCCCAGAAACTGGGCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaacttggggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcg                                                                                                                                      <  2:255992‑M1/251‑97 (MQ=255)                                                                                                         gATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcggg                                                                                                                                    >  1:355084‑M1/1‑153 (MQ=255)                                                                                                                                   ggTATGCCATTGTGTGGTTACCTTCTGTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTTTCTGA‑TCTTTTGTGTGCAcccgcccgcttcgtcgtacgatattcagcataaacgttgggatcatactgaggctccccggcagatgcctcctgtaggcgtgtggggcgacctgcgggaactaaatcccctttccagctataaa                                                                                                          <  1:58556‑M1/251‑125 (MQ=255)                                                                                                                                        gCCATTGTGTGGTTACCTTCTTTTCCGTCCAGAACCTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGACCGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaa                                                                                                     <  2:346807‑M1/251‑130 (MQ=255)                                                                                                                                                          tcttTTTCGTCCAGAAACTGGTCAAGTATCTGAGGCGCAGATGCGCCAGCACCAATCCGCGCCAGAACGGCAGCCGACAGGCCGTCTCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgcccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaacgatcgg                                                                                    >  2:350927‑M1/1‑104 (MQ=255)                                                                                                                                                          tcttTTTCGGCCAGAAACGGGTCAAGTATCTGAGGAGCAGATGCGCCTGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgccctcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaa                                                                                           <  1:350927‑M1/243‑140 (MQ=255)                                                                                                                                                                  ggccagAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGGGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgc                                                                           <  2:374697‑M1/249‑156 (MQ=255)                                                                                                                                                                        aaaCTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttt                                                                     <  2:32278‑M1/251‑162 (MQ=255)                                                                                                                                                                                    tATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaacc                                                         >  2:177698‑M1/1‑78 (MQ=255)                                                                                                                                                                                           ggagCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatac                                                  <  1:503547‑M1/251‑181 (MQ=255)                                                                                                                                                                                           ggagCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgttgttgaaccgaaatac                                                  <  1:161317‑M1/251‑181 (MQ=255)                                                                                                                                                                                            gagCAGATGCGCCAGCACCAATCAGCGCCGGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTGGTGTTCAcccgcccgtttcgtcttacgatattcatcataaactttgggatcatactgaagctacccgccagatgccccctgtaatcgcatcgcgcgaccatcgggaactaaatcacctttccagctataaagcgaagccaaacgaatacctgctgcttgttcaagttttgtttttgaaccgaaataca                                                 <  1:345369‑M1/251‑182 (MQ=255)                                                                                                                                                                                             agCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattaatcataacctttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccatttccagct                                                                                                               >  1:364741‑M1/1‑69 (MQ=255)                                                                                                                                                                                             agCAGATGCGCCAGCACCAATCAGAGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaataca                                                 >  1:167392‑M1/1‑69 (MQ=255)                                                                                                                                                                                               cAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaaa                                               >  2:82711‑M1/1‑67 (MQ=255)                                                                                                                                                                                                           cACCAATCAGCGCCCGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGAGTTCAacagcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaaaagagcgtcagttt                                  >  2:215738‑M1/1‑55 (MQ=255)                                                                                                                                                                                                             ccAATCAGCGCCAGAACGGCAGCCGACAGGCCTTATCTGATT‑TTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctagaaagcgaagccaaacgaatacctgctgcttgtgcaagtttggtttttgaaccgaaatacaaaggagcggcagtttta                                <  2:503199‑M1/251‑199 (MQ=255)                                                                                                                                                                                                                              ggCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaaaagagcgtcagttttaagcatttaaaacacctt               >  2:454708‑M1/1‑36 (MQ=255)                                                                                                                                                                                                                                   ccGACAGGCCGTATCTGA‑TCTTTTGTGTTCCcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagccccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaaaagagcgtcagttttaagcatttaaaacacctttattg          <  1:303882‑M1/251‑221 (MQ=255)                                                                                                                                                                                                                                          gccttatctGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaaaagagcgtcagtgttaagcatttaaaacacctttattgttagtca   <  2:167392‑M1/247‑228 (MQ=255)                                                                                                                                                                                                                                            cGTATCTGA‑TCTTTTGTGCTCAcccgcacgcttcgtattacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgagcgaccttcgggaactaaatcccctttccagctataaagcgacgccaaacgaatacctgctgcttgtgcacgttttgtttttgaaccgaaatacaaaagagcgtcagttttaagcatttaaaacacctttattgttagtctt  >  2:334089‑M1/1‑22 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                AGCGCCTTATCACGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑T‑TTTTGTGTTCATAGATATTTATGATGAGGACGCTCGTGCTTATTGGCAGGATTTTCAATCTTAAAGGAGTACTGATGCTGCAGATAAGACTCAACTTTTTCTGACAATTTTTCTGCTACTTCCAGGAAGACTTGCCGGACGCTCCTTCTGGCTGCTGCCTCATAAAACTCCAGCGCAGCTCCTTCAACACGGTCCATGGCGACATCCAGGTCAAAAATTTCACCGTCAAAGCGTTCTTT  >  NC_000913/1639801‑1640284

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