Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,640,983	C→T	14.7%	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	14.7%	111.8 / 4.3	57	F41F (TTC→TTT)	dtpB	dipeptide and tripeptide permease B
Reads supporting (aligned to +/- strand):  ref base C (9/39);  new base T (7/2);  total (16/41)
Fisher's exact test for biased strand distribution p-value = 1.10e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTAATCTATTCACCGCATCAATATTAAGAATCTCTGACAGATGTAAACTTTTTCGCGCGTTATCCCTTACGCGTTCATACTTTTCAGGATGGTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGC‑TGCAGCAACC‑TCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTC‑GTT‑AAACA‑GC‑TTGGAT‑TCTCG‑CAAGAGCAGGCT‑TTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCT‑ATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTGGCGATTGGCTACTTCATGACCGGCATGTCGCTACTTAAGCCTGACCTGATTTTCATCGCCCTGGGGACTATCGCTGTCGGTAACGGC  >  NC_000913/3640750‑3641220                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                      ttAATCTATTCACCGCATCAATATTAAGAATCTCTGACAGATGTAAACTTTTTCGCGCGTTATCCCTTACGCGTTCATACTTTTCAGGATGGTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATAC‑TGCAGCAACC‑TCGCCCATTTTTAATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCCGGGCGTTATGGC‑GTGTTTTTTC‑Gttta                                                                                                                                                                                                                                                 >  1:899009/1‑237 (MQ=255)    aTCTATTCACCGCATCAATATTAAGAATCTCTGACAGATGTAAACTTTTTCGCGCGTTATCCCTTACGCGTTCATACTTTTCAGGATGGTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGC‑TGCAGCAACC‑TCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGGTTTTTTTT‑GTT‑aaaa                                                                                                                                                                                                                                              >  2:666607/1‑238 (MQ=255)       tATTCACCGCATCAATATTAAGAATCTCTGACAGATGTAAACTTTTTCGCGCGTTATCCCTTACGCGTTCATACTTTTCAGGATGGTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGC‑TGCAGCAACC‑TCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTCTGTT‑AAACaca                                                                                                                                                                                                                                           >  2:969594/1‑237 (MQ=255)         ttCACCGCATCAATATTAAGAATCTCTGACAGATGTAAACTTTTTCGCGCGTTATCCCTTACGCGTTCATACTTTTCAGGATGGTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGC‑TGCAGCAACC‑TCGCCCATTTTTCATGATCTTTTTTTTCGAGTTATGGGAGCGCTTCGGCTACTACGGCGTGCAGGGCGTAATGGC‑GGGTATCTTT‑Gttgaccaacctg                                                                                                                                                                                                                                         >  2:918073/1‑229 (MQ=255)         ttCACCGCATCAATATTAAGAATCTCTGACAGATGTAAACTTTTTCGCGCGTTATCCCTTACGCGTTCATACTTTTCAGGATGGTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGC‑TGCAGCAACC‑TCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTT‑TTT‑TAACA‑AC‑ttc                                                                                                                                                                                                                                      >  2:977198/1‑238 (MQ=255)         ttCACCGCATCAATATTAAGAATCTCTGACAGATGTAAACTTTTTCGCGCGTTATCCCTTACGCGTTCATACTTTTCAGGATGGTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGC‑TGCAGCAACC‑TCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTT‑GTT‑aaaaaacctg                                                                                                                                                                                                                                        >  2:529230/1‑232 (MQ=255)         ttCACCGCATCAATATTAAGAATCTCTGACAGATGTAAACTTTTTCGCGCGTTATCCCTTACGCGTTCATACTTTTCAGGATGGTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGC‑TGCAGCAACC‑TCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTAGGCTATTACTGCGTGAAGGGCGTACTGTC‑GGGTTTTTTC‑GTT‑aaaacaacag                                                                                                                                                                                                                                        >  2:348769/1‑232 (MQ=255)              cGCATCAATATTAAGAATCTCTGACAGATGTAAACTTTTTCGCGCGTTATCCCTTACGCGTTCATACTTTTCAGGATGGTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGC‑TGCAGCAACC‑TCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTCAAGTGCGAACTGGT‑GGTTTTTTTC‑Gtcttatagagcgcagtt                                                                                                                                                                                                                                    >  2:887205/1‑223 (MQ=255)                 aTCAATATTAAGAATCTCTGACAGATGTAAACTTTTTCGCGCGTTATCCCTTACGCGTTCATACTTTTCAGGATGGTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGC‑TGCAGCAACC‑TCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTC‑GTT‑aaa                                                                                                                                                                                                                                               <  2:878796/224‑1 (MQ=255)                        tttacaaTCTCTTACGGTTGTAAGTTTTTTCGCGCCTTTACCCTTTCGCTTTCATACTTTTCAGGATTGTATTGGAAGGTTAAGGAATATGAATTCAACACCCCCCATGTGGATGC‑TGCAGCAACC‑TAGCCCATTTTCCATGATCTTTTTTGTCGAGTTATGGGAGCTATTCGTCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTC‑GTT‑AAACA‑GC‑TTGGAT‑TTTCG‑CCAgcgc                                                                                                                                                                                                                     <  1:938417/234‑1 (MQ=255)                        ttAAGAATCTCTGACAGATGTAAACTTTTTCGCGCGTTATCCCTTACGCGTTCATACTTTTCAGGATGGTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGC‑TGCAGCAACC‑TCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGCCGTGCAGGGCGAACTGGC‑GGGGTTTTTC‑GTT‑TAAAA‑CC‑TTGGGT‑tctaccaaccac                                                                                                                                                                                                                      >  2:571168/1‑230 (MQ=255)                           aGAATCTCTGACAGATGTAAACTTTTTCGCGCGTTATCCCTTACGCGTTCATACTTTTCAGGATGGTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGC‑TGCAGCAACC‑TCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGGTTTTTTT‑GTT‑GAAAA‑AC‑ATGGGT‑TTTTG‑CCAAAGCAgg                                                                                                                                                                                                                  >  2:630340/1‑239 (MQ=255)                           aGAATCTCTGACAGATGTAAACTTTTTCGCGCGTTATCCCTTACGCGTTCATACTTTTCAGGATGGTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGC‑TGCAGCAACC‑TCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGAAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTC‑TTT‑TAATACGC‑TTggtgtttttctttttaac                                                                                                                                                                                                                     >  1:878796/1‑223 (MQ=255)                                         aTGTAAACTTTTTCGCGCGTTATCCCTTACGCGTTCATACTTTTCAGGATGGTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGC‑TGCAGCAACC‑TCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTC‑GTT‑TAACA‑AC‑TTTGGT‑TTTCG‑CAAAAGACGGAG‑TTTGTCAccctt                                                                                                                                                                                                   >  2:263422/1‑235 (MQ=255)                                                   ttttGCGTGTTTTCTCTTACGTTTTTAAACTTTTCAGGATGGTATTGGAATGTTAATAAATATGAATACAACAACACCCATGGGGATGC‑TGCAGCAACC‑TCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTC‑GTT‑AAACA‑GC‑TTGGAT‑TCTCG‑CAAGAGCAGGCT‑TTTGTCACTTTTGGTGCTTTtg                                                                                                                                                                                         <  1:1202589/239‑1 (MQ=255)                                                   tttCGCGCGTTATCCCTTACGCGTTCATACTTTTCAGGATGGTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGC‑TGCAGCAACC‑TCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTC‑GTT‑AAACA‑GC‑TTGGAT‑TCTCG‑CAAGAGCAGGCT‑TTTGTCACTTTTGGTGCTTTtg                                                                                                                                                                                         <  1:294070/239‑1 (MQ=255)                                                        acggtttgcccTTACGCGTTGCAACTTTTCAGGGTGGTATTGGGAGGTTAATAAATATGAATACACAAACCCCCATGGGGATGC‑TGCAGCAACC‑TCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTC‑GTT‑AAACA‑GC‑TTGGAT‑TCTCG‑CAAGAGCAGGCT‑TTTGTCACTTTTGGTGCTTTTGCTgcg                                                                                                                                                                                    <  1:1050608/231‑1 (MQ=255)                                                                tcctttcggcttCCTATTTTTCAGGATGGGTTTGGAATGTTAATAAATATGAATACAACAACCCCCTGGGGGATGC‑TGCAGGCACC‑TCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTTCTGGC‑GGTTTTCTTC‑GTT‑AAACA‑GC‑TTGGAT‑TCTCG‑CAAGAGCAGGCT‑TTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTa                                                                                                                                                                            <  1:1120185/229‑1 (MQ=255)                                                                  cTTACGCGTTCATACTTTTCAGGATGGTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGC‑TGCAGCAACC‑TCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTC‑GTT‑AAACA‑GC‑TTGGAT‑TCTCG‑CAAGAGCAGGCT‑TTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATg                                                                                                                                                                          <  2:592243/239‑1 (MQ=255)                                                                     aCGCGGTCATACTTTTCAGGATGGTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGC‑TGCAGCAAAC‑TCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGGTACTACGGCGTGTCGGGCGTACTGGC‑GGGTTTTTTC‑TTT‑TACAA‑GC‑CTGGGG‑TGTCGCCAAACACAGG‑T‑TTTGTCACTTTTGTTGGtttttgttcgtttgtgcttgag                                                                                                                                                                        >  1:893618/1‑220 (MQ=255)                                                                            cATACTTTTCAGGATGGTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGC‑TGCAGCAACC‑TCGCCCATTTTTCATGATATTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGTCGTGCAGGGCGTACTGGC‑GGTTTTTTTT‑GTT‑TAACA‑ACATTGTGT‑TTTCG‑CAAGAGCAGCTT‑TTTTTCATTTTTTG‑GGTTTTGTTTCCTTGGGCTTTGGCCTCCtttt                                                                                                                                                                >  1:1269462/1‑238 (MQ=255)                                                                               aCTTTTCAGGATGGTATTGGAAGGTTAATAAATATGAATACAACAACCCCCATGGGGATGC‑TGCAGCAACC‑TCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTC‑GTT‑AAACA‑GC‑TTGGAT‑TCTCG‑CAAGAGCAGGCT‑TTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCAt                                                                                                                                                             <  2:1126571/239‑1 (MQ=255)                                                                                cTTTTCAGGATGGTTTTGGAAGGTTAAAAAATTTGAATACAACCACCCCCATGGGGATGC‑TGCAGCAACC‑TCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTC‑GTT‑AAACA‑GC‑TTGGAT‑TCTCG‑CAAGAGCAGGCT‑TTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCAtt                                                                                                                                                            <  1:725239/239‑1 (MQ=255)                                                                                 ttttCAGGATGGTATTGGAAGGTTAATAAATATGAATTCAACAACACCCATGGGGATGC‑TGCAGCAACC‑TCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTC‑GTT‑AAACA‑GC‑TTGGAT‑TCTCG‑CAAGAGCAGGCT‑TTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCAtt                                                                                                                                                            <  2:487893/238‑1 (MQ=255)                                                                                               ttGGAAGGTTAATAAATAGTAATTCAACCACTCCCATGGGGATGC‑TGCAGCAACC‑TCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTTCGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTC‑GTT‑AAACA‑GC‑TTGGAT‑TCTCG‑CAAGAGCAGGCT‑TTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCT‑ATGTCGGc                                                                                                                                            <  2:44214/239‑1 (MQ=255)                                                                                                   gaGGTTAATAAATATTATTACAACAACCCCCATGGGGATGT‑TGCAGTAACCTTCGCCC‑TTTTTCAGGGTCTTTTTTGTCGAGTTATGCGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTC‑GTT‑AAACA‑GC‑TTGGAT‑TTTCG‑CAAGAGCAGGCC‑TTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCT‑ATGTCGGCGacc                                                                                                                                        <  1:516470/238‑1 (MQ=255)                                                                                                    aGGTTAATAAATATGAATACAACAACACCCATGGGGATGC‑TGCAGCAACC‑TCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTC‑GTT‑AAACA‑GC‑TTGGAT‑TCTCG‑CAAGAGCAGGCT‑TTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCT‑ATGTCGGCGacc                                                                                                                                        <  2:381125/238‑1 (MQ=255)                                                                                                       ttAATAAATATGAATACAACAACACCCATGGGGATGC‑TGCAGCAACC‑TCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTT‑GTT‑TAAAA‑CC‑TTTGATGTCT‑G‑CCAAAGCAGGGT‑TTTTGTACTTTTTTTGGCTTTTCTTGGGTGGTCGATGGTCTCACTTTCCTTtgcggcggtttctgccgcgcac                                                                                                                                      >  2:1279607/1‑222 (MQ=255)                                                                                                          ttaaattTTAATACAACAACCCTCATGGGGATGC‑TGCAGCAACC‑TCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTC‑GTT‑AAACA‑GC‑TTTGAT‑TCTCG‑CAAGAGCAGGCT‑TTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCT‑ATGTCGGCGACCACctgct                                                                                                                                 <  1:586893/238‑1 (MQ=255)                                                                                                                      acaatcgcacCCCTGGGGGTGCGTGC‑GCAACC‑TCCCCCATTTTTCAGTATCTTTTTTGTCGATTTATGGGTGCGATTCGGCTACTACGGCGTGCAGTGCGTTCTGGC‑GGTTTTCTTC‑GTT‑AAACA‑GC‑TTGGAT‑TCTCG‑CAAGAGCAGGCT‑TTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCT‑ATGTCGGCGACCACCTGCTGGGGACCAAACg                                                                                                                     <  1:832278/232‑1 (MQ=255)                                                                                                                      acaacatgccccTTGGGGAGGC‑TGCAGCAACC‑GCGCCCATTTTTCATGATCTTGTTTGGCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTC‑GTT‑AAACA‑GC‑TTGGAT‑TCTCG‑CAAGAGCAGGCT‑TTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCT‑ATGTCGGCGACCCCCTGCTGGGGACCAACCg                                                                                                                     <  1:1135640/230‑1 (MQ=255)                                                                                                                          caacaCCCATGGGGATGC‑TGCAGCAACC‑TCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTC‑GTT‑AAACA‑GC‑TTGGAT‑TCTCG‑CAAGAGCAGGCT‑TTTGTCACTTTTGGTGCTTTTGCTgcg                                                                                                                                                                                    <  2:1053474/173‑1 (MQ=255)                                                                                                                          caacaCCCATGGGGATGC‑TGCAGCAACC‑TCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGGTTTTTTC‑TTTAAAACA‑G‑‑TTTGAT‑TTTTG‑CCAAAGCAGGGC‑TTTTTTTCTTTTTGTGGTTTTTCTGCGCTGGGCTATTGGCTCCTTTTCATTTGCGGGTCTTTTTGGGGGCCCCCTCCTGGGGGCCAAAcccc                                                                                                                   >  2:260292/1‑235 (MQ=255)                                                                                                                               cccATGGGGATGC‑CGCAGCAACC‑TCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTC‑GTT‑AAACA‑GC‑TTGGAT‑TCTCG‑CAAGAGCAGGCT‑TTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCT‑ATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGt                                                                                                            <  2:857356/239‑1 (MQ=255)                                                                                                                                    ggggATGC‑TGCAGCAACC‑CTGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTC‑GTT‑AAACA‑GC‑TTGGAT‑TCTCG‑CAAGAGCAGGCT‑TTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCT‑ATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTg                                                                                                       <  2:579155/239‑1 (MQ=255)                                                                                                                                      ggATGC‑TGCAGCAACC‑TCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTC‑GTT‑AAACA‑GC‑TTGGAT‑TCTCG‑CAAGAGCAGGCT‑TTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCT‑ATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTgg                                                                                                      <  2:883745/238‑1 (MQ=255)                                                                                                                                                  gtaCC‑TCCCCCATTTTTCTTTTTTTTTTTTGTCGAGTTATGGGGGCGATTCGGCTTCTTCGGCGTGCAGGGCGTACTGGC‑GGTTTGCTTC‑GTT‑AAACA‑GC‑TTGGAT‑TCTCG‑CAAGAGCAGGCT‑TTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCT‑ATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCt                                                                                          <  1:343699/237‑1 (MQ=255)                                                                                                                                                     cc‑TCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTC‑GTT‑AAACA‑GC‑TTGGAt‑t                                                                                                                                                                                                                                 <  2:365155/104‑1 (MQ=255)                                                                                                                                                        tcttccctttttcgttgttttttttGTCGATTTAGGTGGGCGATTGGGGTTCTCCGGCTGGCATGGCGTACTGGC‑TGTTTTCTTC‑GTT‑AAACA‑GC‑GTGGAT‑TCTCG‑CAAGAGCAGGCT‑TTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTTTATGGCCTCAGTTCCATTGGCGGCT‑ATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTGGCGa                                                                                     <  1:666607/220‑1 (MQ=255)                                                                                                                                                        tCGTCCTTTTTTCATGGTTTTTTTTGTCGAGTTGTGGGAGCGATTCGGCTACTCCGGCGTGCAGGGCGTACTGGC‑GGTTTTCGTC‑GTT‑AAACA‑GC‑TTGGAT‑TCTCG‑CAAGAGCAGGCT‑TTTGTCACTTTTGGTGCTTTGGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCT‑ATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTGGCGa                                                                                     <  1:147749/239‑1 (MQ=255)                                                                                                                                                             tttttttcgttgtctgtgtttgtgAGTTTGGGGACCGATTCGGGCGTTACGGCGTGCCGGGCCTACGGGC‑AGTTTTTGTT‑GGT‑AAACA‑GC‑TTCGAT‑TCTCG‑CAACATCAGGCTGTTTGTC‑CTGTTCGTGCTTTTGCTGCGCTTGTCTTTGTCATCAGTTCCATTGGCGGTT‑ATGTCGACGACCACCTGCTGGGGACCAAACGCAACATTGTTCTTGGAGCACTTGTGCTGGCGATTGGc                                                                                <  1:887205/216‑1 (MQ=255)                                                                                                                                                             cttttttCATTTTTGTTTTTGTCGGTGTATGGGAGCGATTCGGAGTCTATGGCGGTCAGGGCGTACTGGC‑TGTTTTCTTC‑GTT‑AAACA‑GT‑TTGGTT‑TCTCG‑CAAGAGCAGGCT‑TTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCAGTTCCATTGGCGGTT‑ATGTCGGCGTCCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTGGCGATTGGc                                                                                <  1:571168/237‑1 (MQ=255)                                                                                                                                                             cATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTC‑GTT‑AAACA‑GC‑TTGGAT‑TCTCG‑CAAGAGCAGGCT‑TTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCT‑ATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTGGCGATTGGc                                                                                <  2:538717/239‑1 (MQ=255)                                                                                                                                                             cATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTC‑GTT‑AAACA‑GC‑TTGGAT‑TCTCG‑CAAGAGCAGGCT‑TTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCT‑ATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTGGCGATTGGc                                                                                <  2:203599/239‑1 (MQ=255)                                                                                                                                                                           ttttttttCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTC‑GTT‑AAACA‑GC‑TTGGAT‑TCTCG‑CAAGAGCAGGCT‑TTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTACGGCCTCATTTCCATTGGCGGCT‑ATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTGGCGATTGGCTACTTCATGACCgg                                                                  <  2:837076/239‑1 (MQ=255)                                                                                                                                                                           ttttttGTGCGGTTATGGGGGGGTTTCCGCCACCACGGGGTGCAGGGCGGTGTGGG‑TGTTTTATTC‑GTT‑ACACT‑GT‑TGGGAA‑TCTTG‑CAAGAGCGAGGT‑TTTGTCTCTTTGGGTGCTTTTGCTGCGCTGGTTTATGGCCTCATTTCCATTTGCGGTT‑ATGTCGGCGGCCACCTGCTGGGGACCAAACGCACCAGTGTTCGTGGAGCACTTGTGCTGGCGATTGGCTAATTCATGACCgg                                                                  <  1:386772/239‑1 (MQ=255)                                                                                                                                                                                  tCGTGTTATGGGGGCTATTCGGCTACGTCGGCGGGCAGGGCGTACTGGC‑GGTTTTCTTC‑GTT‑AAACA‑GC‑TTGGAT‑TCTCG‑CAAGAGCAGGCT‑TTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCT‑ATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTGGCGATTGGCTACTTCATGACCGGCATGTCg                                                           <  1:260292/239‑1 (MQ=255)                                                                                                                                                                                             gggtggttCGGCTCTTACGGCGTGCAGGGCTTACTGGC‑GTTTTTCTTC‑GTT‑AAACA‑GT‑TTGGAT‑TCGCG‑CAAGAGCAGGCT‑TTTTTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCT‑ATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTGGCGATTGGCTACTTCATGACCGGCATGTCGCTACTTAAGc                                                 <  1:103897/232‑1 (MQ=255)                                                                                                                                                                                                cGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTC‑GTT‑AAACA‑GC‑TTGGAT‑TCTCG‑CAAGAGCAGGCT‑TTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCT‑ATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTGGCGATTGGCTACTTCATGACCGGCATGTCGCTACTTAAGcctg                                              <  2:92089/238‑1 (MQ=255)                                                                                                                                                                                                 gcttctgcttcttgggggtggtgggcgggcTGGG‑GGTTTTCTTC‑GTT‑AAACC‑GT‑TTGGTG‑TCTCG‑CAAGGGCCGGTT‑TTTGTCACTTTTGGTGCTTTTGCTGGGCTGGTCTGGGGCCTCATTTCCATTGGCGGCT‑ATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTGGCGATTGGCTACTTCATGACCGGCATGTCGCTACTTAAGcctga                                             <  1:1057335/216‑1 (MQ=255)                                                                                                                                                                                                           ctacGGCGTGCAGGGCTGACTGGC‑GGTTTTCTTC‑GTT‑AAACA‑GC‑TTGGAT‑TCTCG‑CAAGAGCAGGCT‑TTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCT‑ATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTGGCGATTGGCTACTTCATGACCGGCATGTCGCTACTTAAGCCTGACCTGATTTGCa                                  <  2:599979/239‑1 (MQ=255)                                                                                                                                                                                                                  gTGCAGGGCGTACTGGC‑GGTTTTCTTC‑GTT‑AAACA‑GC‑TTGGAT‑TCTCG‑CAAGAGCAGGCT‑TTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCT‑ATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTGGCGATTGGCTACTTCATGACCGGCATGTCGCTACTTAAGCCTGACCTGATTTTCATCGccc                            <  2:1088136/238‑1 (MQ=255)                                                                                                                                                                                                                     gcggggttttttgggttttttCTTT‑GTT‑AAACG‑GC‑TGGGGT‑TTTCG‑CAAGGGCGGGCT‑TTTGGTACTTTTGGTGCTTTTGCCGCGCGGGTCTATGGCCTCATTTCCATTGGCGGCT‑ATGTCGGCGACCGCCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTGGCGATTGGCTACTTCATGACCGGCATGTCGCTACTTAAGCCTGACCTGATTTTCATCGCCCTg                          <  1:630340/221‑1 (MQ=255)                                                                                                                                                                                                                      aGGGCGTACTGGC‑GGTTTTCTTC‑GTT‑AAACA‑GC‑TTGGAT‑TCTCG‑CAAGAGCAGGCT‑TTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCT‑ATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTGGCGATTGGCTACTTCATGACCGGCATGTCGCTACTTAAGCCTGACCTGATTTTCATCGCCCTggg                        <  2:212887/238‑1 (MQ=255)                                                                                                                                                                                                                       agggcttgtttgtttttttCTTC‑GTT‑AAACA‑GC‑TGGGGT‑TTCCG‑AAGTGGCAGGCT‑TTTTGCATTGTTGGTGCTTTTGCTGCGCTGGTCTAGGGCCGCATTTCCATTGGCGGTT‑ATGGCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTGCTTGGAGCACTTGTGCTGGCGATTGGCTACTTCATGACCGGCATGTCGCTACTTAAGCCTGACCTGATTTTCATCGCCCTggg                        <  1:479896/223‑1 (MQ=255)                                                                                                                                                                                                                                  c‑GGTTTTCTGC‑GTT‑AAACA‑GC‑TTGGAT‑TCTCG‑CAAGAGCAGGCT‑TTTGTCCCTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCT‑ATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTGGCGATTGGCTACTTCATGACCGGCATGTCGCTACTTAAGCCTGACCTGATTTTCATCGCCCTGGGGACTATCGCTGt            <  2:251286/238‑1 (MQ=255)                                                                                                                                                                                                                                       tttCTTC‑GTT‑AAACA‑GC‑TTGGAT‑TCTCG‑CAAGAGCAGGCT‑TTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCT‑ATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTGGCGATTGGCTACTTCATGACCGGCATGTCGCTACTTAAGCCTGACCTGATTTTCATCGCCCTGGGGACTATCGCTGTCGGTa       <  2:1122914/239‑1 (MQ=255)                                                                                                                                                                                                                                            tc‑GTT‑AAATA‑GC‑TTGGAT‑TCTCG‑CAAGAGCAGGTT‑TTTGTCATTTTTGGTGCTTTTGCGGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCT‑ATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCGCTTGTGCTGGCGATTGGCTACTTCATGACCGGCATGTCGCTACTTAAGCCTGACCTGATTTTCATCGCCCTGGGGACTATCGCTGTCGGTAACGGc  <  2:319439/239‑1 (MQ=255)                                                                                                                                                                                                                                            tc‑GTT‑AAACA‑GC‑TTTGAT‑TCTCG‑CAAGAGCAGGCG‑TTTGGCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCT‑ATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTGGCGATTGGCTTCTTCATGACCGGCATGTCGCTACTTAAGCCTGACCTGATTTTCATCGCCCTGGGGACTATCGCTGTCGGTAACGGc  <  2:150793/239‑1 (MQ=255)                                                                                                                                                                                                                                            tc‑GTT‑AAACA‑GC‑TTGGAT‑TCTCG‑CAAGAGCATGCT‑TTTGTCTCTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCT‑ATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTGGCGATTGGCTACTTCATGACCGGCATGTCGCTACTTAAGCCTGACCTGATTTTCATCGCCCTGGGGACTATCGCTGTCGGTAACGGc  <  2:899009/239‑1 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                      TTAATCTATTCACCGCATCAATATTAAGAATCTCTGACAGATGTAAACTTTTTCGCGCGTTATCCCTTACGCGTTCATACTTTTCAGGATGGTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGC‑TGCAGCAACC‑TCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTC‑GTT‑AAACA‑GC‑TTGGAT‑TCTCG‑CAAGAGCAGGCT‑TTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCT‑ATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTGGCGATTGGCTACTTCATGACCGGCATGTCGCTACTTAAGCCTGACCTGATTTTCATCGCCCTGGGGACTATCGCTGTCGGTAACGGC  >  NC_000913/3640750‑3641220

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