Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	309,985	0	T	A	79.5%	104.5 / 18.8	44	G14G (GGA→GGT)	ecpB	putative fimbrial chaperone EcpB
Reads supporting (aligned to +/- strand):  ref base T (3/6);  new base A (18/17);  total (21/23)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGTACGCTAACGTCGCCGCTCCAGATGCCTTCCGGTAGAGTGCTGTAATCTGTTACTGCGGTGGTACCATTGGTGGTACCGCTGATGATGGAGAAGGTGAAACCATCCTGTGCGGTGGT  >  NC_000913/309747‑310227                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                   ggATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTACCAAATCATCCGGCTTGTCCATTGAGATAACATGCCCGTCGGCAAGCGGTGAAAAGAGCCGTTCGAGAAGGATATTGATAAGGGGAACCCTGGCGGTACGGTTTTTGATAGTTTTGCTCAACGTGCTGCTGTCACTGTTCAAAAAACATGATATATCGCCGACATCCACGGCCTGGGCCGGAGATATACCGGAAAAcac                                                                                                                                                                                                                                         >  1:74490/1‑249 (MQ=255)          cGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGCAATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGGGCgaa                                                                                                                                                                                                                                >  2:129034/1‑249 (MQ=255)          cGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCgaga                                                                                                                                                                                                                               <  2:103629/251‑1 (MQ=255)           gCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGagag                                                                                                                                                                                                                              <  1:62048/251‑1 (MQ=255)                           agcagcAAGCTGGCGGGAGTGAGTAGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCttttt                                                                                                                                                                                                              <  1:295470/251‑1 (MQ=255)                              agcaAGCTGGCGGGAGTGAGTAGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGGGCTTTTTCAt                                                                                                                                                                                                           >  1:373990/1‑251 (MQ=255)                                   gCTGGCGGGAGTGAGTAGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACaa                                                                                                                                                                                                      <  1:88330/251‑1 (MQ=255)                                      ggCGGGAGTGAGTAGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTATGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTc                                                                                                                                                                                                   <  2:370080/251‑1 (MQ=255)                                       gCGGGAGTGAGTAGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTcc                                                                                                                                                                                                  <  2:86348/251‑1 (MQ=255)                                          ggAGTGAGTAGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCTAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTc                                                                                                                                                                                               <  1:129034/251‑1 (MQ=255)                                           gagtgagtAGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGCAAGGCAGCAGGGGCTTTTTCATAACCATTCCGTcc                                                                                                                                                                                              >  2:317770/1‑251 (MQ=255)                                            agtgagtAGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCCAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGGGCTTTTTCATAACAATTCCGNcc                                                                                                                                                                                              >  2:440915/1‑250 (MQ=255)                                               gagtAGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAgg                                                                                                                                                                                           >  1:417981/1‑250 (MQ=255)                                                   aGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATaa                                                                                                                                                                                      <  1:62222/251‑1 (MQ=255)                                                             ccGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAgggggg                                                                                                                                                                            <  2:417981/251‑1 (MQ=255)                                                                gCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTa                                                                                                                                                                         <  2:191187/251‑1 (MQ=255)                                                                 cTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTAc                                                                                                                                                                        >  1:438620/1‑251 (MQ=255)                                                                        aTTNAGATAACCTGCCCGTCGTCAAGCGGTGAAGATAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTTCTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACAGCCATCGCCTGGGCCGTAGATATACCGGAAAACAGCATAGCGAGATGCGGAATGTGCTTTGTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTActgctgct                                                                                                                                                                 <  2:452481/251‑8 (MQ=255)                                                                           gagaTAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTgg                                                                                                                                                               <  1:293095/250‑1 (MQ=255)                                                                             aaTAACCTGGCCGTCCTCAAGCGGTGAAGAGAGCCGTTCGAGGCGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGGCACTGTTCATAAACGATGATATATCGCCGACATCCAGCCGCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTAc                                                                                                                                                            <  1:163678/250‑1 (MQ=255)                                                                               tAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCa                                                                                                                                                                                                                           >  1:113319/1‑186 (MQ=255)                                                                               tAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCa                                                                                                                                                                                                                           <  2:113319/186‑1 (MQ=255)                                                                                  ccTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATc                                                                                                                                                         >  1:77740/1‑249 (MQ=255)                                                                                     gCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCag                                                                                                                                                       >  2:352419/1‑248 (MQ=255)                                                                                               aaGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGCGAGATTAACTgg                                                                                                                                          >  1:459189/1‑251 (MQ=255)                                                                                                 gCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTc                                                                                                                                        <  2:373990/251‑1 (MQ=255)                                                                                                 gCGGTGAAGAGAGCAGTTCGAGACGGATATTGATAAGGCGACCACTGTTGGTACTGTTTTTGATCGTTTTGCTCAGCGGGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATAACGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTc                                                                                                                                        <  2:7238/251‑1 (MQ=255)                                                                                                  cGGTGAAGATAGCCGTTCTAGACGGATATTGATAAGGCTACCACTGTCGGTACTGTTTTTGATCCTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTcc                                                                                                                                       <  2:399455/251‑1 (MQ=255)                                                                                                   ggTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGGGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTcc                                                                                                                                       >  1:283821/1‑250 (MQ=255)                                                                                                    gTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAg                                                                                                                                     <  1:374334/251‑1 (MQ=255)                                                                                                              ccGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGa                                                                                                                           <  1:144384/251‑1 (MQ=255)                                                                                                                                                  ggTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGTACGCTAACGTCGCCGCTCCAGATGCCTTCCg                                                                                       <  2:411085/251‑1 (MQ=255)                                                                                                                                                          gtntGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGTACGCTAACGTCGCCGCTCCAGATGCCTTCCGGTAGAGt                                                                                <  2:438620/249‑1 (MQ=255)                                                                                                                                                             tGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGTACGCTAACGTCGCCGCTCCAGATGCCTTCCGGTAGAGTGCTg                                                                            >  2:463985/1‑251 (MQ=255)                                                                                                                                                                   ttttGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGCGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGTACGCTAACGTCGCCGCTCCAGATGCCTTCCGGTAGAGTGCTGTAATCt                                                                      >  1:277688/1‑251 (MQ=255)                                                                                                                                                                       gCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGTACGCTAACGTCGCCGCTCCAGATGCCTTCCGGTAGAGTGCTGTAATCTGtt                                                                   >  1:106158/1‑250 (MQ=255)                                                                                                                                                                                         aCTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGTACGCTAACGTCGCCGCTCCAGATGCCTTCCGGTAGAGTGCTGTAATCTGTTACTGCGGTGGTACCATtgg                                                >  1:205544/1‑251 (MQ=255)                                                                                                                                                                                           tGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGTACGCTAACGTCGCCGCTCCAGATGCCTTCCGGTAGAGTGCTGTAATCTGTTACTGCGGTGGTACCATtggtg                                              >  2:5340/1‑251 (MQ=255)                                                                                                                                                                                                 tAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGTACGCTAACGTCGCCGCTCCAGATGCCTTCCGGTAGAGTGCTGTAATCTGTTACTGCGGTGGTACCATTGGTGGTACCg                                        >  1:284035/1‑251 (MQ=255)                                                                                                                                                                                                      gatgatATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGAAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGTACGCTAACGTCGCCGCTCCAGATGCCTTCCGGTAGAGTGCTGTAATCTGTTACTGCGGTGGTACCATTGGGGGTACCGCtgat                                   >  1:409218/1‑251 (MQ=255)                                                                                                                                                                                                                    gacctCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGTACGCTAACGTCGCCGCTCCAGATGCCTTCCGGTAGAGTGCTGTAATCTGTTACTGCGGTGGTACCATTGGTGGTACCGCTGATGATGGAGAAGGTGa                     <  2:74490/247‑1 (MQ=255)                                                                                                                                                                                                                     aCATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGTACGCTAACGTCGCCGCTCCAGATGCCTTCCGGTAGAGTGCTGTAATCTGTTACTGCGGTGGTACCATTGGTGGTACCGCTGATGATGGAGAAGGTGa                     >  1:128966/1‑250 (MQ=255)                                                                                                                                                                                                                         ccAGCGCCGGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGTACGCTAACGTCGCCGCTCCAGATGCCTTCCGGTAGAGTGCTGTAATCTGTTACTGCGGTGGTACCATTGGTGGTACCGCTGATGATGGAGAAGGTGAAAcc                 <  2:310321/250‑1 (MQ=255)                                                                                                                                                                                                                                       ggAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGTACGCTAACGTCGCCGCTCCAGATGCCTTCCGGTAGAGTGCTGTAATCTGTTACTGCGGTGGTACCATTGGTGGTACCGCTGATGATGGAGAAGGTGAAACCATCCTGTGCggtggt  >  2:300800/1‑251 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                   GGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGTACGCTAACGTCGCCGCTCCAGATGCCTTCCGGTAGAGTGCTGTAATCTGTTACTGCGGTGGTACCATTGGTGGTACCGCTGATGATGGAGAAGGTGAAACCATCCTGTGCGGTGGT  >  NC_000913/309747‑310227

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