Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,604,344	0	C	T	48.1%	37.9 / 29.8	52	S62F (TCT→TTT)	bglJ	bgl operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base C (10/17);  new base T (25/0);  total (35/17)
Fisher's exact test for biased strand distribution p-value = 4.34e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGC‑GTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGG  >  NC_000913/4604194‑4604490                                                                                                                                                       |                                                                                                                                                    cACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttc                                                                                                                                                     <  2:661321/151‑1 (MQ=255)  aCTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttct                                                                                                                                                    <  2:405139/151‑1 (MQ=255)   cTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttttt                                                                                                                                                   >  1:63735/1‑151 (MQ=255)    tCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAtttttttttt                                                                                                                                                  >  1:440512/1‑151 (MQ=255)     cATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttctt                                                                                                                                                   <  1:435297/149‑1 (MQ=255)      aTAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGTGGCGGTtttttttttttttt                                                                                                                                                >  2:349734/1‑151 (MQ=255)         gAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCtttttttttttttt                                                                                                                                              >  1:597547/1‑148 (MQ=255)           aaaaTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTGATTTTTTTTTTTTTTGcc                                                                                                                                           >  1:3857/1‑151 (MQ=255)           aaaaTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGcc                                                                                                                                           >  1:471732/1‑151 (MQ=255)             aaTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTCTTTTTCTGCCAt                                                                                                                                         >  2:121427/1‑151 (MQ=255)              aTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATg                                                                                                                                        <  2:478770/151‑1 (MQ=255)                gCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATga                                                                                                                                       <  1:37810/150‑1 (MQ=255)                     aTGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATgagga                                                                                                                                    >  2:974132/1‑146 (MQ=255)                       gAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGTg                                                                                                                                 >  1:369049/1‑149 (MQ=255)                        agtagtATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAgc                                                                                                                              <  1:99317/151‑1 (MQ=255)                        agtagtATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCTTTTTTTCTTTTTCTTCCCTGAAGAGTg                                                                                                                                 >  2:820552/1‑148 (MQ=255)                           agtaTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTGCCATTAGAAGTGAgcg                                                                                                                             >  2:953145/1‑149 (MQ=255)                            gtaTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGTTAACCCTTATAAGCTCCATACCTATACCAGTCAGGTGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCaa                                                                                                                          <  2:107605/151‑2 (MQ=255)                                   tATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCTTTTTTTTTTTTTTTTCCATGAGAAGGGAGCGCAGAGAggg                                                                                                                    >  1:385615/1‑150 (MQ=255)                                     ttGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTGCCATGAGAAGTGGGCGCAGAGAGGGGtt                                                                                                                 >  1:1012325/1‑151 (MQ=255)                                         gagTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGGAGGGAGGGCAGGGGGGGGTTATCt                                                                                                             >  2:25706/1‑151 (MQ=255)                                          agTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTTCCATGCGAAGTGAGCGCCGAGAGGGGTTAtttt                                                                                                            >  2:940408/1‑151 (MQ=255)                                          agTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCtt                                                                                                            <  1:205597/151‑1 (MQ=255)                                               aTTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCctg                                                                                                       <  1:25741/151‑1 (MQ=255)                                                ttCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGGAGTGAGCGCGGAGAGGGATTTTCTTTcc                                                                                                         >  2:691638/1‑148 (MQ=255)                                                 tCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCCTGAGAAGTGAGCGCCGAGAGGGGTTATCTTTCctgac                                                                                                     >  2:439418/1‑151 (MQ=255)                                                       aGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACt                                                                                               <  2:917075/151‑1 (MQ=255)                                                        gTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCCTGTGAAGGGAGCGCAGAGGGGGATTATCTTGCCTGACTGAAc                                                                                                >  2:822726/1‑149 (MQ=255)                                                          ttGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGTGGGCGCAGAGAGGGGTTATCTTGCCTGACTGAACTgg                                                                                             >  1:496995/1‑150 (MQ=255)                                                            gCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCg                                                                                           <  2:115707/150‑1 (MQ=255)                                                                         ataAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTTCCATGAGAAGTGAGGGCCGGGGGGGGTTATCTTTCCTGACTTAACTGGCGATTAAGTTTTcgc                                                                              >  2:69348/1‑150 (MQ=255)                                                                         ataAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTTTTCCATGAGAAGTGGGGGCAGGGGGGGGTTATCTTTCCTGACTGAACTGGCGATTTAGTTTCcgc                                                                              >  1:517763/1‑150 (MQ=255)                                                                            aGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCTTTTTTTTTTTTTCTTCCCTGAGAAGTGAGCGCAGAGAGGGGTTATTTTGCCTGACTGAACTGGCGATTAAGTTTCCCC‑GGAc                                                                         >  1:776642/1‑151 (MQ=255)                                                                                  aTACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGC‑GTACCCGGCg                                                                   <  1:728829/151‑1 (MQ=255)                                                                                    aCCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGC‑GTACCCGGCGt                                                                  <  1:569731/150‑1 (MQ=255)                                                                                       tataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCATGAGAAGGGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGC‑GTACCCGGCGttt                                                                >  1:594420/1‑149 (MQ=255)                                                                                            cAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTGCACTTATAAGTGGGGGCAGCGGGGGATTATCTTGCCTGACGTAAAAGGCGATTTAATTTTCGC‑GTACACGGCGGTTTGtt                                                            >  2:827837/1‑148 (MQ=255)                                                                                             aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTTGGAAGTGAGGGCCGAGGGGGGTTATCTTTCCTGACTGAACTTGGGGTTTAGTTTCCGC‑GGACCCCGCGTTTTGTtttt                                                         >  2:936655/1‑150 (MQ=255)                                                                                              gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGC‑GTACCCGGCGTTTAGTTATTGc                                                       <  1:64231/151‑1 (MQ=255)                                                                                               tCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCTATGAGAAGTGAGCGCAGAGAGGGGTTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGC‑GTACCCGGCGTTTAGTTAtt                                                         >  1:610741/1‑148 (MQ=255)                                                                                                cAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTTCCATGAGAAGGGAGGGCAGAGGGGGGTTATCCTTCCTGACTGAACTGGGGATTAAGTTTTCGC‑GTACCCGGCGTTTTGTTATTGCgg                                                     >  1:520225/1‑151 (MQ=255)                                                                                                   gAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGC‑GTACCCGGCGTTTAGTTATTGCGgatg                                                  <  2:184260/151‑1 (MQ=255)                                                                                                      tCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGGTTCCGC‑GTACCCGGCGTTTAGTTATTGCGggtg                                                  >  1:682110/1‑148 (MQ=255)                                                                                                         tttCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGGGAAGTGAGCGCAGAGAGGGGTTATCTTGCCTGACTGTACTGGCGATTAAGTTTTCGC‑GTACCCGGCGGTTAGTTATTGCGggtgat                                                >  2:443455/1‑147 (MQ=255)                                                                                                            cAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGC‑GTACCCGGCGTTTAGTTATTGCGGATGATGATAtt                                          <  1:898317/150‑1 (MQ=255)                                                                                                             agtaTGCCATGTCGCGGATCTCGTTTGCGGCGTTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAATTTTCCGC‑GTACCCGGCGTTTAGTTATTGCGGATGATGATAtt                                          <  2:9014/146‑1 (MQ=255)                                                                                                             aGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCTGTTTTTTTTTTTTTGCCATGAGGAGTGAGCGCAGAGAGGGGTTATTTTGCCTGAATTAAATGGGGATTAAGTTTTCGC‑GGACCCGGGGTTTTGTTATTGAGGATGATTATAtt                                          >  1:740492/1‑149 (MQ=255)                                                                                                               gATGCCATGTCGCGGATCTCGTTTGCGGCGGTTATTTTTTTTTTTTTTTCCATGAGAAGTGAGCGCAGAGGGGGGTTATCTTTTCTGACCTAAATGGGGATTAAGTTTTCGCGGTACCC‑GCGTTTTGTTTTTTCGGATGATGAtttt                                          >  2:10414/1‑147 (MQ=255)                                                                                                                  gCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGC‑GTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCt                                    <  1:686801/150‑1 (MQ=255)                                                                                                                       gTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGC‑GTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCt                               <  2:784004/150‑1 (MQ=255)                                                                                                                           cgGATCTCGTTTGCGGCGGTCCTTTTTTCTTTTTCTGCCCTGAGAAGTGAGGGCAGAGAGGGATTATCTTGCCTTACTGAACTGGCGGTTTAGTTTCCGC‑GGACCCGGCGGTTAGTTATTGCGGATGATGATATTGAAGATCGGCTggt                            >  2:555618/1‑147 (MQ=255)                                                                                                                                  cGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGGGAGGGCAGAGAGGGGTTATCTTTTCTGACTGAACTGGCGATTAAGTTTCCGC‑GTACCCGGGGGTTAGTTTTTGCGGGTGATGATATTGAAGCTCGGCTGGTTTGTTCAtt                   >  2:710548/1‑151 (MQ=255)                                                                                                                                  cGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGC‑GTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCAt                    <  1:98780/150‑1 (MQ=255)                                                                                                                                   gTTTGCGGCGGGCCTTTTTTTTTTTTCTGCCATGAGAAGTGAGCGGAGAGGGGGGTTATCTTTCCTGTATTGACTGGCGATTAAGTTTCCGC‑GTGCCCGGGGGTTTGTTATTTCGGATGATGGTATTGAAGCTCGGCTGGTTTGTTTAtt                   >  2:245683/1‑150 (MQ=255)                                                                                                                                       gcggcgGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGC‑GTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCg               >  1:100233/1‑150 (MQ=255)                                                                                                                                        cggcggTCATTTTTTCTTTTTTTGCCCTGAGAAGTGAGCGCCGAGAGGGATTTTCTTGCCTTTCTGAACTGGCGATTAAGTTTTCGC‑GTACCCGGCGGTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCCTTGTCGcc             >  2:933774/1‑151 (MQ=255)                                                                                                                                            ggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGC‑GTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATc          >  2:438721/1‑150 (MQ=255)                                                                                                                                               cATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGC‑GTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCg      >  2:825842/1‑151 (MQ=255)                                                                                                                                                aTTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGC‑GTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGc     >  1:638300/1‑151 (MQ=255)                                                                                                                                                   ttttctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGC‑GTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTg   <  2:123707/150‑1 (MQ=255)                                                                                                                                                   ttttctttttatGCCATGAGAAGTGAGCGCAGTGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGC‑GTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGCTTGGTTCATTGTCGCCATCACCGCTgg  >  2:921498/1‑151 (MQ=255)                                                                                                                                                    tttctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGC‑GTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTgg  >  1:496521/1‑150 (MQ=255)                                                                                                                                                    tttctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGC‑GTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTgg  <  1:848717/150‑1 (MQ=255)                                                                                                                                                    tttctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGC‑GTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATAACCGCTgg  >  2:11118/1‑150 (MQ=255)                                                                                                                                                       |                                                                                                                                                    CACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGC‑GTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGG  >  NC_000913/4604194‑4604490

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