Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	581,416	0	A	G	20.0%	97.9 / 24.3	57	intergenic (+37/‑118)	nohD/tfaD	putative DNA‑packaging protein NohD/putative tail fiber assembly protein TfaD
Reads supporting (aligned to +/- strand):  ref base A (25/19);  new base G (3/8);  total (28/29)
Fisher's exact test for biased strand distribution p-value = 1.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.03e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTATTCTCGACGGGATCCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCA  >  NC_000913/581166‑581665                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                          cGCATTCTGTACGTTCGTGATGTCGCGGATCGCATGTGAAAATGCCAGTATTCTCGACGGGATCCCCCTGTCGGTGCAGCGGCGTGTTCCGGAACTGGAAAACCGACATGTTGATGTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATGTTGTCCGCGCCTGGCTTCGCTCa                                                                                                                                                                                                                                                           <  2:294056/251‑1 (MQ=17)   cATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTATTCTCGACGGGATCCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACt                                                                                                                                                                                                                                                         >  1:568884/1‑251 (MQ=255)                          cgGATCGCAGGTGAAATTGCCAGTATTCTCGACGGGATCCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGAccg                                                                                                                                                                                                                                  <  1:100642/251‑1 (MQ=255)                                       aaaTTGCCAGTATTCTCGACGGGATCCCCCTGCCGGTGCAGCGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGc                                                                                                                                                                                                                     <  1:383139/251‑1 (MQ=255)                                       aaaTTGCCAGTATTCTCGACGGGATCCCCCGGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGc                                                                                                                                                                                                                     <  1:408442/251‑1 (MQ=255)                                               aGTATTCTCGACGGGCTCCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCt                                                                                                                                                                                                              >  1:263574/1‑250 (MQ=255)                                                        gatggtcacccccTGGCGGTGCAGGGGCGTTTCTCGGAACTGGAAACCCGACATGTTTATTCCCTGAAACGGGATATCACCAAAGCCATGAACAAAGCAGCCGCGCTGGAGGAACTGATACCGGGGTTGCTGAGTAAATATACCGAACAGCCAGGTTAACAGGCGCCGGCATTTTGTCCGCGCCGGGTTTAGCTCCCTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTAtc                                                                                                                                                                                                    <  2:203248/243‑1 (MQ=255)                                                          cGGGATCCCCCTGTCGGTGCAGCGGCGTTTTCAGGCACTGGAAAACCGCCATGTTGATTTCCTGCAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTACCCGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCCCAGACCGCCGTTTGCTGGGCGGAGGCTAATTCCTAtctc                                                                                                                                                                                                  >  2:216443/1‑251 (MQ=21)                                                           gggATACCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGCGCCACAGACCGCCGGTGAATGGGCGGTTGCTAATTACTATCTcc                                                                                                                                                                                                 >  2:367745/1‑251 (MQ=255)                                                            ggATCCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGACACGGGATATCATCAAAGCCATGCACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACATGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCAACAGCCCGCCGTTGAATGGGCGGATGCTAATTACTATCTccc                                                                                                                                                                                                >  2:333200/1‑251 (MQ=255)                                                               tCCCCCCGTCGGAGCAGCGGCGTTTTCCGGCACTGGAAAACCGACATGTTGATTTCCTGCAACGGGATATCATCAAAGCCATGAACAACGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGCACCATCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCCCCGTTGAATGGGCGGATGGTAATTACTATCTCCCGa                                                                                                                                                                                              >  2:544884/1‑250 (MQ=21)                                                                        cGGTGCAGCGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGAATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGc                                                                                                                                                                                    >  2:24226/1‑251 (MQ=255)                                                                            gCAGCGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATAc                                                                                                                                                                                >  1:220469/1‑251 (MQ=255)                                                                               gcggcgTTTCCCGGAGCTGGAAAACCGTCATAGTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTGCAGCACTCTGCGCGCGCGGGTCTACTTTCGCTGTACCGACCTGAGCCGCAAACGGCGGTTGAATGGGCAGACGATAATTACTATCTCCCCAAAGAGTCGGCCTACCAg                                                                                                                                                                             <  2:384541/251‑1 (MQ=255)                                                                                 ggcgTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAgg                                                                                                                                                                            >  1:224444/1‑250 (MQ=255)                                                                                     ttttCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAggg                                                                                                                                                                       >  2:177609/1‑251 (MQ=255)                                                                                     ttttCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACACACCGCCGTTGAATGGGCGGATGCTATTTACTATCTCCCGCAAGAATCCGCATACCAGGAAggg                                                                                                                                                                       >  2:27996/1‑251 (MQ=255)                                                                                         ccGGAACTGGAAAACCGAAATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAAAAAAGAAGCCGCGCTGGATGAACTGATACCGGGGTGGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTGCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCt                                                                                                                                                                   <  2:133690/251‑1 (MQ=255)                                                                                         ccGGAACGGGAAAACCGAAAGGTTGATTTCTTGAAACGGGATATCATCAAAGCCAGGAACAAACCAGCCGCGCGGGATGATCTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGTTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTTTTCAGGCCGGAGCCACAGCCCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCt                                                                                                                                                                   <  2:251778/251‑1 (MQ=255)                                                                                          cGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTg                                                                                                                                                                  >  1:524596/1‑251 (MQ=255)                                                                                             aaCTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTggg                                                                                                                                                                >  1:76978/1‑250 (MQ=255)                                                                                                    aaaCCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTg                                                                                                                                                        >  1:324933/1‑251 (MQ=255)                                                                                                    aaaCCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTg                                                                                                                                                        >  1:434620/1‑251 (MQ=255)                                                                                                            aTGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCGCTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCa                                                                                                                                                <  1:244455/251‑1 (MQ=255)                                                                                                                      ccTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATACCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACACACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCcat                                                                                                                                       <  1:27996/250‑1 (MQ=255)                                                                                                                        tGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGc                                                                                                                                                                                                               <  1:39339/176‑1 (MQ=255)                                                                                                                        tGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGc                                                                                                                                                                                                               >  2:39339/1‑176 (MQ=255)                                                                                                                            gaaggacatcattAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTGCAGCACTCTGCGCGCGCGGGTCTACTTTCGCTGTACCGACCTGAGCCGCAAACGGCGGTTGAATGGGCAGACGATAATTACTATCTCCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAAt                                                                                                                                <  1:455931/244‑1 (MQ=255)                                                                                                                              gggATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGc                                                                                                                              >  1:539282/1‑251 (MQ=255)                                                                                                                                gacatcattAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTGCAGCACTCTGCGCGCGCGGGTCTACTTTCGCTGTACCGACCTGAGCCGCAAACGGCGGTTGAATGGGCAGACGATAATTACTATCTCCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCaaca                                                                                                                                             <  1:198669/231‑4 (MQ=255)                                                                                                                                    tcatcaAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACCGCCCTTTCAGCGGGCCATCATGAATGCGATggg                                                                                                                        >  2:278557/1‑251 (MQ=255)                                                                                                                                     catcaAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGc                                                                                                                       <  1:3855/251‑1 (MQ=255)                                                                                                                                        aaaaGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGc                                                                                                                    <  1:122953/250‑1 (MQ=255)                                                                                                                                           aGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGAc                                                                                                                 >  1:262506/1‑251 (MQ=255)                                                                                                                                              cATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTGCAGCACTCTGCGCGCGCGGGTCTACTTTCGCTGTACCGACCTGAGCCGCAAACGGCGGTTGAATGGGCAGACGATAATTACTATCTCCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGGAAAGTTGCCGTTTCAACGAGCGATCATGAATGCgtagggaaaagatta                                                                                                               >  1:63661/1‑236 (MQ=255)                                                                                                                                                 gAACAAAGCAGCTGCGCTGGATGAAATAATCCCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTGCAGCACTCTGCGCGCGCGGGTCTACTTTCGCTGTACCGACCTGAGCCGCAAACGGCGGTTGAATGGGCAGACGATAATTACTATCTCCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAAAGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACAt                                                                                                            >  1:134821/1‑250 (MQ=255)                                                                                                                                                    cAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGGTGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCgg                                                                                                        >  1:502153/1‑250 (MQ=255)                                                                                                                                                     aaaGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGCATGCGATGGGCAGCGACTACATCCGTg                                                                                                       >  1:276972/1‑251 (MQ=255)                                                                                                                                                              gccctggagtaaccgATACCGGGGTTTCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCCCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAAtg                                                                                              <  1:208742/237‑1 (MQ=255)                                                                                                                                                                gcTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTGCAGCACTCTGCGCGCGCGGGTCTACTTTCGCTGTACCGACCTGAGCCGCAAACGGCGGTTGAATGGGCAGACGATAATTACTATCTCCCCAAAGAGTCGGCCTACCAGGGAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGGAACGATTACATACGTGAGGtcaacgtt                                                                                            >  1:388898/1‑244 (MQ=255)                                                                                                                                                                gcTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTGCAGCACTCTGCGCGCGCGGGTCTACTTTCGCTGTACCGACCTGAGCCGCAAACGGCGGTTGAATGGGCAGACGATAATTACTATCTCCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGtcaacgtt                                                                                            <  1:35177/251‑8 (MQ=255)                                                                                                                                                                  tGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTGCAGCACTCTGCGCGCGCGGGTCTACTTTCGCTGTACCGACCTGAGCCGCAAACGGCGGTTGAATGGGCAGACGATAATTACTATCTCCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGtcaacgtt                                                                                            <  1:565931/249‑8 (MQ=255)                                                                                                                                                                            aaTACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTGCAGCACTCTGCGCGCGCGGGTCTACTTTCGCTGTACCGACCTGAGCCGCAAACGGCGGTTGAATGGGCAGACGATAATTACTATCTCCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGcc                                                                                <  1:534053/250‑1 (MQ=255)                                                                                                                                                                            aaTACCGGGGTTGCTGAGTGAATATATCGAACAGGCAGGTTAAGGGGCTGCAGCACTCTGCGCGCGCGGGTCTACTTTCGCTGTACCGACCTGAGCCGCAAACGGCGGTTGAATGGGCAGACGATAATTACTATCTCCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGcc                                                                                <  1:265898/250‑1 (MQ=255)                                                                                                                                                                                  ggggTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTc                                                                          >  1:80876/1‑251 (MQ=255)                                                                                                                                                                                  ggggTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGAGAATGTGGTGAAGTCTGCCCGTGTc                                                                          >  1:82900/1‑251 (MQ=255)                                                                                                                                                                                     gtagCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCGCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGGGAATGTGGTGAAGTCTGCCCGTGTCGGt                                                                       <  1:216443/248‑1 (MQ=255)                                                                                                                                                                                           gAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCGGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTcc                                                                 <  1:367745/251‑1 (MQ=255)                                                                                                                                                                                           gAGTGAATATAGCGAACAGTCAGGTTAACAGGTGGCGGATTTTTTTTCGCGCGGGCCTGCGCTCTCGGTCCAGGCCGTAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTAtta                                                                  <  2:134404/250‑2 (MQ=255)                                                                                                                                                                                              tGAATATATCGATCAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTCCAGGCCGGAGCCACAGGCCGCCGTTGAAGGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATAc                                                                                                                                                                                <  2:79792/137‑1 (MQ=255)                                                                                                                                                                                              tGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGCTGCTAATTACTATCTCCCGAAAGAATCCGCATAc                                                                                                                                                                                >  1:79792/1‑137 (MQ=255)                                                                                                                                                                                                                  gTTAACAGGCGGCGGCAGTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCGGGAGCCACAGGCCGACGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTGCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGTTGGGTGTTTATGc                                          <  2:276972/251‑1 (MQ=255)                                                                                                                                                                                                                     aaCAGGCGCCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTa                                       <  2:263574/251‑1 (MQ=255)                                                                                                                                                                                                                        aGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACtt                                    <  1:177609/251‑1 (MQ=255)                                                                                                                                                                                                                            tGCGGAATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGCAAGAATCCGCATACCAGGAAGGGAGCTGGGAAACACTGCCCTTTCAGCGGGCAATCATGAATGCGATGGGAAGCGCCTACATCCGTGAGGCGAATGTGGTGCAGCCTGCACGCGTCGGTTATTCCAAATGGCTGCTGGGGGTTTGTGCCtcttttttt                                >  2:515152/1‑243 (MQ=255)                                                                                                                                                                                                                             gctgatgttggcccGCGCTGGGCTCCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATag                               <  2:539282/239‑1 (MQ=255)                                                                                                                                                                                                                                         ccGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAagcagc                   <  1:488693/251‑1 (MQ=255)                                                                                                                                                                                                                                                          cACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCa  <  1:514181/251‑1 (MQ=255)                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                          CGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTATTCTCGACGGGATCCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCA  >  NC_000913/581166‑581665

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