Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,421,183	A→C	15.7%	T402P (ACC→CCC)	ulaA →	L‑ascorbate‑specific enzyme IIC permease component of PTS

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,421,183	0	A	C	15.7%	176.3 / 2.2	86	T402P (ACC→CCC)	ulaA	L‑ascorbate‑specific enzyme IIC permease component of PTS
Reads supporting (aligned to +/- strand):  ref base A (45/20);  new base C (2/15);  total (47/39)
Fisher's exact test for biased strand distribution p-value = 3.98e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGAACGCCGTGGTCTGGGGCTTTATGTGGGGCACCATCGGTCAGCTGATTGCGGTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTT‑CGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAA‑GATTTGTCTGGTGATGGGGATGATCGA‑AATCTTTGGTTGCG‑TCT‑GGGCGGTGAAACTCACCGGTATGAGTGCCTGG‑ATGGGCATGGCGGACTGGTCG‑ATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATC‑GGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGTCTGCTTAATAAGGAGTTTTGATTATGACCGTACGTATTCTGG  >  NC_000913/4420947‑4421411                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                        cGAACGCCGTGGTCTGGGGCTTTATGTGGGGCACCATCGGTCAGCTGATTGCGGTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTT‑CGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAA‑GATTTGTCTGGTGATGGGGATGATCGA‑AATCTTTGGTTGCG‑TCT‑GGGCGGTGAAACTCAc                                                                                                                                                                                                                                        <  1:1073261/238‑1 (MQ=255) cGAACGCCGTGGTCTGGGGCTTTATGTGGGGCACCATCGGTCAGCTGATTGCGGTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCGGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTT‑CGCTAACCACGTCGGCGGCTGGCGTGCGGCGCTGAA‑GATTTGTCTGGTGATGGGGATGATCGA‑AATCTTTGGTTGCG‑TCT‑GGGCGGTGAAACTCAc                                                                                                                                                                                                                                        <  2:389388/238‑1 (MQ=255)  gAACGCCGTGGTCTGGGGCTTTATGTGGGTCACCATCGGTCAGCTGATTGCGGTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTT‑CGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAA‑GATTTGTCTGGTGATGGGGATGATCGA‑AATCTTTGGTTGCG‑TCT‑GGGCGGTGAAACTCACCg                                                                                                                                                                                                                                      <  1:951166/239‑1 (MQ=255)    aCGCCGTGGTCTGGGGCTTTATGTGGGGCACCATCGGTCAGCTGATTGCGGTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTT‑CGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAA‑GATTTGTCTGGTGATGGGGATGATCGA‑AATCTTTGGTTGCG‑TCT‑GGGCGGTGAAACTCACCgg                                                                                                                                                                                                                                     <  2:1442036/238‑1 (MQ=255)      gCCGTGGTCTGGGGCTTTATGTGGGGCACCATCGGTCAGCTGATTGCGGTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTT‑CGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAA‑GATTTGTCTGGTGATGGGGATGATCGA‑AATCTTTGGTTGCG‑TCT‑GGGCGGTGAAACTCACCGGTAt                                                                                                                                                                                                                                  <  2:510571/239‑1 (MQ=255)             tCTGGGGCTTTATGTGGGGCACCATCGGTCAGCTGATTGCGGTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTT‑CGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAA‑GATTTGTCTGGTGATGGGGATGATCGA‑AATCTTTGGTTGCG‑TCT‑GGGCGGTGAAACTCACCGGTATGAGTGc                                                                                                                                                                                                                            <  1:1166344/238‑1 (MQ=255)                      ttATGTGGGGCACCATCGGTCAGCTGATTGCGGTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGCCGTGTT‑CGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAA‑GGTTTTTCTGGTGATGGGGATGGTCGA‑AATCTTTGGTTGCG‑TCT‑GGGTGGTGAAACTCCACGGTATGAGTGCCTGG‑ATggtc                                                                                                                                                                                                                 >  2:950120/1‑237 (MQ=255)                       tATGTGGGGCACCCTCGGTCAGCTGGTTGCGGTTGCCATCATGGTAGCCTGCGGCTCCTCGATCCTGATTACTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGTCTCGTT‑CGCTAACCACTTCGACTGCTGGCGGGCGGTGCTGAA‑GCTTTGTGTGTGGATGGGGATGCTACA‑AATACTTGGATGGG‑TCT‑GGGTGGGGACAATCACCGGTATGAGCGccgtgctgggtg                                                                                                                                                                                                                 >  2:432246/1‑229 (MQ=255)                            ggggCACCATCGGTCAGCTGATTGCGGTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTT‑CGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAA‑GATTTGTCTGGTGATGGGGATGATCGA‑AATCTTTGGTTGCG‑TCT‑GGGCGGTGAAACTCACCGGTATGAGTGCCTGG‑ATGGGCATGGCg                                                                                                                                                                                                           <  1:606648/239‑1 (MQ=255)                             gggCACCATCGGTCAGCTGATTGCGGTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTT‑CGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAA‑GATTTGTCTGGTGATGGGGATGATCGA‑AATCTTTGGTTGCG‑TCT‑GGGCGGTGAAACTCACCGGTATGAGTGCCTGG‑ATGGGCATGGCgg                                                                                                                                                                                                          >  1:127794/1‑239 (MQ=255)                                   cATCGGTCAGCTGATTGCGGTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTT‑CGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAA‑GATTTGTCTGGTGATGGGGATGATCTA‑AATCTTTGGTTGCG‑TCT‑GGGCGGTGAAACTCACCGGTATGAGTGCCTGG‑ATGGGCATGGCGGCCTgg                                                                                                                                                                                                     >  1:27933/1‑238 (MQ=255)                                    aTCGGTCAGCTGATTGCGGTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTT‑CGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAA‑GATTTGTCTGGTGATGGGGATGATCGA‑AATCTTTGGTTGCG‑TCT‑GGGCGGTGAAACTCACCGGTATGAGTGCCTGG‑ATGGGCATGGCGGACTGGTc                                                                                                                                                                                                   <  1:1276995/239‑1 (MQ=255)                                        gTCAGCTGATTGCGGTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTT‑CGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAA‑GATTTGTCTGGTGATGGGTATGATCGA‑AATCTTTGGTTGCG‑TCT‑GGGCGGTGAAACTCACCGGTATGAGTGc                                                                                                                                                                                                                            >  1:1201428/1‑211 (MQ=255)                                        gTCAGCTGATTGCGGTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTT‑CGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAA‑GATTTGTCTGGTGATGGGTATGATCGA‑AATCTTTGGTTGCG‑TCT‑GGGCGGTGAAACTCACCGGTATGAGTGc                                                                                                                                                                                                                            <  2:1201428/211‑1 (MQ=255)                                        gTCAGCTGATTGCGGTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTT‑CGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAA‑GATTTGTCTGGTGATGGGGATGATCGA‑AGTCTTTGGTTGCG‑TCT‑GGGCGGTGAAATTCACCGGTATGAGTGCCTGG‑ATGGGCATGGCGGACTGGTGG‑Att                                                                                                                                                                                              >  1:2262/1‑239 (MQ=255)                                                       ttGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTT‑CGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAA‑GATTTGTCTGGTGATGGGGATGATCGA‑AATCTTTGGTTGCG‑TCT‑GGGCGGTGAAACTCACCGGTATGAGTGCCCGG‑ATGGGCATGGCGGACTGGTCG‑ATTCTGGCACCGCCgatg                                                                                                                                                                               >  1:1362320/1‑239 (MQ=255)                                                          gCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTT‑CGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAA‑GGTTTGTCTGGTGATGGGGATGATCGA‑AATCTTTGGTTGCG‑TCT‑GGGCGGTGAAACTCACCGGTATGAGTGCCTGG‑ATGGGCATGGCGGACTGGTCG‑ATTCTGGCACCGCCgatgat                                                                                                                                                                             >  1:1128544/1‑238 (MQ=255)                                                                tGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTT‑CGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAA‑GATTTGTCTGGTGATGGGGATGATCGA‑AATCTTTGGTTGCG‑TCT‑GGGCGGTGAAACTCACCGGTATGAGTGCCTGG‑ATGGGCATGGCGGACTGGGCG‑ATTCTGGCACCGCCGATGGTGCAAg                                                                                                                                                                        >  2:514913/1‑237 (MQ=255)                                                                   tCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGGT‑CGCTAACCACTTCGGCGGATGGCGTGCGGCGCTGAA‑GATTTGTCTGGTGATGGGGGTGGTCGA‑AATCTTTGGTTGCG‑TCT‑GGGCGGTGACACTCACCGGTATGGGTGCCAGG‑GTTGGCCTGGCGGACTGGCTG‑GTTCTGGCACAGCCGCTGATGCAAGGct                                                                                                                                                                     >  2:1194661/1‑237 (MQ=255)                                                                          cGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTT‑CGCTGACCACTTCGGCGGCTGGCGTGCGGCGCTGAA‑GATTTGTCTGGTGATGGGGATGATCGA‑AATCTTTGGTTGCG‑TCT‑GGGCGGTGGAACTCACCGGTATGAGTGCCTGG‑ATGGGCATGGCGGGCTGGTCG‑ATTCTGGCACCGCCGATGATGCAAGGCTTCTtctc                                                                                                                                                              >  2:1288639/1‑237 (MQ=255)                                                                                 tccaTCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTT‑CGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAA‑GATTTGTCTGGTGATAGGGATGATCGA‑AATCTTTGGTTGCG‑TCT‑GGGCGGTGAAACTCACCGGTATGAGTGCCTGG‑ATGGGCATGGCGGACTGGTCG‑ATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATC‑GGTa                                                                                                                                                     >  1:94118/4‑238 (MQ=255)                                                                                              ttCCTGGCTTTATCCCGATTTTCTTCTCTAACGCCACCCTCGGCGTGTT‑CGCTACCCTCTCCGGCGGCGGGCGTGGGGCGCTGAA‑GTTTTGTGTGGTGATGGGGATGATCGA‑AATCTTTGGTTGCG‑TCT‑GGGCGGTGAACCTCCCCGGTATGAGTGCCTGG‑ATGGGCATGGCGGTCTGGTCG‑ATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATC‑GGTATCGCCTTTATGGcc                                                                                                                                       <  1:152190/239‑1 (MQ=255)                                                                                               tCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTT‑CGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAA‑GATTTGTCTGGTGATGGGGATGATCGA‑AATCTTTGGTTGCG‑TCT‑GGGCGGTGAAACTCACCGGTATGAGTGCCTGG‑ATGGGCATGGCGGACTGGTCG‑ATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATC‑GGTATCGCCTTTATGGCCg                                                                                                                                      >  1:841934/1‑239 (MQ=255)                                                                                               tCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTT‑CGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAA‑GATTTGTCTGGTGATGGGGATGATCGA‑AATCTTTGGTTGCG‑TCT‑GGGCGGTGAAACTCACCGGTATGAGTGCCTGG‑ATGGGCATGGCGGACTGGTCG‑ATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATC‑GGTATCGCCTTTATGGCCg                                                                                                                                      >  1:1103176/1‑239 (MQ=255)                                                                                                    gCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTT‑CGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAA‑GATTTGCCTGGTGATGGGGATGATCGA‑AATCTTTGGTTGCG‑TCC‑GGGCGGTGAAACTCACCGGTATGAGTGCCTGG‑ATGGGCATGGCGGACTGGTCG‑ATTCTGGCACCGCCGATGCATCAAGGCGTCTTATCCATC‑GGGATCGCCGTTCTGGCCGtcatc                                                                                                                                 >  2:273234/1‑239 (MQ=255)                                                                                                      tttACCCCGAGTTTCTTCTCTACCCCCACCATCGGCGTGTC‑CGTTACCCCCTGCGGCGGCGGGCGGGCGGCGTTGAA‑GATTTGTTGGGTGATGGGGATGATCGA‑ATTTTTTTTTTGCG‑TCT‑GGGCGGCGAAACTCGCCGGTAGGAGTGCCTGG‑AGGGTCATGGCGGTCTGGTCG‑ATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATC‑GGTATCGCCTTTATGGCCGTcatcat                                                                                                                               <  2:858204/239‑1 (MQ=255)                                                                                                        tATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTT‑CGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAA‑GATTTGTCTGGTGATGGGGATGATCGA‑AATCTTTGGTTGCG‑TCT‑GGGCGGTGAAACTCACCGGTATGAGTGCCTGG‑ATGGGCATGGCGGACTGGTCG‑ATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATC‑GGTATCGCCTTTATGGCCGgcacacttg                                                                                                                             >  2:1165399/1‑233 (MQ=255)                                                                                                         aTCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTT‑CGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAA‑GATTTGTCTGGTGATGGGGATGATCGA‑AATCTTTGGTTGCG‑TCT‑GGGCGGTGAAACTCACCGGTATGAGTGCCTGG‑ATGGGCATGGCGGACTGGTCG‑ATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATC‑GGTATCGCCTTTATGGCCGTCATCATTGt                                                                                                                            >  1:180850/1‑239 (MQ=255)                                                                                                         aTCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTT‑CGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAA‑GATTTGTCTGGTGATGGGGATGATCGA‑AATCTTTGGTTGCG‑TCT‑GGGCGGTGAAACTCACCGGTATGAGTGCCTGG‑ATGGGCATGGCGGACTGGTCG‑ATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATC‑GGTATCGCCTTTATGGCCGTCATCATTGt                                                                                                                            >  1:1047527/1‑239 (MQ=255)                                                                                                              gATGTTCTTCTCTAACGCCACCATCGGCGTGTT‑CGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAA‑GATTTGTCTGGTGATGGGGATGATCGA‑AATCTTTGGTTGCG‑TCT‑GGGCGGTGAAACTCACCGGTATGAGTGCCTGG‑ATGGGCATg                                                                                                                                                                                                              <  1:539453/154‑1 (MQ=255)                                                                                                              gATGTTCTTCTCTAACGCCACCATCGGCGTGTT‑CGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAA‑GATTTGTCTGGTGATGGGGATGATCGA‑AATCTTTGGTTGCG‑TCT‑GGGCGGTGAAACTCACCGGTATGAGTGCCTGG‑ATGGGCATg                                                                                                                                                                                                              >  2:539453/1‑154 (MQ=255)                                                                                                              gATGTTCTTCTCTAACGCCACCATCGGCGTGTT‑CGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAA‑GATTTGTCTGGTGATGGGGATGATCGA‑AATCTTTGGTTGCG‑TCT‑GGGCGGCGAAACTCACCGGTATGAGTGCCTGG‑ATGGGCATGGCGGACTGGTCG‑ATTCTGGCACCGCCGCTGATGCAAGGCTTCTTCTCCATC‑GGTATCGCCTTTATGGCCGTCATCATTGTAAtt                                                                                                                        >  1:425737/1‑238 (MQ=255)                                                                                                               aTGTTCTTCTCTAACGCCACCATCGGCGTGTT‑CGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAA‑GATTTGTCTGGTGATGGGGATGATCGA‑AATCTTTGGTTGCG‑TCT‑GGGCGGTGAAACTCACCGGTATGAGTGCCTGG‑ATGGGCATGGCGGACTGGTCG‑ATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATC‑GGTATCGCCTTTATGGCCGTCATCATTGTAATTGc                                                                                                                      >  1:634121/1‑239 (MQ=255)                                                                                                                ttgtctcccctACCGCCACCATCGGC‑TGTTCCTCTAACCTGGCGGGCGGTTGGCGTGGGGCGCTAAG‑GTTTTGTTTGGTGTGGGGGATGATCTA‑ATTTTTTGG‑TGCGTTGG‑GGGCGTTGAAACCCCCCGGTATGGGTGCCTGG‑ATGGGCTGGGTGGACTGGTCG‑ATTTGGGCACCGCCGATGATGCAAGGCTTCTTCTCCATC‑GGTATCGCCTTTATGGCCGTCATCATTGTAATTGCa                                                                                                                     <  2:993258/230‑1 (MQ=255)                                                                                                                tgtgttcctcaAACCCCACTAGCGGCGTTTT‑CGCTCACCTCTTCGGCGGTGGGCTGGCGGCGTTGAA‑GTTTTGTCTGGGGATGGGGAGGATCAA‑AATTTTTGGTTGCG‑TCG‑GGGCGGTGAACCTCCCCGGTATGAGTGCCTGG‑ATGGGCTTGGCGGACTGGTCT‑ATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATC‑GGTATCGCCTTTATGGCCGTCATCATTGTAATTGCa                                                                                                                     <  2:407825/232‑1 (MQ=255)                                                                                                                 tttcttcACTACCTCCACGATGGGGGGGTT‑CGCTAACCGCTGCGCCGGCGGGCGTGCGGCGCTGAA‑GGTTTGTCGGGGGATGGGGATGATAGA‑ATTCTTTGGTTGCG‑TCG‑GGGGGGTAAAACTCCCCGGTATGAGGCCGTGG‑AGGGGCATGGCGGCTTGGTCG‑ATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATC‑GGTATCGCCTTTATGGCCGTCATCATTGTAATTGCa                                                                                                                     <  2:1025157/237‑1 (MQ=255)                                                                                                                 gTTCTTCTCTAACCCCCCCAGGGGATTTTG‑CGCAAACCCCTTCGGCGGTTGGGGGGCGGCGTTGAA‑GATTTGTCGGGTGAGGGGGATGATAAA‑AACTTTTGGTTGCG‑TCT‑GGGGGGTAAACCTCCCCGTTAGGAGTGCTTGG‑TTGGGC‑TGGCGGGATTGTCGTTTTCTGTCACCGCCGATGATGCAAGGCTTCTTCTCCATC‑GGTATCGCCTTTATGGCCGTCATCATTGTAATTGCa                                                                                                                     <  2:380172/238‑1 (MQ=255)                                                                                                                    cttctCTAACGCCACCATCGGCGTGTT‑CGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAA‑GATTTGTCTGGTGATGGGGATGATCGA‑AATCTTTGGTTGCG‑TCT‑GGGCGGTGAAACTCACCGGTATGAGTGCCTGG‑ATGGGCATGGCGGACTGGTCG‑ATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATC‑GGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTg                                                                                                                  >  1:955380/1‑238 (MQ=255)                                                                                                                    cttctCTAACGCCACCATCGGCGTGTT‑CGCTAACCACTTCGGCCGCTGGCCTGCGGCGCTGAA‑GATTTGTCTGGTGATGGGGATGATTGA‑AATCTTTGGTTGCG‑TCT‑GGGCGGTGTAACTCCCCGGTATGAGAGCCTGG‑GTGGGCATGGCGGACTGGTCG‑ATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATC‑GGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTg                                                                                                                  >  1:661189/1‑238 (MQ=255)                                                                                                                          tAACGCCACCATCGGCGTGTT‑CGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAA‑GATTTGTCTGGTGATGGGGATGATCGA‑AATCTTTGGTTGCG‑TCT‑GGGCGGTGAAACTCACCGGTATGAGTGCCTGG‑ATGGGCATGGCGGGCTGGTCG‑ATTCTGGCACCGCCCATGATGCCAGGCTTCTTCTCCACC‑GGTCTCGCCTTTATGGCCGTCATCATTGTAATTGCCCCGGCTtata                                                                                                           >  2:495076/1‑239 (MQ=255)                                                                                                                          tAACGCCACCATCGGCGTGTT‑CGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAA‑GATTTGTCTGGTGATGGGGATGATCGA‑AATCTTTGGTTGCG‑TCT‑GGGCGGTGAAACTCACCGGTATGAGTGCCTGG‑ACGGGCATGGCGGACTGGTCG‑ATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATC‑GGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTtata                                                                                                           >  1:674678/1‑239 (MQ=255)                                                                                                                            aCGCCACCATCGGCGTGTT‑CGCTAACCACTTCGGCGGCTTGCGTGCGGCGCTGAA‑GATTTGTCTTGTGATAGGGATGATCGA‑AATCTTTGGTTGCG‑TCT‑AGGCGGGGGAACTCACCGGTATGAGTGCCTGG‑TTGGGCATGGAGGACTGGTCG‑ATTCAAGCACCGCCGATAATGCAAGGCTTCTTCTCCCTC‑GGATTCTCCTTTATGGCCTTCGTTCTTTGAATTGCACTCGCttgaatg                                                                                                         >  2:250479/1‑234 (MQ=255)                                                                                                                            aCGCCACCATCGGCGTGTT‑CGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAA‑GATTTGTCTGGTGATGGGGATGATCGA‑AATCTTTGGTTGCG‑TCT‑GGGCGGTGAAACTCACCGGTATGAGTGCCTGG‑ATGGGCATGGCGGACTGGTCG‑ATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATC‑GGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATg                                                                                                         >  1:210955/1‑239 (MQ=255)                                                                                                                                 accaTCGGCGTGTT‑CGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAA‑GATTTGTCTGGTGATGGGGATGATCGA‑AATCTTTGGTTGCG‑TCT‑GGGCGGTGAAACTCACCGGTATGAGTGCCTGG‑ATGGGCATGGCGGCCTGGTCG‑ATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATC‑GGTATCGCCTTTATGGCCGGCATCATTGTAATTGCACTGTCTTATTTGttctt                                                                                                    >  2:658104/1‑239 (MQ=255)                                                                                                                                 accaTCGGCGTGTT‑CGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAA‑GATTTGTCTGGTGATGGGGATGATCGA‑AATCTTTGGTTGCG‑TCT‑GGGCGGTGAAACTCACCGGTATGAGTGCCTGG‑ATGGGCATGGCGGACTGGTCG‑ATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATC‑GGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGttct                                                                                                     >  1:530638/1‑238 (MQ=255)                                                                                                                                             tt‑CGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAA‑GATTTGTCTGTTGATGGGGATGATCGA‑AATCTTTGGTTGCG‑TCT‑GGGCGGTGAAACTCACCGGGATGAGTGCCTGG‑GTGGGCATGGCGGGCTGGTCG‑ATTCTTGCCCCGCCGGTGATGCAAGGCCTCTTCTCCATC‑GGTATCGCCTTTATGGCCGTCACCCTTGTAATTTCCCTGGCTTATTTTTTCTTCGCTTGCcgcgc                                                                                        >  2:739527/1‑239 (MQ=255)                                                                                                                                             tt‑CCCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAA‑GATTTTTCTGGTGATGGGGATGATCGG‑AATCTTTGGTTGCG‑TCT‑GGGCGGTGAAACTCACCGGTATGAGTTCCTGG‑GTGGGCCTGGCGGGCTGGTCG‑ATTCTTGCCCCGCCGCTGATTCAAGGCTTCTTCTCCAAC‑GGGATCCCCTTTTTGGCCGGCCTCATTGTCATTTCCCTGGGCTTTTTGttctttcgttgcccgg                                                                                         >  2:1105453/1‑227 (MQ=255)                                                                                                                                                    accatcTTGGGGGGGTGGCGTGGGGCG‑TGAAGGTTTTGTGTGGTGGTGGGGATAATCAA‑ATTTTTTGGTGGGG‑TCT‑GGGGGGTAAACCCCCCCGTTAGGAGTGCCTGG‑TGGGGCGTGGCGGAGTGGTCG‑ATTCGGGCACCGCCGATGATGCAAGGCTTCTTCTCCATC‑GGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTgcg                                                                                  <  2:27933/234‑1 (MQ=255)                                                                                                                                                    aaCCACTGCGGCGGCTGGCGTGCGGCGTGGAA‑GATTTGTCTGGTGATGGGGATGATAGA‑AATCTTTGGTTGCG‑TCG‑GGGCGGTGAAACTCACCGGTATGAGTGCCGGG‑ATGGGCAGGGCGGACTGGTCG‑ATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATC‑GGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTgcg                                                                                  <  1:948820/239‑1 (MQ=255)                                                                                                                                                     acctcgtcgtgggcgggcgggcgGCGTTGAG‑GATTTGGCTGGTGAGGGGGATG‑TCGAGAATCTGGGGTTGCG‑TCT‑GGGCGGTAAACCTCGCCGGAGTGAGGTCGTGG‑AGGGGCATGGCGGAGTGTGCT‑ATCTTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATC‑GGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTgcgc                                                                                 <  2:406506/228‑1 (MQ=255)                                                                                                                                                        aCTTCGGCGGCTGGCGTGCGGCGCTGAA‑GATTTGTCTGGTGATGGGGATGATCGA‑AATCTTTGGTTGCG‑TCT‑GGGCGGTGAAACTCACCGGTATGAGTGCCTGG‑ATGGGCATGGCGGACTGGTCG‑ATTCTGGCACCGCCGATGAAGCAAGGCTTCTTCTCCATC‑GGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCa                                                                              >  1:435643/1‑239 (MQ=255)                                                                                                                                                             ggcggcTGGCGTGCGGCGCTGAA‑GATTTGTCTGGTGATGGGGATGATCGA‑AATCTTTGGTTGCG‑TCT‑GGGCGGTGAAACTCACCGGTATGAGTGCCTGG‑ATGGGCATGGCGGACTGGTCG‑ATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATC‑GGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCagaaga                                                                         >  1:328754/1‑239 (MQ=255)                                                                                                                                                                 gcTGGCGTGCGGCGCTGAA‑GATTTGTCTGGTGATGGGGATGATCGA‑AATCTTTGGTTGCG‑TCT‑GGGCGGTGAAACTCACCGGTATGAGTGCCTGG‑ATGGGCATGGCGGACTGGTCG‑ATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATC‑GGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGAt                                                                     >  1:763887/1‑239 (MQ=255)                                                                                                                                                                             cgtggAA‑GTTTTGTTGGGGGAGGGGGATGTTGAA‑ATTTTTTGTTTGGG‑TTG‑GGGGGGTAAACCCCACGGTTGTGGGGGCGTGG‑AGGGGAGGGGCGTAGTGGGCG‑TTTGGGGCACCGCCGATGATGCAAGGCTTCTTCTCCATC‑GGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAaaca                                                          <  2:1157736/234‑1 (MQ=255)                                                                                                                                                                             cggtGAG‑GATTTGTCGGGTGATGGGCGTTATCGA‑ATTCTTCGGTTGAT‑TAG‑GGGCGGTGAACCTCACCGGTACGAGT‑CCTGGCGTGGGCATGGCGGCGTGGTCG‑TTTCGGGCGCCGCCGATGATGCAAGGCTTCTTCTCCATC‑GGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGCAGAAGATGCAGAAAaaca                                                          <  2:690293/235‑1 (MQ=255)                                                                                                                                                                              ggtagaggATTTGGCGGGGGAGGGGGAGGATCAA‑AATTTTGGGTGGCG‑TCG‑GGGGGGTGAAACTCCCCGGTATGAGTGCCGGG‑AGGGGGGGGGCGGATGGGTCG‑ATTCGGGCACCGCCGATGATGCAAGGCTTCTTCTCCATC‑GGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAaaca                                                          <  1:1205967/231‑1 (MQ=255)                                                                                                                                                                                  aa‑GTTTTGTTGGGGGAGGGGGATGATC‑A‑AATTTTTGGGTGGG‑TTTGGGGCGTTAAACCTCGCCGGTAGGGGTGCCGGG‑ATGGGCAGGGCGGACTGGTCG‑TTTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATC‑GGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTgg                                                     <  2:137582/238‑1 (MQ=255)                                                                                                                                                                                         tGTCTGGGGAGGGGGAGGAGCGA‑AATCTTGGGTTGCG‑TCT‑GGGGGGTGAAACTCCCCGGTATGAGTGCGTGG‑AGGGGCATGGCGGACTGGTCG‑TTTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATC‑GGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACa                                              <  1:457578/239‑1 (MQ=255)                                                                                                                                                                                           ttggggggtGGGGTTGAGCAA‑AATTTTTGGTTGGG‑GCG‑GGGGGGGGAACCTCCCGGGTATGAGGGCGTGG‑ATGGGCATGGGGGTCTGGTCG‑ATTTGGGCACCGCCGATGATGCAAGGCTTCTTCTCCATC‑GGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGt                                            <  2:722293/231‑1 (MQ=255)                                                                                                                                                                                            tgtggggtGGGAAAGATAGA‑AGTTTGTGGTTGCG‑TGT‑GGGGGGGGAAACTCGCGGGTATGAGGGCTGGG‑GTGGGTAGGGCTGGGTGGTCT‑ATTGTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATC‑GGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGt                                            <  2:335211/231‑1 (MQ=255)                                                                                                                                                                                            ggggggtgggggATGAGCAA‑AGTTTTTGGTTGCT‑TCG‑GGGGGTTAAACCTCCCCGGTATGGGGGCTGGG‑GGGGGGGGGGGGGACTGGTCT‑ATTCGGGCCCCGCCGATGATGCAAGGCTTCTTCTCCATC‑GGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGt                                            <  2:841934/230‑1 (MQ=255)                                                                                                                                                                                             tGGTGATGGGGATGATCGA‑AATCTTTGGTTGCG‑TCT‑GGGCGGTGAAACTCAc                                                                                                                                                                                                                                        <  2:713393/52‑1 (MQ=255)                                                                                                                                                                                             tGGTGATGGGGATGATCGA‑AATCTTTGGTTGCG‑TCT‑GGGCGGTGAAACTCAc                                                                                                                                                                                                                                        >  1:713393/1‑52 (MQ=255)                                                                                                                                                                                             tGGTGATGGGGATGATCGA‑AATCTTTGGTTGCG‑TCT‑GGGCGGTGAAACTCACCGGTATGAGTGCCTGG‑ATGGGCATGGCGGACTGGTCG‑ATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATC‑GGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAaacaac                                                        >  2:169643/1‑225 (MQ=255)                                                                                                                                                                                             tGGTGATGGGGATGATCGA‑AATCTTTGGTTGCG‑TCT‑GGGCGGTGAAACTCACCGGTATGAGTGCCTGG‑ATGGGCATGGCGGACTGGTCG‑ATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATC‑GGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAaacaac                                                        <  1:169643/225‑1 (MQ=255)                                                                                                                                                                                               gTGATGGGGATGATCGA‑AATCTTTGGTTGCG‑TCT‑GGGCGGTGAAACTCACCGGTATGAGTGCCTGG‑ATGGGCATGGCGGACTGGTCG‑ATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATC‑GGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGTCt                                          >  2:722135/1‑237 (MQ=255)                                                                                                                                                                                                tGATGGGGATGATCGA‑AATCTTTGGTTGCG‑TCT‑GGGCGGTGAAACTCACCGGTATGAGTGCCTGG‑ATGGGCATGGCGGACTGGTCG‑ATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATC‑GGGATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGACGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGTCTGCt                                       >  2:1416610/1‑239 (MQ=255)                                                                                                                                                                                                         tgatCGA‑ATTCTTTGGTGGGG‑TGG‑GGGTGGTAAAACCCACGGGTAGGGGTGCCTGG‑GTGGGCAGGGCGGCTTGGTCT‑ATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATC‑GGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGTCTGCTTAATAAgg                               <  2:1100627/238‑1 (MQ=255)                                                                                                                                                                                                         tgatCGA‑AATCTTTGGTTGCG‑TCT‑GGGCGGTGAAACTCACCGGTATGAGTGCCTGG‑ATGGGCATGGCGGACTGGTCG‑ATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATC‑GGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGTCTGCTTAATAAGGa                              >  1:1288730/1‑239 (MQ=255)                                                                                                                                                                                                          gatCGA‑AATCTTTGGTTGCG‑TCT‑GGGCGGTGAAACTCACCGGTATGAGTGCCTGG‑ATGGGCATGGCGGACTGGTCG‑ATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATC‑GGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGTCTGCTTAATAAGGAg                             >  1:388572/1‑239 (MQ=255)                                                                                                                                                                                                            tcaa‑aaTCTTGGGTTGCG‑TCT‑GGGGGGTAAAACCCACGGATATGGGTGCCTGG‑AGGGGCTTGGCGGACTGGTCG‑TTTCGGGCACCGCCGATGATGCAAGGCTTCTTCTCCATC‑GGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGTCTGCTTAATAAGGAGtt                           <  1:950120/236‑1 (MQ=255)                                                                                                                                                                                                            tCGA‑AATTTTGGGTTGCG‑TCT‑GGGGGGTAAAACCCCCCGGTATGAGTGCCTGG‑AGGGGCATGGGGGACTGGTCG‑ATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATC‑GGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGTCTGCTTAATAAGGAGtt                           <  1:514913/239‑1 (MQ=255)                                                                                                                                                                                                             cGA‑ATTCTGTGGTGTGG‑TTT‑GGGCGGTGAAACTCACGTGGATGGGGCCCGGG‑ATGGGCATGGCGGACTGGTCG‑ATTCTGGATCCGCCGATGATGCAAGGCTTCTTCTCCATC‑GGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGTCTGCTTAATAAGGAGtt                           <  1:1082207/238‑1 (MQ=255)                                                                                                                                                                                                             cGA‑AATTTTGGTTTGCT‑TCG‑GGGGGGTGAAACTCACCGGTATGAGTGCCGGG‑AGGGGCATGGCGGACTGGTCG‑ATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGTCTGCTTAATAAGGAGtt                           <  1:1056608/239‑1 (MQ=255)                                                                                                                                                                                                               gaattgtttGGTTGGG‑TCG‑GGGGGGTGAACCTCCCCGGGTTGAGGGCCTGG‑ATGGGTGTGGGGGACTGGTGG‑TTTCGGGCACCGCCGATGATGCAAGGCTTCTTCTCCATC‑GGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGTCTGCTTAATAAGGAGtt                           <  2:634121/231‑1 (MQ=255)                                                                                                                                                                                                                 aaTCTTTGGTTGCG‑TCT‑GGGCGGTGAAACTCACCGGTATGAGTGCCTGG‑ATGGGCATGGCGGACTGGTCG‑ATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATC‑GGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGTCTGCTTAATAAGGAGTTTTGa                       >  2:1129694/1‑239 (MQ=255)                                                                                                                                                                                                                      ttGGTTGCG‑TCT‑GGGCGGTGAAACTCACCGGTATGAGTGCCTGG‑ATGGGCATGGCGGACTGGTCG‑ATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATC‑GGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGTCTGCTTAATAAGGAGTTTTGATTATg                  >  1:701954/1‑239 (MQ=255)                                                                                                                                                                                                                      ttGGTTGCG‑TCT‑GGGCGGTGAAACTCACCGGTATGAGTGCCTGG‑ATGGGCATGGCGGACTGGTCG‑ATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATC‑GGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGTCTGCTTAATAAGGAGTTTTGATTATg                  >  1:110383/1‑239 (MQ=255)                                                                                                                                                                                                                       tGGTTGCG‑TCT‑GGGCGGTGAAACTCACCGGTATGAGTGCCTGG‑ATGGGCATGGCGGACTGGTCG‑ATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATC‑GGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGTCTGCTTAATAAGGAg                             >  1:1334641/1‑227 (MQ=255)                                                                                                                                                                                                                       tGGTTGCG‑TCT‑GGGCGGTGAAACTCACCGGTATGAGTGCCTGG‑ATGGGCATGGCGGACTGGTCG‑ATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATC‑GGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGTCTGCTTAATAAGGAGTTTTGATTATg                  >  1:1230529/1‑238 (MQ=255)                                                                                                                                                                                                                         tttGCG‑TCT‑GGGGGGTAAAACCCCCGGGTATGAGGGCCTGG‑AGGGGCTTGGCGGATTGTGCG‑TTTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATC‑GGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGTCTGCTTAATAAGGAGTTTTGATTATGAc                <  1:625161/237‑1 (MQ=255)                                                                                                                                                                                                                         ggtGCG‑TCG‑GGGCGGGAAACCGCCCGGGTAGGGGTCGCTGT‑GTGGGCTGGGCGTCCTGGTCT‑ATTCTGGCCCCGCCGATGATGCAAGGCTTCTTCTCCATC‑GGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGTCTGCTTAATAAGGAGTTTTGATTATGAc                <  2:127794/236‑1 (MQ=255)                                                                                                                                                                                                                            gCG‑TCG‑GGCCGGTAAAACCCACCTGTATGAGTCCCTGT‑ATGGGCTTGGCGGCCTTGTAT‑TTTCGGGATCCGCCGATGATGCAAGGCTTCTTCTCCATC‑GGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGTCTGCTTAATAAGGAGTTTTGATTATGACcgta            <  2:435643/239‑1 (MQ=255)                                                                                                                                                                                                                                         gTGAAACTCCCCGGTATGAGTGCGTGG‑ATGGGCATGGCGGACTGGTCG‑ATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATC‑GGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGTCTGCTTAATAAGGAGTTTTGATTATGACCGTACGTATTCTgg  <  1:937599/238‑1 (MQ=255)                                                                                                                                                                                                                                          tGAAACTCACCGGTATGAGTGCCTGG‑ATGGGCATGGCGGACTGGTCG‑ATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATC‑GGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGTCTGCTTAATAAGGAGTTTTGATTATGACCGTACGTATTCTgg  >  1:1138145/1‑237 (MQ=255)                                                                                                                                                                                                                                          tGAAACTCACCGGTATGAGTGCCCGG‑ATGGGCATGGCGGACTGGTCG‑ATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATC‑GGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGTCTGCTTAATAAGGAGTTTTGATTATGACCGTACGTATTCTg   >  1:962922/1‑236 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                        CGAACGCCGTGGTCTGGGGCTTTATGTGGGGCACCATCGGTCAGCTGATTGCGGTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTT‑CGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAA‑GATTTGTCTGGTGATGGGGATGATCGA‑AATCTTTGGTTGCG‑TCT‑GGGCGGTGAAACTCACCGGTATGAGTGCCTGG‑ATGGGCATGGCGGACTGGTCG‑ATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATC‑GGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGTCTGCTTAATAAGGAGTTTTGATTATGACCGTACGTATTCTGG  >  NC_000913/4420947‑4421411

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