Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	802,118	0	A	C	38.2%	77.5 / 12.9	56	T78P (ACC→CCC)	ybhI	putative DASS family tricarboxylate or dicarboxylate transporter
Reads supporting (aligned to +/- strand):  ref base A (24/10);  new base C (0/21);  total (25/31)
Fisher's exact test for biased strand distribution p-value = 7.66e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAAGCGAACTGGCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTT‑CCCGGAACCT‑GTCGTA‑CTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTT‑AAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGCTGGGTCTGG  >  NC_000913/801969‑802268                                                                                                                                                          |                                                                                                                                                       tAAGCGAACTGGCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTT‑CCCGGAACCT‑GTCGTA‑CTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTT‑aaaa                                                                                                                                                        <  2:552647/150‑1 (MQ=255)    gCGAACTGGCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTT‑CCCGGAACCT‑GTCGTA‑CTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTT‑AAAaccac                                                                                                                                                    <  1:153977/151‑1 (MQ=255)      gAACTGGCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTT‑CCCGGAACCT‑GTCGTA‑CTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTT‑AAAACCAccg                                                                                                                                                  <  2:236829/151‑1 (MQ=255)       aaCTGGCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTT‑CCCGGAACCT‑GTCGTA‑CTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTT‑AAAACCAccgc                                                                                                                                                 <  1:149592/151‑1 (MQ=255)       aaCTGGCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTT‑CCCGGAACCT‑GTCGTA‑CTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTT‑AAAACCAccgc                                                                                                                                                 >  2:685857/1‑151 (MQ=255)       aaCTGGCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTT‑CCCGGAACCT‑GTCGTA‑CTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTT‑AAAACCAccgc                                                                                                                                                 >  1:300346/1‑151 (MQ=255)         cTGGCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTT‑CCCGGAACCT‑GTCGTA‑CTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTT‑AAAACCAccgccg                                                                                                                                               >  2:269602/1‑151 (MQ=255)               tGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTT‑CCCGGAACCT‑GTCGTA‑CTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTT‑AAAACCACCGCCGTATTa                                                                                                                                          >  1:490687/1‑150 (MQ=255)                  gTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTT‑CCCGGAACCT‑GTCGTA‑CTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTT‑AAAACCACCGCCGTATTAAGCg                                                                                                                                      >  1:973459/1‑151 (MQ=255)                          tGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTT‑CCCGGAACCT‑GTCGTA‑CTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTT‑AAAACCACCGCCGTATTAAGCGGTTActct                                                                                                                              >  1:372380/1‑151 (MQ=255)                           ggTTTTTCCCTGGCGGCCATTGTGGGGTTGGTTATCAAGCTTTC‑CCCGGAACCT‑GTCGTA‑CTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTACCCGCCGGTGCGTTT‑AAACCCACCGCCGTATTAAGCGGTTActct                                                                                                                              <  1:404933/150‑1 (MQ=255)                             tttttACCTGGCGCCCTTTGTGGGGCTGTTTATCAAGCCTTC‑CCCGAACCCT‑GTTTTA‑TTTTTAATTGCCGTTCCTGCCTCAGCGGTGTTTGTCGGTAATTTATCCGACGGTGTTTTA‑AAAACCACCGCCGTATTAAGCGGTTACTCTTc                                                                                                                            <  2:535536/148‑1 (MQ=255)                             tATTTACCTGGCGGACATTGTGGGGCTGGTTATCAAGCCTTT‑CCCGGAACCT‑GTCGTA‑CTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTAGGTAACTTATCCGACGGTGCGTTT‑AAAACCACCGCCGTATTAAGCGGTTACTCTTc                                                                                                                            >  2:1015171/1‑150 (MQ=255)                               tttccgtGGCGGCCTTTGTGGGGCTGGTTATCAGCCCTTT‑CCCGGAACCT‑GTCGTC‑CTTTTAATTGCCGTTGCTCCCTCAATGGTGTTGTTCGGTAATTTACCCGACGGGGGTTTT‑AAACCCACCGCCGTATTAAGCGGTTActc                                                                                                                               <  2:446861/139‑1 (MQ=255)                                ttcccTGGCGCCCTTGGGGGGGCTGGTTATCAGGCCTTT‑CCCGGACCCT‑GTCGTA‑CTTTAATTTGCCTTTGCGACCTCAAGGGTGGTGGTCGGTATCTTATCCGACGGTGTTTTA‑AAACCCACCGCCGTATTAAGCGGTTACTCTTCAgg                                                                                                                         <  1:194516/147‑1 (MQ=255)                                   ccgggcggcCTTTGTGGGGCTTGTTATCAAGCCTTT‑CCCGAAACCT‑GTCGTA‑CTTTTATTTGCCGTTGCTGCCTCAAGGGTGGTGGTCGGTAATTTATCCGACGGTGGTTTG‑AAGCCCACCGCCGTATTAAGCGGTTACTCTTCAgg                                                                                                                         <  2:185398/144‑1 (MQ=255)                                      ggcgtcCGTGGTGGGTCTGGTTATC‑AGCCTTT‑CCCCGAACCC‑GGTGTACCTGTTATTTGCCTTTGCTCCCCCATGGGTGGTGGTCGTTACTTTCCCCGCCGGTGGTTTT‑AAACCCGCCGCCGTATTAAGCGGTTACTCTTCAGGTaccacc                                                                                                                  <  2:866683/151‑1 (MQ=255)                                      ggcggcCATTGTGGGGCTGGTTATCAAGCCTTT‑CCCGGAACCT‑GTCGTA‑CTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTGACTTATCCGACGGTGCGTTT‑AAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTaccacc                                                                                                                  >  2:272757/1‑151 (MQ=255)                                      ggcggcAATTGTGGGGGTGGTTATCATTCTTCT‑CCC‑GAACCTGGCCTTT‑TTTTTAATGGCGTTTGCGGCCGAAATGGTGGTGTTGGGAACTTTACCTGACGGTCCTTTT‑AAAACCTCCGCCGTATTAAGCGGTTACTCTTCAGGTacca                                                                                                                    <  2:90380/149‑1 (MQ=255)                                       gcggcCATTGTGGGGTTGGTTATCAAGCCTTC‑CCCGGACCCT‑GTCGTA‑CTTTTAATTGCCGTTGCTCCCTCAATGGTGGTGGTCGGTAATTTATCCGACGGTGCTTTA‑AAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACct                                                                                                                 <  1:652544/151‑1 (MQ=255)                                         ggcCTTTGTGGGGCTGGTTATCAAGCCTTT‑CCCGGACCCT‑GTCGTA‑CTGTTAATTGCCGTTGCTGCCTAAAGGGTGGTGGTGGGTAACTTACCCGACGGTGCTTTT‑AAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACctg                                                                                                                <  1:68544/150‑1 (MQ=255)                                         ggcCATTGTGGGGCTGGTTATCAAGCCTTT‑CCCGGAACCT‑GTCGTA‑CTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTT‑AAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACctgg                                                                                                               >  2:478525/1‑151 (MQ=255)                                         ggcCATTGTGGGGCTGGGTATCAAGCCTTT‑CCCGGAACCT‑GTCGTA‑CTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTT‑AAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACctgg                                                                                                               >  2:321579/1‑151 (MQ=255)                                          gcCATTGTGGGGCTGGTTATCAAGCCTTT‑CCCGGAACCT‑GTCGTA‑CTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTT‑AAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACctggc                                                                                                              >  1:596856/1‑151 (MQ=255)                                              tggtgGGGTTGGTTAGCAAGCCTTT‑CCCGGAACCT‑TTCGTA‑CTTTTAATTGCCGTGTCTCCCTCAAGGGGGGTGGTCGGTAATTTATCCGACGGTGCGTTT‑AAAACAACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACctggct                                                                                                             <  2:36618/146‑1 (MQ=255)                                                gtgGGGCTGGTTAACAAGCCTTT‑CCCGGACCCT‑GTCGTG‑CTTTTAATTGCCGTTGCTGCCTCAAGGGTGGTGGTCGGTATCTTCTCCGCCGGTTCGTTT‑AAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGgt                                                                                                          <  2:266765/149‑1 (MQ=255)                                                 tgGGGCTGGTTATCAAGCTTTT‑CCCGGAACCT‑GTCGTA‑CTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTT‑AAAACCAACGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGtt                                                                                                       >  2:89599/1‑151 (MQ=255)                                                  ggggTTGGTTATC‑AGCCTTTCCCCGGACCTT‑GTGTTT‑CTTTTATTTGCCGTTGCTCCCCCAATGGTGGGGGTCGGTAATTTCCCCTGCGGTGCTTTT‑AAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGgtg                                                                                                         <  2:370860/148‑1 (MQ=255)                                                  ggggCTGGTTATCAAGCTTTT‑CCGGGAACCT‑GTCGTA‑TTGTTAATTGCCGTTGCTGCCTAAATGGTGGTGGTGGTTAACTTATCCGAGGGTGCTTTA‑AAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGtt                                                                                                       <  1:989990/150‑1 (MQ=255)                                                                aaccTTT‑CCCGGATCCG‑GTCCTT‑CTTTTAATTCCCGTTTCTGCCTCAATGGGTGTGGGCGGTAACTTCTCCTGCTGTGCTTTT‑AAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTa                                                                                            <  2:757391/145‑1 (MQ=255)                                                                aGCCTTT‑CCCGGAACCT‑GTCGTA‑CTGTTAATTGCCGTTGCTGCCTCAATGGGGGTGGTGGGTAACTTATCCGACGGTGCGTTT‑AAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCtt                                                                                        <  1:88343/151‑1 (MQ=255)                                                                   ctttccccGGACCCT‑GTCGTA‑CTGTTAATTGCCGTTGCTGCCTAAATGGTGGTGGTCGTTAATTTATCCGACGGTGCTTTT‑AAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTa                                                                                       <  1:146703/145‑1 (MQ=255)                                                                      t‑CCCGGAACCT‑GTCGTA‑CTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTT‑AAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAgcgc                                                                                  >  2:424339/1‑151 (MQ=255)                                                                        cccGGAACCT‑GTCGTA‑CTGTTATTTGCCTTTGCTGCCTCAATGGTGGTGGTGGGTATCTTATCCGCCGGTGCGTTT‑AAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGcgcc                                                                                 <  2:372034/151‑1 (MQ=255)                                                                        cccGGAACCT‑GTCGTA‑CTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTT‑AAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGcgcc                                                                                 >  2:325736/1‑151 (MQ=255)                                                                                cT‑GTCGTA‑CTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTT‑AAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTtgt                                                                         >  1:191100/1‑151 (MQ=255)                                                                                cT‑GTCCTA‑ATGTTAATTGCCGTTGCTGCCTCAATAGTGGTGGTCGGTAAATTATCCGATGGTGCGTTT‑AAATCCACCGCCGCATTAAGCGGTTACTCTTTAGGAACAACCTGACTGTTGTTCTCGGCGTTTACCTTAAGCGCCGCATTtg                                                                          >  2:792850/1‑150 (MQ=255)                                                                                   gTCGTA‑CTGTTAATTGCGTTGGGTGCCTCAATGGGGGTGGTCGGAAACTTATCCGACGGTGCTTTT‑AAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGa                                                                       <  1:980086/151‑1 (MQ=255)                                                                                       tA‑CTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTT‑AAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGaccac                                                                   >  1:577164/1‑151 (MQ=255)                                                                                                 tggcggTTGCGCCCTCATTGGGGGTGGT‑GGTTATTTTGCCGGAGGGGGGTTTAAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTTCCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCTCATTTGTGACCACCGGttt                                                             <  2:92197/143‑1 (MQ=255)                                                                                                  tGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTT‑AAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAg                                                           >  2:626757/1‑149 (MQ=255)                                                                                                     cGTAGCTCCCTCAAGGGGGTGGGTGGGTAACTTATCCGAGGGGGGTTTA‑AAAACCCCCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGt                                                         <  1:857279/148‑1 (MQ=255)                                                                                                           tgcCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTT‑AAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAAc                                                     <  1:499199/146‑1 (MQ=255)                                                                                                           tgcCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTT‑AAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTAtt                                                >  1:396689/1‑151 (MQ=255)                                                                                                                     gtggtggtCGGTAATTTATCCGACGTTGTGTTT‑AAAACCTCCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATct                                        <  1:809572/149‑1 (MQ=255)                                                                                                                             cGGTAACTTATCCGACGGTGCGTTT‑AAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGt                              >  1:385699/1‑151 (MQ=255)                                                                                                                              ggTAACTTATCCGACGGTGCGTTT‑AAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGt                              >  2:800704/1‑150 (MQ=255)                                                                                                                                taactaacCCGAGGGTGCCTTT‑AAACCTGCCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTa                                           <  1:214728/135‑1 (MQ=255)                                                                                                                                tAACTTATCCGACGGTGCGTTT‑AAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTaa                            >  1:986318/1‑150 (MQ=255)                                                                                                                                    ttacctttcgtgccGTTT‑AACCCCACCGCCGTATTAAGCGGTTACTCTTCAGTTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCTTTTGTGACCACCGGTTTAGCTAAACGTATTGCCGATCTGCTGATTGGTAAAATc                        <  2:894232/137‑1 (MQ=255)                                                                                                                                      acccGAGGGTGCTTTA‑AAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCg                       <  1:971419/147‑1 (MQ=255)                                                                                                                                      aTCCGACGGTGCGTTT‑AAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCgg                      <  1:761263/150‑1 (MQ=255)                                                                                                                                          ggcGGTGCTTTT‑AAACCCGCCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAaca                 <  2:418996/149‑1 (MQ=255)                                                                                                                                             ggTGCGTTT‑AAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAAcacca              <  1:275853/151‑1 (MQ=255)                                                                                                                                                 ctttt‑AAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACg            <  2:558704/147‑1 (MQ=255)                                                                                                                                                 cGTTT‑AAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGct          >  1:454284/1‑151 (MQ=255)                                                                                                                                                  gTTT‑AAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGctg         >  2:439824/1‑151 (MQ=255)                                                                                                                                                    tt‑AAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGctggg       >  1:506185/1‑151 (MQ=255)                                                                                                                                                       aaaaCCCCCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGctggg       <  1:628796/149‑1 (MQ=255)                                                                                                                                                       aaaaCCACCGCCGTATTAAGCGGATACTCTTCAGGTACCACCTGCCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATCTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGCTTGGCAAAATCGGGAACACAACGcagggtc     >  2:686141/1‑151 (MQ=255)                                                                                                                                                          accaccGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGctgggtctgg  >  2:842562/1‑151 (MQ=255)                                                                                                                                                          |                                                                                                                                                       TAAGCGAACTGGCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTT‑CCCGGAACCT‑GTCGTA‑CTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTT‑AAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGCTGGGTCTGG  >  NC_000913/801969‑802268

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