Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,429,225	0	C	T	53.6%	9.2 / 92.2	70	S59S (AGC→AGT)	stfR	putative prophage side tail fiber protein StfR
Reads supporting (aligned to +/- strand):  ref base C (21/11);  new base T (26/11);  total (48/22)
Fisher's exact test for biased strand distribution p-value = 7.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.68e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGTGGCACAGAACACGGCAGCCGCGAAGAAGTCAGCCAGTGATGCCAGCACATCAGCCCGTGAGGCGGC  >  NC_000913/1428982‑1429473                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                         tGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGGTACAGCATGGACGTTGAGAACGGTCAGTACAGCGTTATTc                                                                                                                                                                                                                                                   >  2:255003/1‑251 (MQ=255)                        gaTGACAGCCCGCCGGTTCAGGCGGGGTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACAGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGTTGAACACTCTGGCCTCAGAAAATCCGGAGGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCccgcc                                                                                                                                                                                                                             <  1:323807/249‑1 (MQ=255)                             cAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGCTTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAACCGTAACAACACCCCGGTGGTGGTGAAAACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTcac                                                                                                                                                                                                                        >  1:323834/1‑250 (MQ=255)                                      ggTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGAc                                                                                                                                                                                                              >  1:86022/1‑251 (MQ=255)                                            ggCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTGCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACATAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGGGGTGGGGAACAGGCTGGCCTCAGAAAATCGGGATGAAGCCGGGCTTTATAGCATGGACGTTGAGTACGGTCATTACAGCGTTACTCTGTTGTTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTAc                                                                                                                                                                                                        <  1:212371/251‑1 (MQ=255)                                                 gCTTTTTTGTGGGGTAAATATGGCAGTAAAGATTTCAGGTTAACGGAAAGACGGCACAGGAAAACCGGTACAGAACGGCAAAATCCAGCTGAAAGCAAAACGTAACAGCCCCACGGTTGTGGTGAACACGCTGGCCTCATAAAATCCGGATGAAGCCGGGCGTTACATCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCgtgt                                                                                                                                                                                                   <  1:24294/251‑1 (MQ=255)                                                         gtgGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGAt                                                                                                                                                                                           <  1:8352/251‑1 (MQ=255)                                                           ggggTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAACCCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTACCAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATtc                                                                                                                                                                                         <  2:63714/251‑1 (MQ=255)                                                              gTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATtcac                                                                                                                                                                                       >  1:190936/1‑248 (MQ=255)                                                                  atatGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGCACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAAccg                                                                                                                                                                                  >  2:177041/1‑250 (MQ=255)                                                                    atGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAAccggg                                                                                                                                                                                <  1:119242/251‑4 (MQ=255)                                                                        cAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACg                                                                                                                                                                            >  2:31399/1‑251 (MQ=255)                                                                          gTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCt                                                                                                                                                                          <  1:19093/251‑1 (MQ=255)                                                                               gATTTCAGGTGTACTGAAAGACGGCACAGGAACACCGGTACAGTACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGGGGAAGGATTACCGCCGTCACATGCCGGGCCCATTACCGTGTATGAAGTTTCTCAACCCGGTACGCTGAATg                                                                                                                                                                     >  2:259220/1‑251 (MQ=255)                                                                                        agtccTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCGGCAGGTTGACGGTTTTCGACCATCGCACGCCGGGACCATCACCGTGTATGCAGATTCACAACAGGGGACGCTGAATGCTTTtctct                                                                                                                                                            >  2:261102/5‑250 (MQ=255)                                                                                           ccTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTtctctgt                                                                                                                                                          <  1:268630/249‑4 (MQ=255)                                                                                             tGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGcc                                                                                                                                                       >  2:281113/1‑251 (MQ=255)                                                                                                     ggCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACgg                                                                                                                                                >  1:336756/1‑250 (MQ=255)                                                                                                      gCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACgg                                                                                                                                                >  2:236010/1‑249 (MQ=255)                                                                                                         caGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGGGgat                                                                                                                                           >  2:192830/1‑251 (MQ=255)                                                                                                             aaaaCCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGatgatg                                                                                                                                       >  2:8547/1‑251 (MQ=255)                                                                                                                 ccGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCg                                                                                                                                   >  1:168910/1‑251 (MQ=255)                                                                                                                 acGGTACAGAACTGCACCATTCAGCTGAAAGCAAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGAAGGTTGACGGTTTTCCAACATCGCACGCCGGGACCATCACCGTGTATGAAGTTTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCg                                                                                                                                   >  1:168911/2‑251 (MQ=255)                                                                                                                  cGGTACAGAACTGCACAATCAAGCTGAAAGCAAAACGAAACCGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGACGCAGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATTACGGAGGATGATGCCCGt                                                                                                                                  >  2:181632/1‑251 (MQ=255)                                                                                                                     tACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCg                                                                                                                               >  1:252208/1‑251 (MQ=255)                                                                                                                       cAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCccggccgg                                                                                                                              >  1:8018/1‑250 (MQ=255)                                                                                                                         gAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGAgg                                                                                                                           >  1:321511/1‑251 (MQ=255)                                                                                                                         gAACTGCACAATCCAGCTGAAAGCAAAACGTAAAAGCACCACGGGGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAgg                                                                                                                           <  2:320419/251‑1 (MQ=255)                                                                                                                          aaCTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAAc                                                                                                                                                                                    >  2:351369/1‑193 (MQ=255)                                                                                                                          aaCTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAAc                                                                                                                                                                                    <  1:351369/193‑1 (MQ=255)                                                                                                                          aaCTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGAggt                                                                                                                          <  2:86532/251‑2 (MQ=255)                                                                                                                              gCACAACCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGAGTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTg                                                                                                                      <  1:157979/251‑1 (MQ=255)                                                                                                                                accattcAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGAggtgctgcg                                                                                                                    >  1:115521/7‑244 (MQ=255)                                                                                                                                   attcAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACACTCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGACTGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGAggtgctgcgtcg                                                                                                                 >  1:129026/4‑241 (MQ=255)                                                                                                                                   attcAGCTGAAAGCCAGACGTAACAGCACCACGGTAGAGGTGAACACGGTGGGATCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGAGGAGTACGGTCAGTACAGTGTCATCCTGCAGGATGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGAggtgctgcgtcg                                                                                                                 >  2:252336/4‑241 (MQ=255)                                                                                                                                   aTCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATTCCGGAGGATGATGCCCGTCCGGAGGCACTGCGCCg                                                                                                                 >  2:192412/1‑251 (MQ=255)                                                                                                                                       aGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGtttt                                                                                                             >  1:164849/1‑251 (MQ=255)                                                                                                                                          tGAAAGCCAGACGTAACAGCACAACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCAACGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGAGGTGCTGCGTCGTCTTGa                                                                                                           >  1:129556/1‑250 (MQ=255)                                                                                                                                           gAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCg                                                                                                                                                                                 <  2:296333/179‑1 (MQ=255)                                                                                                                                           gAAAGCAAAACGTAACAGCACCACGGTGGTGGCGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCg                                                                                                                                                                                 >  1:296333/1‑179 (MQ=255)                                                                                                                                                       tAACAGCACCACGGTGGGGGTGAACACGGTGGTCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGCCGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGCTGCCCGGCCGGAGGTGCTGCGTCGTCTTGAACTGATGGTGGa                                                                                               >  2:91613/1‑249 (MQ=255)                                                                                                                                                               ccacGGGGGGGGTGAACCGGCTGGCCTCAGAAAATCCGGAGGAAGCGTGGCGTTACAGCAGGGACGTTGATTACGGTCAGTACAGCGTTATTCGGTGGGGGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGgcg                                                                                     <  2:318751/251‑1 (MQ=255)                                                                                                                                                                  cGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGTCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGAGGTGCTGCGTCGTCTTGAACTGCTGGTGGAAGAGGTGGCGCGt                                                                                  >  2:132350/1‑251 (MQ=255)                                                                                                                                                                   ggtggtggtGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGACGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGt                                                                                  >  2:76729/1‑250 (MQ=255)                                                                                                                                                                        tggtgAACACGCTGGCATAAGACAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGGATGAAGATTCTCAACCCGGCACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGAGGAAGAGGTGGCGCGTAAagcg                                                                            >  2:215578/1‑251 (MQ=255)                                                                                                                                                                         ggtgAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGAGGTGCTGCGTCGTCTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGt                                                                           <  1:281113/251‑1 (MQ=255)                                                                                                                                                                            gactacGCTGGCGTGAGAAAATCCGGAGGAGGCGGGGCGTTACAGCATGGAGGTGGAGTACGCTCAGTACAGCGTTATGCTGTTGGTGGAAGGATTCCCGCCGGCACATGCCGGGACCATTACCGTGTATGAAGATTCGCAACCCGTTACGCTGAATGATTTTCTCGGTGCCATGACGGAGTATGAGGCCCGGCCGGAGGGACTGCGCCGGTTTGAACTGATGGTGGAAGAGGGGGCGAGTAACGCGTcc                                                                         <  2:20417/246‑1 (MQ=255)                                                                                                                                                                             aaCACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGTATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGAGGTGCTGCGTCGTCTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCcgt                                                                       <  1:233291/251‑2 (MQ=255)                                                                                                                                                                             aaCACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGAGGTGCTGCGTCGGCTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCcgt                                                                       <  2:190936/251‑2 (MQ=255)                                                                                                                                                                                 cggtgggcacAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGAACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGCGGATGATGCCCGGCCAGCGGTGCTGCGTCGTCTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCgtggtg                                                                   >  2:105976/10‑251 (MQ=255)                                                                                                                                                                                   gtgggcTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGACATCCTGCAGGATGACGGTTTTCCACCATCGCACGACGGGACCATCACCGTGTATGAAGATTCACAAACGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGAAGGTGCTGCGTCGTCTTGACCTGTTGGTGGAAGAGGTGGCGCGTAACGCGTCcgtgggggc                                                                 >  2:190062/6‑244 (MQ=255)                                                                                                                                                                                   cTGGCATCAGCAAATTCGGGTGAAGCCGGGCCTTACAGCATGGACGCTGAGCACGGTCAGTACCGCGTCATTCAGTTCGTGGAAGGATTCCCGCCGTCACATGCCGGGAACATTACCGTGTATGAAGATTCTCACCCCAGTACGCTGAATGATTTTCTCGGCGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGAAACGCGTCCGCGGTGGc                                                                 >  2:55603/1‑251 (MQ=255)                                                                                                                                                                                    tGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGTGGc                                                                 >  2:8412/1‑250 (MQ=255)                                                                                                                                                                                             aaaaTCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCAAATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGCATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGAGGAAGAGGTGGCGCGTAACGCGTCCGCGGTGGCACAGAACACg                                                       >  1:148090/1‑251 (MQ=255)                                                                                                                                                                                                   cGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGTGGCACAGAACACGGCAGcc                                                 >  1:190553/1‑251 (MQ=255)                                                                                                                                                                                                   cGGATGAAGACGGGCGTTACAGCATGGATGTGGAGAACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCAACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGACTGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGAGGTGCTGCGTCGTCTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGTGGTGGCACAGAGTACGGCAgac                                                 >  2:176546/1‑249 (MQ=255)                                                                                                                                                                                                    ggaggAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGAGGTGCTGCGTCGTCTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGTGGTGGCACAGAGTACGGCAgacg                                                <  2:276128/247‑4 (MQ=255)                                                                                                                                                                                                    ggATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGAGGTGCTGCGTCGTCTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGTGGTGGCACAGAGTACGGCAgacg                                                >  1:106096/1‑248 (MQ=255)                                                                                                                                                                                                             cGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACTGTTTTACACCATCGCACGCCGGGACCATCACCTTTTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGAGGTGCTGCGTCGTCTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTACGTGGTGGCACAGAGTACGGCATACGCgaagaaat                                       <  1:287257/251‑3 (MQ=255)                                                                                                                                                                                                                gCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGAGGTGCTGCGTCGTCTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGTGGTGGCACAGAGTACGGCAGACGCGAAGAAATCAg                                    <  2:245726/251‑1 (MQ=255)                                                                                                                                                                                                                    tACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGAGGTGCTGCGTCGTCTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGTGGTGGCACAGAGTACGGCAGACGCGAAGAAATCAGccgg                                >  1:203506/1‑249 (MQ=255)                                                                                                                                                                                                                          aTGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTCTGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGAGGTGCTGCGTCGTCTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGTGGTGGCACAGAGTACGGCAGACGCGAAGAAATCAGCCGGCGATGc                          >  2:31913/1‑251 (MQ=255)                                                                                                                                                                                                                           tggatgtggAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGAGGTGCTGCGTCGTCTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGTGGTGGCACAGAGTACGGCAGACGCGAAGCAATCAGCCGGCGATGcc                         >  2:318646/9‑251 (MQ=255)                                                                                                                                                                                                                             gatgtggAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGAGGTGCTGCGTCGTCTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGTGGTGGCACAGAGTACGGCAGACGCGCAGAAATCAGCCGGCGATGcc                         >  1:326441/7‑249 (MQ=255)                                                                                                                                                                                                                                    aGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGAGGTGCTGCGTCGTCTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGTGGTGGCACAGAGTACGGCAGACGCGAAGAAATCAGCCGGCGATGCCAgtgcatc                 >  2:335518/1‑244 (MQ=255)                                                                                                                                                                                                                                     gTACGGTCAGTCCAGAGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGCCCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGAGGACGAGGTGGCGCGTACCGCGTCCGCGGTGGCACAGAACACGGCAGCCGCGAAGCAGTCAGCCAGTGATGCCAGCACATc                 >  2:216535/1‑249 (MQ=255)                                                                                                                                                                                                                                         ggTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGGGGCACAGAACACGGCAGCCGCGAAGAAGTCAGCCAGGGATGCCAGCACATCAGCCCg           >  2:198789/1‑251 (MQ=255)                                                                                                                                                                                                                                            cAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGTGGCACAGAACACGGCAGCCGCGAAGAAGTCAGCCAGTGATGCCAGCACATCAGCCCGTGa        <  1:291443/251‑1 (MQ=255)                                                                                                                                                                                                                                                aCAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGCCGGCGGATGATGCCCGTCCGGCGGCACTGCGCCGTTTTGAACTGATGGTGGCAGAGGTGGCGCGTAACGCGTCCGCGGTGGCACAAAACACGGCAGCCGCGAAGAAGTCAGCCAGTGATGCCAGCACATCAACCCGTGAggcg    >  1:344114/1‑251 (MQ=255)                                                                                                                                                                                                                                                  aGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCCCTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGTGGCACAGAACACGGCAGCCGCGCAGAAGTCAGCCAGTGATGCCAGCACATCAACCCGTGAggcggc  >  1:113556/1‑251 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                         TGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGTGGCACAGAACACGGCAGCCGCGAAGAAGTCAGCCAGTGATGCCAGCACATCAGCCCGTGAGGCGGC  >  NC_000913/1428982‑1429473

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