Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	786,855	0	A	C	39.3%	155.4 / 15.4	87	G247G (GGT→GGG)	gpmA	phosphoglyceromutase 1
Reads supporting (aligned to +/- strand):  ref base A (45/6);  new base C (0/33);  total (45/42)
Fisher's exact test for biased strand distribution p-value = 1.33e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGG‑TTT‑GCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATG  >  NC_000913/786710‑787005                                                                                                                                                  |                                                                                                                                                         agagTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTg                                                                                                                                                     >  1:622598/1‑151 (MQ=255)   agTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGG‑t                                                                                                                                                  <  2:637841/151‑1 (MQ=255)        cGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGG‑TTT‑Gc                                                                                                                                             >  2:806679/1‑150 (MQ=255)        cGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGG‑TTT‑GCa                                                                                                                                            <  2:771728/151‑1 (MQ=255)            tCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGG‑TTT‑GCAAccg                                                                                                                                        >  2:1170337/1‑151 (MQ=255)              cGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGG‑TTT‑GCAAccgcc                                                                                                                                      >  1:1211095/1‑151 (MQ=255)                 tGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGG‑TTT‑GCAACCGCCgct                                                                                                                                   >  1:50014/1‑151 (MQ=255)                  aaaTCCCCACAGGCCTGATAAGTCGCGCGGCGTCCCATCAGGCAATGTCCCCCATTGTAAACAAAAAAAAACCCGACTCTCTTGCGGTCGGTTTTCTCTTTTAAAACAAATGACGTTTTCTCCTCTTTCCCCTGT‑TTG‑GAAACCGCCgctg                                                                                                                                  <  2:1153242/150‑1 (MQ=255)                    aTCCCCAGAGGCATAATAAGTCGCGCGGGGTCGCATCAGGCATTGTGCTCAATTTTTAGCAAAAAAAAGGCCGACTCTCGTGCAGTCGCTTTTCTCATTTTAACCGAATGCCTTTTACTTCCTTTTCCCCTGG‑TTT‑GCAACCGCCgctgc                                                                                                                                 <  2:573405/150‑1 (MQ=255)                     tCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGG‑TTT‑GCAACCGCCGCTGCtt                                                                                                                               >  2:234954/1‑151 (MQ=255)                      ccccacAGGCCTAATAAGCCGCGCGGCGTCGCATCAGGCATTGTGCTCCATTTTTAGCAAAAAAAAGGCCGCCTCACTTGCGGTGGGTTTTCTCATTTTAAACAAATGACGTTTCCTTCGTTTTCCCCGGG‑TTT‑GCAACCGCCGCTGCtt                                                                                                                               <  2:758158/148‑1 (MQ=255)                           caGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGTTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTCCCCTGG‑TTT‑GCAACCGCCGCTGCTTTCGCt                                                                                                                          <  1:316974/150‑1 (MQ=255)                               ccTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGG‑TTT‑GCAACCGCCGCTGCTTTCGCTGCGAt                                                                                                                     >  2:891314/1‑151 (MQ=255)                               ccTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGG‑TTT‑GCAACCGCCGCTGCTTTCGCTGCGAt                                                                                                                     >  1:515956/1‑151 (MQ=255)                               ccTGATAAGTCGCACCCCGTAGCACCAGGCATTGGTCTCCCTGGTTAGCAAAAAAAAACCCGACTAACTGGACTTCGGCTTTCCCATTTAAACCGATTGTCTTTTCTTCCGCTTT‑CCCTGGTTTT‑GCAACCGCCGCTGCTTTCGCTGCg                                                                                                                       <  2:237644/149‑1 (MQ=255)                                    tAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGG‑TTT‑GCAACCGCCGCTGCTTTCGCTGCGATCTCGt                                                                                                                >  1:21640/1‑151 (MQ=255)                                    tAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGG‑TTT‑GCAACCGCCGCTGCTTTCGCTGCGATCTCGt                                                                                                                >  2:230162/1‑151 (MQ=255)                                           cgcAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGG‑TTT‑GCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCAtt                                                                                                         >  2:27764/1‑151 (MQ=255)                                              aGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGG‑TTT‑GCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTAcc                                                                                                      >  2:468798/1‑151 (MQ=255)                                                   cgaaTCAGGAAATGTCCCCCTTTGTTACCAACAAAAAAGCCGACCCTCTTGCGGTCGGCTTTCTCTTTTTAAACGAAGGACGTTTCCTCCGCTTTCCCCTGT‑TTT‑GCATCCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTa                                                                                                        <  2:1176343/141‑1 (MQ=255)                                                   cGCTTCAGGCATTTTCCCCCTTTGTTAGAAAAAAAAACGCCGCCTCCTTTGCATTCGTTTTTCTCTTTTAAACCAAATGCCGTTTCCTTCGCTTTCCCCTGT‑TTT‑GCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGat                                                                                                 <  2:783327/151‑1 (MQ=255)                                                   cGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGCCTCCCTTGCAGTCGGTTTTCTCATTTTAAACGAATGCCTTTTCCTCCTCTTTCCCCTGG‑TTG‑GAAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGa                                                                                                  <  2:252200/150‑1 (MQ=255)                                                      aTCAGTCAATGTGATCCATTGTAACAAACAAAAAGCCCGTCTCCCTTGCGGTCGGCTTTCTCATTTAAAACGAATGCCGTTTCCTTCGCTTCTCCCGGG‑TTT‑GAAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGat                                                                                                 <  2:844617/148‑1 (MQ=255)                                                      aTCAGGCAATGTGCTCCATTGTTAGCAAAAAAAAAGCCGCCTCACGTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTCCCCTGG‑TTT‑GCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGataat                                                                                              <  1:1057904/151‑1 (MQ=255)                                                       tCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGG‑TTT‑GCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGataat                                                                                              >  1:711940/1‑150 (MQ=255)                                                         aggaattgtgCCCCTTTGTTAGCAACAAAAAAGCCGACTCACTTGCGGCCGGTTTCCTCATTTTAAACGAAGGACTTTTACTTCGCTTTCCCCTGG‑TTG‑GCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGataata                                                                                             <  2:727574/143‑1 (MQ=255)                                                          ggCAATGTGCTCCATTGTTAGCAAAAAAAAAGCCGCCTCACTTGCAGTCGGTTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTCCCCTGG‑TTT‑GCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGc                                                                                           <  1:86814/150‑1 (MQ=255)                                                          gcaattgtcCTCAATTGTTAAAAAAAAAAAAGCCGACTCACTTGCAGTCGGTTTTCTCATTTAAAACGAATGACGTTTACTTCGCTTTCCCCTGT‑TTT‑GCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAg                                                                                            <  2:668302/144‑1 (MQ=255)                                                               tgtgCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGG‑TTT‑GCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTc                                                                                      >  1:1099936/1‑150 (MQ=255)                                                                 tgCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGG‑TTT‑GCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGc                                                                                   >  2:259927/1‑151 (MQ=255)                                                                  gCTCCATTGTTAGAAACAAAAAAGCCGACTCCCTTGCAGTCGGTTTTCTCTTTTTAACCGAATGACGTTTACTTCGTTTTCCCCTGG‑TTT‑GCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGtt                                                                                        <  2:457140/145‑1 (MQ=255)                                                                   cTCCATTTTTAGAAAAAAAAAAGCCGACTCCCTTGCAGTCGGTTTTCTCATTTAAAACGAATGCCTTTTCCTTCGCTTTCCCTTGT‑TTT‑GAAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCa                                                                                     <  2:425371/147‑1 (MQ=255)                                                                     ccATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGTTTTCTCATTTAAAACGAATGACGTTTACTTCGCTTTCCCCTGG‑TTT‑GCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGtt                                                                               <  1:984025/151‑1 (MQ=255)                                                                     ccATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGG‑TTT‑GCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGtt                                                                               >  2:820971/1‑151 (MQ=255)                                                                      cATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGG‑TTT‑GCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGttt                                                                              >  2:910508/1‑151 (MQ=255)                                                                      cATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGG‑TTT‑GCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGttt                                                                              >  1:187264/1‑151 (MQ=255)                                                                      cATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGG‑TTT‑GCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGttt                                                                              >  1:441717/1‑151 (MQ=255)                                                                           ttAGAAACAAAAAAGCCGCCTCCCTTGCAGTCGGCTTCCTCATTTAAAACGAATGACGTTTACTTCCTTTTCCCCTGG‑TTT‑GCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGtt                                                                               <  1:1170418/145‑1 (MQ=255)                                                                              gCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGG‑TTT‑GCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATtct                                                                      >  2:1208104/1‑151 (MQ=255)                                                                                aaaaaaaaaGCCGACTCCCTTGCAGTCGGTTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTCCCCTGG‑TTT‑GCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATtct                                                                      <  1:520260/146‑1 (MQ=255)                                                                                aaaaaaaaaCCCGACTCCCTTGCAGTCGGCTTTCTCATTTAAAACGAATGACGTTTACTTCGCTTTCCCCTGG‑TTT‑GCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATtct                                                                      <  1:263648/146‑1 (MQ=255)                                                                                 aacaaaaaCCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGCCTTTTACTTCGTTTTCCCCTGT‑TTT‑GCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCtcgt                                                                   <  1:1010416/148‑1 (MQ=255)                                                                                  aaaaaaaCCCGCCTCCCTTGCAGCCGGTTTTCTCTTTTAAAACGAATGCCGTTTCCTTCGCTTTCCCTTGT‑TTT‑GAAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCtcgtc                                                                  <  2:492855/150‑1 (MQ=255)                                                                                   aaaaaaGCCGCCTCACTTGCGGTGGGTTTTCTCATTTAAAACGAATGCCGTTTCCTTCGCTTTCCCCTGG‑TTT‑GCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCtcgtcg                                                                 <  2:1135261/151‑1 (MQ=255)                                                                                    aaaaaGACGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGG‑TTT‑GCAAACGCCGCTGATTTCGCTGCGATCACGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGCAATTCtcgtcg                                                                 >  2:525933/1‑150 (MQ=255)                                                                                     aaaaCCCGCCTCACTTGCGGTCGTTTTTCTCATTTAAAACGAATGACGTTTCCTTCGCTTTCCCCTGT‑TTT‑GCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGaa                                                               <  1:1229440/151‑1 (MQ=255)                                                                                      aaaGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGG‑TTT‑GCAACCGCCGCTGCTTTCGCTGCGATCTTGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAAc                                                              >  2:869591/1‑151 (MQ=255)                                                                                       aaGCCGACTCACTTGCAGTCGGTTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTCCCCTGG‑TTT‑GCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACt                                                             <  1:508215/151‑1 (MQ=255)                                                                                           cGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGG‑TTT‑GCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCAt                                                          >  2:911754/1‑150 (MQ=255)                                                                                            gACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGG‑TTT‑GCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATa                                                         >  2:649340/1‑150 (MQ=255)                                                                                              cTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGG‑TTT‑GCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATacac                                                      >  2:526815/1‑151 (MQ=255)                                                                                               tCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGG‑TTT‑GCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACAcc                                                     >  2:327663/1‑151 (MQ=255)                                                                                                   ttGCTGTCGTCTTTCTAATTTTAAACAAATGACTTTTCCTTCGCTTTCCCCTGT‑TTT‑GCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGc                                                  <  2:552098/150‑1 (MQ=255)                                                                                                     gCAGTCGGTTTTCTCATTTAAAACGAATGACTTTTCCTTCGCTTTCCCCTGG‑TTG‑GCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCgg                                                <  1:1029880/150‑1 (MQ=255)                                                                                                      cAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGG‑TTT‑GCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCgg                                                >  2:1166551/1‑149 (MQ=255)                                                                                                       aGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGG‑TTT‑GCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCAc                                             >  1:312719/1‑151 (MQ=255)                                                                                                           ggttttCTCATTTTAAACAAATGACGTTTACTTCGCTTTCCCCTGG‑TTT‑GCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATCCTCGTCGAACTCATACACCAGCGGCACGc                                           <  1:784761/146‑1 (MQ=255)                                                                                                            gCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGG‑TTT‑GCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGcc                                          >  2:634265/1‑149 (MQ=255)                                                                                                             cTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGG‑TTT‑GCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGt                                       >  2:702104/1‑151 (MQ=255)                                                                                                              tttCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGG‑TTT‑GCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTc                                      >  2:1231055/1‑151 (MQ=255)                                                                                                               ttCTCATTTAAAACGAATGCCTTTTACTCCGTTTTCCCCTGG‑TTT‑CCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTc                                      <  1:295065/150‑1 (MQ=255)                                                                                                                     tttaAAACGAGTGACGTTTACTTCGTTTTCCCTTGG‑TTA‑AAAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAATCGGGat                                 <  1:82302/149‑1 (MQ=255)                                                                                                                     tttaAAACAAATGACGTTTCCTCCGTTTTCCCCTGT‑TTT‑GAAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGa                                  <  2:311893/148‑1 (MQ=255)                                                                                                                       ttAAACGAATGACGTTTACTTCGCTTTACCCTGG‑TTT‑GCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATAtt                              >  2:262868/1‑150 (MQ=255)                                                                                                                         aagcGAATGACGTATACTCACCTTCACCGTGT‑TTA‑GACACCGCCGCTGCTTTCGCTGCGATCTCGTCATCATTACCCAGATAATAGCGTTTCAGCGGATTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGGCGGGata                                <  2:1267209/143‑1 (MQ=255)                                                                                                                         aaaCAAATGCCTTTTCCTTCGTTTTACCCTGG‑TTT‑CCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAg                           <  2:1215086/151‑1 (MQ=255)                                                                                                                          aacaattGACGTTTCCTTCGCTTTCCCTTGG‑TTT‑CCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGc                          <  2:1172148/145‑1 (MQ=255)                                                                                                                          aaCGAATGACGTTTACTTCGCTTTACCCTGG‑TTT‑GCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGc                          >  1:1014975/1‑151 (MQ=255)                                                                                                                           aCGAATGGCTTTTACTCCGCTTTCCCTGGG‑TTT‑GCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCt                         <  2:541524/151‑1 (MQ=255)                                                                                                                               aTGACGTTTACTTCGCTTTACCCTGG‑TTT‑GCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCa                       >  2:578699/1‑149 (MQ=255)                                                                                                                                tGACGTTTCCTTCTCTTTACCCTGT‑TTT‑GCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAg                     <  2:61990/150‑1 (MQ=255)                                                                                                                                tGACGTTTACTTCGCTTTACCCTGG‑TTT‑GCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCa                       >  2:1040859/1‑148 (MQ=255)                                                                                                                                  aCGTTTACTTCGCTTTACCCTGG‑TTT‑GCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAAt                  >  2:1076438/1‑151 (MQ=255)                                                                                                                                    ttttCCTTCGTTTTCCCCTGG‑TTT‑CAAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGa                    <  1:1138077/146‑1 (MQ=255)                                                                                                                                    gTTTACTTCGCTTTACCCTGG‑TTT‑GCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCGGTCGGGATATTAAGCTCAAGAAtct                >  2:1176553/1‑151 (MQ=255)                                                                                                                                     tttACTTCGCTTTACCCTGG‑TTT‑GCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAAtctc               >  1:343862/1‑151 (MQ=255)                                                                                                                                       tACTTCGCTTTACCCTGG‑TTT‑GCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCtctt             >  1:477813/1‑151 (MQ=255)                                                                                                                                        aCTTCGCTTTACCCTGG‑CTT‑GCAACCACCGCTGCTTTCGCTGGGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAAGTCACACACAAGCGGCACGCCAGTCGGGATATCAAGCTCAAGAATCtcttc            >  2:1072261/1‑151 (MQ=255)                                                                                                                                         ccttccTTTTCCCTTG‑TTTAGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATctct              <  2:204693/144‑1 (MQ=255)                                                                                                                                           tCGCTTTACCCTGG‑TTT‑GCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCttcttc         >  1:1090366/1‑151 (MQ=255)                                                                                                                                             gCTTTACCCTCT‑TTT‑CAATCCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGc       <  2:1211012/151‑1 (MQ=255)                                                                                                                                               tttCCCCTGT‑TTT‑GCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTcttctt          <  2:659809/146‑1 (MQ=255)                                                                                                                                               tttCCCCTGT‑TTT‑GAAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCg        <  2:171251/148‑1 (MQ=255)                                                                                                                                               tttACCCTGG‑TTT‑GCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTc     >  1:688405/1‑151 (MQ=255)                                                                                                                                                ttACCCTGG‑TTT‑GCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCCCGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCt      >  2:714156/1‑149 (MQ=255)                                                                                                                                                 tACCCTGG‑TTT‑GCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCAt   >  1:554753/1‑151 (MQ=255)                                                                                                                                                 tACCCTGG‑TTT‑GCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCAt   >  1:1193354/1‑151 (MQ=255)                                                                                                                                                  aCCCTGG‑TTT‑GCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATg  >  1:283197/1‑151 (MQ=255)                                                                                                                                                  |                                                                                                                                                         AGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGG‑TTT‑GCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATG  >  NC_000913/786710‑787005

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