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

TAAGCGAACTGGCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCC‑TTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAAC‑TTATCCGACGGT‑GCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGCTGGG  >  NC_000913/801969‑802263                                                                                                                                                         |                                                                                                                                                   tAAGCGAACTGGCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATCGTGGGGCTGGTTATCAAGCC‑TTTCCCGGAACCTGTCGTACTGTTTGTTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAAC‑TTATCCGACGGC‑GCGTTTaaaa                                                                                                                                                    >  2:839855/1‑150 (MQ=255)    gCGAACTGGCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCC‑TTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAAC‑TTATCCGACGGT‑GCGTTTAAAa‑ccac                                                                                                                                               <  1:1035151/151‑1 (MQ=255)         cTGGCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCC‑TTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAAC‑TTATCCGACGGT‑GCGTTTAAAA‑CCAccgccg                                                                                                                                          >  2:953948/1‑151 (MQ=255)               tGGGTGTTTTTTGGTTTTTACCTGGCGGCCTTTGTGGGGCTGGTTATCAAGCT‑TTTCCCGGACCCTGTCGTACTGTTAATTGCGGTTGCTCCCTCATTGGTGGTGGTGGGTAAT‑TTATCCGACGGT‑GCGTTTAAAC‑CCACCGCCGTAtt                                                                                                                                      <  2:362573/149‑1 (MQ=255)                      tttttGGTATTTACCTGGCGGCCTTTGTGGGGCTGTTTATCAAGCC‑TTTCCCGGACCCTGTCGTATTTTTATTTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAAC‑TTACCCGAGGGT‑GCGTTTAAAC‑CCACCGCCGTATTAAGCGGtt                                                                                                                              <  2:1063317/150‑1 (MQ=255)                       ttttGGTATTTCCCTGGCGGCCATTGTGGGGCTGGTTATCAAGCC‑TTTCCCGGACCCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAAC‑TTATCCGAGGGT‑GCGTTAAAAC‑CCACCGCCGTATTAAGCGGTTa                                                                                                                             <  1:767436/150‑1 (MQ=255)                        tttGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCC‑TTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAAC‑TTATCCGACGGT‑GCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTAct                                                                                                                           >  2:721675/1‑151 (MQ=255)                             tATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCC‑TTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAAC‑TTATCCGACGGT‑GCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTc                                                                                                                       >  2:587630/1‑150 (MQ=255)                              aTTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCC‑TTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAAC‑TTATCCGACGGT‑GCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTc                                                                                                                       >  1:731407/1‑149 (MQ=255)                              aTTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCC‑TTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAAC‑TTATCCGACGGT‑GCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTc                                                                                                                       >  1:68983/1‑149 (MQ=255)                              aTTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCC‑TTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAAC‑TTATCCGACGGT‑GCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTc                                                                                                                       >  1:105724/1‑149 (MQ=255)                               tttCCCTGGCGGCCATTGTGGGGCTGGTTATCAAGCC‑TTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAAC‑TTATCCGACGGT‑GCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAg                                                                                                                     <  1:413139/150‑1 (MQ=255)                               tttCCCTGGCGCCCTTTGTGGGGGGGGTTATCAGGCC‑TTTCCCGGAACCTGTCGTACTGTTAATTGCGGTTGCTGCCTCAATGGTGGGGTTGGGTACC‑TTACCCGACGGT‑GCTTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAg                                                                                                                     <  2:1068849/150‑1 (MQ=255)                                   ccTGGCGGCCATTGTGGGGCTGGTTATCAAGCT‑TTTCCCGGAACCTGTCGTACTGTTAATTGCGGTTGCTGCCTCAATGGTGGTGGTCGGTAAC‑TTATCCGCCGGG‑GCGTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTac                                                                                                                 <  1:209085/150‑1 (MQ=255)                                         ggcCTTTGTGGGGCTGTTTATAAGGCC‑TTTCCGGGAACCTGTCGTACTGTTAATTGCCTTTGCTGCCTCAATGGTGGTGGTCGGTACC‑TTACCCGACGGT‑GCTTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACct                                                                                                            <  1:890592/149‑1 (MQ=255)                                            cATTGTGGGGCTGGTTATCAAGCC‑TTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAAC‑TTATCCGACGGT‑GCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACctggctg                                                                                                       >  2:580364/1‑151 (MQ=255)                                              ttGTGGGGCTGGTTATCAAGCT‑TTTCCCGGAACCTGTCGTACTTTTAATTGCCGTTGCTGCCTAAATGGTGGTGGTCGGTAAC‑TTATCCGACGGT‑GGGTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACctggctgg                                                                                                      <  1:477322/150‑1 (MQ=255)                                               tgtgGGGCTGGTTATCAAGCC‑TTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAAC‑TTATCCGACGGT‑GCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGgtg                                                                                                    >  1:229124/1‑151 (MQ=255)                                               tgtgGGGCTGGTTATCAAGCC‑TTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAAC‑TTATCCGACGGT‑GCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGgt                                                                                                     >  1:932368/1‑150 (MQ=255)                                                  ggggCTGGTTATCAAGCC‑TTCCCCGGAACCTGTCTTATTGTTAATGCCCGTTGCTGCCTCAATGGGGGGGTTGGGTATT‑TTATCCGACGGT‑GCTTTTAAAC‑CCCCCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGtt                                                                                                  <  1:868237/150‑1 (MQ=255)                                                     gCTGGTTATCAAGCC‑TTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAAC‑TTATCCGACGGT‑GCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCg                                                                                              >  2:845495/1‑151 (MQ=255)                                                        ggTTATCAAGCT‑TTTCCCGGAACCTGTCGTACTTTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTATC‑TTATCCGACGGT‑GCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGc                                                                                            <  1:607258/150‑1 (MQ=255)                                                        ggTTATCAAGCC‑TTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAAC‑TTATCCGACGGT‑GCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCg                                                                                           >  1:295230/1‑151 (MQ=255)                                                          ttATCAAGCT‑TTTCCGGGACCCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAAT‑TTTTCCGACGGT‑GCGTTAAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCgg                                                                                             <  1:313389/147‑1 (MQ=255)                                                           ttctaaGCCTTTTCCCGGACCTTGTGTTATTGTTAATTGCCGTTGCTCTTTTAATTGTGGTTGTTGGTAAT‑TTCCCCGACGGT‑CCGTTTAACC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTtctc                                                                                               <  2:524484/141‑1 (MQ=255)                                                                   cttttCCCGGAACCTGTCGTACTTTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAAC‑TTATCCGACGGT‑GTTTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTa                                                                                  <  1:893799/148‑1 (MQ=255)                                                                     tttCCCGGACCCTGTCGTACTGTTAATTGCCGTTGCTCCCTAAATGGTGGTGGTGGGTAAT‑TTATCCGAGGGT‑GCTTTAAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTa                                                                                  <  2:125891/148‑1 (MQ=255)                                                                      ttCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAAC‑TTATCCGACGGT‑GCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAgcg                                                                              >  2:1052664/1‑151 (MQ=255)                                                                      tccgtgaacccTGTCGTACTTTTAATTCCCGTTCCTGCATACTTGGGGGGGTTCTGTCAC‑TTACCCGACTGT‑GGGTTCAAAC‑CCTCCGCCGTATTAAGCGGTTACTCTTCAGGTGCCACCTGGCTGTTGTTCTCGGCGTTTACCTTAAgc                                                                               <  2:39153/141‑1 (MQ=255)                                                                       tCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAAC‑TTATCCGACGGT‑GCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAgcgc                                                                             >  2:947802/1‑151 (MQ=255)                                                                              aCCTGTCGTACTGTTAATTGCCGTTGCTGCCTCACTGGTGGTGGTCGGTAAC‑TTATCCGACGGT‑GCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCAttt                                                                      >  2:17812/1‑151 (MQ=255)                                                                                 tGTCGTTTTGTAAATGCCCGTTGCTGCCTCATGGGGGGTGTTCGGTATC‑TTATCGGCGGGG‑GCGTTTAAAT‑CCGCCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTAGTGTTCTCGGCGTTTACCTTAAGCGCCGCAttt                                                                      <  2:160955/148‑1 (MQ=255)                                                                                  gTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAAC‑TTATCCGACGGT‑GCGTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTtgt                                                                    <  1:542696/149‑1 (MQ=255)                                                                                       aCTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAAC‑TTATCCGACGGT‑GCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGaccac                                                              >  2:952002/1‑150 (MQ=255)                                                                                           ttATTTGCCTTTGCTGCCTCAAGGGTGGTGGTCGGTAAT‑TTATCCGAGGGT‑GCGTTAAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGt                                                          <  2:284163/150‑1 (MQ=255)                                                                                           ttAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAAC‑TTATCCGACGGT‑GCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGtt                                                         >  1:878145/1‑151 (MQ=255)                                                                                           taaaTGGCCGTGGCTGCCTCAATGGTGGGGGT‑GGTAACTTTATCCGACGGT‑GCTTTAAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCgg                                                           <  1:250302/147‑1 (MQ=255)                                                                                             aaTTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAAC‑TTATCCGACGGT‑GCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGttt                                                        >  2:69865/1‑150 (MQ=255)                                                                                             aaTTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAAC‑TTACCCGCCGGT‑GCGTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGtt                                                         <  1:21878/149‑1 (MQ=255)                                                                                                  ccgtgtctgcTTCATTGGTGTTGCTCGGTAGT‑TTGTCCGCCGGT‑GGGTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTa                                                   <  2:1055475/144‑1 (MQ=255)                                                                                                         tgcCTCAATGGTGGTGGTCGGTAAC‑TTATCCGCCGGT‑GCTTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACg                                               <  1:132676/147‑1 (MQ=255)                                                                                                           ccggaaTGTTGGTGTTCGTTAAC‑TCAT‑CGGCGGTGGCTTTTAAGC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTg                                          <  2:328551/146‑1 (MQ=255)                                                                                                               aaTGGTGGTGGTCGGTAAC‑TTATCCGACGGT‑GCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTAt                                     >  1:288510/1‑151 (MQ=255)                                                                                                               aaTGGTGGGGGTCGTTAAC‑TTATCCGACGGT‑GCGTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTAt                                     <  1:217254/151‑1 (MQ=255)                                                                                                                 gggtggtggtCGGTATC‑TTACCCGCCGTT‑GCTTTAAAAA‑CCGCCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATct                                   <  2:232726/150‑1 (MQ=255)                                                                                                                    tggtggtCGGTAAC‑TTATCCGACGGT‑GCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATctgct                                >  1:798054/1‑151 (MQ=255)                                                                                                                    tggtggtCGGTAAC‑TTATCCGACGGT‑GCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATctgc                                 >  2:257875/1‑150 (MQ=255)                                                                                                                       tggtCGGTATC‑TTTCCCGACGGT‑GCGTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGa                              <  2:315501/150‑1 (MQ=255)                                                                                                                          tCGGTAAC‑TTATCCGACGGT‑GCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTgg                          >  2:471491/1‑151 (MQ=255)                                                                                                                           cGGTAAC‑TTATCCGACGGT‑GCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGt                         >  1:201764/1‑151 (MQ=255)                                                                                                                               aaC‑TTATCCGACGTT‑GCGTTTAAAC‑CCCCCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCTGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTaaa                      <  2:148840/150‑1 (MQ=255)                                                                                                                                   ttATCCGACGGT‑GCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGc                                         >  2:419968/1‑128 (MQ=255)                                                                                                                                   ttACCCGACGGT‑GCGTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGc                                         <  1:419994/128‑1 (MQ=255)                                                                                                                                       ccGACGGT‑GCTTTAAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCg                  <  1:953997/147‑1 (MQ=255)                                                                                                                                       ccGACGGT‑GCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTa               >  1:1045874/1‑150 (MQ=255)                                                                                                                                         gACGGT‑GCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAac             >  2:778835/1‑150 (MQ=255)                                                                                                                                           gcgggggctttAAAATCCCCCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGACTGGTGTTCTCTGCGTTTGCCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAaca            <  2:437761/143‑1 (MQ=255)                                                                                                                                           atggttcttttAAAC‑CCGCCGCCGTATTAAGCGGTTACTCTTCCGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAac             <  2:903278/141‑1 (MQ=255)                                                                                                                                             ggtggttttAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGt                <  1:940609/138‑1 (MQ=255)                                                                                                                                             gT‑GCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAAcaccac        >  1:492222/1‑151 (MQ=255)                                                                                                                                             gT‑GCGTTTAAAA‑CCACCGCCGGATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAAcaccac        >  2:173036/1‑151 (MQ=255)                                                                                                                                              t‑GCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACg       >  2:263897/1‑151 (MQ=255)                                                                                                                                              gtccGTTTAAAC‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAAcacca         <  2:292139/147‑1 (MQ=255)                                                                                                                                                  gTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGctg    >  1:77905/1‑151 (MQ=255)                                                                                                                                                  gTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATCTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGctg    >  1:918044/1‑151 (MQ=255)                                                                                                                                                     tAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGctggg  >  2:186108/1‑150 (MQ=255)                                                                                                                                                        aa‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTAt                                     <  1:850047/112‑1 (MQ=255)                                                                                                                                                        aa‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTAt                                     >  2:849999/1‑112 (MQ=255)                                                                                                                                                         |                                                                                                                                                   TAAGCGAACTGGCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCC‑TTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAAC‑TTATCCGACGGT‑GCGTTTAAAA‑CCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGCTGGG  >  NC_000913/801969‑802263

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