Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,305,863	0	A	G	77.2%	91.1 / 24.3	59	N245N (AAT→AAC)	yjcS	linear primary‑alkylsulfatase
Reads supporting (aligned to +/- strand):  ref base A (8/5);  new base G (29/15);  total (38/21)
Fisher's exact test for biased strand distribution p-value = 7.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GGGTTTTCGCGCCGCGCAGAGTGTAGAAGTTGTGCAGGGTATGCGTGGCGTTCTCGGCGGTACACAGGGCTTTCAGGGCCGGAATATAGAAGTGCATTTCGGCTGGCGCTTCGCTACCTGGGGTCATCAGGAAGTCAAACTCCAGGCCGTCGATAATCATCTTCTCGCCAGTTCTGACGATAGTTTTCGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCA‑AGGCCATTACCCACATTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTG  >  NC_000913/4305616‑4306102                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                gggTTTTCGCGCCGCGCAGAGTGTAGAAGTTGTGCAGGGTATGCGTGGCGTTCTCGGCGGTACACAGGGCTTTCAGGGCCGGAATATAGAAGTGCATTTCGGCTGGCGCGTCGCTACCTGGGGTCATCAGGAAGTCAAACTCCAGGCCGTCGATAATCATCTTCTCGCCAGTTCTGACGATAGTTTTCGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCA‑AGGCCATTACCCACATTg                                                                                                                                                                                                                                              <  2:413534/251‑1 (MQ=255)   gTTTTCGCGCCGCGTAGAGTGTAGAAGTTGTGCAGGGTATGCGTGGCGTTCTCGGCGGTACACAGGGCTTTCAGGGCCGGAATATAGAAGTGCATTTCGGCTGGCGCTTCGCTACCTGGGGTCATCAGGAAGTCAAACTCCAGGCCGTCGATAATCATCTTCTCGCCAGTTCGGACGATAGTTTTAGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCA‑AGGCCATTACCCACGTTGcc                                                                                                                                                                                                                                            >  1:285730/1‑251 (MQ=255)        cgcgccgcgtAGCGTGTAGAAGTTGTGCAGGGCATGCGTGGCGTTCTCGGCGGTCAACAGGGCTTTCAGGGCCGGAATATAGAAGTGCATTTCGGCTGGCGCTTCGCTACCTGGGGTCATCAGGAAGTCAAACTCCAGGCCGTCGATAATCATCTTCTCGCCAGTTCGGACGATAGTTTTAGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCA‑AGGCCATTACCCACGTTGCCTTgcg                                                                                                                                                                                                                                       >  2:160152/1‑251 (MQ=255)           gccgcgtAGAGTGTAGAAGTTGTGCAGGGTATGCGTGGCGTTCTCGGCGGTACACAGGGCTTTCAGGGCCGGAATATAGAAGTGCATTTCGGCTGGCGCTTCGCTACCTGGGGTCATCAGGAAGTCAAACTCCAGGCCGTCGATAATCATCTTCTCGCCAGTTCGGACGATAGTTTTAGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCA‑AGGCCATTACCCACGTTGCCTTGCGCGt                                                                                                                                                                                                                                    >  1:225825/1‑251 (MQ=255)             cgcgtaaagTGTAGAAGTTGTGCAGGGTATGCGTGGCGTTCTCGGCGGTACACAGGGCTTTCAGGGCCGGAATATAGAAGTGCATTTCGGCTGGCGCTTCGCTACCTGGGGTCATCAGGAAGGCAAAATCCAGGCCGTCGATAATCATCTTCTCGCCAGTTCGGACGATAGTTTTAGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCA‑AGGCCCTTACCCAAGTTGCCTTGCGCGTtg                                                                                                                                                                                                                                  >  2:320902/8‑251 (MQ=255)                         agaagttggtcggggTATGCGTGGCGTTCTCGGCGGTACACAGGGCTTTCAGGGCCGGAATATAGAATTTCATTTCTGCTGGCGCTTCGCTACCTGGGGTCATCAGGAAGTCAAACTCCAGGCCGGCGATAATCATCTTCTCGCCAGTTCGGACGATAGTTTTAGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCA‑AGGCCATTACCCACGTTGCCTTGCGCGTTGTGCGGCAGTAAc                                                                                                                                                                                                                      <  2:325155/239‑1 (MQ=255)                                gtgCAGGGTATGCGTGGCGTTCTCGGCGGTACACAGGGCTTTCAGGGCCGGAATATAGAAGTGCATTTCGGCTGGCGCTTCGCTACCTGGGGTCATCAGGAAGTCAAACTCCAGGCCGTCGATAATCATCTTCTCGCCAGTTCGGACGATAGTTTTAGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCA‑AGGCCATTACCCACGTTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGt                                                                                                                                                                                                               >  1:221507/1‑251 (MQ=255)                                 tgCAGGGTATGCGTGGCGTTCTCGGCGGTACACAGGGCTTTCAGGGCCGGAATATAGAAGTGCATTTCGGCTGGCGCTTCGCTACCTGGGGTCATCAGGAAGTCAAACTCCAGGCCGTCGATAATCATCTTCTCGCCAGTTCTGACGATAGTTTTCGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCA‑AGGCCATTACCCCCATTTCCTTGCGCGTTGTGCGGCAGTAACAGACCGt                                                                                                                                                                                                               >  2:72630/1‑250 (MQ=255)                                         aTGCGTGGCGTTCTCGGCGGTACACAGGGCTTTCAGGGCCGGAATATAGAAGTGCATTTCGGCTGGCGCTTCGCTACCTGGGGTCATCAGGAAGTCAAACTCCAGGCCGTCGATAATCATCTTCTCGCCAGTTCGGACGATAGTTTTAGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCA‑AGGCCATTACCCACGTTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTaga                                                                                                                                                                                                      <  1:26483/251‑1 (MQ=255)                                            cGTGGCGTTCTCGGCGGTACACAGGGCTTTCAGGGCCGGAATCTAGAAGTGCATTTAGGCTGGCGCTTCGCTACCTGGGGTCATCAGGAAGTCAAACTCCAAGACGTCGATAATCATCTTCTCGCCTGTTCGGACGATAGTTTTAGTCGGTGCATTAATGCTCGGGTCGCACGTTGCCAGCGTCACGCCA‑AGGCCATTACCCACGTTGCCTTGCGCGTTGTGCGGAAGTAACCAACCGTAAGAGTagag                                                                                                                                                                                                     >  1:285846/1‑249 (MQ=255)                                                  gTTCTCGGCGGAACACAGGGCTTTAAGGGCCGGAATCTAGAAGTGCATTTCGGCTGGCGCTTCGCTACATAGGGTCATCCGGCAGTCAAACTCCAGGCCGTCGATAATCATCTTCTCGCCAGTACGGACGATCGTTTTAGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCA‑AGGCCATTAACAACGTTGCCTTGCGCGTTGTGCGGCAGTCACAGAACGTAAGAGAAAAGCGCAcggc                                                                                                                                                                                             >  2:366730/1‑251 (MQ=255)                                                       cggcggTACACAGGGCTTTCAGGGCCGGAATATAGAAGTGCATTTCGGCTGGCGCTTCGCTACCTGGGGTCATCAGGAAGTCAAACTCCAGGCCGTCGATAATCATCTTCTCGCCAGTTCTGACGATAGTTTTCGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCA‑AGGCCATTACCCACATTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCt                                                                                                                                                                                         <  1:368885/250‑1 (MQ=255)                                                                 caGGGCTTTCAGGGCCGGAATATAGAAGTGCATTTCGGCTGGCGCTTCGCTACCTGGGGTCATCAGGAAGTCAAACTCCAGGCCGTCGATAATCATCTTCTCGCCAGTTCGGACGATAGTTTTAGTCGGTGCAATAATGCTCGTGTCGCCCGTTGCCAGCGTCACGCCA‑AGGCCATTACCCACGTTGCCTTGCGCGTGGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTAc                                                                                                                                                                              <  1:28044/251‑1 (MQ=255)                                                                  aGGGCTTTCAGGGCCGGAATATAGAAGTGCATTTCGGCTGGCGCTTCGCTACCTGGGGTCATCAGGAAGTCAAACTCCAGGCCGTCGATAATCATCTTCTCGCCAGTTCTGACGATAGTTTTCGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCA‑AGGCCATTACCCACATTGCCTTGCGCGTTGTGCGGCAGTAACAGACAGTACGAGTAGAGCGCACGGCACCTCATGCTGTTAcc                                                                                                                                                                             >  2:308651/1‑251 (MQ=255)                                                                             ggCCGGAATATAGAAGTGCATTTCGGCTGGCGCTTCGCTACATGGGGTCATCAGGAAGTCAAACTCCAGGCAGTCGATAATCATCTTCTCGCCAGTTCTGACGATAGTTTTCGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCA‑AGGCCATTACCCACATTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGt                                                                                                                                                                  >  1:141410/1‑251 (MQ=255)                                                                                   aaTATAGAAGTGCATTTCGGCTGGCGCTTCGCTACCTGGGGTCATCAGGAAGTCAAACTCCAGGCCGTCGATAATCATCTTCTCGCCAGTTCGGACGATAGTTTTAGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCA‑AGGCCATTACCCACGTTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGc                                                                                                                                                            >  2:19495/1‑251 (MQ=255)                                                                                      ataGAAGTGCATTTCGGCTGGAGCTTCGCTACCTGGGGTCATCAGGAAGTCAAACTCCAGGCCGTCGATAATCATCTTCTCGCCAGTTCGGACGATAGTTTTAGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCA‑AGGCCATTACCCACGTTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGCGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTg                                                                                                                                                          >  2:360340/1‑250 (MQ=255)                                                                                        aGAAGTGCATTTCGGCTGGCGCTTCGCTACCTGGGGTCATCAGGAAGTCAAACTCCAGGCCGTCGATAATCATCTTCTCGCCAGTTCGGACGATAGTTTTAGTCGGTGCAATACTGCTCGGGTCGCCCGTTGCCAGCGTCACGCCA‑AGGCCATTACCCACGTTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATg                                                                                                                                                       >  1:413673/1‑251 (MQ=255)                                                                                             tGCATTTCGGCTGGCGCTTCGCTACCTGGGGTCATCAGGAAGTCAAACTCCAGGCCGTCGATAATCATCTTCTCGCCAGTTCGGACGATCGTTTTAGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCA‑AGGCCATTACCCACGTTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATTATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTc                                                                                                                                                  >  1:184713/1‑251 (MQ=255)                                                                                             ggCATTTCGGCTGGCGCTTCGCTCCCTGGGGTCAGCAGGAAGTCAAACTCCAGCCCGTCGATAATCATCTTCTCGCCAGTGCGGACGATAGTTTTAGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCTCGCCA‑AGGCCATTACCCACGTTGCCTTGCGCTTTGTGCGGCAGAAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCAGGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTc                                                                                                                                                  <  1:363691/250‑1 (MQ=255)                                                                                                aTTTCGGCGGGCGCTTCGCTCCCGGGGGTCACCAGGAAGTCAACCTCCAGGCCGTCGATAATCATCTTCTCGCCAGTTCGGACGATAGTTTTAGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCCCGCCA‑AGGCCATTCCCCACGTTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTc                                                                                                                                               <  1:328488/251‑1 (MQ=255)                                                                                                 tttCGGCTGGCGCTTCGCTACCTGGGGTCATCAGGAAGTCAAACTCCAGGCCGTCGATAATCATCTTCTCGCCAGTTCTGACGATAGTTTTCGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCA‑AGGCCATTACCCACATTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTcc                                                                                                                                              >  1:216856/1‑251 (MQ=255)                                                                                                       cTGGCGCTTCGCTACCTGGGGTCATCAGGAAGTCAAACTCCAGGCCGTCGATAATCATCTTCTCGCCAGTTCGGACGATAGTTTTAGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCA‑AGGCCATTACCCACGTTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATaa                                                                                                                                          >  1:38780/1‑249 (MQ=255)                                                                                                               tCGCTCCCGGGGGCCATAAGAAAGTCAAACCCCAGGCCGTATATAATCTCCTTCTCGCAAGTTCGTAGGATAGTTTTAGTCGGTGCAAGAAGATTCGTGTCGTCCGTTGCCAGCTTCCCTCAA‑AGGCCATTCCCCCGGTTGCTTGGCGCGTTGTGCGGCATTAACAGACCGTAAGGTTAGAGCCCACGGCGGCCCATGATTTTTCCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCa                                                                                                                                  <  1:406889/249‑1 (MQ=255)                                                                                                                   tACCTGGGGTCATCAGGAAGTCAAACTCCAGGCCGTCGCTAATCATCTTCTCGCCAGTTCGGACGATAGTTTTAGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCA‑AGGCCATTACCCACGTTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGc                                                                                                                             >  2:371726/1‑250 (MQ=255)                                                                                                                           gTCATCAGGAAGTCAAACTCCAGGCCGTCGATAATCATCTTCTCGCCAGTTCGGACGATAGTTTTAGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCA‑AGGCCATTACCCACGTTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTAcc                                                                                                                                                                             <  1:352792/194‑1 (MQ=255)                                                                                                                           gTCATCAGGAAGTCAAACTCCAGGCCGTCGATAATCATCTTCTCGCCAGTTAGGACGCTAGTTTTAGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCA‑AGGCCATTACCCACGTTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGAAGAGCGCACGGCGGCCCATGATGTTAcc                                                                                                                                                                             >  2:352792/1‑194 (MQ=255)                                                                                                                                caGGAAGTCAAACTCCAGGCCGTCGATAATCATCTTCTCGCCAGTTCGGACGATAGTTTTAGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCA‑AGGCCATTACCCACGTTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACtt                                                                                                               >  1:243343/1‑251 (MQ=255)                                                                                                                                  ggAAGTCAAACTCCAGGCCGTCGATAATCATCTTCTCGCCAGTTCTGACGATAGTTTTCGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCA‑AGGCCATTACCCACATTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTg                                                                                                             >  1:225690/1‑251 (MQ=255)                                                                                                                                  ggAAGCCAAACTCCAGGCCGTCGATAATCTTCTTCTCCCCAGTTCTGACGATAGTTTTCGCCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCA‑AGGCCATTACCCACATTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTCCCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTg                                                                                                             <  1:72630/251‑1 (MQ=255)                                                                                                                                    aagtaaacccccAGGCCTTAGATATTCATCTTCTCGCCAGTTCTGACGATAGTTTTCTGCGGTGCAATAATGCCCGGGTCGCCCGTTGCCAGCGTCACGCCA‑AGGCCATTACCCACATTGCCTGTCGCGTGGTGCGGCAGTAACAGCCCGTAAGTGTAGAGCGCCGGGGGGCTCATGAGGTTCCCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGc                                                                                                            <  2:351607/240‑1 (MQ=255)                                                                                                                                    aaGTAAAACTCCAGGCCGTCGATAATCATCTTCCCGCCAGTTCGGACGATAGTTTTAGCCGGTGAAATAATGTTCGGGTCGCCCGTGGCCAGCGTCACGCAA‑AGGCCATTCCCCACGTTGCCTTGCGCGTTGGGCGGCAGTAACAGACCGTAAGAGTGGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCCCCTGAACTTTGcc                                                                                                           <  1:372839/251‑1 (MQ=255)                                                                                                                                     aGTCACACTCCAGGCCGTCGATAATCATCTTCTCGCCAGTTCGGACGATAGTTTTAGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCA‑AGGCCATTACCCACGTTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCg                                                                                                          >  2:416570/1‑251 (MQ=255)                                                                                                                                              ccAGGCCGTCGATAATCATCTTCTCGCCAGTTGGGACGATAGTTTTAGTCGTGGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCA‑AGGCCATTCCCCACGTTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCTTGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTaa                                                                                                   <  1:334667/249‑1 (MQ=255)                                                                                                                                              ccAGGCCGTCGATAATCATCTTCCCGCCAGTTCGGACGATAGTTTTAGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCA‑AGGCCATTACCCACGTTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACa                                                                                                 <  1:19495/251‑1 (MQ=255)                                                                                                                                                 ggCCGCCGATAATCTTCTTCTCGCCATTTCGGACGATAGTTTTAGTCGGTGCAATAAGGCTGGGGTCGCCCGTTGCCAGCGCCACTCCA‑AGGCCATTACCCACGTTGCTTTGCGCGTGGTGCGGCAGTAACAGACCGTAAGAGTAGGGCGCACGGCGCCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCg                                                                                              <  1:360340/251‑1 (MQ=255)                                                                                                                                                       cGATAATCATCTTCTCGCCAGTTCGGACGATAGTTTTAGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCA‑AGGCCATTACCCACGTTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGACTTc                                                                                         >  1:204108/1‑250 (MQ=255)                                                                                                                                                        gATAATCATCTTCTAGCCAGTTCGGCCGATAGTTTTAGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCA‑AGGCCATTACCCACGTTGCATTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCCTGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCgaattc                                                                                         >  2:51254/1‑244 (MQ=255)                                                                                                                                                         aTAATCATCATCTCGCCAGTTCTGACGATAGTTTTCGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCA‑AGGCCATTACCCACATTGCCTTGCGCGTTGTGCGGCAGTAACAGCCCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTc                                                                                         >  1:132781/1‑248 (MQ=255)                                                                                                                                                                tcttctCGCCAGTTCGGACGATAGTTTTAGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCA‑AGGCCATTACCCACGTTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGACTTCAGAGATAATg                                                                               >  1:383126/1‑251 (MQ=255)                                                                                                                                                                    ctcGCAAGTTCTGACGATAGTTTTCGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCA‑AGGCCATTACCCACATTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCtt                                                                           >  1:303289/1‑251 (MQ=255)                                                                                                                                                                                gACGATAGTTTTAGTCGGTGCAATAATGCTCGGGTCGACCGTTGCCAGCGTCACGCCA‑AGGCCATTACCCACGTTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACCTCGCCTTCAGAGATAATGCCTTTCACGCCACCAt                                                               >  1:198363/1‑251 (MQ=255)                                                                                                                                                                                 aCGATAGTTTTAGTCGGTGCAATAATGCTCGGGCCGCCCGTTGCCAGCGTCACGCCA‑AGGCCTTTACCCACGTTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCAGGGCGGCTCAGGTTTTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGACTTCAGAGATAATGCCTTTCACGCCACCATa                                                              <  1:371726/251‑1 (MQ=255)                                                                                                                                                                                  cGATAGTTTTAGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCA‑AGGCCATTACCCACGTTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACg                                                                                                                                                                   >  1:293460/1‑149 (MQ=255)                                                                                                                                                                                     tAGTTTTAGTCGGTGCAATAATGATCGGGTCGCCCGTTGCCAGAGTCACGACA‑AGACCATTACCCACGTTGCCTTACGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGCGCGCACGGCGGCACATGATGTTACCCGCAAGACCGGTTTCGCTGATAGCTTCGTCCCTAAAGACTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAAAATCGACTTCAGAGATAAAGCCTTTCACGgcaccagagagg                                                          >  2:95362/1‑245 (MQ=255)                                                                                                                                                                                             gtcgGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCA‑AGGCCATTACCCACGTTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGACTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGtgtg                                                  >  2:96980/1‑251 (MQ=255)                                                                                                                                                                                                          tGCTCGGGTCGCCCGTTGCCAGCGTCACGCCA‑AGGCCATTACCCACGTTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGCTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGCCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAgtgt                                       >  1:124436/1‑249 (MQ=255)                                                                                                                                                                                                            cccGGGTCGTCGGTTGCCAGCGTCA‑GCCACAGGCCATTACCCACAGTGCCTTGCGCTTTGTGCGCGAGAAACATCCCTTAAGAGTAGAGCGCACGGGGGCTAATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGAt                                   <  2:225690/249‑1 (MQ=255)                                                                                                                                                                                                               gggTCGCCCGTTGCCAGCGTCACGCCA‑AGGCCATTACCCACGTTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGACTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAAc                                >  1:421942/1‑251 (MQ=255)                                                                                                                                                                                                                      ccGTTGCCAGCGTCCCGCCA‑AGGCCATTCCCCTCGTTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGGGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGACTTCAGAGATAATGCCTTTCACGCCGCCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGcaaca                         <  1:160152/251‑1 (MQ=255)                                                                                                                                                                                                                       cGTTGCCAGCGTCACGCCA‑AGGCCATTACCCACATTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAgt                                         >  1:158036/1‑234 (MQ=255)                                                                                                                                                                                                                            ccAGCGTCACGCCA‑AGGCCATTCCCCACGTTGCTTTGCGCGTTGTGCGGCAGAAACAGCCCGTAAGAGTAGAGCGCCCGGCGGCTCATGATGTTCCCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGACTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGt                   <  1:51254/251‑1 (MQ=255)                                                                                                                                                                                                                                   cACGCCA‑AGGCAATTACCCAAGTTCCCTTGCGCGTTGTGCGGCAGACACAGACCGTAAGAGTAGAGCGCAAGGCGGCTCATGATATTACCCGCCAGCACGCTTTCGCTGATGGCTTCGTCCATAAAACCTGCAGGAGCAATCACCTGAACTTTGCCGGATTTAACATCGACTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTAGGAGCGGCTGGGAGTGTAGATAACGGCAACAATCGGGTTTTGCg            >  2:68197/1‑251 (MQ=255)                                                                                                                                                                                                                                    aCGCCA‑AGGCCATTACCCACCTTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCCGCACGTTTTCGCTGATGGCTTCGTCCATACAGCCTGCAGGGGCCATCACCTGAACTTTGCCGGATTTAACCTCGGCTTCAGAGATAATGCCTTTCACGCCACCATCGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCgg           >  2:207917/1‑251 (MQ=255)                                                                                                                                                                                                                                         catggCCATTACCCACGTTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGACTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATATAGGCTACAAACGGATTTTGCGGAca        >  2:349739/4‑249 (MQ=255)                                                                                                                                                                                                                                             ccaaTTACCCCCTTTGCCTTGCGTGTTGGGCGGCTGTAAGAGACGGAAAGAGTAGGGCGCGCGGCGGCTCAGGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGGCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGACTGCAGAGATAATGCCTGTCACGCCACCGTAGTGTTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCt   <  2:413673/248‑1 (MQ=255)                                                                                                                                                                                                                                              ccATTACCCACGTTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGAGCAATCACCTGAACTTTGCCGGATTTAACATCGACTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTg  >  1:191193/1‑251 (MQ=255)                                                                                                                                                                                                                                               cATTACCCACGTTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGACTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCtg                                              >  2:117140/1‑206 (MQ=255)                                                                                                                                                                                                                                               cATTACCCACGTTGCCTGGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTCCCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGACTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTATGGCtg                                              <  1:117140/206‑1 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                GGGTTTTCGCGCCGCGCAGAGTGTAGAAGTTGTGCAGGGTATGCGTGGCGTTCTCGGCGGTACACAGGGCTTTCAGGGCCGGAATATAGAAGTGCATTTCGGCTGGCGCTTCGCTACCTGGGGTCATCAGGAAGTCAAACTCCAGGCCGTCGATAATCATCTTCTCGCCAGTTCTGACGATAGTTTTCGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCA‑AGGCCATTACCCACATTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTG  >  NC_000913/4305616‑4306102

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