Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,252,101	0	T	C	28.1%	107.6 / 48.6	65	intergenic (+1627/‑405)	malM/ubiC	maltose regulon periplasmic protein/chorismate lyase
Reads supporting (aligned to +/- strand):  ref base T (27/19);  new base C (10/8);  total (37/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TCCAACTTATCGATCAGGCCTTACGATCAACGGATGATATGATCAAAGCAAATGTTTGGCATCTTTATAAAGAGTGGATTCGTAGTGATGATGTTTCACCTATATTTATAGAAACTGAAGATAATTTAAGAACCTTTAACACGAATGAATTAACACGAAACGATAATATCTTTATCCTGTTCTCCTCAGTCGATGATGGGCCAGTTATGGTGGTAAGCTCCCAGCGCTTACATGATATGTTGAATCCTACAAAAGATACCAATTGGAATTCCACGTATATCTACAAATCCAGACATGAGATGTTGCCTGTTAATCTTACTCAGGAAACACTTTTCAGCTCCAAATCTCATGGTAAATATGCGCTTTTCCCCATTTTTACTGCGAGTTGGCGAGCTCATCGTATAATGAATAAGGGTGTTTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGC  >  NC_000913/4251862‑4252349                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                         tCCAACTTATCGATCACGCCTTACGATCAACGGATGATATGATCAAAGCAAATGTTTGGCATCTTTATAAAGAGTGGATTCGTAGTGATGATGTTTCACCTATATTTATAGAAACTGAAGATAATTTAAGAACCTTTAACACGAATGAATTAACACGAAACGATAATATCTTTATCCTGTTCTCCTCAGTCGATGATGGGCCAGTTATGGTGGTAAGCTCCCAGCGCTTACATGATATGCTGAATCCTAc                                                                                                                                                                                                                                                >  1:148350/1‑250 (MQ=255)    aaCTTATCGATCAGGCCTTACGATCAACGGATGATATGATCAAAGCAAATGTTTGGCATCTTTATAAAGAGTGGATTCGTAGTGATGATGTTTCACCTATATTTATAGAAACTGAAGATAATTTAAGAACCTTTAACACGAATGAATTAACACGAAACGATAATATCTTTATCCTGTTCTCCTCAGTCGATGATGGGCCAGTTATGGTGGTAAGCTCCCAGCGCTTACATGATATGCTGAATCCTACaaaa                                                                                                                                                                                                                                            <  2:198369/251‑1 (MQ=255)        tATCGATCAGGCCTTACGATCAACGGATGATATGATCAAAGCAAATGTTTGGCATCTTTATAAAGAGTGGATTCGTAGTGATGATGTTTCACCTATATTTATAGAAACTGAAGATAATTTAAGAACCTTTAACACGAATGAATTAACACGAAACGATAATATCTTTATCCTGTTCTCCTCAGTCGATGATGGGCCAGTTATGGTGGTAAGCTCCCAGCGCTTACATGATATGTTGAATCCTACAAAAGAt                                                                                                                                                                                                                                         <  2:314790/250‑1 (MQ=255)        tATCGATCAGGCCTTACGATCAACGGATGATATGATCAAAGCAAATGTTTGGCATCTTTATAAAGAGTGGATTCGTAGTGATGATGTTTCACCTATATTTATAGAAACTGAAGATAATTTAAGAACCTTTAACACGAATGAATTAACACGAAACGATAATATCTTTATCCTGTTCTCCTCAGTCGATGATGGGCCAGTTATGGTGGTAAGCTCCCAGCGCTTACATGATATGTTGAATCCTACAAAAGAt                                                                                                                                                                                                                                         >  1:314790/1‑250 (MQ=255)         aTCGATTAGGCCTTACGATCAACGGATGATATGATCAAAGCAAATGTTTGGCATCTTTATAAAGAGTGGCTTCGTAGAGATGATGTTTCACCTATATTTATAGAAACTGAAGATAATTTAAGAACCTTTAACACGAATGAATTAACACGAAACGATAATATCTTTATCCTGTTCTCCTCAGTCGATGATGGGCCAGTTATGGTGGTAAGCTCCCAGCGCTTCCATGATATGCTGCATCCTACAAAAGATAc                                                                                                                                                                                                                                       >  2:238280/1‑251 (MQ=255)          tCGATCAGGCCTTACGATCAACGGATGATATGATCAAAGCAAATGTTTGGCATCTTTATAAAGAGTGGATTCGTAGTGATGATGTTTCACCTATATTTATAGAAACTGAAGATAATTTAAGAACCTTTAACACGAATGAATTAACACGAAACGATAATATCTTTATCCTGTTCTCCTCAGTCGATGATGGGCCAGTTATGGTGGTAAGCTCCCAGCGCTTACATGATCTGTTGAATCCTACAAAAGATAcc                                                                                                                                                                                                                                      >  1:430080/1‑251 (MQ=255)                       aCGATCAACGGATGATATGATCAAAGCAAATGTTTGGCATCTTTATAAAGAGTGGATTCGTAGTGATGATTTTTCTCCTATATTTATAGAAACTGAAGATAATTTAAGAACCTTTAACACGAATGAATTAACACGAAACGATAATATCTTTATCCTGTTCTCCTCAGTCGATGATGGGCCAGTTATGGTGGTAAGCTCCCAGCGCTTACATTATATGTTGAAGCCTACAAAAGAGACCAATGGGAATTcc                                                                                                                                                                                                                          <  1:569161/250‑1 (MQ=255)                        cGATCAACGGATGATATGATCAACGCAAATGTTTGGCATCTTTATAAAGAGTGGATTCGTAGTGATGATGTTTCACCTAGATTTATAGAAACTGAAGATAATTTAAGAACCTTTAACACTAATGAATTAACACGAAACGATAATATCTTTATCCTGTTCTCCTCAGTCGATGATGGGCCAGTTATGGTGGTAAGCTCTCAGCGCTTACATGATATGTTGAATCCTACAAAAGATACCAATTGGAATTCCAc                                                                                                                                                                                                                        <  2:490458/251‑1 (MQ=255)                          aTCAACGGATGCTATAATCAAAGCAAATGTTTGGCATCTTTATAAAGAGTGGATTCGTAGTGATGATGTTTCACCTATATTTATAGAAACTGAAGATAATTTAAGAACCTTTAACACGAATGAATTAACACGAAACGATAATATCTTTATCCTGTTCTCCTCAGTCGATGATGGGCCAGTTATGGTGGTAAGCTCCCAGCGCTTACATGATATGCTGAATCCTACAAAAGATACCAATTGGAATTCCACGt                                                                                                                                                                                                                      >  2:46219/1‑251 (MQ=255)                             aaCGGATGATATGATCAAAGCAAATGTTTGGCATCTTTATAAAGAGTGGATTCGTAGTGATGATGTTTCACCTATATTTATAGAAACTGAAGATAATTTAAGAACCTTTAACACGAATGAATTAACACGAAACGATAATATCTTTATCCTGTTCTCCTCAGTCGATGATGGGCCAGTTATGGTGGTAAGCTCCCAGCGCTTACATGATATGCTGAATCCTACAAAAGATACCAATTGGAATTCCACGtata                                                                                                                                                                                                                   <  1:478980/251‑1 (MQ=255)                                     atatgatCAAAGCAAATGTTTGGCATCTTTATAAAGAGTGGATTCGTAGTGATGATGTTTCACCTATATTTATAGAAACTGAAGATAATTTAAGAACCTTTAACACGAATGAATTAACACGAAACGATAATATCTTTATCCTGTTCTCCTCAGTCGATGATGGGCCAGTTATGGTGGTAAGCTCCCAGCGCTTACATGATATGTTGAATCCTACAAAAGATACCAATTGGAATTCCACGTATATCTACaaa                                                                                                                                                                                                           <  1:581516/251‑1 (MQ=255)                                     atatgatCAAAGCAAATGTTTGGCATCTTTATAAAGAGTGGATTCGTAGAGATGATGTTTCACCTATATTTATAGAAACTGACGATAATTTAAGAACCTTTAACACGACTGAATTAACACGAAACGATAATATCTTTATCCTGTTCTCCTCAGTCGATGATGGGCCAGTTATGGTGGTAAGCTCCCAGCGCTTACATGATATGTTGAATCCTACAAAAGATACCAATTGGAATTCCACGTATATCTACaa                                                                                                                                                                                                            >  1:386583/1‑250 (MQ=255)                                            cAAAGCAAATGTTTGGCATCTTTATAAAGAGTGGATTCGTAGTGATGATGTTTCACCTATATTTATAGAAACTGAAGATAATTTAAGAACCTTTAACACGAATGAATTAACACGAAACGATAATATCTTTATCCTGTTCTCCTCAGTCGATGATGGGCCAGTTATGGTGGTAAGCTCCCAGCGCTTACATGATATGTTGAATCCTACAAAAGATACCAATTGGAATTCCACGTATATCTACAAATCCAGAc                                                                                                                                                                                                    <  2:504374/251‑1 (MQ=255)                                             aaaGCAAATGTTTGGCATCTTTATAAAGAGTGGATTCGTAGTGATGATGTTTCACCTATATTTATAGAAACTGAAGATAATTTAAGAACCTTTAACACGAATGAATTAACACGAAACGATAATATCTTTATCCTGTTCTCCTCAGTCGATGATGGGCCAGTTATGGTGGTAAGCTCCCAGCGCTTACATGATATGCTGAATCCTACAAAAGATACCAATTGGAATTCCACGTATATCTACAAATCCAGAc                                                                                                                                                                                                    >  1:47221/1‑250 (MQ=255)                                                  aaaTGTTTGGCATCTTTATAAAGAGTGGATTCGTAGTGATGATGTTTCACCTATATTTATAGAAACTGAAGATAATTTAAGAACCTTTAACACGAATGAATTAACACGAAACGATAATATCTTTATCCTGTTCTCCTCAGTCGATGATGGGCCAGTTATGGTGGTAAGCTCCCAGCGCTTACATGATATGCTGAATCCTACAAAAGATACCAATTGGAATTCCACGTATATCTACAAATCCAGACATg                                                                                                                                                                                                 >  2:629935/1‑248 (MQ=255)                                                   aaTGTTTGGCATCTTTATAAAGAGTGGATTCGTAGTGATGATGTTTCACCTATATTTATAGAAACTGCAGATAATTTAAGAACCTTTAACACGAATGAATTAACACGAAACGATAATATCTTTATCCTGTTCTCCTCAGTCGATGATGGGCCAGTTATGGTGGTAGGCTCCCAGCGCTTACATGATATGTTGAATCCTACAAAAGATACCCATTGGAATTCAACGTATATCTACAAATCCAGACATgatat                                                                                                                                                                                             >  2:39532/1‑248 (MQ=255)                                                      gTTTGGCATCTTTATAAAGAGCGGATTCGTAGTGATGATGTTTCACCTATATTTCTAGAAACTGAAGATAATTTAAGAACCTTTAACACGAATGAATTAACACGAAACGATAATATCTTTATCCTGTTCTCCTCAGTCGATGATGGGCCAGTTATGGTGGTAAGCTCCCAGCGCTTACATGATATGTTGAATCCTACAAAAGATACCAATTGGAATTCCACGTATATCTACAAATCCAGACATGAGATGtt                                                                                                                                                                                          >  2:331368/1‑251 (MQ=255)                                                                 ttATAAAGAGTGGATTCGTAGTGATGATGTTTCACCTATATTTATAGAAACTGAAGATAATTTAAGAACCTTTAACACGAATGAATTAACACGAAACGATAATATCTTTATCCTGTTCTCCTCAGTCGATGATGGGCCAGTTATGGTGGTAAGCTCCCAGCGCTTACATGATATGTTGAATCCTACAAAAGATACCAATTGGAATTCCACGTATATCTACAAATCCAGACATGAGATGTTGCCTGTTa                                                                                                                                                                                  <  1:311833/248‑1 (MQ=255)                                                                   ataAAGAGTGGATTCGTAGTGATGATGTTTCACCTATATTTATAGAAACTGAAGATAATTTAAGAACCTTTAACACGAATGAATTAACACGAAACGATAATATCTTTATCCTGTTCTCCTCAGTCGATGATGGGCCAGTTATGGTGGTAAGCTCCCAGCGCTTACATGATATGCTGAATCCTACAAAAGATACCAATTGGAATTCCACGTATAGCTACAAATCCAGACATGAGATGTTGCCTGTTAATCtt                                                                                                                                                                             <  1:46219/251‑1 (MQ=255)                                                                     aaaGAGTGGATTCGTAGTGATGATGTTTCACCTATATTTATAGAAACTGAAGATAATTTAAGAACCTTTAACACGAATGAATTAACACGAAACGATAATATCTTTATCCTGTTCTCCTCAGTCGATGATGGGCCAGTTATGGTGGTAAGCTCCCAGCGCTTACATGATATGTTGAATCCTACAAAAGATACCAATTGGAATTCCACGTATATCTACAAATCCAGACATGAGATGTTGCCTGTTAATCTTAc                                                                                                                                                                           <  2:520443/251‑1 (MQ=255)                                                                                       gatgatGTTTCACCTATATTTATAGACACTGAAGATAATTTAAGAACCTTTAACACGAATGAATTAACACGAAACGATAATATCTTTATACTGTTCTCCTCAGTCGATGATGGGCCAGTTATGGTGGTAAGCTCCCAGCGCTTACATGATATGTTGAATCCTACAAAAAATACCAATTGGAATTCCACGTATATCTACACATCCAGACATGAGATGTTGCCTGTTAATCTTACTCAGGAAACACTTTTc                                                                                                                                                           >  2:578304/1‑249 (MQ=255)                                                                                        atgatgTTTCACCTATATTTATAGAAACTGAAGATAATTTAAGAACCTTTAACACGAATGAATTAACACGAAACGATAATATCTTTATCCTGTTCTCCTCAGTCGATGATGGGCCAGTTATGGTGGTAAGCTCCCAGCGCTTACATGATATGTTGAATCCTACAAAAGATACCAATTGGAATTCCACGTATATCTACAAATCCAGACATGAGATGTTGCCTGTTAATCTTACTCAGGAAACACTTTTCAGc                                                                                                                                                        <  2:26134/251‑1 (MQ=255)                                                                                           atgTTTCACCTATATTTATAGAAACTGAAGATAATTTAAGAACCTTTAACACGAATGAATTAACACGAAACGATAATATCTTTATCCTGTTCTCCTCAGTCGATGATGGGCCAGTTATGGTGGTAAGCTCCCAGCGCTTACATGATATGTTGAATCCTACAAAAGATACCAATTGGAATTCCACGTATATCTACAAATCCAGACATGAGATGTTGCCTGTTAATCTTACTCAGGAAACACTTTTCAGCTcc                                                                                                                                                     <  1:234260/251‑1 (MQ=255)                                                                                           atgTTTCACCTATATTTATAGAAACTGAAGATAATTTAAGAACCTTTAACACGAATGAATTAACACGAAACGATAATATCTTTATCCTGTTCTCCTCAGTCGATGATGGGCCAGTTATGGTGGTAAGCTCCCAGCGCTTACATGATATGTTGAATCCTACAAAAGATACCAATTGGAATTCCACGTATATCTACAAATCCAGACATGAGATGTTGCCTGTTAATCTTACTCAGGAAACACTTTTCAGCTcc                                                                                                                                                     >  2:234260/1‑251 (MQ=255)                                                                                                  aCCTATATTTATAGAAACTGAAGATAATTTAAGAACCTTTAACACGAATGAATTAACACGAAACGATAATATCTTTATCCTGTTCTCCTCAGTCGATGATGGGCCAGTTATGGTGGTAAGCTCCCAGCGCTTACATGATATGTTGAATCCTACAAAAGATACCAATTGGAATTCCACGTATATCTACAAATCCAGACATGAGATGTTGCCTGTTAATCTTACTCAGGAAACACTTTTCAGCTCCAAAtctc                                                                                                                                              >  2:519472/1‑251 (MQ=255)                                                                                                    cTATATTTATAGAAACTGAAGATAATTTAAGAACCTTTAACACGAATGAATTAACACGAAACGATAATATCTTTATCCTGTTCTCCTCAGTCGATGATGGGCCAGTTATGGTGGTAAGCTCCCAGCGCTTACATGATATGCTGAATCCTACAAAAGATACCAATTGGAATTCCACGTATATCTACAAATCCAGACATGAGATGTTGCCTGTTAATCTTACTCAGGAAACACTTTTCAGCTCCAAATCTCAt                                                                                                                                            >  2:339682/1‑251 (MQ=255)                                                                                                          ttATAGAAACTGAAGATAATTTAAGAACCTTTAACACGAATGAATTAACACGAAACGATAATATCATTATCCTGTTCTCCTCAGTCGATGATGGGCCAGTTATGGTGGTAAGCTCCCAGAGCTTACATGATATGCTGAATCCTACAAAAGATACCAATTGGAATACCACGTATATCTACAAATCCAGACATGAGATGTTGCCTGTTAAACTTACACATTAAACACTTTTCAGCTCCAAATCTCATGGTaa                                                                                                                                       <  2:467788/250‑1 (MQ=255)                                                                                                           tataGAAACTGAAGATAATTTAAGAACCTTTAACACGAATGAATTAACACGAAACGATAATATCTTTATCCTGTTCTCCTCAGTCGATGATGGGCCAGTTATGGTGGTAAGCTCCCAGCGCTTACATGATATGCTGAATCCTACAAAAGATACCAATTGGAATTCCACGTATATCTACAAATCCAGACATGAGATGTTGCCTGTTAATCTTACTCAGGAAACACTTTTCAGCGCCAAATCTCATGGTaaa                                                                                                                                      <  2:148350/250‑1 (MQ=255)                                                                                                                   cTGAAGATAATTTAAGAACCTTTAACACGAATGAATTAACACGAAACGATAATATCTTTATCCTGTTCTCCTCAGTCGATGATGGGCCAGTTATGGTGGTAAGCTCCCAGCGCTTACATGATATGTTGAATCCTACAAAAGATACCAATTGGAATTCCACGTATATCTACAAATCCAGACATGAGATGTTGCCTGTTAATCTTACTCAGGAAACACTTTTCAGCTCCAAATCTCATGGTAAATATGCGCtt                                                                                                                             >  1:245555/1‑251 (MQ=255)                                                                                                                          tAATTTAAGAACCTTTACCACGAATGAATTAACACGAAACGATAATATCTTTATCCTGTTCTCCTCAGTCGATGATGGGCCAGTTATGGTGGTAAGCTCCAAGCGCTTACATGATATGTTGAATCCTACAAAAGATACCAATTGGAATTCCACGTATATCTACAAATCCAGACATGAGATGTTGCCTGTTAATCTTACTCAGGAAACACTTTTCAGCTCCAAATCTCATGGTAAATATGCGCTTTTcccc                                                                                                                       >  2:335951/1‑250 (MQ=255)                                                                                                                           aaTTTAAGAACCTTTAACACGAATGAATTAACACGAAACGATAATATCTTTATCCTGTTCTCCTCAGTCGATGATGGGCCAGTTATGGTGGTAAGCTCCCAGCGCTTACATGATATGTTGAATCCTACAAAAGATACCAATTGGAATTCCACGTATATCTACAAATCCAGACATGAGATGTTGCCTGTTAATCTTACTCAGGAAACACTTTTCAGCTCCAAATCTCATGGTAAATATGCGCTTTTCCCCAt                                                                                                                     >  1:516773/1‑251 (MQ=255)                                                                                                                                 aGAACCTTTAACACGAATGAATTAACACGAAACGATAATATCTTTATCCTGTTCTCCTCAGTCGATGATGGGCCAGTTATGGTGGTAAGCTCCCAGCGCTTACATGATATGTTGAATCCTACAAAAGATACCAATTGGAATTCCACGTATATCTACAAATCCAGACATGAGATGTTGCCTGTTAATCTTACTCAGGAAACACTTTTCAGCTCCAAATCTCATGGTAAATATGCGCTTTTCCCCATTTTTAc                                                                                                               >  1:129779/1‑251 (MQ=255)                                                                                                                                   aaCCTTTAACACGAATGAATTAACACGAAACGATAATATCATTATCCTGTTCTCCTCAGTCGATGATGGGCCAGTTATGGCGGTAAGCTCCCAGCGCTTACATGATATGTTGCATCCTACAAAAGATACCAATTGGAATTCCACGTATATCTACAAATCCAGACATGAGATGTTGCCTGTTAATCTTACTCAGGAAACACTTTTCAGCTCCAAATCTCATGGTAAATATGCGCTTTTCCCCATTTTTACTg                                                                                                             >  2:285507/1‑251 (MQ=255)                                                                                                                                          aaCACGAATGAATTAACACGAAACGATAATATCTTTATCCTGTTCTCCTCAGTCGATGATGGGCCAGTTATGGTGGTAAGCTCCCAGCGCTTACATGATATGTTGAATCCTACAAAAGATACCAATTGGAATTCCACGTATATCTACAAATCCAGACATGAGATGTTGCCTGTTAATCTTACTCAGGAAACACTTTTCAGCTCCAAATCTCATGGTAAATATGCGCTTTTCCCCATTTTTACTGCGAGTTg                                                                                                      >  1:578057/1‑251 (MQ=255)                                                                                                                                          aaCACGAATGAATTAACACGAAACGATAATATCTTTATCCTGTTCTCCTCAGTCGATGATGGGCCAGTTATGGTGGTAAGCTCCCAGCGCTTACATGATATGTTGAATCCTACAAAAGATACCAATTGGAATTCCACGTATATCTACAAATCCAGACATGAGATGTTGCCTGTTAATCTTACTCAGGAAACACTTTTCAGCTCCAAATCTCATGGTAAATATGCGCTTTTCCCCATTTTTACTGCGAGTTg                                                                                                      <  2:526235/251‑1 (MQ=255)                                                                                                                                          aaCACGAATGAATTAACACGAAACGATAATATCTTTATCCTGTTCTCCTCAGTCGATGATGGGCCAGTTATGGTGGTAAGCTCCCAGCGCTTACATGATATGCTGAATCCTACAAAAGATACCAATTGGAATTCCACGTATATCTACAAATCCAGACATGAGATGTTGCCTGTTAATCTTACTCAGGAAACACTTTTCAGCTCCAAATCTCATGGTAAATATGCGCTTTTCCCCATTTTTACTGCGAGTTg                                                                                                      <  1:145303/251‑1 (MQ=255)                                                                                                                                              cGAATGAATTAACACGAAACGATAATATCTTTATCCTGTTCTCCTCAGTCGATGATGGGCCAGTTATGGTGGTAAGCTCCCAGCGCTTACATGATATGTTGAATCCTACAAAAGATACCAATTGGAATTCCACGTATATCTACAAATCCAGACATGAGATGTTGCCTGTTAATCTTACTCAGGAAACACTTTTCAGCTCCAAATCTCATGGTAAATATGCGCTTTTCCCCATTTTTACTGCGAGTTGGCGa                                                                                                  <  2:516773/251‑1 (MQ=255)                                                                                                                                               gaatgaatTAACACGAAACGATAATATCTTTATCCTGTTCTCCTCAGTCGATGATGGGCCAGTTATGGTGGTAAGCTCCCAGCGCTTACATGATATGTTGAATCCTACAAAAGCTACCAATTGGAATTCCACGTCTATCTACAAATCCAGACATGAGATGTTGCCTGTTAATCTTACTCAGGAAACACTTTTCAGCTCCAAATCTCATGGTAAATATGCGCTTTTCCCCATTTTTACTGCGAGTTGGCg                                                                                                   >  1:105091/1‑249 (MQ=255)                                                                                                                                                       tAACACGAAACGATAATATCTTTATCCTGTTCTCCTCAGTCGATGATGGGCCAGTTATGGTGGTAAGCTCCCAGCGCTTACATGATATGGTGAATCCTACAAAAGATACCAATTGGAATTCCACGTATATCTACAAATCCAGACATGAGATGTTGCCTGTTAATCTTACTCAGGAAACACTTGTCAGCTCCAAATCTCATGGTAAATATGCGCTTTTCCCCAGTTTTACTGCGAGTTGGCGAGCTCATCg                                                                                          <  2:386583/250‑1 (MQ=255)                                                                                                                                                            ggAAACGAGAATATCGTTATCCGGTTCTCCTCAGTCGATGATGGGCCAGTTATGGTGGTAAGCTCCCAGCGCTTACATGATATGTTGAATCCTACAAAAGATACCAATTGGAATTCCACGTATATCTACAAATCCAGACATGAGATGTTGCCTGTTAATCTTACTCAGGAAACACTTTTCAGCTCCAAATCTCATGGTAAATATGCGCTTTTCCCCATTTTTACTGCGAGTTGGCGAGCTCATCGTATAAt                                                                                    <  2:245555/250‑1 (MQ=255)                                                                                                                                                             gAAACGATAATATCTTTATCCTGTTCTCCTCAGTCGATGATGGGCCAGTTATGGTGGTAAGCTCCCAGCGCTTACATGATATGCTGAATCCTACAAAAGATACCAATTGGAATTCCGCGTATATCTACAAATCCAGACATGAGATGTTGCCTGTTAATCTTACTCAGGAAACACTTTTCAGCTCCAAATCTCATGGTAAATATGCGCTTTTCCCCATTTTTACTGCGAGTTGGCGAGCTCATCGTATAATg                                                                                   >  1:162646/1‑251 (MQ=255)                                                                                                                                                                   ataataTCTTTATCCTGTTCTCCTCAGTCGATGATGGGCCAGTTATGGTGGTAAGCTCCCAGCGCTTACATGATATGTTGAATCCTACAAAAGATACCAATTGGAATTCCACGTATATCTACAAATCCAGACATGAGATGTTGCCTGTTAATCTTACTCAGGAAACACTTTTCAGCTCCAAATCTCATGGTAAATATGCGCTTTTCCCCATTTT‑ACTGCGAGTTGGCGAGCTCATCGTATAATGAATAAgg                                                                            >  2:221588/1‑251 (MQ=255)                                                                                                                                                                            ttATCCTGTTCTCCTCAGTCGATGATGGGCCAGTTATGGTGGTAAGCTCCCAGCGCTTACATGATATGTTGAATCCTACAAAAGATACCAATTGGAATTCCACGTATATCTACAAATCCAGACATGAGATGTTGCCTGTTAATCTTACTCAGGAAACACTTTTCAGCTCCAAATCTCATGGTAAATATGCGCTTTTCCCCATTTTTACTGCGAGTTGGCGAGCTCATCGTATAATGAATAAGGGTGttt                                                                      <  1:39532/249‑1 (MQ=255)                                                                                                                                                                            ttATCCTGTTCTCCTCAGTCGATGATGGGCCAGTTATGGTGGTAAGCTCCCAGCGCTTACATGATATGTTGAATCCTACAAAAGATACCAATTGGAATTCCACGTATATCTACAAATCCAGACATGAGATGTTGCCTGTTAATCTTACTCAGGAAACACTTTTCAGCTCCAAATCTCATGGTAAATATGCGCTTTTCCCCATTTTTACTGCGAGTTGGCGAGCTCATCGTATAATGAATAAGGGTGttt                                                                      <  2:131140/249‑1 (MQ=255)                                                                                                                                                                                   gTTCTCCTCAGTCGATGATGGGCCAGTTATGGTGGTAAGCTCCCAGCGCTTACATGATATGCTGAATCCTACAAAAGATACCAATTGGAATTCCACGTATATCTACAAATCCAGACATGAGATGTTGCCTGTTAATCTTACTCAGGAAACACTTTTCAGCTCCAAATCTCATGGTAAATATGCGCTTTTCCCCATTTTTACTGCGAGTTGGCGAGCTCATCGTATAATGAATAAGGGTGTTTAAGTAAAgg                                                             <  2:28710/251‑1 (MQ=255)                                                                                                                                                                                      ctcctcAGTCGATGATGGGCCAGTTATGGTGGTAAGCTCCCAGCGCTTACATGATATGTTGAATCCTACAAAAGATACCAATTAGAATTCCACGTATATCTACAAATCCAGACATGAGATGTTGCCTGTTAATCTTACTCAGGAAACACTTTTCAGCTCCAAATCTCATGGTAAATATGCGCTTTTCCCCATTTTTACTGCGAGTTGGCGAGCTCATCGTATATTGATTAAGGGTGTTTAAGTAATGGaa                                                           >  2:310757/1‑250 (MQ=255)                                                                                                                                                                                       tcctcAGTCGATGATGGGCCAGTTATGGTGGTAAGCTCCCAGCGCTTACATGATATGTTGAATCCTACAAAAGATACCAATTGGAATTCCACGTATATCTACAAATCCAGACATGAGATGTTGCCTGTTAATCTTACTCAGGAAACACTTTTCAGCTCCAAATCTCATGGTAAATATGCGCTTTTCCCCATTTTTACTGCGAGTTGGCGAGCTCATCGTATAATGAATAAGGGTGTTTAAGTAAAGGaaa                                                          >  2:146505/1‑250 (MQ=255)                                                                                                                                                                                        cctcAGTCGATGATGGGCCAGTTATGGTGGTAAGCTCCCAGCGCTTACATGATATGTTGAATCCTACAAAAGATACCAATTGGAATTCCACGTATATCTACAAATCCAGACATGAGATGTTGCCTGTTAATCTTACTCAGGAAACACTTTTCAGCTCCAAATCTCATGGTAAATATGCGCTTTTCCCCATTTTTACTGCGAGTTGGCGAGCTCATCGTATAATGAATAAGGGTGTTTAAGTAAAGGaaa                                                          <  2:105091/249‑1 (MQ=255)                                                                                                                                                                                        cctcAGTCGATGATGGGCCAGTTATGGTGGTAAGCTCCCAGCGCTTACATGATATGCTGAATCCTACAAAAGATACCAATTGGAATTCCACGTATATCTACAAATCCAGACATGAGATGTTGCCTGTTAATCTTACTCAGGAAACACTTTTCAGCTCCAAATCTCATGGTAAATATGCGCTTTTCCCCATTTTTACTGCGAGTTGGCGAGCTCATCGTATAATGAATAAGGGTGTTTAAGTAAAGGAAAAc                                                        >  2:138223/1‑251 (MQ=255)                                                                                                                                                                                            aGTCGATGATGGGCCAGTTATGGTGGTAAGCTCCCAGCGCTTACATGATATGTTGAATCCTACAAAAGATACCAATTGGAATTCCACGTATATCTACAAATCCAGACATGAGATGTTGCCTGTTAATCTTACTCAGGAAACACTTTTCAGCTCCAAATCTCATGGTAAATATGCGCTTTTCCCCATTTTTACTGCGAGTTGGCGAGCTCATCGTATAATGAATAAGGGTGTTTAAGTAAAGGAAACCATc                                                     >  2:658250/1‑250 (MQ=255)                                                                                                                                                                                               cgaggatgGGCCAGTTATGGTGGTAAGCGCCCAGCGCTTACATGATATGTTGAATCCTACAAAAGATACCAATTGGAATTCCACGTATATCTACAAATCCAGACATGAGATGTTGCCTGTTAATCTTACTCAGGAAACACTTTTCAGCTCCAAATCTCATGGTAAATATGCGCTGTTCCCCATTTTTACTGCGAGTTGGCGAGCTCATCGTATAATGAATAAGGGTGTTTAAGTAAAGGAAAActgcaccg                                                 <  2:363355/247‑8 (MQ=255)                                                                                                                                                                                               cGATGATGGGCCAGTTATGGAGGTAAGCTCCCAGCGCTTACATGATATGTTGAATCCTACAAACGCTACCAATTGGAATTCCACGTATATCTACAAATCCAGACATGAGATGTTGCCTGTTAATCTTACTCAGGAAACACTTTTCAGCTCCAAATCTCATGGTAACTATGCGCTTTTCCCCATTTTTACTGCGAGTTGGCGAGCTCATCGTCTAATGAATAAGGGTGTTTAAGTAAAGGAAAACTTCACCg                                                 >  1:363355/1‑251 (MQ=255)                                                                                                                                                                                                gatgatGGGCCAGTTATGGTGGTAAGCTCCCAGCGCTTACATGATATGTTGAATCCTACAAAAGATACCAATTGGAATTCCACGTATATCTACAAATCCAGACATGAGATGTTGCCTGTTAATCTTACTCAGGAAACACTTTTCAGCTCCAAATCTCATGGTAAATATGCGCTTTTCCCCATTTTTACTGCGAGTTGGCGAGCTCATCGTATAATGAAt                                                                                >  2:157464/1‑219 (MQ=255)                                                                                                                                                                                                gatgatGGGCCAGTTATGGTGGTAAGCTCCCAGCGCTTACATGATATGTTGAATCCTACAAAAGATACCAATTGGAATTCCACGTATATCTACAAATCCAGACATGAGATGTTGCCTGTTAATCTTACTCAGGAAACACTTTTCAGCTCCAAATCTCATGGTAAATATGCGCTTTTCCCCATTTTTACTGCGAGTTGGCGAGCTCATCGTATAATGAAt                                                                                <  1:157464/219‑1 (MQ=255)                                                                                                                                                                                                gatgatGGGCCAGTTATGGTGGTAAGCTCCCAGCGCTTACATGATATGCTGAATCCTACAAAAGATACCAATTGGAATTCCACGTATATCTACAAATCCAGACATGAGATGTTGCCTGTTAATCTTACTCAGGAAACACTTTTCAGCTCCAAATCTCATGGTAAATATGCGCTTTTCCCCATTTTTACTGCGAGTTGGCGAGCTCATCGTATAATGAATAAGGGTGTTTAAGTAAAGGAAAACATCACCGt                                                <  2:258788/251‑1 (MQ=255)                                                                                                                                                                                                 atgatgGGCCAGTTATGGTGGTAAGCTCCCAGCGCTTACATGATATGTTGAATCCTACAAAAGATACCAATTGGAATTCCACGTATATCTACAAATCCAGACATGAGATGTTGCCTGTTAATCTTACTCAGGAAACACTTTTCAGCTCCAAATCTCATGGTAAATATGCGCTTTTCCCCATTTTTACTGCGAGTTGGCGAGCTCATCGTATAATGAATAAGGGTGTTTAAGTAAAGGAAAACATCAcc                                                  >  2:442679/1‑248 (MQ=255)                                                                                                                                                                                                 atgatgGGCCAGTTATGGTGGTAAGCTCCCAGCGCTTACATGATATGTTGAATCCTACAAAAGATACCAATTGGAATTCCACGTATATCTACAAATCCAGACATGAGATGTTGCCTGTTAATCTTACTCAGGAAACACTTTTCAGCTCCAAATCTCATGGTAAATATGCGCTTTTCCCCATTTTTACTGCGAGTTGGCGAGCTCATCGTATAATGAATAAGGGTGTTTAAGTAAAGGAAAACATCAcc                                                  <  1:442679/248‑1 (MQ=255)                                                                                                                                                                                                 atgatgGGCCAGTTATGGTGGTAAGCTCCCAGCGCTTACATGATATGCTGAATCCTACAAAAGATACCAATTGGAATTCCACGTATATCTACAAATCCAGACATGAGATGTTGCCTGTTAATCTTACTCAGGAAACACTTTTCAGCTCCAAATCTCATGGCAAATATGCGCTTTTCCCCATTTTTACTGCGCGTTGGCGAGCTCATCGTATAATGAATAAGGGTGTTTAAGTAAAGGAAAACATCACCGtt                                               >  1:37418/1‑251 (MQ=255)                                                                                                                                                                                                      gggCCAGTTATGGTGGTAAGCTCCCAGCGCTTACATGATATGTTGAATCCTACAAAAGATACCAATTGGAATTCCACGTATATCTACAAATCCAGACATGAGATGTTGCCTGTTAATCTTACTCAGGAAACACTTTTCAGCTCCAAATCTCATGGTAAATATGCGCTTTTCCCCATTTTTACTGCGAGTTGGCGAGCTCATCGTATAATGAATAAGGGTGTTTAAGTAAAGGAAAACATCACCGTTCCTgg                                          >  1:502710/1‑251 (MQ=255)                                                                                                                                                                                                             ttATGGTGGTAAGCTCCCAGCGCTTACATGATATGTTGAATCCTACAAAAGATACCAATTGGAATTCCACGTATATCTACAAATCCAGACATGAGATGTTGCCTGTTAATCTTACTCAGGAAACACTTTTCAGCTCCAAATCTCATGGTAAATATGCGCTTTTCCCCATTTTTACTGCGAGTTGGCGAGCTCATCGTATAATGAATAAGGGTGTTTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTg                                   >  2:74128/1‑251 (MQ=255)                                                                                                                                                                                                                     gtAAGCTCCCAGCGCTTACATGATATGCTGAATCCTACAAAAGATACCAATTGGAATTCCACGTATATCTACAAATCAAGACATGAGATGTTGCCTGTTAATCTTACTCAGGAAACACTTTTCAGCTCCAAATCTCATGGTAAATATGCGCTTTTCCCCATTTTTACTGCGAGTTGGCGAGCTCATCGTATAATGAATAAGGGTGTTTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTg                            >  1:710982/1‑250 (MQ=255)                                                                                                                                                                                                                            cccAGCGCTTACATGATATGTTGAATCCTACAAAAGATACCAATTGGAATTCCACGTATATCTACAAATCCAGACATGAGATGTTGCCTGTTAATCTTACTCAGGAAACACTTTTCAGCTCCAAATCTCATGGTAAATATGCGCTTTTCCCCATTTTTACTGCGAGTTGGCGAGCTCATCGTATAATGAATAAGGGTGTTTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCt                    <  2:402422/251‑1 (MQ=255)                                                                                                                                                                                                                                         tGATATGTTGAAGCCTACAAAAGATACCAATTGGAATTCCACGTATATCTACAAATCCAGACATGAGATGTTGCCTGTTAATCTTACTCAGGAAACACTTTTCAGCTCCAAATCTCATGGTAAATATGCGCTTTTCCCCATTTTTACTGCGAGTTGGCGAGCTCATCGTATAATGAATAAGGGTGTTTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCg       >  1:272998/1‑251 (MQ=255)                                                                                                                                                                                                                                             atGTTGAATCCTACAAAAGATACCAATTGGAATTCCACGTATATCTACAAATCCAGACATGAGATGTTGCCTGTTAATCTTACTCAGGAAACACTTTTCAGCTCCAAATCTCATGGTAAATATGCGCTTTTCCCCATTTTTACTGCGAGTTGGCGAGCTCATCGTATAATGAATAAGGGTGTTTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGcc     >  1:58271/1‑249 (MQ=255)                                                                                                                                                                                                                                              tGTTGAATCCTACAAAAGATACCAATTGGAATTCCACGTATATCTACAAATCCAGACATGAGATGTTGCCTGTTAATCTTACTCAGGAAACACTTCTCAGCTCCAAATCTCATGGTAAATATGCGCTTTTCCCCATTTTTACTGCGAGTTGGCGAGCTCATCGTATAATGAATAAGGGTGTTTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGc  >  1:660553/1‑251 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                         TCCAACTTATCGATCAGGCCTTACGATCAACGGATGATATGATCAAAGCAAATGTTTGGCATCTTTATAAAGAGTGGATTCGTAGTGATGATGTTTCACCTATATTTATAGAAACTGAAGATAATTTAAGAACCTTTAACACGAATGAATTAACACGAAACGATAATATCTTTATCCTGTTCTCCTCAGTCGATGATGGGCCAGTTATGGTGGTAAGCTCCCAGCGCTTACATGATATGTTGAATCCTACAAAAGATACCAATTGGAATTCCACGTATATCTACAAATCCAGACATGAGATGTTGCCTGTTAATCTTACTCAGGAAACACTTTTCAGCTCCAAATCTCATGGTAAATATGCGCTTTTCCCCATTTTTACTGCGAGTTGGCGAGCTCATCGTATAATGAATAAGGGTGTTTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGC  >  NC_000913/4251862‑4252349

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