Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	802,118	0	A	C	44.8%	91.6 / 13.0	69	T78P (ACC→CCC)	ybhI	putative DASS family tricarboxylate or dicarboxylate transporter
Reads supporting (aligned to +/- strand):  ref base A (24/13);  new base C (0/30);  total (24/45)
Fisher's exact test for biased strand distribution p-value = 3.74e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGAACTGGCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACCTG‑TCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAA‑CTTATCCG‑ACGGTGCGTTTAAAA‑CCA‑CCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGCTGGGTC  >  NC_000913/801972‑802265                                                                                                                                                      |                                                                                                                                                      gCGAACTGGCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACCTG‑TCGTACTGTTAATTGCCGTTGCTGCCTCAATGGGGGCGGTCGGTAA‑CTTATCCG‑ACGGTGCGTTTAAAa‑cca‑c                                                                                                                                                 >  1:1178485/1‑151 (MQ=255)  cGAACTGGCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACCTG‑TCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAA‑CTTATCCG‑ACGGTGCGTTTAAAa‑cca‑cc                                                                                                                                                <  1:381982/151‑1 (MQ=255)   gAACTGGCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACCTG‑TCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAA‑CTTATCCG‑ACGGTGCGTTTAAAA‑CCA‑ccg                                                                                                                                               <  1:616501/151‑1 (MQ=255)             gggTGTTTTTGGGTTTTTCCCGGGGGGCTTTTGTGGGGCTGTTTATCAAGCCTTTCCCGGAACCTC‑TCGTACTGTTAATTGCCGTTGCTCCCTAAATGGTGGTGGTCGGTAC‑CTTACCCG‑GCGGTGCTTTTAAAC‑CCA‑CCGCCGTATTAAg                                                                                                                                     <  1:201813/151‑1 (MQ=255)                 gCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACCTG‑TCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAA‑CTTATCCG‑ACGGTGCGTTTAAAA‑CCA‑CCGCCGTATTAAGCGGt                                                                                                                                 >  2:944943/1‑151 (MQ=255)                   tttttGGTTTTTCCCGGGCGGCTATTGTGGGGCTGTTTATCAAGCTTTCCCCGGACCCTG‑TCGTACTTTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAA‑CTTACCCG‑ACGGTGCGTTAAAAC‑CCA‑CCGCCGTATTAAGCgg                                                                                                                                  <  1:94654/148‑1 (MQ=255)                      tgggTATTTACCTGGCGGCCATGGGGGGGTTGGTTATCAAGCTTTTCCCGGACCCTG‑TCGTACTGTTAATTGCCGTTGCTGCCTAAATGGTGGGGGTCGGTAA‑CTTATCCG‑ACGGTGCGTTTAAAA‑CCA‑CCGCCGTATTAAGCGGTTAct                                                                                                                             <  2:287542/148‑1 (MQ=255)                        ggTTTTTACCTGGCGGCCTTTGTGGGGCTGGTTATCAACCCTTTCCCGGACCCTG‑TCGTCTTTTTATTTGCGTTTGCTGCCTCAATGTTGGGGGTGGGTAA‑TTTATCCG‑ACGGTGCTTTTAAAC‑CCA‑CCGCCGTATTAAGCGGTTActct                                                                                                                           <  2:771810/150‑1 (MQ=255)                               cccgggcgccCTTTGGGGGTCGGTTTATCAGGCCTTCCCCGAAACCTG‑TCGTTCTTTTATTCCCCGTTGCTCCCTCAAGGGTGGTGGTGGGTAA‑CTTTCCCG‑ACGGGGCTTTAAAAC‑CCA‑CCGCCGTATTAAGCGGTTACTCTTCAGGTa                                                                                                                    <  1:768439/146‑1 (MQ=255)                                cggggtggcCATTGGGGGGCTGGTTATC‑ACCCTTTCCCCGGCGCTGCGGGTCTTGTTAATTCCCTTTGCTCCCTAATGGGTGGGGGCGGGTAC‑TTTTTCCG‑CCTGGGCTTTTAAAA‑CCC‑CCGCCGTATTAAGCGGTTACTCTTCAGGTac                                                                                                                   <  2:780619/144‑1 (MQ=255)                                ccgggcgccCATTGTGGGGCTGTTTATCAGGCCTTCCCCGGAACCGG‑TCTTACTTTTAATTGCCGTTGCTGCCTAAAGGGGGGGGGTCGGTAC‑CTTACCCG‑AGGGGGTTTTTAAAC‑CCA‑CCGCCGTATTAAGCGGTTACTCTTCAGGTa                                                                                                                    <  1:172327/146‑1 (MQ=255)                                 cTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACCTG‑TCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAA‑CTTATCCG‑ACGGTGCGTTTAAAA‑CCA‑CCGCCGTATTAAGCGGTTACTCTTCAGGTacc                                                                                                                  >  1:378036/1‑150 (MQ=255)                                  tGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACCTG‑TCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAA‑CTTATCCG‑ACGGTGCGTTTAAAA‑CCA‑CCGCCGTATTAAGCGGTTACTCTTCAGGTaccac                                                                                                                >  2:1053578/1‑151 (MQ=255)                                   gggggcCATTGTGGGGCTGGTTATCAAGCCTTCCCCGGACCTTT‑TCGTACTGTTAATTGCCGTTGCTGCCTAAATGGTGGTGGTCGGTAA‑CTTATCCG‑AGGGTGCTTTAAAAC‑CCA‑CCGCCGTATTAAGCGGTTACTCTTCAGGTaccacc                                                                                                               <  2:992280/148‑1 (MQ=255)                                        cctttGGGGGGTTGTTTATCAGGCCTTCCCCGGAACCCT‑TCGTACTGTTATTTGCGGTTGCCCCCAAAATGGTGGTGGGCGGTAA‑TTTATCGG‑ACGGTGCTTTAAAAC‑CCA‑CCGCCGTATTAAGCGGTTACTCTTCAGGTACCACct                                                                                                              <  1:776793/144‑1 (MQ=255)                                           ttGTGGGGTTGGTTTTCAAGCCTTGCCCGGACCTTG‑TCTTACTGTAAATGGCCGTGGCTCCCTGATGGTTGGTGGGCGTTAC‑CGTCCCCG‑CCGTGGGGTTAAAAC‑CCT‑CCGCCGTATTAAGCGGTTACTCTTCAGGTACCACctggctgg                                                                                                        <  2:640412/150‑1 (MQ=255)                                           ttGGGGGGTTGTTTATCAAGCTTTCCCCGGACCTTG‑TCGTACTTTTAATTGCCTTTGCTGCCTAAAGGGTGGTGGTCGGTAT‑CTTACCCG‑AGGGTG‑GTTTAAAACCCA‑CCGCCGTATTAAGCGGTTACTCTTCAGGTACCACctggctg                                                                                                         <  2:101769/149‑1 (MQ=255)                                            tgtgGGGCTGGTTATCAAGCCTTTCCCGGAACCTG‑TCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAA‑CTTATCCG‑ACGGTGCGTTTAAAA‑CCA‑CCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGgtg                                                                                                      >  2:1070707/1‑151 (MQ=255)                                                gggCTGGTTATCAAGCCTTTCCCGGAACCTG‑TCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAA‑CTTATCCG‑ACGGTGCGTTTAAAA‑CCA‑CCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTtc                                                                                                   >  1:1002278/1‑150 (MQ=255)                                                 ggCTGGTTATCAAGCCTTTCCCGAAACCTG‑TCGTATTTTTAATTGCCGTTGCTGCCTCAATGGTGGGGGTCGGTAA‑CTTATCCG‑ACGGGGCTTTTAAAC‑CCA‑CCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTtctc                                                                                                 <  2:622116/151‑1 (MQ=255)                                                    tGGTTATCAAGCCTTTCCCGGAACCTG‑TCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAA‑CTTATCCG‑ACGGTGCGTTTAAAA‑CCA‑CCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGc                                                                                              >  2:393097/1‑151 (MQ=255)                                                     gtttATCACCCCTTTCCGGGACCTTT‑TCGTACTTTTAATTGCGTTTCCTGCCTCAATGGTGGTGGTCGATAT‑TTTACCCG‑ACGTTGTGTTTAAAA‑CCA‑CCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCg                                                                                             <  1:586996/149‑1 (MQ=255)                                                      gTTATCAAGCCTTTCCCGGAACCTG‑TCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAA‑CTTATCCG‑ACGGTGCGTTTAAAA‑CCA‑CCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGt                                                                                            >  2:349162/1‑151 (MQ=255)                                                       ttATCAAGCTTTTCCCGGAACCTT‑TCTTACTTTTAATTGCCGTTGCTGCCTCAAGGGTGGTGGGCGGTAA‑TTTACCCG‑GGGGGGCGTTAAAAC‑CCA‑CCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGc                                                                                              <  1:1085835/148‑1 (MQ=255)                                                       ttATCAAGCCTTTCCCGGAACCTG‑TCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAA‑CTTATCCG‑ACGGTGCGTTTAAAA‑CCA‑CCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGtt                                                                                           >  2:1065150/1‑151 (MQ=255)                                                        tATCAAGCCTTTCCCGGAACCTG‑TCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAA‑CTTATCCG‑ACGGTGCGTTTAAAA‑CCA‑CCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGttt                                                                                          >  1:390813/1‑151 (MQ=255)                                                         aTCACGCCTTTCCCGGAACCTG‑TCGTACTGTTAATTGCCGTTGCTGCATCAATGGTGGCGGTCGGAAA‑CTTATCCG‑ACGGTGCGTTTAAAA‑CCA‑CCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGt                                                                                            >  1:446268/1‑148 (MQ=255)                                                            aaccttttCCGGGAACCTG‑CCTTTTTTTTATTTGCCGTGGCGGCCCAATTGGTGGTGGTGGGTAT‑TTTACCCG‑TCGGTGCTTTTAAAC‑CCA‑CCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTa                                                                                         <  2:479237/143‑1 (MQ=255)                                                              cccTTCCCGGGAGCCTG‑TCGTCCTGTTAATTGCCGTTGGTCCCCCAAGGGGGGGGGTCGGTAT‑TTTATCCG‑ACGGTGCTTTAAACC‑CCA‑CCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTa                                                                                    <  2:576499/150‑1 (MQ=255)                                                                 ttcccgggcaCCTG‑TCTTATTTTTATTTGCCGTTGTTCCCGCAAGGGGGGTGGTCGGTAA‑CTTATCCG‑GGGGTGCTTTAAAAC‑CCA‑CCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAg                                                                                  <  1:428067/141‑1 (MQ=255)                                                                    cccGGAACCTG‑TCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAA‑CTTATCCG‑ACGGTGCGTTTAAAA‑CCA‑CCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGcgcc                                                                              >  2:936629/1‑151 (MQ=255)                                                                        gAACCTG‑TCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAA‑CTTATCCG‑ACGGTGCGTTTAAAA‑CCA‑CCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCAt                                                                          >  2:840123/1‑151 (MQ=255)                                                                             tG‑TCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAA‑CTTATCCG‑ACGGTGCGTTTAAAA‑CCA‑CCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTtgtg                                                                     >  2:525034/1‑151 (MQ=255)                                                                              g‑TCGTACTGTTAATTGCCGTTGCTCCCTCATTGGTGGTGGTCGGTAA‑CTTACCCG‑AGGGTGCGTTAAACC‑CCA‑CCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTtgtg                                                                     <  1:771743/150‑1 (MQ=255)                                                                                tcttccttttATTTCCCGTTGCTCCCTAATTGGTGGTGGTC‑GTTA‑CTTTTCCGAAGGGTGCTTTAAAAC‑CCA‑CCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCAtt                                                                         <  1:260716/137‑1 (MQ=255)                                                                                 cGTACTTTTAATTGCCGTTGCCGCCTCAAGGGTGGTGGTGGTTAA‑CTTACCCG‑CGGGTGCTTTAAACA‑CCA‑CCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTtgtg                                                                     <  2:508076/148‑1 (MQ=255)                                                                                  gTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAA‑CTTATCCG‑ACGGTGCGTTTAAAA‑CCA‑CCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGacc                                                                  >  1:797586/1‑150 (MQ=255)                                                                                      tGTTAATTGCCGTTGCTGCCTAAATGGTGGTGGTCGGTAA‑CTTATCCG‑ACGGTGCGTTAAAAC‑CCA‑CCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGaccac                                                                <  2:928361/148‑1 (MQ=255)                                                                                         tAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAA‑CTTATCCG‑ACGGTGCGTTTAAAA‑CCA‑CCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGttt                                                          >  2:502008/1‑151 (MQ=255)                                                                                          aattgcgttgccttccTCAGTGGGGGGGGTCGTTAT‑CTTCTCCG‑ACGGTTATTTTAAAA‑CCT‑CCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGtt                                                           <  2:1099782/135‑1 (MQ=255)                                                                                              gCCGTTGCTGCCTAAATGGTGGTGGTCGGTAT‑CTTATCCG‑ACGGTGCGTTAAAAA‑CCA‑CCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAgg                                                       <  2:882077/149‑1 (MQ=255)                                                                                                cGTTGCTGCCTCAATGGTGGTGGTCGGTAA‑CTTATCCG‑ACGGTGCGTTTAAAA‑CCA‑CCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGG‑GTTCTCGGCGTTTACCTTAAGCACCGCATTTGTGACCACCGGGTTAGGGaaa                                                   >  2:685612/1‑150 (MQ=255)                                                                                                  ttGCTGCCTCAATGGTGGTGGTCGGTAA‑CTTATCCG‑ACGGTGCGTTTAAAA‑CCA‑CCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACg                                                 >  2:855634/1‑151 (MQ=255)                                                                                                  ttGCTGCCTAATGGGTGGTGGTCGTTAA‑CTTATCCG‑ACGGTGTGTTTAAAC‑CCA‑CCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTa                                                         <  1:591155/143‑1 (MQ=255)                                                                                                   tgctgcCTCAATGGTGGTGGTCGGTAA‑CTTATCCG‑ACGGTGCGTTTAAAA‑CCA‑CCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGCAAACGt                                                >  1:1068535/1‑151 (MQ=255)                                                                                                      tgcCTTAGTGGTGGGGGTGGTTAA‑TTTATCCG‑GCGGTG‑GTTTAAAACCCA‑CCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAAc                                                  <  1:1017854/146‑1 (MQ=255)                                                                                                      tgcCTAAATGGTGGGGTTTAGTAT‑TTTATCCG‑ACGGTTCTTTTAAAC‑CCT‑CCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTaa                                                    <  1:438937/144‑1 (MQ=255)                                                                                                        ccTCAGTGGTGGTGGT‑GGTAACCTTACCCG‑GCGGTGTTTTAAAAA‑CCG‑CCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGt                                                <  1:979078/146‑1 (MQ=255)                                                                                                        ccTAAAGGGTGGGGGTCGTTAA‑CTTAGCCG‑AGGGTGCGTTAAAAC‑CCA‑CCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTAtt                                             <  2:252574/149‑1 (MQ=255)                                                                                                                gtggtggtCGGTAA‑CTTATCCG‑ACGGTGCGTTTAAAA‑CCA‑CCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATctgc                                   >  2:937259/1‑151 (MQ=255)                                                                                                                  ggtgttCGGTAT‑CTTACCCG‑CGGGTGCGTTAAAAA‑CCT‑CCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATctgct                                  <  2:1002327/150‑1 (MQ=255)                                                                                                                     ggggggcaa‑CTTAGCGG‑GCGGTGCTTTAAAAG‑CCC‑CCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGAt                               <  2:306352/143‑1 (MQ=255)                                                                                                                      ggtgggtacttttgccG‑ACGTTGCGTTTAAAC‑CCA‑CCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGAtt                              <  2:871686/137‑1 (MQ=255)                                                                                                                         ggTAA‑CTTATCCG‑ACGGTGCGTTTAAAA‑CCA‑CCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGt                           >  2:717974/1‑150 (MQ=255)                                                                                                                           tAA‑GTTATCCA‑AGGGCGC‑TTTAAAA‑CCAGCCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTaa                         <  1:390338/150‑1 (MQ=255)                                                                                                                             actttccccG‑AGGGTGCTTTAAAAC‑CCA‑CCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTaa                         <  1:318768/142‑1 (MQ=255)                                                                                                                                ttATCCG‑CCGGTGCGTTTAAAC‑CCA‑CCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAAt                      <  2:220711/149‑1 (MQ=255)                                                                                                                                ttATCCG‑ACGGTGCTTTAAAAC‑CCC‑CCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATc                     <  2:509412/150‑1 (MQ=255)                                                                                                                                 tacccg‑ccgGGGCTTTAAAAC‑CCA‑CCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCg                    <  1:970518/147‑1 (MQ=255)                                                                                                                                  aTCCG‑ACGGTGCGTTTAAAA‑CCA‑CCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGt                  >  2:42112/1‑151 (MQ=255)                                                                                                                                   cccG‑ACGGTGCGTTAAAAC‑CCA‑CCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTa                 <  1:33371/150‑1 (MQ=255)                                                                                                                                    ccG‑AGGTTGCTTTAAAAA‑CCA‑CCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTaa                <  1:1169255/151‑1 (MQ=255)                                                                                                                                    ccG‑ACGGTGCGTTAAAAC‑CCA‑CCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCgg                   <  2:722308/148‑1 (MQ=255)                                                                                                                                     cG‑ACGGTGCGTTTAAAA‑CCA‑CCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGt                  >  2:702593/1‑148 (MQ=255)                                                                                                                                      g‑ACGTTGTGTATAAAC‑CCA‑CCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAac               <  2:770726/150‑1 (MQ=255)                                                                                                                                      g‑ACGGTGCGTTTAAAA‑CCA‑CCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAac               >  2:87698/1‑150 (MQ=255)                                                                                                                                          ggTGCGTTAAAAG‑CCG‑CCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATc                     <  1:72366/141‑1 (MQ=255)                                                                                                                                           gggCGTTTAAAC‑CCA‑CCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAAcaccac          <  1:373693/149‑1 (MQ=255)                                                                                                                                                tttAAAC‑CCA‑CCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACg         <  1:197739/147‑1 (MQ=255)                                                                                                                                                tttAAAC‑CCA‑CCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGctg      <  1:447556/150‑1 (MQ=255)                                                                                                                                                  tAAAA‑CCA‑CCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGctggg    >  1:109116/1‑150 (MQ=255)                                                                                                                                                   aaac‑CCA‑CCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGctgggtc  <  2:1194583/151‑1 (MQ=255)                                                                                                                                                    aaa‑CCA‑CCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGctgggtc  >  2:1087879/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                      GCGAACTGGCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACCTG‑TCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAA‑CTTATCCG‑ACGGTGCGTTTAAAA‑CCA‑CCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGCTGGGTC  >  NC_000913/801972‑802265

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