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.4%	135.9 / ‑4.4	39	coding (13/534 nt)	rrrQ	putative lysozyme
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (18/20);  total (18/21)
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

TGCCAGCGCCTTATCACGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑T‑TTTTGTGTTCATAGATATTTATGATGAGGACGCTCGTGCTTATTGGCAGGATTTTCAATCTTAAAGGAGTACTGATGCTGCAGATAAGACTCAACTTTTTCTGACAATTTTTCTGCTACTTCCAGGAAGACTTGCCGGACGCTCCTTCTGGCTGCTGCCTCATAAAACTCCAGCGCAGCTCCTTCAACACGGTCCATGGCGACATCCAGGTCAAAAATTTCACCGTCAAAGCGTTCTTT  >  NC_000913/1639797‑1640284                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                ggCCAGCGCCTTATCACGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTGTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑tct                                                                                                                                                                                                                                                <  2:95048/250‑3 (MQ=255)                            gTTGACCTGGTCGCATTTTTCCTGCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAAGCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtctt                                                                                                                                                                                                                     <  1:198256‑M1/251‑18 (MQ=255)                               gACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacg                                                                                                                                                                                                                  >  1:208176‑M1/1‑231 (MQ=255)                                                   tCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCCCCACTGGGGCACCCCGACAGATGTTCCATATGCCAGAACCATCGCGGTATGCCATGGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCGcccgcccgcttcgtcttacgatattcatcataaactttg                                                                                                                                                                                               <  2:298119‑M1/250‑40 (MQ=255)                                                      gcaGTTGCATATTGGGAAAAACGTTTTTTCCATCCTCCGCTGTGGCGCCCCGACAGATGTTCCATATGCCAGATCCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTGTCTGAGGAGCAGATTCGCCATCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCTcccgcccgcttcgtcttacgatattcatcataaactttgggat                                                                                                                                                                                           <  2:313763‑M1/250‑44 (MQ=255)                                                       cAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatc                                                                                                                                                                                          >  2:278387‑M1/1‑207 (MQ=255)                                                         gTTTCATATTGGGACAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCAATTGTGTGGGTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGCCAGGCCGTATCTGA‑TCTTTTGTGTTAAcccgcccgcttcgtcttacgatattcatcataaactttgggatcat                                                                                                                                                                                        >  2:306375‑M1/1‑205 (MQ=255)                                                             cATATTGGGAAAAACGGTTTTTCCATCCACAACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggcgcatactg                                                                                                                                                                                    >  2:57460‑M1/1‑201 (MQ=255)                                                                   gggAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctc                                                                                                                                                                              >  1:167590‑M1/1‑195 (MQ=255)                                                                                               gtgGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCGTTTTCGGCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaacttttggatcatactgaagctccccgccagatgcctcctgtagacgca                                                                                                                                                    <  1:378117‑M1/249‑83 (MQ=255)                                                                                                          aCAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTTTTTGGTTACCTTCTTTTGCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCCGCACCAATCAGCGCCAGAACGGCAGCCGTCAGGCGTTATTTGAGT‑TTTTGTGTTCCcccgtccgcttcgtcttacgatattcatcataaactttgggctcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttc                                                                                                                                       <  2:367874‑M1/251‑96 (MQ=255)                                                                                                          aCAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttc                                                                                                                                       <  1:105089‑M1/251‑96 (MQ=255)                                                                                                                ggTCCATATGCCAGAACCATCGCGGTATGCCATTGTTTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGGATCTGAGGAGCAGATGCGCCAGCTCCAATCAGCGCCAGAACGGCAGCCAACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataatctttgggatcatactgaatctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaac                                                                                                                                 <  1:469508‑M1/251‑102 (MQ=255)                                                                                                                            aGAACCATCGCGGGATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGACGTATCTGA‑TCTTTTGGGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagatccccgccagatgcctcctgaagacgcatcgcgcgcccttcgggaactaactccccttt                                                                                                                     >  2:209446‑M1/1‑138 (MQ=255)                                                                                                                                 catagcggtatGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGAAGCGCGAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGACGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgcaagatgcctcctgtagacgcatcgcgcgaccttcgggagctaaatcccctttccagc                                                                                                                >  1:450608‑M1/5‑133 (MQ=255)                                                                                                                                      cgGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctata                                                                                                            >  1:276929‑M1/1‑128 (MQ=255)                                                                                                                                          aTGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTAGCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagc                                                                                                        <  1:30829‑M1/250‑127 (MQ=255)                                                                                                                                            gcccttgtgtGGTTCCCTTCTATTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCTCCAGCCCCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttactatattcatcataaactttgggatcatactgaagctccccgccagatgccgcctgtatacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaa                                                                                                     <  1:138267‑M1/247‑130 (MQ=255)                                                                                                                                             ccATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaag                                                                                                    <  1:440894‑M1/251‑131 (MQ=255)                                                                                                                                             ccATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcga                                                                                                      >  2:459919‑M1/1‑121 (MQ=255)                                                                                                                                                       gTTACCTTCTTTTTCGTCCAGAAACGGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaat                                                                                          <  1:439759‑M1/251‑141 (MQ=255)                                                                                                                                                         tACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaata                                                                                         <  1:459919‑M1/250‑142 (MQ=255)                                                                                                                                                                  tttCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgctt                                                                               >  2:150047‑M1/1‑100 (MQ=255)                                                                                                                                                                    ccttccataatctgaccaagtatctgaggacccGCTGCGCCAGCACCAATCATCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGGGTTCAcccgcccgcttcgtcttacgatattcatcataaactttggtatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgt                                                                             <  2:41802‑M1/235‑154 (MQ=255)                                                                                                                                                                        ccAGAAACTGGTCAAGTACCTGAGGAGCAGATGCGCCAGCACCAACCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaa                                                                         <  1:209446‑M1/251‑158 (MQ=255)                                                                                                                                                                                 ggTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgttt                                                                >  2:193131‑M1/1‑85 (MQ=255)                                                                                                                                                                                    cAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgttttt                                                              <  1:150047‑M1/250‑169 (MQ=255)                                                                                                                                                                                     aaGTATCTGAGGGGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttga                                                            <  1:148586‑M1/251‑171 (MQ=255)                                                                                                                                                                                     aaGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgactacctgctgcttgtgcaagttttgtttttg                                                             >  2:364783‑M1/1‑81 (MQ=255)                                                                                                                                                                                       gTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcatgttttgtttttgaac                                                          <  1:121180‑M1/251‑173 (MQ=255)                                                                                                                                                                                            tgaagagcagaTGCGCCAACACCAATCAGCTCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgacactccacgccagatgcctcctgtagacgcatcgcgcgaccttcgggtactaaatcccctttccagctctcaagcgaagccaaacgaaaacctgctgcttgtgcaagttttgtttttgaacagaa                                                      >  2:370007‑M1/5‑74 (MQ=255)                                                                                                                                                                                              gggggcagaGGCCCCAGCCCCCATCAGCGCACGAACGGCAGCGGCCAGGCCGTATCTGA‑TCTTTTTTGTTCCcccgcccgcttcgtcttacgatattcaacataaactgtgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaata                                                   <  2:194460‑M1/247‑180 (MQ=255)                                                                                                                                                                                                    aGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTCTCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagcaccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgccgctagagcacgctttgtttttgaaccgaactacaaa                                               >  2:458985‑M1/1‑66 (MQ=255)                                                                                                                                                                                                        gcgcCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaaaagagc                                          <  1:306375‑M1/250‑189 (MQ=255)                                                                                                                                                                                                                aCCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaaaagagcgtcagtttt                                 >  1:430194‑M1/1‑54 (MQ=255)                                                                                                                                                                                                                     ccaGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaaaagagcgtcagttttaagca                            <  2:167590‑M1/250‑203 (MQ=255)                                                                                                                                                                                                                       aGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgagatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaaaagagcgtccgttttaagcatt                          >  1:129150‑M1/1‑47 (MQ=255)                                                                                                                                                                                                                                         gACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaaaagagcgtcagttttaagcatttaaaacacctttattgtt        <  1:57460‑M1/251‑223 (MQ=255)                                                                                                                                                                                                                                               ccGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaaaagagcgtcagttttaagcatttaaaacacctttattgttagtcat  >  1:205204‑M1/1‑23 (MQ=255)                                                                                                                                                                                                                                               ccGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaaaagagcgtcagttttaagcatttaaaacacctttattgttagtcat  >  2:293857‑M1/1‑23 (MQ=255)                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                TGCCAGCGCCTTATCACGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑T‑TTTTGTGTTCATAGATATTTATGATGAGGACGCTCGTGCTTATTGGCAGGATTTTCAATCTTAAAGGAGTACTGATGCTGCAGATAAGACTCAACTTTTTCTGACAATTTTTCTGCTACTTCCAGGAAGACTTGCCGGACGCTCCTTCTGGCTGCTGCCTCATAAAACTCCAGCGCAGCTCCTTCAACACGGTCCATGGCGACATCCAGGTCAAAAATTTCACCGTCAAAGCGTTCTTT  >  NC_000913/1639797‑1640284

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