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	95.9%	171.8 / NA	50	I5V (ATC→GTC)	rrrQ	putative lysozyme
Reads supporting (aligned to +/- strand):  ref base . (0/1);  major base C (18/29);  minor base A (1/1);  total (19/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCT‑GAT‑TTTTGTGTTCATAGATATTTATGATGAGGACGCTCGTGCTTATTGGCAGGATTTTCAATCTTAAAGGAGTACTGATGCTGCAGATAAGACTCAACTTTTTCTGACAATTTTTCTGCTACTTCCAGGAAGACTTGCCGGACGCTCCTTCTGGCTGCTGCCTCATAAAACTCCAGCGCAGCTCCTTCAACACGGTCCATGGCGACATCCAGGTCAAAAATTTCACCGTCAAAGCGTTCT  >  NC_000913/1639811‑1640282                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                              cccgcTCAATGGCGTTGATCTGGTCGATTTTTTCCTTCGACAGTTTCATATTGGGAAAATCGGTTTTTCCATCCACCACTGGTGCACCCCGACAGTTGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCGTCTTGTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCACGCTCCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCT‑GATCTGTTGTGTTCAccc                                                                                                                                                                                                                                 <  2:297965‑M1/248‑4 (MQ=255)  ccgCTCAATGGCGTTGACCGGGGCGCATTTTTACTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCT‑GATCTTTTGTGTTCAcccgc                                                                                                                                                                                                                               <  1:488608‑M1/250‑6 (MQ=255)       cAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCT‑GATCTTTTGTGTTCAcccgcccgct                                                                                                                                                                                                                          <  1:88419‑M1/251‑11 (MQ=255)       cAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCT‑GATCTTTTGTGTTCAcccgcccg                                                                                                                                                                                                                            <  1:103541‑M1/249‑9 (MQ=255)                            aTTTTTCCTTCGACAGTTCCATAGTGGGAAAAACGGTTTTTCCACCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCTCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCT‑GATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcat                                                                                                                                                                                                     <  2:150563‑M1/251‑32 (MQ=255)                             tttttCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGACGTATCT‑GATCTTTTGTGGTTAcccgcccgcttcgtcttacgatattcatcat                                                                                                                                                                                                     >  1:404510‑M1/1‑219 (MQ=255)                                  catgagactgtttcatagtGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATAGGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGAGGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCT‑GATCTTTTGTGTTCAcccgcacgcttcgtcttacgatattcatcataaact                                                                                                                                                                                                <  2:107453‑M1/241‑37 (MQ=255)                                            tttCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGTTGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCT‑GATATTTGGTGTTCAcccgcccgcttcgtcttacgatattcatcatgaaatttgggatcat                                                                                                                                                                                      >  1:6293‑M1/1‑204 (MQ=255)                                                atatTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCT‑GATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactg                                                                                                                                                                                  >  1:551386‑M1/1‑200 (MQ=255)                                                   ttGGGAAAAACGGTTTTTCCATCCACCCCTGTGGCACCCCGACAGCTGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGAATCTGAGGAGCAGATGCGCCAGAACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATAT‑GATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgcagc                                                                                                                                                                              >  2:64480‑M1/1‑197 (MQ=255)                                                     gggAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGGCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCT‑GATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcagcataaactttgggatcatactgaagctc                                                                                                                                                                            <  2:551386‑M1/251‑57 (MQ=255)                                                          aaaCGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCACGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGCAAGCCGACAGGCCGTATCT‑GATCTTTTGTGGTCAcccgcccgcttcgtcttacgatattcatcntaaactttggtatcatactgaagatccccgc                                                                                                                                                                       >  2:275806‑M1/1‑190 (MQ=255)                                                           aaCGGTTTGTCCATCCACCACGGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTTTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCT‑GATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgcc                                                                                                                                                                      <  1:411994‑M1/251‑63 (MQ=255)                                                                  tttCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCT‑GATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcc                                                                                                                                                               >  1:334969‑M1/1‑182 (MQ=255)                                                                  tttCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGGGGTTCCCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACCGGCCGTATCT‑GATCTTTTGTGTTCAcccgaccgcttcgacttacaatattcatcataaactttgggatcatactgcagctcccagacatatgcc                                                                                                                                                               >  2:94570‑M1/1‑182 (MQ=255)                                                                     ccatccaCCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCT‑GATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcc                                                                                                                                                            <  1:341576‑M1/251‑73 (MQ=255)                                                                       atccaCCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCT‑GATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctg                                                                                                                                                          >  1:312950‑M1/1‑177 (MQ=255)                                                                           accacTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCT‑GATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagaagcctcctgtat                                                                                                                                                       >  2:112426‑M1/1‑173 (MQ=255)                                                                                 gtgGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCT‑GATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatc                                                                                                                                                <  1:429916‑M1/251‑85 (MQ=255)                                                                                           gACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCT‑GATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgacctt                                                                                                                                      >  2:466633‑M1/1‑157 (MQ=255)                                                                                               gATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCT‑GATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcggg                                                                                                                                  <  1:267459‑M1/251‑99 (MQ=255)                                                                                                   gTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCCCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCT‑GATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaact                                                                                                                              <  1:206665‑M1/251‑103 (MQ=255)                                                                                                         atGCCAGAACCACCGCGGTAGGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAACCAGCGCCAGAACGGCAGCCGACAGGCCGTATCT‑GATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcc                                                                                                                        <  2:91421‑M1/251‑109 (MQ=255)                                                                                                                                           tACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCT‑GATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatac                                                                                      >  1:327448‑M1/1‑109 (MQ=255)                                                                                                                                                 ttttttcGTCACGAAACGGGTAAAGTATCTAGGGAGCAGATGCTCCAGCACCAACCAGCGCCAGACCGCCAGCCTCCAGGCCGTATCT‑GATCTTTGGTGTTCCcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccatacgaatacctgct                                                                                 <  2:475110‑M1/249‑148 (MQ=255)                                                                                                                                                    ttacgtcCGGAACCTGGTCAAGTACCGGAGGAGCAGATGGGCCAGCCCCATTCAGCGCCAGAACGCCGCCCGCAAGGCCGTATCT‑GATATTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcagactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgct                                                                              <  2:312950‑M1/247‑151 (MQ=255)                                                                                                                                                     ttCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCT‑GATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgctt                                                                             <  1:83405‑M1/250‑152 (MQ=255)                                                                                                                                                      tCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCCGAACGGCAGCCGACAGGCCGTATCT‑GATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcctcgcgcgaccttcgggaactaaatcccctttccagctctaaagcgaagccaaacgaatacctgctgcttgt                                                                           >  1:461172‑M1/1‑98 (MQ=255)                                                                                                                                                       cGTCCAGAAACGGGGCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCGGCGTCAGAACGGCAGCCGACAGGCCGTATCT‑GATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattaatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgt                                                                           <  1:275806‑M1/250‑154 (MQ=255)                                                                                                                                                              aaaCTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCT‑GATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttt                                                                   <  1:137991‑M1/251‑162 (MQ=255)                                                                                                                                                                      acattatcgaaggagcagttgcgccACAACCAACCAGCGCACGAAGGGCTGCCTACAGGCCGTATCT‑GATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgttttt                                                            <  2:532501‑M1/240‑169 (MQ=255)                                                                                                                                                                           aTCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCT‑GATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccg                                                      >  1:198243‑M1/1‑77 (MQ=255)                                                                                                                                                                             cTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCT‑GATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccga                                                     >  2:408805‑M1/1‑75 (MQ=255)                                                                                                                                                                               gaggggcagatgctccctacccaatcagcgcctgaacggcacCCGTCAGGCCTTATTT‑GATCTTTTGTGTTCCccctcccgcttcgtcttacgatattcatcataaactttgggatcatactgaaactccccgccatatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgatccgaaa                                                   <  2:402163‑M1/230‑178 (MQ=255)                                                                                                                                                                                  gagCAGATGCGCCAGCCCCAATCAGCGCCAGACCGGCAGCCGACAGGCCGTATCT‑GATCTTTTTTGTTCCcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaataca                                               <  2:111690‑M1/251‑182 (MQ=255)                                                                                                                                                                                           cgcCCGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCT‑GATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaaaagagcgt                                      <  2:550006‑M1/251‑191 (MQ=255)                                                                                                                                                                                           cgcCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCT‑GATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaaaagagcgt                                      <  1:436798‑M1/251‑191 (MQ=255)                                                                                                                                                                                                     gcacagcGCCAGAACGGCGGCAGACAGCCCGGACCT‑GATCTTTTGTGTTCCcccgcccgcttcgtcttacgatattcatcataaacgtcgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctgtccagctataaagcgatgacaaacgaagacctgctgcttgtgcaagtattgtttttgaaccgaaatacaaaagagcatcagttttaag                            <  2:188787‑M1/248‑201 (MQ=255)                                                                                                                                                                                                       tCAGCGCCAGAACGGCAGCCGACAGGCCGTATCT‑GATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttganccgaaatacaaaagagcgtcagt                                  >  2:265743‑M1/1‑49 (MQ=255)                                                                                                                                                                                                       tCAGCGCCAGAACGGCAGCCGACAGGCCGTATCT‑GATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaaaagagcgtcagt                                  <  1:265743‑M1/243‑195 (MQ=255)                                                                                                                                                                                                        cAGCGCCAGACCGGCAGCCGACAGGCCGTATCT‑GATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaaaagagcgtcagttttaagcat                         <  1:77156‑M1/251‑204 (MQ=255)                                                                                                                                                                                                            gcCAGAACGGCAGCCGACAGGCCGTATCT‑GATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaaaagagcgtcagttttaagcatttaa                     <  1:36647‑M1/251‑208 (MQ=255)                                                                                                                                                                                                               aGAACGGCAGCCGACAGGCCGTATCT‑GATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaaaagagcgtcagttttaagcatttaaaac                  >  1:353942‑M1/1‑41 (MQ=255)                                                                                                                                                                                                                  aCGGCAGCCGACAGGCCGTATCT‑GATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaaaagagcgtcagttttatgcatttaaaccacc               >  1:238782‑M1/1‑38 (MQ=255)                                                                                                                                                                                                                    ggCAGCCGACAGGCCGTATCT‑GATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagtt                                                                    <  1:557004‑M1/196‑161 (MQ=255)                                                                                                                                                                                                                    ggCAGCCGACAGGCCGTATCT‑GATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagtt                                                                    >  2:557004‑M1/1‑36 (MQ=255)                                                                                                                                                                                                                      cgcccgacTGGCCGTATCTACAC‑TTTTGAGTTCAcccccccgcttcgtcttacgatattcatcataaactttgggaacattcggaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttatgaaccgaaatacaaaagagcgtcagttttaagcatttaaaacaccctta           <  2:1970‑M1/248‑218 (MQ=255)                                                                                                                                                                                                                       aGCCGACAGGCCGTATCT‑GATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaaaagagcgtcagttttaagcatttaaaacacctttat          >  1:505728‑M1/1‑33 (MQ=255)                                                                                                                                                                                                                          cGACAGGCCGTATCT‑GATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaaaagagcgtcagttttaagcatttaaaacacctttattgt       <  1:316168‑M1/251‑222 (MQ=255)                                                                                                                                                                                                                               ggCCGTATCT‑GATCTTTTGTGTTCAcccgcccgcttcgtcttgcgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaaaagagcgtcagttttaagcatttaaaacacctttattgttagtc  <  1:2025‑M1/251‑227 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                              CACGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCT‑GAT‑TTTTGTGTTCATAGATATTTATGATGAGGACGCTCGTGCTTATTGGCAGGATTTTCAATCTTAAAGGAGTACTGATGCTGCAGATAAGACTCAACTTTTTCTGACAATTTTTCTGCTACTTCCAGGAAGACTTGCCGGACGCTCCTTCTGGCTGCTGCCTCATAAAACTCCAGCGCAGCTCCTTCAACACGGTCCATGGCGACATCCAGGTCAAAAATTTCACCGTCAAAGCGTTCT  >  NC_000913/1639811‑1640282

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