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

AACTGGCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAG‑CCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGT‑GCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGCTGGGTCTGG  >  NC_000913/801975‑802268                                                                                                                                                  |                                                                                                                                                        aaCTGGCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAG‑CCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGT‑GCGTTTAAAA‑CCAccgc                                                                                                                                                 <  1:164937/151‑1 (MQ=255) aaCTGGCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAG‑CCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGT‑GCGTTTAAAA‑CCAccgc                                                                                                                                                 <  1:443444/151‑1 (MQ=255)     ggCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAG‑CCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGT‑GCGTTTAAAA‑CCACCGCCGt                                                                                                                                              >  1:88113/1‑150 (MQ=255)     ggCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAG‑CCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGT‑GCGTTTAAAA‑CCACCGCCGt                                                                                                                                              >  2:56481/1‑150 (MQ=255)                  tttGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAG‑CCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGT‑GCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTAct                                                                                                                                >  2:186175/1‑151 (MQ=255)                   ttGGTATTTACCTGGCGGCCATTGTGGGGCTGTTTATCAAG‑CCTTTCCCGGACCCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTACCCGACGGT‑GCGTTTAAAC‑CCACCGCCGTATTAAGCGGTTAct                                                                                                                                <  1:341856/150‑1 (MQ=255)                         tttCCCTGGCGGCCATTGGGGGGCTGGTTATCAAG‑CCTTTCCGGGACCCTGTCGTACTGTTAATTGCCGTTGCGGCCTAAATGGTGGTGGTCGTTAACTTATCCGAGGGT‑GCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCa                                                                                                                           <  1:779752/149‑1 (MQ=255)                         tttACCTGGCGGCCATTGTGGGGCTGGTTATCAAG‑CCTTTCCCGGAACCTGTCGTACTTTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGT‑GCGTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCa                                                                                                                           <  2:926757/149‑1 (MQ=255)                           tACCTGGCGGCCATTGTGGGGCTGGTTATCAAG‑CCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGT‑GCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGt                                                                                                                        >  2:293976/1‑150 (MQ=255)                            acttGGCGGCCATGGTGGGTTTGGTT‑CCAAGCCCTTTCCCGGACCCGTTCGTTTTGTTCACTCCTGTTGTTACCAAAATGTTGTTGGGC‑AAAACTTTACCACAGGTCGCGTTAAACA‑CCACCGCCGTATTAAGCGGTTACTCTTCAgc                                                                                                                         <  1:214909/145‑2 (MQ=255)                             ccgggcggcATTTGTGGGGTTGGTTTTCAAC‑CCTTTCCCGGACCCTGTCGTATTGTTAATTGCGGTTGCTGCTTCAAGGGTGGTGGGCGGTAATTTACCCGACGGT‑GCTTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTa                                                                                                                       <  2:93077/146‑1 (MQ=255)                              ctggcgcactgtgtgGGGTGGGTTCTCAAT‑CCTTTCCGGAACCCTGTCGTCCTGTTAATTGCCGTTGCTGCTACAATGGTGGGGGTCGTTAACGTATCCGCCGGT‑GCTTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGt                                                                                                                        <  2:280478/136‑1 (MQ=255)                               tGGCGGCCATTGTGGGGCTGGTTATCAAG‑CCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGT‑GCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTaccac                                                                                                                   >  2:378938/1‑151 (MQ=255)                                ggcggcCTTTGTGGGGCTGGTTATCAAG‑CCTTTCCCGGAACCTGTCTGATGGTTAATTGCCGTTGCTGCCTCAATGGTGGTGTTCGGTAACTTATCCGACGGG‑GCTTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTaccacc                                                                                                                  <  1:546302/151‑1 (MQ=255)                                   ggcCATTGTGGGGCTGGTTATCAAG‑CCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGT‑GCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACctgg                                                                                                               >  2:834920/1‑151 (MQ=255)                                    gcCTTTGTGGGGCTGGTTATAAAG‑CCTTTCCCGGACCCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGT‑GCGTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACct                                                                                                                 <  2:107700/148‑1 (MQ=255)                                      cATTGTGGGGCTGGTTATCAAG‑CCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGT‑GCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACctggctg                                                                                                            >  2:621343/1‑151 (MQ=255)                                      cATTGTGGGGCTGGTTATCAAG‑CCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGAGGTGGTCGGTAACTTATCCGACGGT‑GCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACctggctg                                                                                                            >  2:443686/1‑151 (MQ=255)                                          gtgGGGCTGGTTATCAAG‑CCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGT‑GCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGgtgt                                                                                                        >  1:853179/1‑151 (MQ=255)                                              ggCTGGTTATCAAG‑CCTTTCCCGGACCCTGTCGTACTGTTAATTGCCGTTGCTCCCTCAATGGTGGTGGTCGGTAATTTACCCGACGGT‑GCGTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTtctc                                                                                                    <  1:959314/151‑1 (MQ=255)                                                         aaG‑CCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGTTCGGTACCTTATCCGAGGGT‑GCGTTAAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTAc                                                                                           <  2:528988/149‑1 (MQ=255)                                                           accttttCCGGGACCTTGTCGTACTTTTAATTGCCGTTGCTGCCTAATTGGTGGTGGTGGTTAATTTATCCGACGGG‑GGGTTTAAAC‑CCTCCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTAc                                                                                           <  1:684554/144‑1 (MQ=255)                                                              cTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGT‑GCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTa                                                                                       >  1:934548/1‑149 (MQ=255)                                                              cTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGT‑GCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTa                                                                                       >  1:148254/1‑149 (MQ=255)                                                                  cccGGAACCTGTCGTACTTTTAATTGCCGTTGCTGCCTCAAGGGTGGTGGTCGGTAACTTATCCGACGGT‑GCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGcgcc                                                                                 <  2:618474/151‑1 (MQ=255)                                                                   ccGGACCCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAAGGGTGGTGGTCGGTACCTTACCCGACGGT‑GCTTTAAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGcgccg                                                                                <  1:277444/151‑1 (MQ=255)                                                                     ggAACGTTTCGTATTGTTAATTGCTGTGGCTCCTTAAATGTGGGTGGTCGTTAATTTTTCCTGCGGT‑GTGTTAAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGcgccgc                                                                               <  2:845965/150‑1 (MQ=255)                                                                       acccTGTCGTACTGTTATTTGCCTTTGCTGCCTCAATGGTGGTGGTGGGTAACTTATCCGCCGGT‑GCGTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCAtt                                                                            <  1:913528/149‑1 (MQ=255)                                                                          cTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGT‑GCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTtg                                                                          >  1:951762/1‑150 (MQ=255)                                                                           tGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGT‑GCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTtg                                                                          >  2:163529/1‑149 (MQ=255)                                                                            gTCGTACTTTTAATTCCCTTTGCTGCCTCAAGGGGGGTGGTCGGTAATTTATCCGACGGT‑CCGTTAAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTtgt                                                                         <  2:423193/149‑1 (MQ=255)                                                                                 aCTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGT‑GCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGaccacc                                                                  >  1:353865/1‑151 (MQ=255)                                                                                    tttAATGGCCTTTGCTCCCTCAATGGTGGTGGTCGTTATTTTACCCGACGGT‑GCGTTAAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCgg                                                                <  1:150118/149‑1 (MQ=255)                                                                                    gTTAATTGCGGTTTCTGCCTAATTGGTGGTGGTCGGTAATTTACCGGGGGGT‑GCGTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCg                                                                 <  1:434390/149‑1 (MQ=255)                                                                                     ttAATTGCCGTTGCTGCCTAAATGGGGGTGGTGGTTAATTTATCCGACGGT‑GCTTTAAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCg                                                                 <  1:623744/148‑1 (MQ=255)                                                                                            cctttGCTGCCTCAAGGGTGGGGGTCGGAAACTGACCCGGGGGT‑GCTTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTaa                                                       <  1:149621/148‑1 (MQ=255)                                                                                             cGTTGCGCCCTAAAGGGTGGTGGTCGTTAACTTCCCCGGCGGT‑GCTTTAAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTaaa                                                      <  2:373951/151‑1 (MQ=255)                                                                                              gTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGT‑GCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAAc                                                     >  1:356612/1‑151 (MQ=255)                                                                                                tgctgcTTAAATGGGGGTGTTCGTTACCTTACCCGACGG‑‑GCCTTTAAAACCCCCCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACg                                                    <  1:729458/150‑1 (MQ=255)                                                                                                 gttgcCTGAATGGTGGTGATCGTTAATTTATCCGGCGGT‑GGGTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGt                                                   <  1:152160/148‑1 (MQ=255)                                                                                                  ctgcCTCAATGGTGGTGTTCGGTAATTTACCCGACGGT‑GCTTTAAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTa                                                  <  1:99654/150‑1 (MQ=255)                                                                                                       tCAATGGTGGTGGTCGGTAACTTATCCGACGGT‑GCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCt                                            >  2:602372/1‑151 (MQ=255)                                                                                                        cAATGGTGGTGGTCGGTAACTTATCCGACGGT‑GCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCt                                            >  1:595502/1‑150 (MQ=255)                                                                                                         agtggtggtggGCTGTAACTTATCCGACGGT‑GCTTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTa                                           <  1:599168/148‑1 (MQ=255)                                                                                                         aaTGGTGGTGGTCGGTAACTTATCCGACGGT‑GCGTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCt                                            <  2:649334/149‑1 (MQ=255)                                                                                                         aaTGGTGGTGGTCGGTAACTTATCCGACGGT‑GCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCt                                            <  2:321872/149‑1 (MQ=255)                                                                                                         aaTGGTGGTGGTCGGTAACTTATCCGACGGT‑GCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTAt                                          >  1:45175/1‑151 (MQ=255)                                                                                                           tggtggtggTCGGTAACTTATCCGACGGT‑GCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATct                                        >  1:977528/1‑151 (MQ=255)                                                                                                           tggtggtggTCGGTAACTTATCCGACGGT‑GCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATct                                        >  2:257104/1‑151 (MQ=255)                                                                                                               ggtggtCGGTAACTTTTCCGCCGGT‑GCGTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATctgctg                                    <  2:157655/151‑1 (MQ=255)                                                                                                               ggtggtCGGTAACTTATCCGACGGT‑GCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATctgctg                                    >  2:966093/1‑151 (MQ=255)                                                                                                               ggtggtCGGTAACTTACCCGACGGT‑GCGTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATctgct                                     <  2:617375/150‑1 (MQ=255)                                                                                                                  ggtCGGTAACTTAGCCGACGGT‑GCTTTAAAAC‑CCCCCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGAtt                                 <  1:375568/151‑1 (MQ=255)                                                                                                                      gttAACTTATCCGAGGGT‑GCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTa                             <  2:603938/149‑1 (MQ=255)                                                                                                                       gTAACTTATCCGACGGT‑GCGTTAAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTaa                            <  2:818/151‑1 (MQ=255)                                                                                                                        tAACTTATCCGACGGT‑GCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGGAAACGTATTGCCTATCTGCTGATTGGTaa                            >  1:942910/1‑150 (MQ=255)                                                                                                                             taccggACGGT‑CCTTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCg                       <  1:253219/145‑1 (MQ=255)                                                                                                                             tATCCGACGGT‑GCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCgg                      >  2:16271/1‑151 (MQ=255)                                                                                                                                      gtggcttttAAAC‑CCTCCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAAcacca              <  1:395010/145‑1 (MQ=255)                                                                                                                                           gTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGct          <  2:33346/150‑1 (MQ=255)                                                                                                                                            tttAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGctgg        >  2:679427/1‑151 (MQ=255)                                                                                                                                              tAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGctgggt      >  1:181742/1‑151 (MQ=255)                                                                                                                                                  a‑ccaccGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGctgggtctgg  >  2:814288/1‑151 (MQ=255)                                                                                                                                                  |                                                                                                                                                        AACTGGCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAG‑CCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGT‑GCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGCTGGGTCTGG  >  NC_000913/801975‑802268

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