Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,429,997	C→T	15.8%	S45F (TCT→TTT)	cycA →	D‑alanine/D‑serine/glycine transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,429,997	0	C	T	15.8%	132.2 / 4.6	65	S45F (TCT→TTT)	cycA	D‑alanine/D‑serine/glycine transporter
Reads supporting (aligned to +/- strand):  ref base C (18/36);  new base T (8/3);  total (26/39)
Fisher's exact test for biased strand distribution p-value = 2.07e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATATCGATGAGTCCCGATACAGATTCGTCGTATCATAGACTGACTAAAGGCCGTAGAGCCTGAACAACACAGACAGGTACAGGAAGAAAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTG‑CCATT‑GG‑TAC‑GGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGA‑TC‑ATTTTCGTTTATATGAT‑CA‑TTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATTGCTGCTTTCGAATCTGG‑AATACAAATCTTTTAGTGACTTCGC‑TTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCT  >  NC_000913/4429770‑4430235                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                     attccGATTAGTCCCGATACAAATTCGTCGTATCATAGACTGAACAAAGGCCGTAGATCCTGAACAACACATACAGGTAAAGGAAGAAAAAACCATGGTAGATCCGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTG‑CCATT‑GG‑TAC‑GGGGTTGTTTATGGGGTCTGGCAAAACGa                                                                                                                                                                                                                                           <  2:270935/235‑1 (MQ=255)       aTGAGTCCCGATACAGATTCGTCGTATCATAGACTGACTAAAGGCCGTAGAGCCTGAACAACACAGACAGGTACAGGAAGAAAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGGG‑CCATTGGG‑AAC‑GGGGTGGTTTATGGGG‑CTGGCaaaaagtatta                                                                                                                                                                                                                                       >  1:359873/1‑230 (MQ=255)            tCCCGATACAGATTCGTCGTATCATAGACTGACTAAAGGCCGTAGAGCCTGAACAACACAGACAGGTACAGGAAGAAAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGGGCCCATT‑GG‑‑AC‑GGGGGTGTTTTTGTGGGTTTGCAAAAAGAATAGcccttcc                                                                                                                                                                                                                                >  2:1029941/1‑234 (MQ=255)             cccGATACAGATTCGTCGTATCATAGACTGACTAAAGGCCGTAGAGCCTGAACAACACAGACAGGTACAGGAAGAAAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTG‑CCATT‑GG‑TAC‑GGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCg                                                                                                                                                                                                                               <  2:259410/239‑1 (MQ=255)                       aTTCGTCGTATCATAGACTGACTAAAGGCCGTAGAGCCTGAACAACACAGACAGGTACAGGAAGAAAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTG‑CCCTT‑GG‑TAC‑GGGGGTGTTTTTTGGGGCTTGCAAAACAATTAGCCTTTCCGGGCCGGCGac                                                                                                                                                                                                                     >  1:1100964/1‑238 (MQ=255)                       aTTCGTCGTATCATAGACTGACTAAAGGCCGTAGAGCCTGAACAACACAGACAGGTACAGGAAGAAAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTAGCTACGGCGCAATCTCACAAACCGACATATTCAGCTAATTGCCATTGGCGGAG‑CGCTT‑TG‑TAT‑GGGGGTGGTTTTTGGGTCGGGGAAAACAAGTAGACCTGCCCGGGCgggggc                                                                                                                                                                                                                     >  2:883964/1‑234 (MQ=255)                            tcgtATCATAGACTGACTAAAGGCCGTAGAGCCTGAACAACACAGACAGGTACAGGAAGAAAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTG‑CCCTT‑GG‑TAC‑GGGGTTGTTTATTGGGTCTTGCCAAACGATTTGCCCTTCCGGGGCGTCGA‑TG‑Attt                                                                                                                                                                                                              >  1:1099334/1‑239 (MQ=255)                            aggttacATAGACTGAAAAAAGGCGGGAGAGCCTGAACAACAAAGACAGGTAAAGGAAAAAAAAAACTGGGTAGTACCGGTAAAAGGCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGTTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTG‑CCATT‑GG‑TAC‑GGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGA‑TC‑Att                                                                                                                                                                                                               <  2:9754/232‑1 (MQ=255)                             aacaTCATAGAATGATAAAAGGCCGTAGAGCCTGAAGAACACAGACAGGTACAGGAAGAAAAAAACATGGTAGAGCAGGTAAAAGTCGTTGCCGATGACCAGGCTCCGGCTTAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTG‑CCATT‑GG‑TAC‑GGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGA‑TC‑Attt                                                                                                                                                                                                              <  2:23695/235‑1 (MQ=255)                                 tCATAGACTGACTAAAGGCCGTAGAGCCTGAACAACACAGACAGGTACAGGAAGAAAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTG‑CCCTT‑GG‑TAA‑GGGGGTGGTTTTGGGGGTTGGGAAAACAATTAACCTTTCCGGGCCGTCGagaatattttgt                                                                                                                                                                                                           >  1:230845/1‑228 (MQ=255)                                     aGACTGACTAAAGGCCGTAGAGCCTGAACAACACAGACAGGTACAGGAAGAAAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTG‑CCATT‑GG‑TAC‑GGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGA‑TC‑ATTTTCGTTtata                                                                                                                                                                                                     <  1:1085544/239‑1 (MQ=255)                                       actgactAAAGGCCGTAGAGCCTGAACAACACAGACAGGTACAGGAAGAAAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTG‑CCATT‑GG‑TAC‑GGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGA‑TC‑ATTTTCGTTTATATg                                                                                                                                                                                                   <  1:249283/239‑1 (MQ=255)                                         tgactAAAGGCCGTAGAGCCTGAACAACACAGACAGGTACAGGAAGAAAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTG‑CCATT‑TG‑TAC‑TGGGGTGTTTATTGGGGCTGGGAAAACGATTTACCCTTCCCGGCCGCCGC‑TC‑ATTTTTGTTTttatctt                                                                                                                                                                                                 >  1:9754/1‑236 (MQ=255)                                           actAAAGGCCGTAGAGCCTGAACAACACAGACAGGTACAGGAAGAAAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTG‑GCATT‑GG‑TAA‑GGGGTTGGTTATGGGGGCTGGCAAAACGATTATCCTTGCCGGGCCGGCGA‑TT‑CTTTTCGGTTatatatgca                                                                                                                                                                                               >  1:145966/1‑234 (MQ=255)                                           actAAAGGCCGTAGAGCCTGAACAACACAGACAGGTACAGGAAGAAAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTG‑CCCTT‑TG‑TAC‑GGGGTTGTTTTTGGGGTTCTGCAAAACAATTAGCCTTACCGGGCCGTCGG‑TC‑GCTTTCTTTTttatattca                                                                                                                                                                                               >  1:292096/1‑234 (MQ=255)                                           actAAAGGCCGTAGAGCCTGAACAACACAGACAGGTACAGGAAGAAAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGGG‑CCCTT‑GG‑GAA‑GGGGGTGTTTTTGGGGTTTGGCAAAAACATTAGGCTTTCCGGGCCGGGCG‑TC‑ATTTTCGGTTATTTgcttt                                                                                                                                                                                               >  1:23695/1‑235 (MQ=255)                                           actAAAGGCCGTAGAGCCTGAACAACACAGACAGGTACAGGAAGAAAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAAACGACATATTCAGCTTATTGCCATTGGCGGGT‑GCATT‑TG‑TACGGGGGTGGTTTTTGGGGTTGTGAAAAACAATGTGCTTTGCCGGGCC‑TGGG‑TC‑ATTTTTGttttctttaaca                                                                                                                                                                                               >  2:1044462/1‑230 (MQ=255)                                            ctAAAGGCCGTAGAGCCTGAACAACACAGACAGGTACAGGAAGAAAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTG‑CCATT‑GG‑TAC‑GGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGA‑TC‑ATTTTCGTTTATATGAT‑CA‑t                                                                                                                                                                                            <  1:883964/239‑1 (MQ=255)                                            aaaaaGGCCGTAGGTCCTGAACAACACAGACAGGTACAGGAAGAAAAAAACAGGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATACTCAGCTTATTGCCATTGGCGGTG‑CCATT‑GG‑TAC‑GGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGA‑TC‑ATTTTCGTTTATATGAT‑CA‑t                                                                                                                                                                                            <  1:737964/237‑1 (MQ=255)                                                     gTAGAGCCTGAACAACACAGACAGGTACAGGAAGAAAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTG‑CCATT‑TG‑TAC‑GGGGGTGCTTATGGGGTCTGGCAAAACGATTTACGTTGCCGGGCCGCCGA‑TA‑CTTTTCGTTTTAATAGT‑CA‑TTTGttttta                                                                                                                                                                                   >  1:751691/1‑237 (MQ=255)                                                     gTAGAGCCTGAACAACACAGACAGGTACAGGAAGAAAAAAACATGGTAGATCAGGCAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGTTACGGCGCAATCTCACCAACCGAGATATTCCGATTATGGCCATTGGCGGTG‑CCACT‑GGTTTC‑GGGGGTGTTT‑TGAGGGCTGGCACAAAGATTTGCCATTCCGGGGCGTGGAGCC‑GTTTTTGTTT‑TTTTAT‑CA‑TTTGttttta                                                                                                                                                                                   >  2:910751/1‑237 (MQ=255)                                                                    cacaGACAGGTACAGGAAGAAAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTG‑GCCTT‑TG‑ATC‑GGGGTTGGTTATTGGGGTTGGCCAAACGGTTTACCCTGGCCGGTCGGGGG‑TC‑CTTTTTGTTTTTATAAT‑TATTTGGTTTTTTGTTTTTTTTAGtat                                                                                                                                                                     >  2:119709/1‑237 (MQ=255)                                                                         acagGTACAGGAAGAAAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTG‑CCATT‑GG‑TAA‑CGGGGTGTTTTTGGGGGCTGGCCAAAAGATTAGCCTTTCCGGGGCGGCGG‑TC‑CTTTTCGTTTTTATGAT‑TA‑TTGGGTTTATGCTCTTTTTTGTGATGCgg                                                                                                                                                                >  1:55174/1‑238 (MQ=255)                                                                           agGTACAGGAAGAAAAAAACATGGGAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTG‑CCATT‑GG‑TAC‑GGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGA‑TC‑ATTTTCGTTTATATGAT‑CA‑TTGGTTTTATGCTCTTTTTCGTGATGCGGGCa                                                                                                                                                             <  1:634803/239‑1 (MQ=255)                                                                              tACAGGAAGAAAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTG‑CCATT‑GG‑TAC‑GGGGTTGTTTATGGGTTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGA‑TC‑ATTTTCGTTTATATGAT‑CA‑TTGGTTTTATGCTCTTTTTCGTGATGCGGGc                                                                                                                                                              <  1:398454/235‑1 (MQ=255)                                                                                       aaaaaaaCATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTG‑CCATT‑GG‑TAC‑GGGGGTGTTTTTGGGGTCTGGCCAAACGATTAGCCTTTCCGGGCCGTCGG‑TC‑ATTTTTGTTTTTTTGAT‑CT‑TTGGGTTTTTTCTCTTTTTTGTGCTTCGGGGAATTGGGGGAttt                                                                                                                                                 >  1:1294739/1‑238 (MQ=255)                                                                                        aaaaaaCATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGGG‑CCCTT‑TG‑TAC‑GGGGGTGTTTTTGGGGACTGGCAAAAAGTCTAGCCCTGCCGGGCGCTGGC‑TCTATTTTC‑TTTTTATTAT‑GA‑GTGGGTTTATTCTCTTTTTTGTTGTGCGGGCACAggggggtaggg                                                                                                                                                >  2:385104/1‑233 (MQ=255)                                                                                                 gggaGATCAGGTAAAAGGCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTG‑CCATT‑GG‑TAC‑GGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGA‑TC‑ATTTTCGTTTATATGAT‑CA‑TTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATTGCTGCTTTCGa                                                                                                                                       <  1:644528/236‑1 (MQ=255)                                                                                                       tCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCCCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTG‑CCATT‑GG‑TAC‑GGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGA‑TC‑ATTTTCGTTTATATGAT‑CA‑TTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATTGCTGCTTTCGAATCTg                                                                                                                                  <  1:1102601/238‑1 (MQ=255)                                                                                                        cAGGTAAAGGTCGTTTCCTATTATCGGGGTCGCGTTGAAACGCTGCTACGGCGCATACTCACAACCCTACCTATTCAGCTTATTGCCATTGGCGGTG‑CCATT‑GG‑TAC‑GGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGA‑TC‑ATTTTCGTTTATATGAT‑CA‑TTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATTGCTGCTTTCGAATCTg                                                                                                                                  <  2:751691/237‑1 (MQ=255)                                                                                                         aGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTG‑CCATT‑GG‑TAA‑GGGGGTGTTTTTGGGGGCTGGCAAAAAGAGTAGCCCTGCCCGGCCGGCGG‑CC‑ATTTTTGGTTTTATGAT‑CA‑TTGGTTTTTTGTCCTTTTTTTTGTTGCGGGCAAAGGGGGGATTGTTGATTTTGaaatttgaa                                                                                                                               >  2:644528/1‑232 (MQ=255)                                                                                                              aaaggctgtGCCGATGAACAGGCCCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTG‑CCATT‑GG‑TAC‑GGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGA‑TC‑ATTTTCGTTTATATGAT‑CA‑TTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATTGCTGCTTTCGAATCTGG‑AATACaa                                                                                                                         <  1:1326915/231‑1 (MQ=255)                                                                                                              aaaGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTG‑CCATT‑GG‑TAC‑GGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGA‑TC‑ATTTTCGTTTATATGAT‑CA‑TTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATTGCTGCTTTCGAATCTGG‑AATACaa                                                                                                                         <  1:1029941/239‑1 (MQ=255)                                                                                                                     ttGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTG‑CCCTT‑TG‑TAC‑GGGGCGGTTTTTTGGGTTTTGCCAAACCGCTTGCCTTGCCCGGCCGCTGA‑CG‑GTTTTTGTTTGTATCAT‑CA‑CTGGTTTTATGCCCCTTTTTGTGATCCGGGATATGCGGGAATTGTGGCTTTTTACTCTGG‑CATACAAACCtttt                                                                                                                  >  2:194983/1‑239 (MQ=255)                                                                                                                                  ggCTCCGGCTGAACAGTCGGTACGGGACAATCTCACAAACCGACATATTCAGTTTATTGCCATTGGCGGTG‑CCATT‑GG‑TAC‑GGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGA‑TC‑ATTTTCGTTTATATGAT‑CA‑TTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATTGCTGCTTTCGAATCTGG‑AATACAAATCTTTTAGTGACTTCGC‑tt                                                                                                    <  1:176564/239‑1 (MQ=255)                                                                                                                                      ccGGTTGAAAAGTCGTTACGGGCCAATCTCACAAACCGGTATATTCAGGTTATTGCCATTGGGGGGG‑CCATT‑GG‑TAC‑GGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGA‑TC‑ATTTTCGTTTATATGAT‑CA‑TTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATTGCTGCTTTCGAATCTGG‑AATACAAATCTTTTAGTGACTTCGC‑TTCCGa                                                                                                <  2:55174/239‑1 (MQ=255)                                                                                                                                       ctgtgtaaCGGTCCCTTCGGCGCAATCTCACAAACCGACATGTTCAGCTTATGTCCATTGGCGGTG‑CCATT‑GG‑TGC‑GGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGA‑TC‑ATTTTCGTTTATATGAT‑CA‑TTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATTGCTGCTTTCGAATCTGG‑AATACAAATCTTTTAGTGACTTCGC‑TTCCGAt                                                                                               <  1:792496/233‑1 (MQ=255)                                                                                                                                            gAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTG‑CCATT‑GG‑TAC‑GGGGTTTTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGA‑TC‑ATTTTCGTTTATATGAT‑CA‑TTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATTGCTGCTTTCGAATCTGG‑AATACAAATCTTTTAGTGACTTCGC‑TTCCGATTTACt                                                                                          <  1:968533/239‑1 (MQ=255)                                                                                                                                             gacagtgtgttcGGCGCAATCCCCCAAACCGGCATATTTAGCTTATTTGCATTGGGGGGT‑CCATT‑GG‑TGC‑GGGGTTGTTTTGGGGGTCTGGGAAAACGATTAGCCTTGGCGGGCCGTCGA‑TC‑ATTTTCGTTTATATGAT‑CA‑TTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATTGCTGCTTTCGAATCTGG‑AATACAAATCTTTTAGTGACTTCGC‑TTCCGATTTACTc                                                                                         <  2:743173/228‑1 (MQ=255)                                                                                                                                              aCAGGCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTG‑CCATT‑GG‑TAC‑GGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGA‑TC‑ATTTTCGTTTATATGAT‑CA‑TTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATTGCTGCTTTCGAATCTGG‑AATACAAATCTTTTAGTGACTTCGC‑TTCCGATTTACTc                                                                                         <  1:1044462/238‑1 (MQ=255)                                                                                                                                                   cGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTG‑CCATT‑GG‑TAC‑GGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGA‑TC‑ATTTTCGTTTATATGAT‑CA‑TTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATTGCTGCTTTCGAATCTGG‑AATACAAATCTTTTAGTGACTTCGC‑TTCCGATTTACTCGGGcc                                                                                    <  1:803210/238‑1 (MQ=255)                                                                                                                                                                tctcACAAACCGACATATTCAGCTTATTGCCATTGGCGGTG‑CCATT‑TG‑TAC‑GGGGGTGTTTATGGGGGCTTGCAAAAAGATTAGCCTTTCCGGGCCGTCGG‑TT‑ATTTTTGTTTATTATAT‑TA‑TTGGTTTTTTTGTCTTTTTTGTGGTTTGGGGCATTGGGGAATTGTTGCTTTTTGATCTGGTGATA‑AAAACTTTTTTTGGTTTCGC‑TTCCGTTTTATTTCGGGCGTggggggggaatt                                                                      >  2:465131/1‑230 (MQ=255)                                                                                                                                                                tctcACAAACCGACATATTCAGCTTATTGCCATTGGCGGTG‑ACATT‑GG‑TAC‑GGGGGTGTTTTTGGGGGTTGGCAAAAAGGTTAAGCTTTGCGGGCCGGCGT‑TC‑ATTTTCGTTTTTATGATCCA‑TTGGTTTT‑TTCTCTTTTTCGTGATGTGGGCCATGGGGGGATTGCTGCGTTTGTATTTGG‑TATACAAAACTTTTTTTGACTTTGCGTTCCG‑TTTACTTGAGCCGGGGGCAGGAGAtt                                                                      >  2:173164/1‑239 (MQ=255)                                                                                                                                                                       aaCCGACATATTCAGCTTATTGCCATTGGCGGTG‑CCATT‑GG‑TAC‑GGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGA‑TC‑ATTTTCGTTTATATGAT‑CA‑TTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATTGCTGCTTTCGAATCTGG‑AATACAAATCTTTTAGTGACTTCGC‑TTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCgg                                                               <  1:1222621/239‑1 (MQ=255)                                                                                                                                                                        aCCGACATATTCAGCTTATTGCCATTGGCGGTG‑CCATT‑GG‑TAC‑GGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGA‑TC‑ATTTTCGTTTATATGAT‑CA‑TTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATTGCTGCTTTCGAATCTGG‑AATACAAATCTTTTAGTGACTTCGC‑TTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGc                                                              <  1:202196/239‑1 (MQ=255)                                                                                                                                                                             tattttccgtttATTGCCTTTGGGGGGG‑CCATT‑GG‑TAC‑GGGGTTGTTTTTGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTAGA‑TC‑ATTTTCGTTTATATGAT‑CA‑TTGGTTTTATGCTCTTTTTCGTGATGCAGGCAATGGGGGAATTGCTGCTTTCGAATCTGG‑AATACAAATCTTTTAGTGACTTCGC‑TTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGa                                                          <  2:600457/228‑1 (MQ=255)                                                                                                                                                                               tatTCAGCTTATTGCCATTGGCGGTG‑CCATT‑GG‑TAC‑GGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGA‑TC‑ATTTTCGTTTATATGAT‑CA‑TTGGTTTTATGCTCTTTTTCGTGAt                                                                                                                                                                    <  1:997468/132‑1 (MQ=255)                                                                                                                                                                                          ttGCCTTTTGCGGTG‑TCATT‑GG‑TGG‑GGGGTTTTTTATGGGGTCTGGAAAAACGATTAGCCTTGCCGGGCCGTCGA‑TC‑ATTTTCGTTTATATGAT‑CA‑TTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATTGCTGCTTTCGAATCTGG‑AATACAAATCTTTTAGTGACTTCGC‑TTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTctgc                                            <  2:359873/239‑1 (MQ=255)                                                                                                                                                                                          ttGCCATTGGCGGTG‑CCATT‑GG‑TAC‑GGGTTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGA‑TC‑ATTTTCGTTTATATGAT‑CA‑TTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATTGCTGCTTTCGAATCTGG‑AATACAAATCTTTTAGTGACTTCGC‑TTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTctgc                                            <  1:610013/239‑1 (MQ=255)                                                                                                                                                                                           ttctgttGGTGGGG‑CCATT‑GT‑GAG‑GGGGTTGTTTATGGGGTTTGGCAAAACTATTAGCCTTGCCGGGCCGTCGA‑TC‑ATTTTCGTTTATATGAT‑TA‑TTGGTTTTATGCTCTGTTTCGTGATGGGGGGAATGGGGGAATTGCTGCTTTCGAATTTGG‑AATAGAAATCTTTTAGTGACTTCGC‑TTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTctgct                                           <  2:292096/234‑1 (MQ=255)                                                                                                                                                                                           tGCCATTGGCGGTG‑CCATT‑GG‑TAC‑GGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGA‑TC‑ATTTTCGTTTATATGAT‑CA‑TTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATTGCTGCTTTCGAATCTGG‑AATACAAATCTTTTAGTGACTTCGC‑TTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTctgct                                           <  1:385104/239‑1 (MQ=255)                                                                                                                                                                                             ccATTGGCGGTG‑CCATT‑GG‑TAC‑GGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGA‑TC‑ATTTTCGTTTATATGAT‑CA‑TTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATTGCTGCTTTCGAATCTGG‑AATACAAATCTTTTAGTGACTTCGC‑TTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGc                                                              <  1:728866/218‑1 (MQ=255)                                                                                                                                                                                                ttGGCGGTC‑CCATT‑GG‑TTC‑GGGGTTGTTTATGGGGACTGGCAAAACGATTAGCCTTGCCGGGCCGTGGA‑TC‑ATTTTCGTTTATATGAT‑CA‑TTGGTTTTATGCTCGTTTTCGTGATGCGGGCAATGGGGGAATTGCTGCTTTCGAATCTGG‑AATACAAATCTTTTAGTGACTTCGC‑TTCCGATTTACTCGGGCCGTGGGCTGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGtt                                      <  1:431681/239‑1 (MQ=255)                                                                                                                                                                                                    cGGTG‑CTATT‑GG‑GAC‑GGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGA‑TC‑ATTTTCGTTTATAGGAT‑CA‑TTGGTTTTATGCCCTTTTTCGTGATGCGGGCAATGGGGGAATTGCTGCTTTCGAATCTGG‑AATACAAATCTTTTAGTGACTTCGC‑TTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTa                                   <  1:910751/238‑1 (MQ=255)                                                                                                                                                                                                    cGGTG‑CCATT‑GG‑TAC‑GGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGA‑TC‑ATTTTCGTTTATATGAT‑CA‑TTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATTGCTGCTTTCGAATCTGG‑AATACAAATCTTTTAGTGACTTCGC‑TTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTa                                   <  1:119709/238‑1 (MQ=255)                                                                                                                                                                                                     ggTG‑CCATT‑GG‑TAC‑GGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGA‑TC‑ATTTTCGTTTATATGAT‑CA‑TTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATTGCTGCTTTCGAATCTGG‑AATACAAATCTTTTAGTGACTTCGC‑TTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTaa                                  <  1:1154227/238‑1 (MQ=255)                                                                                                                                                                                                     ctgtcccATG‑GG‑TAC‑GGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGA‑TC‑ATTTTCGTTTATATTAT‑CA‑TTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATTGCTGCTTTCGAATCTGG‑AATACAAATCTTTTAGTGACTTCGC‑TTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTa                                   <  1:5238/233‑1 (MQ=255)                                                                                                                                                                                                                     c‑GGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGA‑TC‑ATTTTCGTTTATATGAT‑CA‑TTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATTGCTGCTTTCGAATCTGG‑AATACAAATCTTTTAGTGACTTCGC‑TTCCGATTTACTCGGGCCGTGGGCAGGTTATTTCCCCGGCTGGACTTACTTGTTTTGCTGGGTTGTAAACGGTATGGgagaa                    >  2:621713/1‑238 (MQ=255)                                                                                                                                                                                                                       ggggTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGA‑TC‑ATTTTCGTTTATATGAT‑CA‑TTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATTGCTGCTTTCGAATCTGG‑AATACAAATCTTTTAGTGACTTCGC‑TTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCGGACg                   >  1:605724/1‑239 (MQ=255)                                                                                                                                                                                                                         ggTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGA‑TC‑ATTTTCGTTTATATGAT‑CA‑TTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATTGCTGCTTTCGAATCTGG‑AATACAAATCTTTTAGTGACTTCGC‑TTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACgtg                 >  1:1186217/1‑239 (MQ=255)                                                                                                                                                                                                                         ggTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGA‑TC‑ATTTTCGTTTATATGAT‑CA‑TTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATTGCTGCTTTCGAATCTGG‑AATACAAATCTTTTAGTGACTTCCC‑TTCCGATTTATTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTCATGGGTCTGCTGGGTTGTCACCGGTATGGCAGAccgg                 >  1:1274760/1‑236 (MQ=255)                                                                                                                                                                                                                             gttTATGGGGTTTGGCAAAACGATTAGCCTTGCCGGGCCGTCGA‑TC‑ATTTTCGTTTATATGAT‑CA‑TCGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATTGCTGCTTTCGAATCTGG‑AATACAAATCTTTTAGTGACTTCGC‑TTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGtggtgg             <  1:490106/239‑1 (MQ=255)                                                                                                                                                                                                                              tttATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGA‑TC‑ATTTTCGTTTATATGAT‑CA‑TTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATTGCTGCTTTCGAATCTGG‑AATACAAATCTTTTAGTGACTTCGC‑TTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGGGGc            >  1:1154529/1‑239 (MQ=255)                                                                                                                                                                                                                                    gggTTTGGCAAAACGATTAGGCTTGCGGGGTCGTCGA‑TC‑ATTTTCGTTTATATGAT‑CA‑TTGGTTTTATGCTCTTTTTCGGGATGCGGGCAATGGGGGAATTGGTGCTTTCGAATCTGG‑AATACAAAGCTTTTAGTGACTTCGC‑TTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCAc      <  2:1294739/239‑1 (MQ=255)                                                                                                                                                                                                                                     ggTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGA‑TC‑ATTTTCGTTTATATGAT‑CA‑TTGGTTTTATGCTCTTGTTCGGGATGCGGGCAATGGGGGAATTGCTGCTTTCGAATCTGG‑AATACAAATCTTTTAGTGACTTCGC‑TTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCAc      <  1:362593/238‑1 (MQ=255)                                                                                                                                                                                                                                        cTGGCAAAACGATTAGCCTTGCCGGGCCGTCGA‑TC‑ATTTTCGTTTATATGAT‑CA‑TTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATTGCTGCTTTCGAATCTGG‑AATACAAATCTTTTAGTGACTTCGC‑TTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTaa                                  >  1:740318/1‑207 (MQ=255)                                                                                                                                                                                                                                        cTGGCAAAACGATTAGCCTTGCCGGGCCGTCGA‑TC‑ATTTTCGTTTATATGAT‑CA‑TTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATTGCTGCTTTCGAATCTGG‑AATACAAATCTTTTAGTGACTTCGC‑TTCCGATTTACTCGGGCCGTGGGCAGGATATTTAACCGGCTGGACTTCCTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCt  >  1:372008/1‑239 (MQ=255)                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                     ATATCGATGAGTCCCGATACAGATTCGTCGTATCATAGACTGACTAAAGGCCGTAGAGCCTGAACAACACAGACAGGTACAGGAAGAAAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTG‑CCATT‑GG‑TAC‑GGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGA‑TC‑ATTTTCGTTTATATGAT‑CA‑TTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATTGCTGCTTTCGAATCTGG‑AATACAAATCTTTTAGTGACTTCGC‑TTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCT  >  NC_000913/4429770‑4430235

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