Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,305,803	0	C	A	72.6%	76.0 / 31.1	62	T265T (ACG→ACT)	yjcS	linear primary‑alkylsulfatase
Reads supporting (aligned to +/- strand):  ref base C (12/5);  new base A (30/15);  total (42/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.00e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CGTCGTTACCCCACATATCCAGCGTTTCGTTCAGATACTCGGTCCACTTGCTGGTGTCGCGGGTTTTCGCGCCGCGCAGAGTGTAGAAGTTGTGCAGGGTATGCGTGGCGTTCTCGGCGGTACACAGGGCTTTCAGGGCCGGAATATAGAAGTGCATTTCGGCTGGCGCTTCGCTACCTGGGGTCATCAGGAAGTCAAACTCCAGGCCGTCGATAATCATCTTCTCGCCAGTTC‑TGACGATAGTTTTCGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCAAGGCCATTACCCACATTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGG  >  NC_000913/4305556‑4306046                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                    ggtcgtTACCCCACATATCCAGCGGTTCGTTCAGATTCTCGGTCCGCTGGCTGGTGTCGCGGGTTTTCGCGCCGCGCAGAGTGTAGAAGTTGTGCATGGTTTGCGTGGCTTTCTCGGCGGTACACAGGGCTTTCAGGGCCGGAATATAGAAGTTCATGTCTGCTGGCGCTGCGCTACCTGGGGTCATCAGGAAGTCAAACTCCAGGCCGTCGATAATCATCTTCTCGTCAGTTC‑TGACGATAGTTTtcgtc                                                                                                                                                                                                                                                  <  2:170095/250‑1 (MQ=255) cgtcgtTACCCCACATATCCAGACTTTCGTTCAGATACTCGGTCCACTTACTGGTGTCGCGGGTTTCCGCGCCGCGTAGAGTGAAGAAGATGTGCAGGGTATGCGTGGCGTTCTCGGAGGGACACAGGGCTTTCAGGGCCGGAATATAGGAGTGCATTTCGGCTGGCGCTTCGCTACATGGGGTCATCAGAAAGTCAAACTCCAGGCCGTCGATAATCATCTTCTCGGCAGTTC‑GGAAGATTGTTTtagtc                                                                                                                                                                                                                                                  >  1:110916/1‑251 (MQ=255)             ccctatCCAGTGTTTCGTTCAGAGACTCGGGCCTCTTGCTGGTGTCGCGGGTTTTCGCGCCGCGTAGAGTGTAGAAGTTGTGCAGGGTATGCGTGGCGTTCGCGGCGGTACACAGGGCTTTCAGGGCCGGAATATAGAAGTGCATTGCGGCTGGCGCTTCGCTACCTGGGGTCATCAGGAAGTCAAACTCCAGGCCGTCGATAATCATCTTCTCGCCAGTTC‑GGACGATAGTTTTAGTCGGTGCaataagt                                                                                                                                                                                                                                      <  2:227516/248‑3 (MQ=255)               ttatCCAGCGTTTCGTTCTTATACTCAGTACACTTGCTATTGTCGCGGGTTTTCGCGCCGCGTAGAGGGTAGAAGTTGTACAGGGTATGCGTGGCGTTCACGGCGGTACACAGGGAGTTCAGAGCCGGAATATAGAAGTGCATTTCGGCTGGCGCTTCGCGAACTGAGGTCATCAAGAAGACAAACTACAGGACCTCGACAATCATCTTCTCGCCAGTTCGTGA‑GAAAATTTTAGTCGGTGCAATAATGCt                                                                                                                                                                                                                                    >  2:114893/2‑251 (MQ=255)                   ccAGCGTTTCGTTCAGATACTCGGTCCACTTGCTGGTGTCGCGGGTTTTCGCGCCGCGTAGAGTGTAGAAGTTGTGCAGGGTATGCGTGGCGTTCTCGGCGGTACACAGGGCTTTCAGGGCCGGAATATAGAAGTGCATTTCGGCTGGCGCTTCGCTACCTGGGGGCATCAGGAAGTCAAACTCCAGGCCGTCGATAATCATCTTCTCGCCAGTTC‑GGACGATAGTTTTAGTCGGTGCAATAATGCTCggg                                                                                                                                                                                                                                >  2:357568/1‑251 (MQ=255)                    cAGCGTTTCGTTCAGATACTCGGTCCACTTGCTGGTGTCGCGGGTTTTCGCGCCGCGCAGAGTGTAGAAGTTGTGCAGGGTATGCGTGGCGTTCTCGGCGGTACACAGGGCTTTCAGGGCCGGAATATAGAAGTGCATTTCGGCTGGCGCTTCGCTACCTGGGGTCATCAGGAAGTCAAACTCCAGGCCGTCGATAATCATCTTCTCGCCAGTTC‑TGACGATAGTTTTCGTCGGTGCAATAATGCTCGGGt                                                                                                                                                                                                                               >  1:248055/1‑251 (MQ=255)                       cGTTTCGTTCAGATACTCGGTCCACTTGCTGGTGTCGCGGGTTTTCGCGCCGCGCAGAGTGTAGAAGTTGTGCAGGGTATGCGTGGCGTTCTCGGCGGTACACAGGGCTTTCAGGGCCGGCATATAGAAGTGCATTTCGGCTGGCGCTTCGCTACCTGGGGTCATCAGGAAGTCAAACTCCAGGCCGTCGATAATCATCTTCTCGCCAGTTC‑TGACGATAGTTTTCGTCGGTGCAATAATGCTCGGGTCGc                                                                                                                                                                                                                            >  1:348337/1‑251 (MQ=255)                        gTTTCGTTCAGGTACTCGGTCCACTGGCTGGTGTCGCGGGTTTTCGCGCCGCGTAGAGTGTAGAAGTTGTGCAGGGTATGCGTGGCGTTCTCGGCGGTACACAGGGCTTTCAGGGCCGGCATATAGACGTGCCTTTCGGCTGGCGCTTCGCTACCTGGGGTCATCAGGAAGTCAAACTCCAGGGCGTCGATAATCATCTTCTCGCCAGTTC‑GGGCGGTAGTTTTACTCGGTGCAATAAAGCTCGGGTCGcc                                                                                                                                                                                                                           >  2:19310/1‑251 (MQ=255)                         tttCGTTCAGATACTCGGTCCACTTGCTGGGGTCGCGGGTTATCGCGACGCGCAGCGCGAAGAAGATGTGCAGGGGATGCGTGGCGTTCTCGGCGGTACCCAGGGCTTTCAGGGCCGGAATATAGCAGAGCATTTCGGCTGGCGCTTCGCTACCTGGGGCCATCAGGAAGTCAAACTCCAGGACGTCGATAATCATCTTTTCGCCAGTTC‑GGACGATAGTTTTAGTCGGTGGAATAATGCTCGAGTCGccc                                                                                                                                                                                                                          >  2:28044/1‑251 (MQ=255)                          ttCGTTCAGATACTCGGTCCACTTGCTGGTGTCGCGGGTTTTCGCGCCGCGTAGAGTGTAGAAGTTGTGCAGGGTATGCGTGGCGTTCTCGGCGGTACACAGGGCTTTCAGGGCCGGAATATAGAAGTGCATTTCGGCTGGCGCTTCGCTACCTGGGGTCATCAGGAAGTCAAACTCCAGGCCGTCGATAATCATCTTCTCGCCAGTTC‑GGACGATAGTTTTAGTCGGTGCAATAATGCTCGGGTCGCCCg                                                                                                                                                                                                                         >  1:160126/1‑251 (MQ=255)                               tCAGATACTCGGTCCACTTGCTGGTGTCGCGGGTTTTCGCGCCGCGTAGAGTGTAGAAGTTGTGCAGGGTATGCGTGGCGTTCTCGGCGGTACACAGGGCTTTCAGGGCCGGAATATAGAAGTGCATTTCGGCTGGCGCTTCGCTACCTGGGGTCCTCAGGAAGTCAAACTCCAGGCCGTCGATAATCATCTTCTCGCCCGTTC‑GGACGATAGTTTTAGTCGGTGCAATAATGCTCGGGTCGCCCTTTGcc                                                                                                                                                                                                                    >  2:328488/1‑251 (MQ=255)                                cAGATACTCGGTCCACTTGCTGGTGTCGCGGGTTTTCGCGCCGCGTAGAGTGTAGAAGTTGTGCAGGGTATGCGTGGCGTTCTCGGCGGTACACAGGGCTTTCAGGGCCGGAATATAGAAGTGCATTTCGGCTGGCGCTTCGCTACCTGGGGTCATCAGGAAGTCAAACTCCAGGCCGGCGGTAATCATCTTCTCGCCCGTTC‑GGACGATAGGTTTAGTCGGTGCAATACTGCTCGGGGCGCCCGTTGccc                                                                                                                                                                                                                   >  2:334667/1‑250 (MQ=255)                                     aCTCGGTCCACTTGCTGGTGTCGCGGGTTTTCGCGCCGCGCAGAGTGTAGAAGTTGTGCAGGGTATGCGTGGCGTTCTCGGCGGTACACAGGGCTTTCAGGGCCGGAATATAGAAGTGCATTTCGGCTGGCGCTTCGCTACCTGGGGTCATCAGGAAGGAAAACTACAGGCCGTCGATAATCATCTTCGTGCCAGTTC‑TGACGGTAGTTTTCGTCGGTGCCATCATGCGCCGGTCGCCCGTTTCCAGCGTc                                                                                                                                                                                                              >  2:6138/1‑251 (MQ=255)                                         ggTCCACTTGCTGGTGTCGCGGGTTTTCGCGCCGCGCAGAGTGTAGAAGTTGTGCAGGGTATGCGTGGCGTTCTCGGCGGTACACAGGGCTTTCAGGGCCGGAATATAGAAGTGCATTTCGGCTGGCGCTTCGCTACCTGGGGTCATCAGGAAGTCAAACTCCAGGCCGTCGATAATCATCTTCTCGCCAGTTC‑TGACGATAGTTTTCGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGc                                                                                                                                                                                                          >  1:351607/1‑251 (MQ=255)                                            ccACTAGCTGGTGACGCGGGGTTTCGAGCCGCGCAGAGTGTAGAAGTTGTGCAGGGTATGCGTGGCGTTCTCGGCGGCACACAGGGCTTTCCGGGCCGGAATATAGAAGTGCATTTCGGCTGGCGCTTCGCTACCTGGGGTCATCAGGAAGTCAAACTCCAGGCCGTCGATAATCATCTTCTCGCCAGTTC‑TGACGCTAGGTTTCGTCGGTGCAATAATGCTCGGGGCGCCCGTTGCCAGCGTCCCGCCa                                                                                                                                                                                                        >  2:386393/1‑250 (MQ=255)                                             cACTTGCTGGTGTCGCGGGTTTTCGCGCCGCGCAGAGTGTAGAAGTTGTGCAGGGTATGCGTGGCGTTCTCGGCGGTACACAGGGCTTTCAGGGCCGGAATATAGAAGTGCATTTCGGCTGGCGCTGCGCTACCTGGGGTCATCAGGAAGTCAAACTCCAGGCCGTCGATAATCATCTTCTCGCCAGTTC‑TGACTATAGTGTTCGTCGGTGCAATAATTCTCGGGTCGCCCGTTGCCAGCGTCACGCCAAg                                                                                                                                                                                                      <  1:37028/251‑1 (MQ=255)                                                     ggTGTCGCGGGTTTTCGCGCCGCGTAGAGTGTAGAAGTTGTGCAGGGTATGCGTGGCGTTCTCGGCGGTACACAGGGCTTTCAGGGCCGGAATATAGAAGTGCATTTCGGCTGGCGCTTCGCTACCTGGGGTCATCAGGAAGTCAAACTCCAGGCCGTCGATAATCATCTTCTCGCCAGTTC‑GGACGATAGTTTTAGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCAAGGCCATTAc                                                                                                                                                                                              <  1:357568/251‑1 (MQ=255)                                                            cgGGTTTTCGCGCCGCGTAGAGTGTAGAAGTTGTGCAGGGTATGCGTGGCGTTCTCGGCGGTACACAGGGCTTTCAGGGCCGGAATATAGAAGTGCATTTCGGCTGGCGCTTCGCTACCTGGGGTCATCAGGAAGTCAAACTCCAGGCCGTCGATAATCATCTTCTCGCCAGTTC‑GGACGATAGTTTTAGTCGGAGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCAAGGCCATTACCcacgtt                                                                                                                                                                                       >  1:291022/1‑248 (MQ=255)                                                            cgGGTTTTCGCGCCGCGTAGAGTGTAGAAGTTGGGCAGGGTATGCTTGGCGTTCCCGGCGGGACACAGGGCTTTCAGGGCCGGAATATAGAAGTGCATTTCGGCTGGCGCTTCGCTACCTGGGGTCATCAGGAAGTCAAACTCCAGGCCGTCGATAATCATCTTCGCGCCAGTTCGTG‑CGATAGTTTTATGCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGTGTCACGCCAAGGCCATTACCcacgtt                                                                                                                                                                                       <  2:74330/251‑4 (MQ=255)                                                             gggTTTTCGCGCCGCGCAGAGTGTAGAAGTTGTGCAGGGTATGCGTGGCGTTCTCGGCGGTACACAGGGCTTTCAGGGCCGGAATATAGAAGTGCATTTCGGCTGGCGCGTCGCTACCTGGGGTCATCAGGAAGTCAAACTCCAGGCCGTCGATAATCATCTTCTCGCCAGTTC‑TGACGATAGTTTTCGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCAAGGCCATTACCCACATTg                                                                                                                                                                                      <  2:413534/251‑1 (MQ=255)                                                               gTTTTCGCGCCGCGTAGAGTGTAGAAGTTGTGCAGGGTATGCGTGGCGTTCTCGGCGGTACACAGGGCTTTCAGGGCCGGAATATAGAAGTGCATTTCGGCTGGCGCTTCGCTACCTGGGGTCATCAGGAAGTCAAACTCCAGGCCGTCGATAATCATCTTCTCGCCAGTTC‑GGACGATAGTTTTAGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCAAGGCCATTACCCACGTTGcc                                                                                                                                                                                    >  1:285730/1‑251 (MQ=255)                                                                    cgcgccgcgtAGCGTGTAGAAGTTGTGCAGGGCATGCGTGGCGTTCTCGGCGGTCAACAGGGCTTTCAGGGCCGGAATATAGAAGTGCATTTCGGCTGGCGCTTCGCTACCTGGGGTCATCAGGAAGTCAAACTCCAGGCCGTCGATAATCATCTTCTCGCCAGTTC‑GGACGATAGTTTTAGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCAAGGCCATTACCCACGTTGCCTTgcg                                                                                                                                                                               >  2:160152/1‑251 (MQ=255)                                                                       gccgcgtAGAGTGTAGAAGTTGTGCAGGGTATGCGTGGCGTTCTCGGCGGTACACAGGGCTTTCAGGGCCGGAATATAGAAGTGCATTTCGGCTGGCGCTTCGCTACCTGGGGTCATCAGGAAGTCAAACTCCAGGCCGTCGATAATCATCTTCTCGCCAGTTC‑GGACGATAGTTTTAGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCAAGGCCATTACCCACGTTGCCTTGCGCGt                                                                                                                                                                            >  1:225825/1‑251 (MQ=255)                                                                         cgcgtaaagTGTAGAAGTTGTGCAGGGTATGCGTGGCGTTCTCGGCGGTACACAGGGCTTTCAGGGCCGGAATATAGAAGTGCATTTCGGCTGGCGCTTCGCTACCTGGGGTCATCAGGAAGGCAAAATCCAGGCCGTCGATAATCATCTTCTCGCCAGTTC‑GGACGATAGTTTTAGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCAAGGCCCTTACCCAAGTTGCCTTGCGCGTtg                                                                                                                                                                          >  2:320902/8‑251 (MQ=255)                                                                                     agaagttggtcggggTATGCGTGGCGTTCTCGGCGGTACACAGGGCTTTCAGGGCCGGAATATAGAATTTCATTTCTGCTGGCGCTTCGCTACCTGGGGTCATCAGGAAGTCAAACTCCAGGCCGGCGATAATCATCTTCTCGCCAGTTC‑GGACGATAGTTTTAGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCAAGGCCATTACCCACGTTGCCTTGCGCGTTGTGCGGCAGTAAc                                                                                                                                                              <  2:325155/239‑1 (MQ=255)                                                                                            gtgCAGGGTATGCGTGGCGTTCTCGGCGGTACACAGGGCTTTCAGGGCCGGAATATAGAAGTGCATTTCGGCTGGCGCTTCGCTACCTGGGGTCATCAGGAAGTCAAACTCCAGGCCGTCGATAATCATCTTCTCGCCAGTTC‑GGACGATAGTTTTAGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCAAGGCCATTACCCACGTTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGt                                                                                                                                                       >  1:221507/1‑251 (MQ=255)                                                                                             tgCAGGGTATGCGTGGCGTTCTCGGCGGTACACAGGGCTTTCAGGGCCGGAATATAGAAGTGCATTTCGGCTGGCGCTTCGCTACCTGGGGTCATCAGGAAGTCAAACTCCAGGCCGTCGATAATCATCTTCTCGCCAGTTC‑TGACGATAGTTTTCGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCAAGGCCATTACCCCCATTTCCTTGCGCGTTGTGCGGCAGTAACAGACCGt                                                                                                                                                       >  2:72630/1‑250 (MQ=255)                                                                                                     aTGCGTGGCGTTCTCGGCGGTACACAGGGCTTTCAGGGCCGGAATATAGAAGTGCATTTCGGCTGGCGCTTCGCTACCTGGGGTCATCAGGAAGTCAAACTCCAGGCCGTCGATAATCATCTTCTCGCCAGTTC‑GGACGATAGTTTTAGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCAAGGCCATTACCCACGTTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTaga                                                                                                                                              <  1:26483/251‑1 (MQ=255)                                                                                                        cGTGGCGTTCTCGGCGGTACACAGGGCTTTCAGGGCCGGAATCTAGAAGTGCATTTAGGCTGGCGCTTCGCTACCTGGGGTCATCAGGAAGTCAAACTCCAAGACGTCGATAATCATCTTCTCGCCTGTTC‑GGACGATAGTTTTAGTCGGTGCATTAATGCTCGGGTCGCACGTTGCCAGCGTCACGCCAAGGCCATTACCCACGTTGCCTTGCGCGTTGTGCGGAAGTAACCAACCGTAAGAGTagag                                                                                                                                             >  1:285846/1‑249 (MQ=255)                                                                                                              gTTCTCGGCGGAACACAGGGCTTTAAGGGCCGGAATCTAGAAGTGCATTTCGGCTGGCGCTTCGCTACATAGGGTCATCCGGCAGTCAAACTCCAGGCCGTCGATAATCATCTTCTCGCCAGTAC‑GGACGATCGTTTTAGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCAAGGCCATTAACAACGTTGCCTTGCGCGTTGTGCGGCAGTCACAGAACGTAAGAGAAAAGCGCAcggc                                                                                                                                     >  2:366730/1‑251 (MQ=255)                                                                                                                   cggcggTACACAGGGCTTTCAGGGCCGGAATATAGAAGTGCATTTCGGCTGGCGCTTCGCTACCTGGGGTCATCAGGAAGTCAAACTCCAGGCCGTCGATAATCATCTTCTCGCCAGTTC‑TGACGATAGTTTTCGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCAAGGCCATTACCCACATTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCt                                                                                                                                 <  1:368885/250‑1 (MQ=255)                                                                                                                             caGGGCTTTCAGGGCCGGAATATAGAAGTGCATTTCGGCTGGCGCTTCGCTACCTGGGGTCATCAGGAAGTCAAACTCCAGGCCGTCGATAATCATCTTCTCGCCAGTTC‑GGACGATAGTTTTAGTCGGTGCAATAATGCTCGTGTCGCCCGTTGCCAGCGTCACGCCAAGGCCATTACCCACGTTGCCTTGCGCGTGGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTAc                                                                                                                      <  1:28044/251‑1 (MQ=255)                                                                                                                              aGGGCTTTCAGGGCCGGAATATAGAAGTGCATTTCGGCTGGCGCTTCGCTACCTGGGGTCATCAGGAAGTCAAACTCCAGGCCGTCGATAATCATCTTCTCGCCAGTTC‑TGACGATAGTTTTCGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCAAGGCCATTACCCACATTGCCTTGCGCGTTGTGCGGCAGTAACAGACAGTACGAGTAGAGCGCACGGCACCTCATGCTGTTAcc                                                                                                                     >  2:308651/1‑251 (MQ=255)                                                                                                                                         ggCCGGAATATAGAAGTGCATTTCGGCTGGCGCTTCGCTACATGGGGTCATCAGGAAGTCAAACTCCAGGCAGTCGATAATCATCTTCTCGCCAGTTC‑TGACGATAGTTTTCGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCAAGGCCATTACCCACATTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGt                                                                                                          >  1:141410/1‑251 (MQ=255)                                                                                                                                               aaTATAGAAGTGCATTTCGGCTGGCGCTTCGCTACCTGGGGTCATCAGGAAGTCAAACTCCAGGCCGTCGATAATCATCTTCTCGCCAGTTC‑GGACGATAGTTTTAGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCAAGGCCATTACCCACGTTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGc                                                                                                    >  2:19495/1‑251 (MQ=255)                                                                                                                                                  ataGAAGTGCATTTCGGCTGGAGCTTCGCTACCTGGGGTCATCAGGAAGTCAAACTCCAGGCCGTCGATAATCATCTTCTCGCCAGTTC‑GGACGATAGTTTTAGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCAAGGCCATTACCCACGTTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGCGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTg                                                                                                  >  2:360340/1‑250 (MQ=255)                                                                                                                                                    aGAAGTGCATTTCGGCTGGCGCTTCGCTACCTGGGGTCATCAGGAAGTCAAACTCCAGGCCGTCGATAATCATCTTCTCGCCAGTTC‑GGACGATAGTTTTAGTCGGTGCAATACTGCTCGGGTCGCCCGTTGCCAGCGTCACGCCAAGGCCATTACCCACGTTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATg                                                                                               >  1:413673/1‑251 (MQ=255)                                                                                                                                                         tGCATTTCGGCTGGCGCTTCGCTACCTGGGGTCATCAGGAAGTCAAACTCCAGGCCGTCGATAATCATCTTCTCGCCAGTTC‑GGACGATCGTTTTAGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCAAGGCCATTACCCACGTTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATTATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTc                                                                                          >  1:184713/1‑251 (MQ=255)                                                                                                                                                         ggCATTTCGGCTGGCGCTTCGCTCCCTGGGGTCAGCAGGAAGTCAAACTCCAGCCCGTCGATAATCATCTTCTCGCCAGTGC‑GGACGATAGTTTTAGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCTCGCCAAGGCCATTACCCACGTTGCCTTGCGCTTTGTGCGGCAGAAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCAGGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTc                                                                                          <  1:363691/250‑1 (MQ=255)                                                                                                                                                            aTTTCGGCGGGCGCTTCGCTCCCGGGGGTCACCAGGAAGTCAACCTCCAGGCCGTCGATAATCATCTTCTCGCCAGTTC‑GGACGATAGTTTTAGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCCCGCCAAGGCCATTCCCCACGTTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTc                                                                                       <  1:328488/251‑1 (MQ=255)                                                                                                                                                             tttCGGCTGGCGCTTCGCTACCTGGGGTCATCAGGAAGTCAAACTCCAGGCCGTCGATAATCATCTTCTCGCCAGTTC‑TGACGATAGTTTTCGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCAAGGCCATTACCCACATTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTcc                                                                                      >  1:216856/1‑251 (MQ=255)                                                                                                                                                                   cTGGCGCTTCGCTACCTGGGGTCATCAGGAAGTCAAACTCCAGGCCGTCGATAATCATCTTCTCGCCAGTTC‑GGACGATAGTTTTAGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCAAGGCCATTACCCACGTTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATaa                                                                                  >  1:38780/1‑249 (MQ=255)                                                                                                                                                                           tCGCTCCCGGGGGCCATAAGAAAGTCAAACCCCAGGCCGTATATAATCTCCTTCTCGCAAGTTC‑GTAGGATAGTTTTAGTCGGTGCAAGAAGATTCGTGTCGTCCGTTGCCAGCTTCCCTCAAAGGCCATTCCCCCGGTTGCTTGGCGCGTTGTGCGGCATTAACAGACCGTAAGGTTAGAGCCCACGGCGGCCCATGATTTTTCCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCa                                                                          <  1:406889/249‑1 (MQ=255)                                                                                                                                                                               tACCTGGGGTCATCAGGAAGTCAAACTCCAGGCCGTCGCTAATCATCTTCTCGCCAGTTC‑GGACGATAGTTTTAGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCAAGGCCATTACCCACGTTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGc                                                                     >  2:371726/1‑250 (MQ=255)                                                                                                                                                                                       gTCATCAGGAAGTCAAACTCCAGGCCGTCGATAATCATCTTCTCGCCAGTTC‑GGACGATAGTTTTAGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCAAGGCCATTACCCACGTTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTAcc                                                                                                                     <  1:352792/194‑1 (MQ=255)                                                                                                                                                                                       gTCATCAGGAAGTCAAACTCCAGGCCGTCGATAATCATCTTCTCGCCAGTTA‑GGACGCTAGTTTTAGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCAAGGCCATTACCCACGTTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGAAGAGCGCACGGCGGCCCATGATGTTAcc                                                                                                                     >  2:352792/1‑194 (MQ=255)                                                                                                                                                                                            caGGAAGTCAAACTCCAGGCCGTCGATAATCATCTTCTCGCCAGTTC‑GGACGATAGTTTTAGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCAAGGCCATTACCCACGTTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACtt                                                       >  1:243343/1‑251 (MQ=255)                                                                                                                                                                                              ggAAGTCAAACTCCAGGCCGTCGATAATCATCTTCTCGCCAGTTC‑TGACGATAGTTTTCGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCAAGGCCATTACCCACATTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTg                                                     >  1:225690/1‑251 (MQ=255)                                                                                                                                                                                              ggAAGCCAAACTCCAGGCCGTCGATAATCTTCTTCTCCCCAGTTC‑TGACGATAGTTTTCGCCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCAAGGCCATTACCCACATTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTCCCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTg                                                     <  1:72630/251‑1 (MQ=255)                                                                                                                                                                                                aagtaaacccccAGGCCTTAGATATTCATCTTCTCGCCAGTTC‑TGACGATAGTTTTCTGCGGTGCAATAATGCCCGGGTCGCCCGTTGCCAGCGTCACGCCAAGGCCATTACCCACATTGCCTGTCGCGTGGTGCGGCAGTAACAGCCCGTAAGTGTAGAGCGCCGGGGGGCTCATGAGGTTCCCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGc                                                    <  2:351607/240‑1 (MQ=255)                                                                                                                                                                                                aaGTAAAACTCCAGGCCGTCGATAATCATCTTCCCGCCAGTTC‑GGACGATAGTTTTAGCCGGTGAAATAATGTTCGGGTCGCCCGTGGCCAGCGTCACGCAAAGGCCATTCCCCACGTTGCCTTGCGCGTTGGGCGGCAGTAACAGACCGTAAGAGTGGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCCCCTGAACTTTGcc                                                   <  1:372839/251‑1 (MQ=255)                                                                                                                                                                                                 aGTCACACTCCAGGCCGTCGATAATCATCTTCTCGCCAGTTC‑GGACGATAGTTTTAGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCAAGGCCATTACCCACGTTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCg                                                  >  2:416570/1‑251 (MQ=255)                                                                                                                                                                                                          ccAGGCCGTCGATAATCATCTTCTCGCCAGTTG‑GGACGATAGTTTTAGTCGTGGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCAAGGCCATTCCCCACGTTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCTTGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTaa                                           <  1:334667/249‑1 (MQ=255)                                                                                                                                                                                                          ccAGGCCGTCGATAATCATCTTCCCGCCAGTTC‑GGACGATAGTTTTAGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCAAGGCCATTACCCACGTTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACa                                         <  1:19495/251‑1 (MQ=255)                                                                                                                                                                                                             ggCCGCCGATAATCTTCTTCTCGCCATTTC‑GGACGATAGTTTTAGTCGGTGCAATAAGGCTGGGGTCGCCCGTTGCCAGCGCCACTCCAAGGCCATTACCCACGTTGCTTTGCGCGTGGTGCGGCAGTAACAGACCGTAAGAGTAGGGCGCACGGCGCCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCg                                      <  1:360340/251‑1 (MQ=255)                                                                                                                                                                                                                   cGATAATCATCTTCTCGCCAGTTC‑GGACGATAGTTTTAGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCAAGGCCATTACCCACGTTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGACTTc                                 >  1:204108/1‑250 (MQ=255)                                                                                                                                                                                                                    gATAATCATCTTCTAGCCAGTTC‑GGCCGATAGTTTTAGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCAAGGCCATTACCCACGTTGCATTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCCTGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCgaattc                                 >  2:51254/1‑244 (MQ=255)                                                                                                                                                                                                                     aTAATCATCATCTCGCCAGTTC‑TGACGATAGTTTTCGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCAAGGCCATTACCCACATTGCCTTGCGCGTTGTGCGGCAGTAACAGCCCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTc                                 >  1:132781/1‑248 (MQ=255)                                                                                                                                                                                                                            tcttctCGCCAGTTC‑GGACGATAGTTTTAGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCAAGGCCATTACCCACGTTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGACTTCAGAGATAATg                       >  1:383126/1‑251 (MQ=255)                                                                                                                                                                                                                                ctcGCAAGTTC‑TGACGATAGTTTTCGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCAAGGCCATTACCCACATTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCtt                   >  1:303289/1‑251 (MQ=255)                                                                                                                                                                                                                                             gACGATAGTTTTAGTCGGTGCAATAATGCTCGGGTCGACCGTTGCCAGCGTCACGCCAAGGCCATTACCCACGTTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACCTCGCCTTCAGAGATAATGCCTTTCACGCCACCAt       >  1:198363/1‑251 (MQ=255)                                                                                                                                                                                                                                              aCGATAGTTTTAGTCGGTGCAATAATGCTCGGGCCGCCCGTTGCCAGCGTCACGCCAAGGCCTTTACCCACGTTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCAGGGCGGCTCAGGTTTTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGACTTCAGAGATAATGCCTTTCACGCCACCATa      <  1:371726/251‑1 (MQ=255)                                                                                                                                                                                                                                               cGATAGTTTTAGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCAAGGCCATTACCCACGTTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACg                                                                                                           >  1:293460/1‑149 (MQ=255)                                                                                                                                                                                                                                                  tAGTTTTAGTCGGTGCAATAATGATCGGGTCGCCCGTTGCCAGAGTCACGACAAGACCATTACCCACGTTGCCTTACGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGCGCGCACGGCGGCACATGATGTTACCCGCAAGACCGGTTTCGCTGATAGCTTCGTCCCTAAAGACTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAAAATCGACTTCAGAGATAAAGCCTTTCACGgcaccagagagg  >  2:95362/1‑245 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                    CGTCGTTACCCCACATATCCAGCGTTTCGTTCAGATACTCGGTCCACTTGCTGGTGTCGCGGGTTTTCGCGCCGCGCAGAGTGTAGAAGTTGTGCAGGGTATGCGTGGCGTTCTCGGCGGTACACAGGGCTTTCAGGGCCGGAATATAGAAGTGCATTTCGGCTGGCGCTTCGCTACCTGGGGTCATCAGGAAGTCAAACTCCAGGCCGTCGATAATCATCTTCTCGCCAGTTC‑TGACGATAGTTTTCGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCAAGGCCATTACCCACATTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGG  >  NC_000913/4305556‑4306046

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