Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	786,855	0	A	C	41.0%	87.0 / 10.2	62	G247G (GGT→GGG)	gpmA	phosphoglyceromutase 1
Reads supporting (aligned to +/- strand):  ref base A (18/18);  new base C (0/25);  total (18/44)
Fisher's exact test for biased strand distribution p-value = 7.93e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAA‑GCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGG‑TTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCG  >  NC_000913/786705‑786999                                                                                                                                                        |                                                                                                                                                  ccGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAA‑GCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTa                                                                                                                                                   <  2:848177/151‑1 (MQ=255)  cGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAA‑GCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTAc                                                                                                                                                  <  2:393309/151‑1 (MQ=255)     cAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAA‑GCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCt                                                                                                                                               <  2:466564/151‑1 (MQ=255)        agcggCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAA‑GCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTAcccgg                                                                                                                                              <  1:430004/146‑3 (MQ=255)           ggCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAA‑GCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGG‑ttt                                                                                                                                         <  2:557664/150‑1 (MQ=255)           ggCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAA‑GCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGG‑ttt                                                                                                                                         <  1:677155/150‑1 (MQ=255)                aTCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAA‑GCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGG‑TTTGCAAcc                                                                                                                                   >  1:790544/1‑151 (MQ=255)                 tCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAA‑GCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGG‑TTTGCAAccg                                                                                                                                  >  2:906443/1‑151 (MQ=255)                   cGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAA‑GCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGG‑TTTGCAAccgcc                                                                                                                                <  1:668719/151‑1 (MQ=255)                    gATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAA‑GCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGG‑TTTGCAAccgcc                                                                                                                                <  1:251792/150‑1 (MQ=255)                        aaTCCCCACAGGCCTGATAAGCCGCGCGGCGTCCCATCAGGCAATGTCCTCTATTGTTAGCAACAAAAAA‑CCCGACTCCCTTGCGGCCGTTTTCCTCATTTAAAACAAATGACTTTTACTTCACTTTCCCCTGG‑TTTGCAACCGCCgctg                                                                                                                            <  1:48301/150‑1 (MQ=255)                           caccacAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAA‑GCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGG‑TTTGCAACCGCCGCTGCttt                                                                                                                        >  2:29801/1‑151 (MQ=255)                                 aGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAA‑GCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGG‑TTTGCAACCGCCGCTGCTTTCGCTGc                                                                                                                  >  2:450113/1‑151 (MQ=255)                                  ggccgaaTAAGCCGCGCAGCGTCCCATCAGGCAATGTCCTCCATTTTTAGCAACAAAAAA‑GCCGCCTCACTTGCGGTCGCTTTTCTCATTTTAAACGAATGACGTTTACTTCG‑TTTTCCCTGGTTTTGCAACCGCCGCTGCTTTCGCTg                                                                                                                   <  1:310302/143‑1 (MQ=255)                                    ccTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAA‑GCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGG‑TTTGCAACCGCCGCTGCTTTCGCTGCGAt                                                                                                               >  1:874605/1‑151 (MQ=255)                                      tGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAA‑GCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGG‑TTTGCAACCGCCGCTGCTTTCGCTGCGAtct                                                                                                             >  2:28035/1‑151 (MQ=255)                                        atatGTCGCGCAGCTTCGCTTCAGGCAGTGTGCTCCTTTGTTACCAACAAAAAC‑GCCGCCTCGCTGGCAGTGGGTTTTCTCTTTTTAAACAATTGACTTTTACTTCCTTTAACCCTGG‑TTTGCAACCGCCGCTGCTTTCGCTGCGAtctc                                                                                                            <  2:943010/150‑1 (MQ=255)                                          aaGTGGCGCACCGTCTCATCAGACAATGTGCTCCTTTGTTAGCAACAAAAAA‑CCCGTCTCATTTGCATTCGTTTTTCTCATTTAAAACAAATGGCGTTTACTTCTCTTTCCCCTGT‑TTGGCAACCGCCGCTGCTTTCGCTGCGAt                                                                                                               <  1:381169/145‑1 (MQ=255)                                          aaGTCGCGCGGCTTCGCATAGGGCATTGTGCTCCATTGTAAGCAAAAAAAAA‑GCCGACTCACTTCCAGTCGGCTTACTCATTTTAAACGAATGCCGTTTCCTTCGCTTTACCCTGT‑TTTGAAACCGCCGCTGCTTTCGCTGCGAtct                                                                                                             <  1:674599/147‑1 (MQ=255)                                          aaGTCGCGCAGCGTGGAATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAA‑GCCGACTCCCTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTCCTTCGCTTTCCCCTGG‑TTTGAAACCGCCGCTGCTTTCGCTGCGAtctc                                                                                                            <  2:435829/148‑1 (MQ=255)                                          aaGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAA‑GCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACTTTTACTTCGCTTTACCCTGG‑TTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGt                                                                                                          <  2:45567/150‑1 (MQ=255)                                                   aGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAA‑GCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGG‑TTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTAcc                                                                                                >  2:325114/1‑151 (MQ=255)                                                        cGCATCGGGCAATGTGCTCCATTTTTAGCAAAAAAAAA‑GCCGACTCCCTTGCAGTCGTTTTTCTCATTTTAAACGAATGCCTTTTCCTTCGCTTTACCCTGT‑TTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTAcc                                                                                                <  1:899373/146‑1 (MQ=255)                                                          cATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAA‑GCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGG‑TTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGata                                                                                          >  1:219089/1‑150 (MQ=255)                                                             cAGGCAATGTCCTCCATTGTTAGCAAAAAAAAA‑CCCGCCTCACTTGCGGTCGGTTTTCTCATTTTAAACGAATGACGTTTCCTTCGCTTTCCCCTGT‑TTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGataata                                                                                       <  2:425729/150‑1 (MQ=255)                                                               ggCAATATCCTCAATTGTTATAAAAAAAAAA‑GCAGACTCCTCTGATGTCTGCTTCTTAATTTTAAACAAATTAGTTTCCCTTCGTTTTCCCCGGG‑TAGGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGctaata                                                                                       <  2:747701/148‑1 (MQ=255)                                                                        cTCCATTGTTAGCAACAAAAAA‑GCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGG‑TTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCgg                                                                           >  1:713501/1‑151 (MQ=255)                                                                             ttGTTAGAAAAAAAAAA‑GCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGG‑TTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGttt                                                                        <  2:192420/149‑1 (MQ=255)                                                                              tttttagaacaaAAAAGGCCGCCTCACTTGCAGTCGGTTTTCTCATTTAAAACAAATGACGTTTTCTTCGCTTTCCCCTGG‑TTGAAATCCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCgg                                                                           <  1:155171/139‑1 (MQ=255)                                                                                  aGCAACAAAAAA‑GCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGG‑TTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATtc                                                                 >  1:102827/1‑151 (MQ=255)                                                                                  aGCAAAAAAAAG‑CCCGACTCCCTTGCAGTGGGTTTTCTCTTTTTAAACGAATGACTTTTACTTCGCTTTCCCCTGG‑TTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAAtt                                                                  <  2:489012/150‑1 (MQ=255)                                                                                   gCAACAAAAAA‑GCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGG‑TTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATtct                                                                >  2:248707/1‑151 (MQ=255)                                                                                      acaaAAAA‑GCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGG‑TTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCtcgt                                                             >  2:391367/1‑151 (MQ=255)                                                                                      aaaaaaaa‑GCCGCCTCCCTGGCAGTGGGCTTTCTCATTTAAAAGGATTGACGTTTCCTCCGCTTTACCCTGT‑TTGGCATCCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCtcgt                                                             <  1:427169/149‑1 (MQ=255)                                                                                       caaAAAA‑GCCGACTCACGTGCAGAGGCTTTTCTCATTTTAAACGAATGGCTTTTACTTCTCTTTCCCATGT‑TTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTTCCCAGATAATAGCGTTTCAGCGGTTTGAAATTCtcgt                                                             <  2:143580/150‑1 (MQ=255)                                                                                        aaaaaa‑GCCGCCTCCCTTGCAGTCGGTTTTCTCATTTTAACCGAATGACTTTTCCTTCGCTTTCCCCTGT‑TTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCtcgt                                                             <  1:333607/149‑1 (MQ=255)                                                                                                ccGACTCACTTGCAGTCGTCTTTCTCATTTTAAACGAAGGACGTTTACTTCGCTTTCCCCGGG‑TTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCAt                                                    <  1:328058/151‑1 (MQ=255)                                                                                                    cctaCTTGCATGTGGTTTTCTATTTTTAAACGAGGGTCTTTTTCTTCGATTACCCCTGG‑TTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCa                                                     <  2:362464/143‑1 (MQ=255)                                                                                                         ttGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGG‑TTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCg                                           >  1:789588/1‑151 (MQ=255)                                                                                                              gtcgttttcctcATTTTAAACAAATAGCGTTTCTTTCGCTTTACCTTGG‑TTGGAAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCg                                           <  1:549347/137‑1 (MQ=255)                                                                                                                 ggttttCTCATTTTAAACGAATGACGTTTACTCCGCTTTCCCCTGT‑TTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGcc                                    <  2:514809/147‑1 (MQ=255)                                                                                                                 ggtttcctaaTTTAAACCGATGGACGTTTTCTCCGTTTTACCTTGT‑TTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCa                                   <  1:44098/142‑1 (MQ=255)                                                                                                                    tttCTCATTTTAAACGAATGACGTTTACTTCGCTTTCCCCTGG‑TTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGcc                                    <  2:417051/147‑1 (MQ=255)                                                                                                                     ttCTCTTTTTAAACGAATGATGTTTCCTCCGCTTTCCCCTTG‑TTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTc                                <  1:112644/150‑1 (MQ=255)                                                                                                                     ttCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGG‑TTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGt                                 >  2:566752/1‑149 (MQ=255)                                                                                                                       ctcATTTTAAACGAATGACGTTTACTTCGCTTTCCCCTGG‑TTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCg                               <  2:585347/149‑1 (MQ=255)                                                                                                                           ttttAAACGATTGACTTTTACTCCGCTTTCCCTTGG‑TTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCg                               <  1:395247/145‑1 (MQ=255)                                                                                                                           ttttAAACGAATGCCGTTTCCTTCGCTTCACCCTGG‑TTTCCACCCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCgg                              <  1:505846/146‑1 (MQ=255)                                                                                                                           ttttAAACGAATGACGTTTACTTCGCTTTACCCTGG‑TTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGatat                         >  1:436557/1‑151 (MQ=255)                                                                                                                            tttAAACGAAGGACTTTTATTTCGCTTTCCCCTGT‑TTTGCACCCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGa                            <  1:905255/147‑1 (MQ=255)                                                                                                                             ttAAACGAATGCCGTTTACTTCGCTTTCCCCTGG‑TTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATAtt                        <  1:749452/150‑1 (MQ=255)                                                                                                                             ttAAACAAACGCCGTTTCTTCCTCTTTCCCCTGG‑CTTCAAACCGCCGCTGCTTTCGCTGCGTTCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGata                          <  1:413263/148‑1 (MQ=255)                                                                                                                              tAAACGAATGACGTTTCCTTCGTTTTCCCCTGG‑TTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTa                       <  1:4180/150‑1 (MQ=255)                                                                                                                              tAAACGAATGACGTTTACTTCGCTTTACCCTGG‑TTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTa                       >  2:137414/1‑150 (MQ=255)                                                                                                                                    aaTGCCGTTTCCTTCGCTTTCCCCTGT‑TTTACAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCa                 <  1:248667/150‑1 (MQ=255)                                                                                                                                    aaTGACGTTTCCTTCCCTTTCCCCTGG‑TTTGAAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCa                 <  1:789065/150‑1 (MQ=255)                                                                                                                                           tttGCTCCGTTTTCCCGTGG‑TTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAAtct          <  2:108254/150‑1 (MQ=255)                                                                                                                                               cTTCGCTTTACCCTGG‑TTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATctct        <  1:676823/148‑1 (MQ=255)                                                                                                                                                ttCCCTTTCCCCTGG‑TTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATctct        <  2:43514/147‑1 (MQ=255)                                                                                                                                                 tCGCTTTACCCTGG‑TTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCttcttc   >  1:109127/1‑151 (MQ=255)                                                                                                                                                  cGCTTTACCCTGG‑TTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCg  >  1:371288/1‑151 (MQ=255)                                                                                                                                                     tttCCCCTGT‑TTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCg  <  1:324643/148‑1 (MQ=255)                                                                                                                                                     tttACCCTGT‑TTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTcttctt    <  2:865892/146‑1 (MQ=255)                                                                                                                                                        |                                                                                                                                                  CCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAA‑GCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGG‑TTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCG  >  NC_000913/786705‑786999

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