Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	786,855	0	A	C	42.9%	167.4 / 15.2	91	G247G (GGT→GGG)	gpmA	phosphoglyceromutase 1
Reads supporting (aligned to +/- strand):  ref base A (42/10);  new base C (0/39);  total (42/49)
Fisher's exact test for biased strand distribution p-value = 1.12e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTC  >  NC_000913/786706‑787002                                                                                                                                                      |                                                                                                                                                    cGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTAc                                                                                                                                                    >  1:305857/1‑151 (MQ=255)  gTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTAc                                                                                                                                                    <  1:165796/150‑1 (MQ=255)     agagTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTg                                                                                                                                                >  1:431515/1‑151 (MQ=255)     agagTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTg                                                                                                                                                <  2:593819/151‑1 (MQ=255)              tATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAAc                                                                                                                                       >  2:8569/1‑151 (MQ=255)               aTCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAAcc                                                                                                                                      >  2:223227/1‑151 (MQ=255)               aTCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAAcc                                                                                                                                      >  2:536571/1‑151 (MQ=255)               aTCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAAc                                                                                                                                       <  1:781255/150‑1 (MQ=255)               aTCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAAc                                                                                                                                       <  1:1779/150‑1 (MQ=255)                  cGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAAccgcc                                                                                                                                   <  1:27167/151‑1 (MQ=255)                       aaTCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCgctgc                                                                                                                              >  2:1049719/1‑151 (MQ=255)                            ccacGGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTTTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGTTTTCTCATTTTAACCGAATGACTTTTACTTCGCTTTACCCTGTTTTGCAACCGCCGCTGCTTTc                                                                                                                          <  2:670110/150‑1 (MQ=255)                             cacaGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGc                                                                                                                        >  1:675986/1‑151 (MQ=255)                              acaGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCt                                                                                                                       >  1:125983/1‑151 (MQ=255)                                 ggCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTTTTAGCAAAAAAAAAGCCGACTCATTGGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGc                                                                                                                     <  1:1051705/150‑1 (MQ=255)                                 ggCCTGAAAAGCCGCGCAGGGGCGCATCGGGCAATGTGCTCCATTGTGAGCAAAAAAAAAGCCGACCCACTTGCAGTCGGCTTTCTCATTTTAAGCGAATGATGTTTACTTCGTTTTACCCGGTTTTGCAACCGCCGCTGCTTTCGCTg                                                                                                                      <  2:452943/149‑1 (MQ=255)                                     tGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGAtct                                                                                                                >  2:820254/1‑151 (MQ=255)                                       aTAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGCCTCACTTGCAGTCGGCTTTCTCATTTTAACCGAATGACGTTTACTTCGCTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGAtctc                                                                                                               <  2:1197448/150‑1 (MQ=255)                                       aTAAGTCGCGCAGCGTCGCAGCGGGAAATGTCCTCCATTTTTAGCAAAAAAAAGGCCGTCTCGCTTGCAGTCGGTTTCCTCATTTAAAACAAATGACGTTTACTCCGCTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGAtctc                                                                                                               <  1:135502/150‑1 (MQ=255)                                       aTAAGCCGCGCAGCGTCGCATCGGGCAATGTGCTCCATTGTTAGAAACAAAAAACCCGACTCACTTGCGGTCGGTTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGAtctc                                                                                                               <  2:481546/150‑1 (MQ=255)                                         aggccgcgcgGCGCCGCATCAGGCATTGTGCTCCATTTTTAGCAACAAAAAAGCCGACTCCCTTGCAGTCGGCTTTCTCATTTTAAACGAATGACTTTTACTTCGCTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTc                                                                                                            <  2:327982/147‑1 (MQ=255)                                            tCGCGCAGCGTCGCATCGGGAATTTTGCTCCATTTTAACAAACAAAAAAGCCGCCTCCCTTGCAGTCGCTTTTCTCATTTTAAACAAATGCCGTTTCCTTCGCTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCa                                                                                                           <  2:1092808/149‑1 (MQ=255)                                              gcgcGGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGTTTTCTCTTTTTAAACGAATGACGTTTACTTCGCTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCa                                                                                                        <  1:1135840/150‑1 (MQ=255)                                              gcgcAGCGCCGCATCAGGAAAGGTGCTCCATTTTTAGAAAAAAAAAAGCCGACTCACTTGCAGTCGGTTTTCTCATTTTAAACGAATGCCTTTTACTTCGCTTTACCCTGGTTTCCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCa                                                                                                        <  1:52736/150‑1 (MQ=255)                                                gcAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCAtt                                                                                                      >  1:917050/1‑150 (MQ=255)                                                      tCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGTTTTCTCATTTTAAACGAATGACGTTTCCTTCGCTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTAcc                                                                                                   <  1:92783/147‑1 (MQ=255)                                                          aTCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGataat                                                                                           >  1:555742/1‑151 (MQ=255)                                                             aGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGc                                                                                        >  2:213890/1‑151 (MQ=255)                                                              ggCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCCCTTGCAGTCGGTTTTCTCATTTAAAACGAATGACGTTTACTTCGCTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGc                                                                                        <  1:559472/150‑1 (MQ=255)                                                               gCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGt                                                                                      >  2:772652/1‑151 (MQ=255)                                                                cAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGtt                                                                                     >  1:581469/1‑151 (MQ=255)                                                                  aTGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTc                                                                                   >  1:338358/1‑151 (MQ=255)                                                                    gtcccccATTGTTAGCAAAAAAAAAGCCGCCTCACTTGCAGTCGTCTTTCTCATTTTAAACGAATGCCGTTTACTTCGCTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTc                                                                                   <  2:1173938/144‑1 (MQ=255)                                                                       cTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCgg                                                                              >  2:281345/1‑151 (MQ=255)                                                                       cTCCATTGTTAACAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAACCGAATGACTTTTACTTCGCTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGc                                                                                <  2:42542/149‑1 (MQ=255)                                                                        tCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGt                                                                             >  2:837362/1‑151 (MQ=255)                                                                         ccATTGTAAGCAACAAAAAAGCCGACTCCTTGGCGGTCGGCTTTCTCTTTTTAAACGATTGCCGTTTACTTCCCTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCgg                                                                              <  2:374439/149‑1 (MQ=255)                                                                           attttaaccaaaaaAAAAGCCGCCTCACTTGCAGTCGCCTTTCTTATTTAAAACAAATGACGTTTATTTCGCTTTACCCTGTTTTGCATCCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGttt                                                                           <  1:802271/142‑1 (MQ=255)                                                                            ttGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGt                                                                             <  1:37774/147‑1 (MQ=255)                                                                            ttGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTg                                                                          >  2:975476/1‑150 (MQ=255)                                                                                    aacaaAAAACCCGACTCACTTGCAGTCGGTTTTCTCATTTTAAACGAATGACTTTTACTTCGCTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATtct                                                                   <  1:496882/149‑1 (MQ=255)                                                                                       aaaaaaGCCGACTCACTTGCAGTCGTCTTTCTCATTTAAAACGAATGACGTTTACTTCTCTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCtcgtc                                                               <  2:1174476/150‑1 (MQ=255)                                                                                          aaaGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAAc                                                           >  2:414774/1‑151 (MQ=255)                                                                                           aaGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGTGGTTTGAAATTCTCGTCGAACt                                                          >  2:1100289/1‑151 (MQ=255)                                                                                           aaGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACt                                                          >  1:441601/1‑151 (MQ=255)                                                                                              ccGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCAt                                                       >  1:436051/1‑151 (MQ=255)                                                                                                  cTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATacac                                                   >  2:950986/1‑151 (MQ=255)                                                                                                  cTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATacac                                                   >  1:498072/1‑151 (MQ=255)                                                                                                   tCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACAcc                                                  >  2:580342/1‑151 (MQ=255)                                                                                                      cTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGc                                               >  1:323292/1‑151 (MQ=255)                                                                                                       ttGCAGTCGTCTTTCTCATTTTAAACGAATGACTTTTACTTCGCTTTCCCCGGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCa                                                 <  1:837101/148‑1 (MQ=255)                                                                                                        tGCAGTCGGCTTTCTCATTTTAAACGAATGACTTTTACTTCGTTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCg                                              <  1:347490/150‑1 (MQ=255)                                                                                                         gCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGc                                            <  1:345938/151‑1 (MQ=255)                                                                                                          cattCGGCTTTCTCATTTTAAACGAATGTCTTTTACTTCGCTTTCCCCTGTTTTGAAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGc                                            <  2:317086/147‑1 (MQ=255)                                                                                                           aGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCAc                                          >  1:133459/1‑151 (MQ=255)                                                                                                           aGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCAc                                          >  1:271617/1‑151 (MQ=255)                                                                                                           aGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCAc                                          >  1:933104/1‑151 (MQ=255)                                                                                                           aGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCAc                                          >  1:704838/1‑151 (MQ=255)                                                                                                               ggCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGcc                                       >  1:576999/1‑150 (MQ=255)                                                                                                                gCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGcc                                       >  2:272630/1‑149 (MQ=255)                                                                                                                  tttCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGt                                    >  2:754681/1‑150 (MQ=255)                                                                                                                  tttCTCATTTAAAACAAATGCCGTTTCCTCCGCTTTCCCCGGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGt                                    <  2:39100/150‑1 (MQ=255)                                                                                                                   ttCTCATTTAAAACGAATGCCGTTTACTTCGCTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTc                                   <  1:1093246/150‑1 (MQ=255)                                                                                                                       cATTTTAAACGAATGACTTTTACTTCGCTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCgg                                 <  2:617164/148‑1 (MQ=255)                                                                                                                        aTTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGat                              >  2:285371/1‑150 (MQ=255)                                                                                                                        aTTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGat                              >  1:376125/1‑150 (MQ=255)                                                                                                                         ttttAAACGAATGACTTTTTCTTCGTTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGata                             <  2:721076/150‑1 (MQ=255)                                                                                                                         ttttAAACGAATGACTTTTACTTCGTTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGata                             <  2:962901/150‑1 (MQ=255)                                                                                                                         ttttAAACGAATGACGTTTACTTCGCTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGata                             <  2:1165568/150‑1 (MQ=255)                                                                                                                         ttttAAACGAATGACGTTTACTTCGCTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGat                              <  1:389153/149‑1 (MQ=255)                                                                                                                         ttttAAACGAATGACGTTTACTTCGCTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGat                              <  1:495668/149‑1 (MQ=255)                                                                                                                         tttaAAACGACGGACTTTTCCTTCGCTTTCCCCGGTTTTGCATCCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGat                              <  2:466731/149‑1 (MQ=255)                                                                                                                             aaccGAATGCCTTTTCCTCCCCTTTCCCCGGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTa                          <  2:683215/146‑1 (MQ=255)                                                                                                                              aaCGAATGACTTTTCCTTCGCTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGc                       <  2:301007/151‑1 (MQ=255)                                                                                                                                cGAATGCCGTTTACTTCGCTTTCCCCTGTTTTGAAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCt                      <  1:1125903/150‑1 (MQ=255)                                                                                                                                 gAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTc                     >  1:413553/1‑150 (MQ=255)                                                                                                                                  aaTGACTTTTACTTCGCTTTCCCCGGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCaa                   <  1:906603/151‑1 (MQ=255)                                                                                                                                     gACGTTTACTTCCCTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGaa                <  2:156049/151‑1 (MQ=255)                                                                                                                                       cttttCCTTCTCTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGa                 <  2:675967/146‑1 (MQ=255)                                                                                                                                        gTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAAtct             >  1:76373/1‑151 (MQ=255)                                                                                                                                         tttACTCCGTTTTCCCCTGGTTTACAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGa                 <  2:1099315/146‑1 (MQ=255)                                                                                                                                          ttACTTCCTTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAAtctc            <  1:678361/150‑1 (MQ=255)                                                                                                                                           ttcTTCGTTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATctct           <  2:793272/148‑1 (MQ=255)                                                                                                                                             cTTCGCTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCtctt          <  2:1092217/149‑1 (MQ=255)                                                                                                                                             cTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCtcttct        >  1:112227/1‑151 (MQ=255)                                                                                                                                                cGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCg     >  2:521746/1‑151 (MQ=255)                                                                                                                                                 gCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGc    >  1:333252/1‑151 (MQ=255)                                                                                                                                                  cTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCt   >  1:685444/1‑151 (MQ=255)                                                                                                                                                   tttCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCt   <  2:924200/150‑1 (MQ=255)                                                                                                                                                   tttCCCCGGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCt   <  2:982900/150‑1 (MQ=255)                                                                                                                                                   tttACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTc  >  2:43823/1‑151 (MQ=255)                                                                                                                                                      |                                                                                                                                                    CGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTC  >  NC_000913/786706‑787002

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