Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,640,983	C→T	20.1%	F41F (TTC→TTT)	dtpB →	dipeptide and tripeptide permease B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,640,983	0	C	T	20.1%	82.9 / 2.9	48	F41F (TTC→TTT)	dtpB	dipeptide and tripeptide permease B
Reads supporting (aligned to +/- strand):  ref base C (5/32);  new base T (7/3);  total (12/36)
Fisher's exact test for biased strand distribution p-value = 1.06e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTAATCTATTCACCGCATCAATATTAAGAATCTCTGACAGATGTAAACTTTTTCGCGCGTTATCCCTTACGCGTTCATACTTTTCAGGATGGTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGC‑GTACTGGC‑GG‑TTTTCTTCGTTA‑AACAGCTTGGATTCTCGCAAGAGCA‑GGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTGGCGATTGGCTACTTCATGACCGGCATGTCGCTACTTAAGCCTGACCTGATTTTCATCGCCCTGGGGACTATCGCTGTCGGTA  >  NC_000913/3640749‑3641215                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                           aTTAATCTATTCACCGCATCAATATTAAGAATCTCTGACAGATGTAAACTTTTTCGCGCGTTATCCCTTACGCGTTCATACTTTTCAGGATGGTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGC‑GTACTGGC‑GG‑TTTTCTTCGTTa                                                                                                                                                                                                                                        <  1:610303/239‑1 (MQ=255)   taattttttcccctctttaaaagtAAGACTCTGTGACGTATGGAAAGTTTTGGGCGAGTTATCCCTTACTCGTTCATCCTTTTCAGGATGGTATTGGAAGGGTAATAAATATGAAGTCAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATATGTCATGATCGTGGTTGTCGAGTTATGGGAGCGATGCGGCTACTACGGCGTGCAGGGC‑GTACTGGC‑GT‑TTTTCTTCGTTa‑aa                                                                                                                                                                                                                                     <  2:49060/216‑1 (MQ=255)          ttCACCGCATCAATATTAAGAATCTCTGACAGATGTAAACTTTTTCGCGCGTTATCCCTTACGCGTTCATACTTTTCAGGATGGTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGTAGCGATTCGGCGACTACGGCGTTCAGGGA‑GTACTGGC‑GG‑TTTTTTTTCTTa‑aaaaacatg                                                                                                                                                                                                                              >  1:1047575/1‑232 (MQ=255)          ttCACCGCATCAATATTAAGAATCTCTGACAGATGTAAACTTTTTCGCGCGTTATCCCTTACGCGTTCATACTTTTCAGGATGGTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGC‑GTACTGGCGGG‑TTTTTTTTGTTa‑aaaaactg                                                                                                                                                                                                                               >  1:924345/1‑233 (MQ=255)          ttCACCGCATCAATATTAAGAATCTCTGACAGATGTAAACTTTTTCGCGCGTTATCCCTTACGCGTTCATACTTTTCAGGATGGTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGC‑GTACTGGC‑GGTTTTTTTTCTTTA‑AACacttt                                                                                                                                                                                                                               >  1:1033194/1‑235 (MQ=255)               cGCATCAATATTAAGAATCTCTGACAGATGTAAACTTTTTCGCGCGTTATCCCTTACGCGTTCATACTTTTCAGGATGGTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGC‑GTACTGGC‑GG‑TGTTCTTCGTTT‑AACAAATTGGAttt                                                                                                                                                                                                                         >  1:1095151/1‑238 (MQ=255)               cGCATCAATATTAAGAATCTCTGACAGATGTAAACTTTTTCGCGCGTTATCCCTTACGCGTTCATACTTTTCAGGATGGTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGC‑GTACTGGC‑GG‑GTTTTTTTCTTT‑AACAGAtttgggt                                                                                                                                                                                                                          >  1:536743/1‑233 (MQ=255)               cGCATCAATATTAAGAATCTCTGACAGATGTAAACTTTTTCGCGCGTTATCCCTTACGCGTTCATACTTTTCAGGATGGTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGAAGGGC‑GTACTGGC‑GG‑GTTTTTTCttttatcaagttgggttt                                                                                                                                                                                                                          >  1:46497/1‑224 (MQ=255)                            aGAATCTCTGACAGATGTAAACTTTTTCGCGCGTTATCCCTTACGCGTTCATACTTTTCAGGATGGTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGC‑GTCCTGGC‑GG‑GTTTTTTTTTTA‑AAAAACATGGATTCTCTCacaaaacgg                                                                                                                                                                                                            >  2:573952/1‑231 (MQ=255)                                tctcTGACAGATGTAAACTTTTTCGCGCGTTATCCCTTACGCGTTCATACTTTTCAGGATGGTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATTGGAGCGATTCGGCTACTACGGCGTGCAGGGG‑GTACTGGC‑GG‑TTTTTTTTTTTACAAAAACTTGG‑GTGTCGCAAGAACA‑ggggtt                                                                                                                                                                                                       >  2:344755/1‑235 (MQ=255)                                 ctctGACAGATGTAAACTTTTTCGCGCGTTATCCCTTACGCGTTCATACTTTTCAGGATGGTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGC‑GTACTGGCGGG‑TTTTTTTTTTTA‑ACAACCTTGGTTTTTCGccaaacaaggggtt                                                                                                                                                                                                        >  1:723925/1‑226 (MQ=255)                                    tGACAGATCGAAACTTTTTCGCGCGTTATCCCTTACGCGTTCATACTTTTCAGGATGGTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGC‑GTACTGGC‑GG‑TTTTCTTCGTTA‑AACAGCTTGGATTCTCGCAAGAGCA‑GGCTTTTGTc                                                                                                                                                                                                   <  1:226640/239‑1 (MQ=255)                                          ttGTTAACTTTTTCGCGCTTTATCCTTTACGCGTTCATCCTTTTCAGGATGGGATGGGAAGGTTAATAAATATGAATCCAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGC‑GTACTGGC‑GG‑TTTTCTTCGTTA‑AACAGCTTGGATTCTCGCAAGAGCA‑GGCTTTTGTCACtttt                                                                                                                                                                                             <  2:944822/238‑1 (MQ=255)                                                         cgcgTTATCCCTTACGCGTTCATACTTTTCAGGGTGGTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGC‑GTACTGGC‑GG‑TTTTCTTCGTTA‑AACAGCTTGGATTCTCGCAAGAGCA‑GGCTTTTGTCACTTTTGGTGCTTTTGCTgcg                                                                                                                                                                              <  1:434977/239‑1 (MQ=255)                                                           cgttttcctgttcgggTTCATTCTTGTCAGGATGGTATTGGAAGGTTAATAAATATGAATCCAACAACCCCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGC‑GTACTGGC‑GG‑TTTTCTTCGTTA‑AACAGCTTGGATTCTCGCAAGAGCA‑GGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCt                                                                                                                                                                            <  2:947321/224‑1 (MQ=255)                                                           cgTTATCCCTTACGCGTTCATTCTTTTCAGGATGGTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGC‑GTACTGGC‑GG‑TTTTCTTCGTTA‑AACAGCTTGGATTCTCGCAAGAGCA‑GGCTTTTGTCACTTTTGGTGCTTTTGCTgcgc                                                                                                                                                                             <  1:354868/238‑1 (MQ=255)                                                            gTTATCCCTTACGGTTTCATACTTTTCAGGATGGTATTGGAAGGTTAATAAATTTGAATACAACAACTCCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGC‑GTACTGGC‑GG‑TTTTCTTCGTTA‑AACAGCTTGGATTCTCGCAAGAGCA‑GGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCt                                                                                                                                                                            <  1:460954/238‑1 (MQ=255)                                                                      aCGCGTTCATACTTTTCAGGATGGTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGC‑GTACTGGC‑GG‑TTTTTTTTGTTT‑AAAACCTTCGGTTTTCCCCAGACCAGGGGTTTTTTCCCTTTTTGGGTttttttgccgcggggctttgg                                                                                                                                                                   >  1:154498/1‑221 (MQ=255)                                                                          gTTCATACTTTTCAGGATGGTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGC‑GTACTGGC‑GGTTTTTTTTCGTTA‑AA‑AACTTTGAGTTTCCCCAAGGGA‑GGCTTTTTTTGCTTTTGTGGCtttttttttgccggttttttggctc                                                                                                                                                               >  1:890560/1‑216 (MQ=255)                                                                                               ttttgAAGGTTAATAAATATTAATACAACAACCCCCTTGGGGTTGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGC‑GTACTGGC‑GG‑TTTTCTTCGTTA‑AACAGCTTGGATTCTCGCAAGAGCA‑GGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCgg                                                                                                                                        <  2:573347/235‑1 (MQ=255)                                                                                                ttGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGC‑GTACTGGC‑GG‑TTTTCTTCGTTA‑AACAGCTTGGATTCTCGCAAGAGCA‑GGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGc                                                                                                                                       <  2:80470/239‑1 (MQ=255)                                                                                                        ttAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGC‑GTACTGGC‑GG‑TTTTTTTCTTTA‑AACAGCCTTGGTTCTCTCCAGAGAC‑GGCCTTTTTTACCTTTTGTGGTTTTTCTGCGCTGGTGTTTTGCCCTCTTTCCATTtgccgccttattgcggtacacc                                                                                                                                >  1:1066453/1‑222 (MQ=255)                                                                                                        ttAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGC‑GTACTGGC‑GG‑TTTTCTTCGTTA‑AACAGCTTGGATTCTCGCAAGAGCA‑GGCTTTTGTCACTTTTGGTGCTTTtg                                                                                                                                                                                   <  1:748450/187‑1 (MQ=255)                                                                                                           ataaattTGAATACAACAACCACTATGGTGAGTCGGCCGCATCCGTGCCCATTTTGCATGATGTTTGTTGTCGAGTTAGGGGAGCGCTTCGGCTACTACGGCGTGCAGGGC‑GTACTGGC‑GG‑TTTTCTTCGTTA‑AACAGCTTGGATTCTCGCAAGAGCA‑GGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCTGCTATGTCGGCGACCCCctgct                                                                                                                            <  2:1047575/239‑1 (MQ=255)                                                                                                                       acaacaacaCCCATGGGGATGCTGCAGCCACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGC‑GTACTGGC‑GG‑TTTTCTTCGTTA‑AACAGCTTGGATTCTCGCAAGAGCA‑GGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACg                                                                                                                <  1:573952/239‑1 (MQ=255)                                                                                                                           caaccCCCATGGGGTTTGTGCCGCAACCTCGCCCTTTTTTCATGGTCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGC‑GTACCGGC‑GG‑TTTTCTTCGTTA‑AACAGCTTGGATTCTCGCAAGAGCA‑GGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCAcc                                                                                                            <  2:376701/239‑1 (MQ=255)                                                                                                                                     ggggtttcttccacaACCTTGCCCATTTTTCATTATCTTTTTGGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGC‑GTACTGGC‑GG‑TTTTCTTCGTTA‑AACAGCTTGGATTCTCGCAAGAGCA‑GGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTg                                                                                                  <  2:890560/226‑1 (MQ=255)                                                                                                                                           ggtgcAGCAACCTGGCCCATTTTTCATGATTTTTTTTGGCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGC‑GTACTGGC‑GG‑TTTTCTTCGTTA‑AACAGCTTGGATTCTCGCAAGAGCA‑GGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCAc                                                                                            <  1:117412/237‑1 (MQ=255)                                                                                                                                                 taaacccctcctgtttttCATGATTGTTTTTTTCGGTTTAGGGGGGAGATTCGGCTACTACGGCGTGAAGGGC‑GTTCTGGC‑GG‑TTTTCTTCGTTA‑AACAGCTTCGATTCGCGCAAGAGCA‑GGCTTTTGTATCTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCCCCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGc                                                                                      <  2:951451/226‑1 (MQ=255)                                                                                                                                                           cccttttttttgtttttttttttCGAGTTATGGGGGTGATTTGGCTTCCTCGGGGTGCAGGGG‑TGACTGGT‑GG‑TTTTCTTCGGTA‑AACAGCTTGGATTCTCGGAAGAGCA‑GGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTGGCGATTgg                                                                            <  2:1097023/224‑1 (MQ=255)                                                                                                                                                            cttttttCTTGTTCTTTTTTGTGGTGTGTTGGGGGCGATTCGGCTACTGCGGCCTGCGGGGA‑GTGCTGGC‑GG‑TTTTCTTCGTTA‑AACAGCTTGGATTCTCGCAAGAGCA‑GGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGTTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTGGCGATTGGc                                                                           <  2:536743/237‑1 (MQ=255)                                                                                                                                                                                 ttgggTTATGGGGGCGTTTTGGCTTCTACGGGGTGGCGGGCGGTATTGG‑‑GG‑TTGTCTTCGTTA‑AACAGCTTGGATTCCCGCAAGAGCA‑GGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTGGCGATTGGCTACTTCATGACCGGCATGTCg                                                      <  2:51499/235‑1 (MQ=255)                                                                                                                                                                                 tgtgtttttgggtgggATTTGGTCACTACGTGGGGCGGGGC‑GTTCTGGC‑GG‑TTTTCTTTGTTA‑AACTGCTTGGATTCTCTCAAGACCA‑TGCTTTTGTCTCTTTTGGTGCTTTTGCTGTGCTGGTCGATGGCCTCATTTCCATAGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTGGCGATTGGCTACTTCATGACCGGCATGTCg                                                      <  2:771540/231‑1 (MQ=255)                                                                                                                                                                                                ttttCGGCTTCTGCGGCTGGCAGGGG‑GTACTGGC‑GG‑TTTTTTTCGTTA‑AACAGCTTGGATTCTCGCAAGAGCA‑GGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTGGCGATTGGCTACTTCATGACCGGCATGTCGCTACTTAAGcctgac                                       <  2:13586/237‑1 (MQ=255)                                                                                                                                                                                                 gtttggCTACTTCGGCTGTCGGGGC‑GTACTGGC‑TG‑TTTTCGTCGTTA‑AACCGCTTGGATTTTGTCAAGAGCA‑GGGTTTTGTCGCTTTGGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCAGTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTGGCGATTGGCTACTTCATGACCGGCATGTCGCTACTTAAGcctgac                                       <  2:154498/234‑1 (MQ=255)                                                                                                                                                                                                    cGGCTACTACGGCGTGCAGGGC‑GTACTGGC‑GG‑TTTTCTTCGTTA‑AACAGCTTGGATTCTCGCAAGAGCA‑GGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTGGCGATTGGCTACTTCATGACCGGCATGTCGCTACTTAAGcctgacctga                                   <  1:334997/239‑1 (MQ=255)                                                                                                                                                                                                     tgcttcctgggggggcggggC‑GTACTGGT‑GG‑TTTTCTTCTGTA‑AACGTGTTGGATTCTCGCAAGCGCA‑GGCTTTTGTCACTTGTGGTTCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATGTGCGGCTATGTCGGCGACCCCCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTGGCGATTGGCTACTTCATGACCGGCATGTCGCTACTTAAGCCTGACCTGAt                                  <  2:1096954/222‑1 (MQ=255)                                                                                                                                                                                                     ggCCACTACGGCGTGCAGGGC‑GTTCTGGC‑GG‑TTTTCTTCGTTA‑AACAGTTTGGATTCGCGCAAGAGCA‑GGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTAGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTGGCGATTGGCTACCTCATGACCGGCATGTCGCTACTTAAGCCTGACCTGAt                                  <  1:512602/239‑1 (MQ=255)                                                                                                                                                                                                          cttggcggTGCAGGGG‑CTACTGGG‑CT‑GTTTCTTCGTTA‑ACCTGCTTGGATTCTCGAGAGGGTA‑GGTTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCTTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTGGCGATTGGCTACTTCATGACCGGCATGTCGCTACTTAAGCCTGACCTGATTTTCa                             <  2:46497/232‑1 (MQ=255)                                                                                                                                                                                                          cttcGGCGTGGCGGGG‑TTTCTTGC‑GG‑TTTTCTTCGTTA‑AACAGGTTGGATTCTCGCAAGAGCA‑GGCTTTTGTCCCTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCAGTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTGGCGATTGGCTACTTCATGACCGGCATGTCGCTACTTAAGCCTGACCTGATTTTCa                             <  2:1066453/236‑1 (MQ=255)                                                                                                                                                                                                            tcGGCGTGCAGGGC‑GTACTGGC‑GG‑TTTTCTTCGTTA‑AACAGCTTGGATTCTCGCAAGAGCA‑GGCTTTTGTCCCTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTGGCGATTGGCTACTTCATGACCGGCATGTCGCTACTTAAGCCTGACCTGATTTTCATc                           <  1:768499/238‑1 (MQ=255)                                                                                                                                                                                                                ggggCAGGGC‑GGACGGGC‑GG‑TTTGCTGCGTTA‑AACAGCTTGGATTCTCGCAAGAGCA‑GGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTGGCGATTGGCTACTTCATGACCGGCATGTCGCTACTTAAGCCTGACCTGATTTTCATCGcc                        <  1:189173/235‑1 (MQ=255)                                                                                                                                                                                                                  gcgggggT‑GTGCTGGT‑GG‑TTTTTTTCGTTA‑AGCTGGTTGGAGTCTCGCAAGAGCA‑GGCTTTTGTCTCTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTGGCGATTGGCTACTTCATGACCGGCATGTCGCTACTTAAGCCTGACCTGATTTTCATCGCCCTg                     <  2:723925/235‑1 (MQ=255)                                                                                                                                                                                                                              tctgtgtg‑TTTTCTTTGTTA‑AACAGCTTGGATTCTCGCAAGGGCA‑GGCTTTTGTCATTTTGGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTGGCGATTGGCTACTTCATGACCGGCATGTCGCTACTTAAGCCTGACCTGATTTTCATCGCCCTGGGGACTATCg           <  2:1033194/232‑1 (MQ=255)                                                                                                                                                                                                                                 gC‑GT‑TTTTTTTTGTTA‑AACGGCTTGGATTCTCGCAGGAGCA‑GGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTGCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAGCGCACCATTGTTCTTGGAGCACTTGTGCTGGCGATTGGCTACTTCATGACCGGCATGTCGCTACTTAAGCCTGACCTGATTTTCATCGCCCTGGGGACTATCGCTGt       <  2:269549/239‑1 (MQ=255)                                                                                                                                                                                                                                  c‑GG‑TTTTCTTCGTTA‑AACAGCTTGGATTCTCGCAAGAGCA‑GGCTTTTGTCCCTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTGGCGATTGGCTACTTCATGACCGGCATGTCGCTACTTAAGCCTGACCTGATTTTCATCGCCCTGGGGACTATCGCTGt       <  1:292496/238‑1 (MQ=255)                                                                                                                                                                                                                                       ttttCTTCGTTA‑AACAGCTTGGATTCTCGCAAGAGCA‑GGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTGGCGATTGGCTACTTCATGACCGGCATGTCg                                                      <  1:221712/188‑1 (MQ=255)                                                                                                                                                                                                                                        tttCTTGCTTA‑AACCGCTTGGATTCTCGCAAGAGCA‑GGCTTTTGTCACTTTTGGGGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTGGCGATTGGCTACTTCATGACCGGCATGTCGCTACTTAAGCCTGACCTGATTTTCATCGCCCTGGGGACTATCGCTGTCGGTa  <  1:344755/239‑1 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                           ATTAATCTATTCACCGCATCAATATTAAGAATCTCTGACAGATGTAAACTTTTTCGCGCGTTATCCCTTACGCGTTCATACTTTTCAGGATGGTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGC‑GTACTGGC‑GG‑TTTTCTTCGTTA‑AACAGCTTGGATTCTCGCAAGAGCA‑GGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTGGCGATTGGCTACTTCATGACCGGCATGTCGCTACTTAAGCCTGACCTGATTTTCATCGCCCTGGGGACTATCGCTGTCGGTA  >  NC_000913/3640749‑3641215

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