Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	602,635	A→G	12.8%	E226G (GAG→GGG)	pheP →	phenylalanine transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	602,635	0	A	G	12.8%	175.7 / 2.1	69	E226G (GAG→GGG)	pheP	phenylalanine transporter
Reads supporting (aligned to +/- strand):  ref base A (15/44);  new base G (6/4);  total (21/48)
Fisher's exact test for biased strand distribution p-value = 5.69e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

TAACCTGGTGAACGTGCGCTTATATGGCGAAACCGAGTTCTGGTTTGCGTTGATTAAAGTGCTGGCAATCATCGGTATGATCGGCTTTGGCCTGTGGCTGCTGTTTTCTGGTCACGGCGGCGAGAAAGCCAGTATCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGAT‑TTTGTCGCTGGC‑GGTAATTATGTTC‑TCCTTCGGCGGT‑CT‑GGA‑GCTGATTGGGATTACTGCCGCTGAAG‑CGCGCGATCCGGAAAAAAGCA‑TTCCAAAAGCGGTAAATCAGGTGGTGTATC‑GCATCCTGCT‑GTTTTACATCGGTTCACT‑GGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGT‑AGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAGCTTCTGCGCTGAACTTCGTCA  >  NC_000913/602405‑602856                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                     tAACCTGGTGAACGTGCGTTTATATGGCGAAACCGAGTTCTGGTTTGCGTTGATTAAAGTGCTGGCAATCATCGGTATGATCGGCTTTGGCCTGTGGCTGCTGTTTTCTGGTCACGGCGGCGAGAAAGCCAGTATCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGAT‑TTTGTCGCTGGC‑GGTAATTATGTTC‑TCCTTCGGCGGT‑CT‑GGA‑GCTGAtt                                                                                                                                                                                                                              <  1:629229/238‑1 (MQ=255)        gTGAACGTGCGCTTATATGGCGAAACCGAGTTCTGGTTTGCGTTGATTAAAGTGCTGGCAATCATCGGTATGATCGGCTTTGGCCTGTGGCTGCTGTTTTCTGGTCACGGCGGCGAGAAAGCCAGTATCGACAACCTCTGGCGCTCCGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGAT‑TTTGTCGCTGGG‑GGTGATTATGTTT‑TCCTTTGGCGGG‑CT‑TGG‑GCCGATTGGGGTTat                                                                                                                                                                                                                      >  1:750702/1‑238 (MQ=255)        gTGAACGTGCGCTTATATGGCGAAACCGAGTTCTGGTTTGCGTTGATTAAAGTGCTGGCAATCATCGGTATGATCGGCTTTGGCCTGTGGCTGCTGTTTTCTGGTCACGGCGGCGAGAAAGCCAGTATCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGGT‑TTTGTCGCTGGC‑GGTAATTTTGTTC‑TCCTTCGGCGGT‑CT‑GGG‑GCTGGTTTGGGtt                                                                                                                                                                                                                        >  1:227248/1‑237 (MQ=255)         tGACCGTGCGCTTTGATGGGGAAACCGAGTTCTGGTTTGCGTTGATTAAAGTGCTGGCAATCATCGGTATGATCGGCTTTGGCCTGTGGCTGCTGTGTGCTGGTCACGGCGGCGAGAAAGCCAGTATCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGAT‑TTTGTCGCTGGC‑GGTAATTATGTTC‑TCCTTCGGCGGT‑CT‑GGA‑GCTGATTGGGATTAc                                                                                                                                                                                                                      <  2:227248/238‑1 (MQ=255)         tGAACGTGCGCTTATATGGCGAAACCGAGTTCTGGTTTGCGTTGATTAAAGTGCTGGCAATCATCGGTATGATCGGCTTTGGCCTGTGGCTGCTGTTTTCTGGTCACGGCGGCGAGAAAGCCAGTATCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGAT‑TTTGTCGCTGGC‑GGTAATTATGTTC‑TCCTTCGGCGGT‑CT‑GGA‑GCTGATTGGGATTAc                                                                                                                                                                                                                      <  1:1091959/238‑1 (MQ=255)         ggAAGGTGCGCTTATATGGCGAAACGGAGTTCTGGTTTGCGTTGATTAAAGTGCTGGCAATCATCGGTATGATCGGCTTTGGCCTGTGGCTGCTGTTTTCTGGTCACGGCGGCGAGAAAGCCAGTATCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGAT‑TTTGTCGCTGGC‑GGTAATTATGTTC‑TCCTTCGGCGGT‑CT‑GGA‑GCTGATTGGGATTAc                                                                                                                                                                                                                      <  1:181967/237‑1 (MQ=255)             cGTGCGCTTATATGGCGAAACCGAGTTCTGGTTTGCGTTGATTAAAGTGCTGGCAATCATCGGTATGATCGGCTTTGGCCTGTGGCTGCTGTTTTCTGGTCACGGCGGCGAGAAAGCCCGTATCGACAACCTCTGGCGCTACGGTGGTTTTTTCGCCACCGGCTGGAATGGGCTGAT‑TTTGTCGCTGGG‑GGGGATTATTTTCTTCCTTCGGGGGT‑CG‑GGA‑GGGGGTTggggtagtggcc                                                                                                                                                                                                                  >  2:679673/1‑230 (MQ=255)                       ttgggCGAAACCGAGTTCTGGTTTGCGTTGATTAAAGTGCTGGCAATCATCGGTATGATCGGCTTTGGCCTGTGGCTGCTGTTTTCTGGTCACGGCGGCGAGAAAGCCAGTATCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGAT‑TTTGTCGCTGGC‑GGTAATTATGTTC‑TCCTTCGGCGGT‑CT‑GGA‑GCTGATTGGGATTACTGCCGCTGAAg‑cgcg                                                                                                                                                                                                      <  1:233297/236‑1 (MQ=255)                         tGGCGAAACCGCGTTCTGGTTTGCGTTGATTAAAGTGCTGGCAATCATCGGTATGATCGGCTTTGGCCTGTGGCTGCTGTTTTCTGGTCACGGCGGCGGGAAACCCTGTATCGACAACCTCTGGCGCTACGGGGGTTTCTTCGCCACCGCCTGGACTGGGCGGATGTTTTTCGC‑GGC‑GGCACTTATAATT‑TCCTCCTGCGGT‑GT‑GGA‑GCCGCGTGGGTGTTCTGCCGCTGAAg‑ccggcg                                                                                                                                                                                                    >  2:757323/1‑234 (MQ=255)                            cGAAACCGAGTTCTGGTTTGCGTTGATTAAAGTGCTGGCAATCATCGGTATGATCGGCTTTGGCCTGTGGCTGCTGTTTTCTGGTCACGGCGGCGAGAAAGCCAGTATCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGAT‑TTTGTCGCTGGC‑GGTAATTATGTTC‑TCCTTCTGCGGT‑CT‑GGA‑GGTGATTTGGGTTACTGCCGCTTAAG‑agcccggt                                                                                                                                                                                                  >  1:758042/1‑236 (MQ=255)                            cGAAACCGAGTTCTGGTTTGCGTTGATTAAAGTGCTGGCAATCATCGGTATGATCGGCTTTGGCCTGTGGCTGCTGTTTTCTGGTCACGGCGGCGAGAAAGCCAGTATCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCCGAT‑TTTGTCGCTGGC‑GGCAATTTTGTTC‑TTCTTCTGCGGT‑CC‑GGAGGTTGTTTGGGTTTACtccggcgaatgagcgcgt                                                                                                                                                                                                    >  1:846656/1‑221 (MQ=255)                                       ttgtttttggtgtGATTAAAGTGCTGGCAATCATGGGGATGGTCGGTTTTGGCCTGTGGCTGCTGTTTGCGGGGCTCGGCGGCGAGAAAGCCAGTATCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGAT‑TTTGTCGCTGGC‑GGTAATTATGTTC‑TCCTTCGGCGGT‑CT‑GGA‑GCTGATTGGGATTACTGCCGCTGAAG‑CGCGCGATCCGGAAAAAAg                                                                                                                                                                                       <  2:947481/226‑1 (MQ=255)                                                   tGATTAAAGGGCTGGCAATCAGCGGTAAGATCGGCTTGGGCCGTTTGCTGCTGTTTTCTGGTCACGGCGGCGAGAAAGCCAGTATCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGAT‑TTTGTCGCTGGC‑GGTAATTATGTTC‑TCCTTCGGCGGT‑CT‑GGA‑GCTGATTGGGATTACTGCCGCTGAAG‑CGCGCGATCCGGAAAAAAGCA‑TTCCAAAAGCg                                                                                                                                                                         <  1:350757/239‑1 (MQ=255)                                                     aTTAAAGTGCTGGCAATCATCGGTATGATCGGCTTTGGCCTGTGGCTGCTGTTTTCTGGTCACGGCGGCGAGAAAGCCAGTATCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGAT‑TTTGTCGCTGGC‑GGTAATTATGTTC‑TCCTTCGGCGGT‑CT‑GGA‑GCTGATTGGGATTACTGCCGCTGAAG‑CGCGCGATCCGGAAAAAAGCA‑TTCCAAAAGCGGt                                                                                                                                                                       <  1:972479/239‑1 (MQ=255)                                                      ttAAAGTGCGGGCAATCATCGGTATGATCGGTTTTGGCCTGTGGCTGCTGTTTTCTGGTCACGGCGGCGAGAAAGCCAGTATCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGAT‑TTTGTCGCTGGC‑GGTAATTATGTTC‑TCCTTCGGCGGT‑CT‑GGA‑GCTGATTGGGATTACTGCCGCTGAAG‑CGCGCGATCCGGAAAAAAGCA‑TTCCAAAAGCGGt                                                                                                                                                                       <  1:770450/238‑1 (MQ=255)                                                             gCTGGCAATCATCGGTATGATCGGCTTTGGCCTGTGGCTGCTGTTTTCTGGTCACGGCGGCGAGAAAGCCAGTATCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGAT‑TTTGTCGCTGGC‑GGTAATTATGTTC‑TCCTTCGGCGGT‑CT‑GGA‑GCTGATTGGGATTACTGCCGCTGAAG‑CGCGCGATCCGGAAAAAAGCA‑TTCCAAAAGCGGTAAATCAg                                                                                                                                                                <  1:980737/238‑1 (MQ=255)                                                             gCTGGCAATCATCGGTATGATCGGCTTTGGCCTGTGGCTGCTGTTTTCTGGTCACGGCGGCGAGAAAGCCAGTATCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGAT‑TTTGTCGCTGGC‑GGTAATTATGTTC‑TCCTTCGGCGGT‑CT‑GGA‑GCTGATTGGGATTACTGCCGCTGAAG‑CGCGCGATCCGGAAAAAAGCA‑TTCCAAAAGCGGTAAATCAg                                                                                                                                                                <  1:179225/238‑1 (MQ=255)                                                              cTGGCAATCATCGGTATGATCGGCTTTGGCCTGTGGCTGCTGTTTTCTGGTCACGGCGGCGAGAAAGCCAGTATCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGAT‑TTTGTCGCTGGC‑GGTAATTATGTTC‑TCCTTCGGCGGT‑CT‑GGA‑GCTGATTGGGATTACTGCCGCTGAAG‑CGCGCGATCCGGAAAAAAGCA‑TTCCAAAAGCGGTAAATCAgg                                                                                                                                                               <  1:1063710/238‑1 (MQ=255)                                                               tGGCAATCATCGGTATGATCGGCTTTGGCCTGTGGCTGCTGTTTTCTGGTCACGGCGGCGAGAAAGCCAGTATCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGAT‑TTTTTCGCTGGC‑GGGAATTTTGTTC‑TCCTTCGGGGGG‑CT‑GGG‑GCGTAGTGGGGTTACTGCCCCCGAAG‑GGCGCGGTCCGGGAAAAAACA‑ATCCaaaaaggggaaaaaagggg                                                                                                                                                             >  2:179225/1‑224 (MQ=255)                                                               tGGCAATCATCGGTATGACCGGCTTTGGCCTGTGGCTGCTGTTTTCGGGCCACGGCGGCGAGAAAGCCAGTATCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGAT‑TTTGTCGCTGGC‑GGTAATTATGTTC‑TCCTTCGGCGGT‑CT‑GGA‑GCTGATTGGGATTACTGCCGCTGAAG‑CGCGCGATCCGGAAAAAAGCA‑TTCCAAAAGCGGTAAATCAggtg                                                                                                                                                             <  1:148866/239‑1 (MQ=255)                                                                ggCAATCATCGGGATGATCGGCTTTGGCCTGTGGCTGCTGTTTTCTGGTCACGGCGGCGAGAAAGCCAGTATCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGAT‑TTTGTCGCTGGC‑GGTAATTATGTTC‑TCCTTCGGCGGT‑CT‑GGA‑GCTGATTGGGATTACTGCCGCTGAAG‑CGCGCGATCCGGAAAAAAGCA‑TTCCAAAAGCGGTAAATCAggtg                                                                                                                                                             <  1:842533/238‑1 (MQ=255)                                                                 gCAATCATCGGTATGATCGGCTTTGGCCTGTGGCTGCTGTTTTCTGGTCACGGCGGCGAGAAAGCCAGTATCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGAT‑TTTGTCGCTGGC‑GGTAATTATGTTC‑TCCTTCGGCGGT‑CT‑GGG‑GCTGATTGGGATTACTGCCGCTGAAG‑CGCGCGATCCGGAAAAAACAA‑TTCCACAAACGGGAAAttcgggggg                                                                                                                                                           >  1:7686/1‑231 (MQ=255)                                                                 gCAATCATCGGTATGATCGGCTTTGGCCTGTGGCTGCTGTTTTCTGGTCACGGCGGCGAGAAAGCCAGGATCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGGT‑TTTGTCGCTGGCGGGTAATTTTTTTC‑‑CCTTCGGCGGTGCT‑GAG‑GGTTATTGGGTTTAGTGCC‑CTGAAGACGCCGGACCCGGAAAAAAGAA‑TCCAAAAAGGGGaaaacaaggggg                                                                                                                                                            >  2:1063710/1‑231 (MQ=255)                                                                  cAATCCTCGGTATGATCGGCTTTGGCCTGTGGCTGCTGTTTTCTGGTCACGGCGGCGAGAAAGCCAGTATCGACAACCTCGGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGAT‑TTTGTCGCTGGC‑GGTAATTATGTTC‑TCCTTCGGCGGT‑CT‑GGA‑GCTGATTGGGATTACTGCCGCTGAAG‑CGCGCGATCCGGAAAAAAGCA‑TTCCAAAAGCGGTAAATCAggtggt                                                                                                                                                           <  1:679673/238‑1 (MQ=255)                                                                    tggtttgttgttttCGGTTTTGTCCTGGGGCGCCGGTTTTCGGTGCCGGGCGGCAAGAAAGCCGTTATCGACAACCCCTGGCGGTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGAT‑TTTGTCGCTGGC‑GGTAATTATGTTC‑TCCTTCGGCGGT‑CT‑GGA‑GCTGATTGGGATTACTGCCGCTGAAG‑CGCGCGATCCGGAAAAAAGCA‑TTCCAAAAGCGGTAAATCAGGTGGTGTa                                                                                                                                                        <  2:120883/226‑1 (MQ=255)                                                                              tGATCGGCTTTGGCCTGTGGCTGCTGTTTTCTGGGCACGGCGGCGAGAAAGCCAGTATCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGAT‑TTTGTCGCTGGC‑GGTAATTATGTTC‑TCCTTCGGCGGT‑CT‑GGA‑GCTGATTGGGATTACTGCCGCTGAAG‑CGCGCGATCCGGAAAAAAGCA‑TTCCAAAAGCGGTAAATCAGGTGGTGTATC‑GCATCctg                                                                                                                                             <  1:1004790/239‑1 (MQ=255)                                                                              tGATCGGCTTTGGCCGGTGGCTGCTGTTTTCTGGTCGCGGCGGCGAAAAAGCCAGTATCGACAACCTCTGGCTCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGAT‑TTTGTCGCCGGC‑GGTAATTATGTTC‑TCCTTCGGCGGT‑CT‑GGA‑GCTGATTGGGATTACTGCCGCTGAAG‑CGCGCGATCCGGAAAAAAGCA‑TTCCAAAAGCGGTAAATCAGGTGGTGTATC‑GCATCctg                                                                                                                                             <  1:757323/239‑1 (MQ=255)                                                                                      tttGGCCTGTGGCTGCTGTTTTCTGGTCACGGCGGCGAGAAAGCCAGTATCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGAT‑TTTGTCGCTGGC‑GGGAATTATGTTC‑TCCTTCGGCGGG‑CT‑TGG‑GCTGATTTGGATTACTGCCGCTGAAG‑CGCGCGGTCCGGGAAAAAACA‑TTCCCAAAACGGTAAATCAGGGGGGGTATC‑GCATCCTGCT‑Ggttt                                                                                                                                     >  1:364151/1‑238 (MQ=255)                                                                                         tgtgtgggggggttggtttCGGGGCCCGGAGGGGAGAAAGCCGGTTACGACACCCGCTGGGTCTCGGGTGGTTTCTTCGGCACGGGATGGAATGGGCTGAT‑TTTGTCGCTGGC‑GGTAATTATGTTT‑TCCTTGGGCGGT‑CT‑GGA‑GCTGATTGGGATTACTGCCGCTGAAG‑CGCGCGATCCGGAAAAAAGCA‑TTCCAAAAGCGGTAAATCAGGTGGTGTATC‑GCATCCTGCG‑GGTTTACAt                                                                                                                                 <  2:1032952/223‑1 (MQ=255)                                                                                         gggcTGTGGCTGCTGTTTTCTGGTCACGGCGGCGGAAAAGCCAGTATCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGAT‑TTTGTCGCTGGC‑GGTAATTATGTTC‑TCCTTCGGCGGT‑CT‑GGA‑GCTGATTGGGATTACTGCCGCTGAAG‑CGCGCGATCCGGAAAAAAGCA‑TTCCAAAAGCGGTAAATCAGGTGGTGTATC‑GCATCCTGCT‑GTTTTACAt                                                                                                                                 <  2:7686/236‑1 (MQ=255)                                                                                              tgtgctgttTTTTTTTGGTCGCGGGGGCGAGAAAGCCAGTATCGACAACCTTTGGCGGTACGGTGGTTTCTTCGGCACCGGCTGGAATGGGCTGAT‑TTTGTCGGTGGC‑GGTAATTATGTTC‑TCCTTCGGCGGT‑CT‑GGA‑GCTGATTGGGATTACTGCCGCTGAAG‑CGCGCGATCCGGAAAAAAGCA‑TTCCAAAAGCGGTAAATCAGGTGGTGTATC‑GCATCCTGCT‑GTTTTACATCGGtt                                                                                                                            <  2:13746/236‑1 (MQ=255)                                                                                              gtggtggctgTTTTCTGGTCACGGCGGCGAGAAAGCCAGTATCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGAT‑TTTGTCGCTGGC‑GGTAATTATGTTC‑TCCTTCGGCGGT‑CT‑GGA‑GCTGATTGGGATTACTGCCGCTGAAG‑CGCGCGATCCGGAAAAAAGCA‑TTCCAAAAGCGGTAAATCAGGTGGTGTATC‑GCATCCTGCT‑GTTTTACATCGGtt                                                                                                                            <  1:495708/233‑1 (MQ=255)                                                                                                   tgctgTTTTCTGGTCACGGCGGCGAGAAAGCCAGTATCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGAT‑TTTGTCGCTGGC‑GGTAATTATGTTC‑TCCTTCGGCGGT‑CT‑GGA‑GCTGATTGGGATTACTGCCGCTGAAG‑CGCGCGATCCGGAAAAAAGCA‑TTCCAAAAGCGGTAAATCAGGTGGTGTATC‑GCATCCTGCT‑GTTTTACATCGGTTCACt                                                                                                                        <  1:711766/238‑1 (MQ=255)                                                                                                        ttttCTGGTCACGGCGGCGAGAAAGCCAGTATCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGAT‑TTTGTCGCTGGC‑GGTAATTATGTTC‑TCCTTCGGCGGT‑CT‑GGA‑GCTGATTGGGATTACTGCCGCTGAAG‑CGCGCGATCCGGAAAAAAGCA‑TTCCAAAAGCGGTAAATCAGGTGGTGTATC‑GCATCctg                                                                                                                                             <  1:7594/213‑1 (MQ=255)                                                                                                              tgttccgggggGGAGAAAGCCCGTTTCGAAAACCTTTGGCCGTTGGGGGGTTTCTTGGCCACGGGCTGGAATGGGCGTAT‑TGTGTCGCTGGC‑GGTAATTATGTTC‑TCCTTCGGCGGT‑CT‑GGA‑GCTGATTGGGATTACTGCCGCTGAAG‑CGCGCGATCCGGAAAAAAGCA‑TTCCAAAAGCGGTAAATCAGGTGGTGTATC‑GCATCCTGCT‑GTTTTACATCGGTTCACT‑GGTGGTTTTACt                                                                                                           <  2:961920/234‑1 (MQ=255)                                                                                                                 cACGGGGGGCAGAAAGCCAGTATCGACAACCTTCGGCGCTACGGTTGTTTCTTCGCCAGCGGCTGGAATGGGCTTAT‑TTTGTCGGTGGC‑GGTAATTATGTTC‑TCATTCGGCGGT‑CT‑GGA‑GCTGATTGGGATTACTGCCGCTGAAG‑CGCGCGATCCGGAAAAAAGCA‑TTCCAAAAGCGGTAAATCAGGTGGTGTATC‑GCATCCTGCT‑GTTTTACATCGGTTCACT‑GGTGGTTTTACTgg                                                                                                         <  1:861376/238‑1 (MQ=255)                                                                                                                 cACGGCGGGGAGAAAGCCAGTATCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCCACGGCTGGAATGGGCTGAT‑TTTGTCGCTGGC‑GGTAATTATGTTC‑TCCTTCGGCGGT‑CT‑GGA‑GCTGATTGGGATTACTGCCGCTGAAG‑CGCGCGATCCGGAAAAAAGCA‑TTCCAAAAGCGGTAAATCAGGTGGTGTATC‑GCATCCTGCT‑GTTTTACATCGGTTCACT‑GGTGGTTTTACTgg                                                                                                         <  1:812992/238‑1 (MQ=255)                                                                                                                  ccgggggggggaAAGCACGTTTCGCCAACCTCTGGCGCTACGGTTGTTTCTTCGCCACCGGCTGGGATGGGCTGAT‑TTTGTCGCTGGC‑GGTAATTATTTTC‑CCCTTCGGCGGT‑TG‑GGA‑GCTGATTGGGATTACTGCCGCTGAAG‑CGCGCGATCCGGAAAAAAGCA‑TTCCAAAAGCGGTAAATCAGGTGGTGTATC‑GCATCCTGCT‑GTTTTACATCGGTTCACT‑GGTGGTTTTGCTGgcg                                                                                                       <  2:1076208/229‑1 (MQ=255)                                                                                                                  ccggcggCGAGAAAGCCAGTATCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGAT‑TTTGTCGCTGGC‑GGTAATTATGTTC‑TCCTTCGGCGGT‑CT‑GGA‑GCTGATTGGGATTACTGCCGCTGAAG‑CGCGCGATCCGGAAAAAAGCA‑TTCCAAAAGCGGTAAATCAGGTGGTGTATC‑GCATCCTGCT‑GTTTTACATCGGTTCACT‑GGTGGTTTTACTGgcg                                                                                                       <  1:563917/238‑1 (MQ=255)                                                                                                                  ccggcgcgaagaAAGCCAGTTTCGACAACCTCGGGCTCTGCGGTGGTTTCCTCTGCAGCGGCTGGAATGGGCTGAT‑TTTGTCGCTGGC‑GGTAATTATGTTC‑TCCTTCGGCGGT‑CT‑GGA‑GCTGATTGGGATTACTGCCGCTGAGG‑CGCGCGATCCGGAAAAAAGCA‑TTCCAAAAGCGGTAAATCAGGTGGTGTATC‑GCATCCTGCT‑GTTTTACATCGGTTCACT‑GGTGGTTTTACTGgc                                                                                                        <  1:470497/229‑1 (MQ=255)                                                                                                                  aCGGGGGCGAGAAAGCCAGTACCGACACCCGCTGGCGTTACGGGGGTTTCTTCGCCACCGGCTGGAATGGGCTGAT‑TTTGTCGCTGGC‑GGTAATTATGTTC‑TCCTTCGGCGGT‑CT‑GGA‑GCTGATTGGGATTACTGCCGCTGAAG‑CGCGCGATCCGGAAAAAAGCA‑TTCCAAAAGCGGTAAATCAGGTGGTGTATC‑GCATCCTGCT‑GTTTTACATCGGTTCACT‑GGTGGTTTTACTGgcg                                                                                                       <  2:636677/239‑1 (MQ=255)                                                                                                                  aCGGCGGCGAGAAAGCCAGTATCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGAT‑TTTGTCGCTGGC‑GGTAATTATGTTC‑TCCTTCGGCGGT‑CT‑GGA‑GCTGATTGGGATTACTGCCGCTGAAG‑CGCGCGATCCGGAAAAAAGCA‑TTCCAAAAGCGGTAAATCAGGTGGTGTATC‑GCATCCTGCT‑GTTTTACATCGGTTCACT‑GGTGGTTTTACTGgcg                                                                                                       <  1:1006827/239‑1 (MQ=255)                                                                                                                    gggggagggaAAACCAGTTTCGACAACCTTGTGCGTTAGGGTGGTTTCTTGGCCACCGGCTGGAATGGGCTGAG‑TTGGTCGCTGGC‑GGTAATTATGTTC‑TCCTTCGGCGGT‑CT‑GGA‑GCTGATTGGGATTACTGCCGCTGAAG‑CGCGCGATCCGGAAAAAAGCA‑TTCCAAAAGCGGTAAATCAGGTGGTGTATC‑GCATCCTGCT‑GTTTTACATCGGTTCACT‑GGTGGTTTTACTGGCGct                                                                                                     <  2:364151/231‑1 (MQ=255)                                                                                                                                  cggttacgccccccTTTGGCGCTGGGGTGGTTTCTTCGCCCCCGGCTGGAATGGGCGGAT‑TTTGTCGCTGGC‑GGTAATTATGTTC‑TCCTGCGGCGGT‑CT‑GGA‑GCTGATTGGGATTACTGCCGCTGAAG‑CGCGCGATCCGGAAAAAAGCA‑TTCCAAAAGCGGTAAATCAGGTGGTGTATC‑GCATCCTGCT‑GTTTTACATCGGTTCACT‑GGTGGTTTTACTGGCGCTCTATCCgtgggtgg                                                                                       <  2:871048/227‑1 (MQ=255)                                                                                                                                   aGTATCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGAT‑TTTGTCGCTGGC‑GGGAATTATGTTC‑TCCTTCGGCGGT‑CT‑GGA‑GCTGACTGGGATTACTGCCGCTGAAA‑CGGGCGATCCGCAGCAAAGCA‑ATCCCAAAAGGGTAAATCAGGAGGGGGATCGGCAT‑CAGCG‑GTTTTACATCGGTTTAAT‑GGTGGGTTTTCTgcggccttctccaggggggga                                                                                      >  2:821334/1‑219 (MQ=255)                                                                                                                                          aCAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGAT‑TTTGTCGCTGGC‑GGTAATTATGTTC‑TCCTTCGGCGGT‑CT‑GGA‑GCTGATTGGGATTACTGCCGCTGAAG‑CGCGCGATCCGGAAAAAAGCA‑TTCCAAAAGCGGTAAATCAGGTGGTGTATC‑GCATCCTGCT‑GTTTTACATCGGTTCACT‑GGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGaaa                                                                               <  1:604110/239‑1 (MQ=255)                                                                                                                                                  tGGCGTTTCGGTGGTTTTTTCGTCTCCGGCTGGAATGGGCTGAT‑TTTGTCGGTGGC‑GGTAATTATGTTC‑TCCTTCGGCGGT‑CT‑GGA‑GCTGATTGGGATTACTGCCGCTGAAG‑CGCGCGATCCGGAAAAAAGCA‑TTCCAAAAGCGGTAAATCAGGTGGTGTATC‑GCATCCTGCT‑GTTTTACATCGGTTCACT‑GGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAg                                                                       <  2:369769/239‑1 (MQ=255)                                                                                                                                                         aCGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGAT‑TTTGTCGCTGGC‑GGTGATTATGTTC‑TCCTTCGGCGGT‑CT‑GGA‑GGTGATTGGGGTTCCTGCCGCTGAAA‑CGGCCGATCCGGAAAAAAGCA‑TTCCCAAAACGGGCAATCCGGAGGGGTATC‑GCCTCCTGGC‑GTTTTTCCACCGTGTACTGGGGGGTTTTATGGACGCTCTACACGTGGG‑GGATGGGAAATCCaaacgaacgccc                                                                >  2:226828/1‑229 (MQ=255)                                                                                                                                                         aCGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGAT‑TTTGTCGCTGGC‑GGTAATTATGTTC‑TCCTTCGGCGGT‑CT‑GGA‑GCTGATTGGGATTACTGCCGCTGAAG‑CGCGCGATCCGGAAAAAAGCA‑TTCCAAAAGCGGTAAATCAg                                                                                                                                                                >  1:633719/1‑146 (MQ=255)                                                                                                                                                         aCGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGAT‑TTTGTCGCTGGC‑GGTAATTATGTTC‑TCCTTCGGCGGT‑CT‑GGA‑GCTGATTGGGATTACTGCCGCTGAAG‑CGCGCGATCCGGAAAAAAGCA‑TTCCAAAAGCGGTAAATCAg                                                                                                                                                                <  2:633719/146‑1 (MQ=255)                                                                                                                                                         aCGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGAT‑TTTGTCGCTGGC‑GGTAATTATGTTC‑TCCTTCGGCGGG‑CT‑GGA‑GCTGGTTTGGGTTACTTCCCCTGAAG‑CGCGCGACCCGGCAAAAAGCA‑TTCCAAAAGCGGTAAATCAGGGGGGGGATC‑GGCTCCTGGTGGTTTTACACCGG‑TCCCT‑GGGGGGTTTTCTGGGGCGCTATTCGTGGGTGGAAGGGGAAATCAACAGT‑AGCCCg                                                               >  2:370391/1‑239 (MQ=255)                                                                                                                                                                    ttcGCCACCGGCTGGAATGGGCTGAT‑TTTGTCGCTGGC‑GGTAATTATGTTC‑TCCTTCGGCGGT‑CT‑GGA‑GCTTGTTGGGGTTACTGCCGCTGACG‑CGCGCGATCCGGAAAAAAGCA‑TTCCAAAAACGGGAAATCAGGGGGTGGATC‑GCATCCTTAT‑GTTTTACATCGGTTCACT‑GGTGGGTTTTTTGGGGCCCTCTCCGTGGGGGGGAGGGAAAACCaaaaaaaaaccgcttgggtgtgt                                                     >  2:677319/1‑218 (MQ=255)                                                                                                                                                                    ttcGCCACCGGCTGGAATGGGCTGAT‑TTTGTCGCTGGC‑GGTAATTATGTTC‑TCCTTCGGCGGT‑CT‑GGA‑GCTGATTGGGATTACTGCCGCTGAAG‑CGCCCGGTCCGGAAAAAAGCA‑TTCCCAAAGCGGTAAATCAGGGGGGGGATC‑GCATCCTGCT‑GGTTTACATCGGTTCACT‑GGTGGGTTTTCTTGCGCCCTTTCCCGGGGGGGGAGGGGAAACCAAAAGTAAGCCGGTTtgggtgat                                                     >  2:779650/1‑233 (MQ=255)                                                                                                                                                                    ttcCCCCCCGGCTGGAATGGGCTGGT‑TTTGGCGCTGGC‑GGTAATTATGTTC‑GCCTTCGGCGGT‑CT‑GGA‑GCTGATTGGGATTACTGCCGCTGAAG‑CGCGCGATCCGGAAAAAAGCA‑TTCCAAAAGCGGTAAATCAGGTGGTGTATC‑GCATCCTGCT‑GTTTTACATCGGTTCACT‑GGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGT‑AGCCCGTTTGTgatgat                                                    <  2:758042/239‑1 (MQ=255)                                                                                                                                                                              tgttgtatgggttttt‑tttGTCGGTGGC‑CGTAATTATGTTC‑CTCTTCGGCGGT‑CT‑GGG‑CCTGATTGGGATTCCTGCCGCTGAAG‑CGCGCGATCCGGAAAAAAGCA‑TTCCAAAAGCGGTAAATCAGGTGGTGTATC‑GCATCCTGCT‑GTTTTACATCGGTTCACT‑GGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGT‑AGCCCGTTTGTGATGATTTTCCATAAt                                          <  2:652118/224‑1 (MQ=255)                                                                                                                                                                              gCTGGAATGGGCTGAT‑TTTGTCGCTGGC‑GGTAATTATGTTC‑TCCTTCGGCGGT‑CT‑GGA‑GCTGATTGGGATTACTGCCGCTGAAG‑CGCGCGATCCGGAAAAAAGCA‑TTCCAAAAGCGGTAAATCAGGTGGTGTATC‑GCATCCTGCT‑GTTTTACATCGGTTCACT‑GGTGGTTTTACTGGCGct                                                                                                     <  1:670186/181‑1 (MQ=255)                                                                                                                                                                                 ggAATGGGCTTAT‑TTTGTCGCTGGC‑GGTAATTATTTTC‑TCCTTCGGCGGT‑CT‑GGA‑GCTGATTGGGATTACTGCCGCTGAAG‑CGCGCGATCCGGAAAAAAGCA‑TTCCAAAAGCGGTAAATCAGGTGGTGTATC‑GCATCCTGCT‑GTTTTACATCGGTTCACT‑GGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGT‑AGCCCGTTTGTGATGATTTTCCATAAtct                                        <  1:581440/238‑1 (MQ=255)                                                                                                                                                                                    aTGGGCTGAT‑TTTGTCGCTGGC‑GGTAATTATGTTC‑TCCTTCGGCGGT‑CT‑GGA‑GCTGATTGGGATTACTGCCGCTGAAG‑CGCGCGATCCGGAAAAAAGCA‑TTCCAAAAGCGGTAAATCAGGTGGTGTATC‑GCATCCTGCT‑GTTTTACATCGGTTCACT‑GGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGT‑AGCCCGTTTGTGATGATTTTCCATAATCTCGa                                     <  1:938146/238‑1 (MQ=255)                                                                                                                                                                                      gggCTGAT‑TTTGTCGGTGGC‑GGTAATTATGTTC‑TCCTTCGGCGGT‑CT‑GGA‑GCTGATTGGGATTACTGCCGCTGAAG‑CGCGCGATCCGGAAAAAAGCA‑TTCCAAAAGCGGTAAATCAGGTGGTGTATC‑GCATCCTGCT‑GTTTTACATCGGTTCACT‑GGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGT‑AGCCCGTTTGTGATGATTTTCCATAATCTCGACAg                                  <  1:608710/239‑1 (MQ=255)                                                                                                                                                                                        ggttttgttttcgccTGGG‑GTTATTTATTTTC‑TCCCTCGTGGTT‑CG‑GGG‑TCTGATTGGGATTACGGCTGCTTAGG‑CGCGGGATCCGG‑AAAAAGCATTTCCAAAAGCGGTAAAGCAGGTGGTGTATC‑GCATCCTGCT‑GTTTTACATCGGTTCACG‑GGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGT‑AGCCCGTTTGTGATGATTTTCCATAATCTCGACAGc                                 <  2:456790/225‑1 (MQ=255)                                                                                                                                                                                        gCTGAT‑TTTGTCGCTGGC‑GGTCATTATGTTC‑TCCTTCGGCGGT‑CT‑GGA‑GCTGATTTGGGTTACTGCCGCTGAAG‑CCCGCAGACCGGAAAAAAGAA‑TTCCAAAAGCGGTAAATCCGGGGGTGCAAC‑GCATACTGCT‑GTTTTATATCGCGTCACT‑GTGGGGTTTATTGGCCCCCTTTCCGGGGGTGGAAGGGAAAACCAACCAA‑AGCCCTGTTGTGTTTATTTTTCATacaaacgacacac                                >  2:438924/1‑227 (MQ=255)                                                                                                                                                                                             t‑tttGTCGCTGGC‑GGTAATTATGTTC‑TCCTTCGGCGGT‑CT‑GGA‑GTTGATTGGGATTTATGCCCATGAAG‑CGCGCGATCCGGAAAAAAACA‑TTCCAAAATCGGTAAATCAAGGGGGGTATC‑CCATCCTTCT‑GTTTTACCTCGGTTACCT‑TGGGGTTTTACCGGGGCTGTATCCTGGGGTGGGAGTGAAATCCCACAGG‑AACCCGTTTGGGTAGGTTTTTCATAAACTAGAACGCaaacgg                           >  2:657421/1‑235 (MQ=255)                                                                                                                                                                                                ttGTCGCTGGC‑GGTAATTATGTTC‑TCCTTCGGCGGT‑CT‑GGA‑GCTGATTGGGATTACTGCCGCTGAAG‑CGCGCGATCCGGAAAAAAGCA‑TTCCAAAAGCGGTAAATCAGGTGGTGTATC‑GCATCCTGCT‑GTTTTACATCGGTTCACT‑GGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGT‑AGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGt                         <  1:756887/239‑1 (MQ=255)                                                                                                                                                                                                      ccgtgtggttcttgtGTTC‑TCCTTTGTCGGT‑CTGGGG‑GCTATTTGTGA‑TACGGCGCCTGAAG‑CGCGCCATCCGGAAAAAAGCA‑TTCCAAAAGGGGTAATTCAGGGTGTGTATC‑GCATCCTGGT‑GTTTTCCATCGGTGCACT‑GGTGGTTTTACTGGCGTTCTAGCCGTGGGTGGAAGTGAAATCCAACAGT‑AGCCCGTTTGTTATGATGTTCCATAATCTCGACAGCAACGTGGTAGCtt                    <  2:863836/225‑1 (MQ=255)                                                                                                                                                                                                                     ggtc‑tcCTTCGGCGGT‑CT‑GGG‑GCTGATTGGGATTGCCGCCGCTGAAG‑CCGGCGATCCGGAAAAAAGCA‑TTCCAAAAGCGGTAAATCAGGTGGTGTATC‑GCATCCTGCT‑GTTTTACATCGGTTCACT‑GGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGT‑AGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAGCTTCTGCGCTGAACTTCg     <  1:754238/237‑1 (MQ=255)                                                                                                                                                                                                                        c‑tcCTTCGGCGGT‑CT‑GGA‑GCTGATTGGGATTACTGCCGCTGAAG‑CGCGCGATCCGGAAAAAAGCA‑TTCCAAAAGCGGTAAATCAGGTGGTGTATC‑TGATCCTGCT‑GTTTTACATCGGTTCACT‑GGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGT‑AGCCCGTTTGTGATGATTTTCCATAATCTCGCCAGCAACGTGGTAGGTTTTGCGCTGAACTTCGTCa  >  1:448322/1‑239 (MQ=255)                                                                                                                                                                                                                        c‑tcCTTCGGCGGT‑CT‑GGA‑GCTGATTGGGATTACTGCCGCTGAAG‑CGCGCGATCCGGAAAAAAGCA‑TTCCAAAAGCGGTAAATCAGGTGGTGTATC‑GCATCCTGCT‑GTTTTACATCGGTTCACT‑GGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATc                                                                             <  2:358377/165‑1 (MQ=255)                                                                                                                                                                                                                        c‑tcCTTCGGCGGT‑CT‑GGA‑GCTGATTGGGATTACTGCCGCTGAAG‑CGCGCGATCCGGAAAAAAGCA‑TTCCAAAAGCGGTAAATCAGGTGGTGTATC‑GCATCCTGCT‑GTTTTACATCGGTTCACT‑GGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATc                                                                             >  1:358377/1‑165 (MQ=255)                                                                                                                                                                                                                        c‑tcCTTCGGCGGT‑CT‑GGA‑GCTGATTGGGATTACTGCCGCTGAAG‑CGCGCGATCCGGAAAAAAGCA‑TTCCAAAAGCGGTAAATCAGGTGGGGTATC‑GCATCCTGCT‑GTTTTACATCGGTTCACT‑GGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCACCGGT‑AGCCCGTTTGGGATGATTTTCCATAATCTCGCAACCAACGGGGGAGTGTCTGCGCGGAACTTGGtaa  >  2:469704/1‑237 (MQ=255)                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                     TAACCTGGTGAACGTGCGCTTATATGGCGAAACCGAGTTCTGGTTTGCGTTGATTAAAGTGCTGGCAATCATCGGTATGATCGGCTTTGGCCTGTGGCTGCTGTTTTCTGGTCACGGCGGCGAGAAAGCCAGTATCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGAT‑TTTGTCGCTGGC‑GGTAATTATGTTC‑TCCTTCGGCGGT‑CT‑GGA‑GCTGATTGGGATTACTGCCGCTGAAG‑CGCGCGATCCGGAAAAAAGCA‑TTCCAAAAGCGGTAAATCAGGTGGTGTATC‑GCATCCTGCT‑GTTTTACATCGGTTCACT‑GGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGT‑AGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAGCTTCTGCGCTGAACTTCGTCA  >  NC_000913/602405‑602856

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