Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,934,874	Δ1 bp	8.7%	coding (640/1317 nt)	fucP →	L‑fucose transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,934,874	0	A	.	8.7%	127.6 / 5.5	46	coding (640/1317 nt)	fucP	L‑fucose transporter
Reads supporting (aligned to +/- strand):  ref base A (23/17);  new base . (2/2);  total (27/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.69e-01

TGTTACGGTATTAGGGCCGGAAAGTAGTGGTCACTTCCGCTTAAATCTTGCGCAAACATTTAACTCGTTTGGCGCAATTATCGCGGTTGTCTTTGGGCAAAGTCTTATTTTGTCTAACGTGCCACATCAATCGCAAGACGTTCTCGATAAAATGTCTCCAGAGCAATTGAGTGCGTATAAACACAGCCTGGTATTATCGGTACAGACACCTTATATGATCATCGTGGCTATCGTGTTACTGGTCGCCCTGCTGATCATGCTGACGAAATTCCCGGCATTGCAGAGTGATAATCACAGTGACGCCAAACAAGGATCGTTCTCCGCATCGCTTTCTCGCCTGGCGCGTATTCGCCACTGGCGCTGGGCGGTATTAGCGCAATTCTGCTATGTCGGCGCACAAACGGCCTGCTGGAGCTATTTGATTCGCTACGCTGTAGAA  >  NC_000913/2934657‑2935095                                                                                                                                                                                                                          |                                                                                                                                                                                                                              ggttcggggagtAGGGGCGAACAGTAGTGGTCACTTCCGGTTAAATATTTCGCGATCATTTAACGCGTTTCGCGCAATTATCGCTGTTGTCTTTGGGCAAAGTCTTATTTTGTCTAACGTTCCACATCAATCGCAAGACGTTCTCGATAAAATGTCTCCAGAGCAATTGAGTGCGTATAAACGCAGCCTGGTATTATCGGTACAGACACCTTATTTGATCATCGTGGc                                                                                                                                                                                                                     <  2:509727/217‑1 (MQ=255)   ttACGGTATTAGGGCCGGAAAGTAGTGGTCACTTCCGTTTAAATCTTGCGCAAACATTTAACTCGTTTGGCGCAATTATCGCGGTTGTCTTTGGGCAAAGTCTTATTTTGTCTAACGTGCCACATCAATCGCAAGACGTTCTCGATAAAATGTCTCCAGAGCAATTGAGTGCGTATAAACACAGCCTGGTATTATCGGTACAGACACCTTATATGATCATCGTGGCTa                                                                                                                                                                                                                   <  1:656840/228‑1 (MQ=255)    tACGGTTTTAGGGCCGGAAAGTAGTGGTCACTTCCGGTTAAATCTTGCGCAAACATTTAACTCGTTTGGCGCAATTATCGCGGTTGTCTTTGGGCAAAGTCTTATTTTGTCTAACGTGCCACATCAATCGCAAGACGTTCTCGATAAAATGTCTCCAGAGCAATTGAGTGCGTATAAACACAGCCTGGTATTATCGGTACAGACACCTTATATGATCATCGTGGCTAt                                                                                                                                                                                                                  <  1:12779/228‑1 (MQ=255)       tgtttttgggCTGGAAATTGGTGGTCGCTTCCGATTAAATGTTGCGAAAACATTTACCTCGTTTGGCGCAATTAGCGCGGTTGTTTTTGGGTAAAGTCTTATTTTGTCTAACGTGCCACATTGATCGCAAGACGATCGCGATAAAATGTTTCCAGAGCAATTGAGTGCGTATAAACACACTCTGGTATTATCGGTACAGACTCCTTATATGATCATCGTGGCTATc                                                                                                                                                                                                                 <  2:62310/219‑1 (MQ=255)         tattatGGCCGGAAAGTAGTGGGCCCTTCCGCTTAAATCTTGCGCAAATATTGACCTCTTTTGGCGCAATTATCGCGGTTGTCTTTGGGCAAAGTCTTATTTTGTCTAACGTGCCACATCAATCGCAAGACGTTCTCTATAAAATGTCTCCAGAGCAATTGAGTGCGTATAAACACAGCCTGGTATTATCGGTACAGACACCTTATATGATCAGCGTGGCTATCgtgt                                                                                                                                                                                                             <  1:616274/228‑1 (MQ=255)                  cGGAAAGTAGTGGTCACTTCCGCTTAAATCTTGCGCAAACATTTAACTCGTTTGGCGCAATTATCGCGGTTGTCTTTGGGCAAAGTCTTATTTTGTCTAACGTGCCACATCAATCGCAAGACGTTCTCGATAAAATGTCTCCAGAGCAATTGAGTGCGTATAAACACAGCCTGGTATTATCGGTACAGACACCTTATATGATCATCGTGGCTATCGTGTTACTGGTCg                                                                                                                                                                                                    <  1:439515/228‑1 (MQ=255)                   ggAAAGTAGTGGTCACTTCCGCTTAAATCTTGCGCAAACATTTAACTCGTTTGGCGCAATTATCGCGGTTGTCTTTGGGCAAAGTCTTATTTTGTCTAACGTTCCACATCAATCGCAAGACGTTCTCGATAAAATGTCTCCACAGCAATTTCGTGCGAATAAACACAGCCTGGAATTATCGTAACAAACACCTTATATGATAATCTTGTCGATCGGTTTACCGGTccc                                                                                                                                                                                                   >  2:489920/1‑226 (MQ=255)                   ggAAAGTAGTGGTCACTTCCGCTTAAATCTTGCGCAAACATTTAACTCGTTTGGCGCAATTATCGCGGTTGTCTTTGGGCAAAGTCTTATTTTGTCTAACGTGCCACATCAATCGCAAGACGTTCTCGATCAAATGTTTCCAGAGCAATTGAGTGCGTATAAACACAGCCTGGTATTATCGGTACAGACACCTTATATGTTCATCGTGCCTATCGTGTTACTGGTc                                                                                                                                                                                                     >  1:697646/1‑226 (MQ=255)                   ggAAAGTAGTGGTCACTTCCGCTTAAATCTTGCGCAAACATTTAACTCGTTTGGCGCAATTATCGCGGTTGTCTTTGGGCAAAGTCTTATTTTGTCTAACGTGCCACATCAATCGCAAGACGTTCTCGATAAAATGTCTCCAGAGCAATTGAGTTCGTATAAACACAGCCTGGTATTATCGGTACAGACACCTTATATGATCATCGTGGCTATCGGGTTACTGGTGGc                                                                                                                                                                                                   >  1:550628/1‑228 (MQ=255)                   ggAAAGTAGTGGTCACTTCCGCTTAAATCTTGCGCAAACATTTAACTCGTTTGGCGCAATTATCGCGGTTGTCTTTGGGCAAAGTCTTATTTTGTCTAACGTGCCACATCAATCGCAAGACGTTCTCGATAAAATGTCTCCAGAGCAATTGAGTGCGTCTAAACACAGCCTGGTATTATCGGTACAGACACCTTATATG‑TCATCGTGGCTATCGCGTTAATgtgcccc                                                                                                                                                                                                  >  2:172052/1‑222 (MQ=255)                     aaaGTAGTGGTCACTTCCGCTTAAATCTTGCGCAAACATTTAACTCGTTTGGCGCAATTATCGCGGTTGTCTTTGGGCAAAGTCTTATTTTGTCTAACGTGCCACATCAATCGCAAGACGTTCTCGATAAAATGTCTCCAGAGCAATTGAGTGCGTATAAACACAGCCTGGTATTATCGGTACAGACACCTTATAGGATAATAGTGGCTAACGCGTTACTGGTCGcc                                                                                                                                                                                                  >  1:123566/1‑227 (MQ=255)                          agtGGTCACTTCCGCTTAAATCTTGCGCAAACATTTAACTCGTTTGGCGCAATTATCGCGGTTGTCTTTGGGCAAAGTCTTATTTTGTCTAACGTGCCACATCAATCGCAAGACGTTCTCGATAAAATGTCTCCAGAGCAATTGAGTGCGGATAAACACAGCCTGGTATTATCGGTCCAGAAACCTTTTTTGATCATCTGGGCTATCGTGTTAATGGTCGACctgctg                                                                                                                                                                                            >  2:159159/1‑228 (MQ=255)                                                   cgcAAACATTTAACTCGTTTGGCGCAATTATCGCGGTTGTCTTTGGGCAAAGTCTTATTTTGTCTAACGTGCCACATCAATCGCAAGACGTTCTCGATAAAATGTCTCCAGAGCAATTGAGTGCGTATAAACACAGCCTGGTATTATCGGTACAGACACCTTATATGATCATCGTGGCTATCGTGTTACTGGTCGCCCTGCTGATCATGCTGACGAAATTCCCGGCAt                                                                                                                                                                   >  1:393288/1‑228 (MQ=255)                                                   cgcAAACATTTAACTCGTTTGGCGCAATTATCGCGGTTGTCTTTGGGCAAAGTCTTATTTTGTCTAACGTGCCACATCAATCGCAAGACGTTCTCGATAAAATGTCTCCAGAGCAATTGAGTGCGTATAAACACAGCCTGGTATTATCGGTACAGACACCTTATATG‑TCATCGTGGCTATCGTGTTACTGGTCGCCCTGCTGATCATGCTTACGAATTTCCCGGCAtt                                                                                                                                                                  >  1:444666/1‑228 (MQ=255)                                                      aaaaATTTAACCGTTTTGGCGCAATTATCCCGGTTGTCTTTGGGCAAAGGTTTATTTTGTCTAACGTGCCACATCAATCGCAAGACGTTCTCGATAAAATGTCTCCAGAGCAATTGAGTGCGTATAAACACAGCCTGGTATTATCGGTACAGACACCTTATATG‑TCATCGTGGCTATCGTGTTACTGGTCGCCCTGCTGATCATGCTGACGAAATTCCCGGCATTGCa                                                                                                                                                               <  2:444666/228‑1 (MQ=255)                                                       aaaaagaaagtCGTTTGGCGCAATTAGCGCGGTTGGCTTTGGGCAAAGTCTTATTTTGTCTAACTTGCCACATCAATCGCAAGACGTTCTCGATAAAATGTCTCCAGAGCAATTGAGTGCGTATAAACACAGCCTGGTATTAGCGGTACAGACACCTTATATGATCATCGTGGCTATCGTGTTACTGGTCGCCCTGCTGATCATGCTGACGAAATTCCCGGCATTGCa                                                                                                                                                               <  2:615692/218‑1 (MQ=255)                                                       aaCATTTAACTCGTTTGGCACAATTATCGCGGTTGTCTTTGGGCAAAGTCTTATTTTGTCTAACGTGCCACATCAATCGCAAGACGTTCTCGATAAAATGTCTCCAGAGCAATTGAGTGCGTATAAACACAGCCTGGTATTATCGGTACAGACACCTTATATGATCATCGTGGCTATCGTGTTACTGGTCGCCCTGCTGATCATGCTGACGAAATTCCCGGCATTGCa                                                                                                                                                               <  1:506267/228‑1 (MQ=255)                                                              taCTCGATTGGCCCAATTACCGCTGTTGTCTTTGGGCAAAGTCTTATTTTGTCGACCGTGCCACCTCAATCCGAAGCCGTTCTCGATAAAATGTCCCCAGAGCAATTGAGTGCGTATAAACACAGCCTGGTATTATCGGTACAGTCTCCTTATATGATCATCGTGGCTATCGTGTTACTGGTCGCCCTGCTGATCATGCTGACGAAATTCCCGGCATTGCAGAGTGAt                                                                                                                                                        <  1:554966/227‑1 (MQ=255)                                                                    tttGGCGCAATTATCGCGGTTGTCTTTGGGCAAAGTCTTATTTTGTCTAACGTGCCACATCAATCGCAAGACGTTCTCGATAAAATGTCTCCAGAGCAATTGAGTGCGTATAAACACAGCCTGGTATTATCGGTACAGACACCTTATATGATCATCGTGGCTATCGTGTTACTGGTCGCCCTGCTGATAATGCTGACGAAATTCCCGGCATTGCAGAGTGATAATcac                                                                                                                                                  >  1:618919/1‑228 (MQ=255)                                                                         cgcAATTATCGCGGTTGTCTTTGGGCAAAGTCTTATTTTGTCTAACGTGCCACATCAATCGCAAGACGTTCTCGATAAAATGTCTCCAGAGCAATTGAGTGCGTATAAACACAGCCTGGTATTATCGGTACAGACACCTTATATGATCATCGTGGCTATCGTGTTACTGGTCGCCCTGCTGATCATGCTGACGAAATTCCCGGCATTGCAGAGTGTTAATcctgtct                                                                                                                                              >  1:615692/1‑221 (MQ=255)                                                                                                     aGTCTTATTTTGTCTAACGTGCCACATCAATCGCAAGACGTTCTCGATAAAATGTCTCCAAAGCAATTGAGTGCGTATAAACACAGCCTGGTATTATCGGTACAGACACCTTATATGATCATCGTGGCTATCGTGTTACTGGTCGCCCTGCTGATCATGCTGACGAAATTCCCGGCATTGCAGAGTGATAATCACAGTGACGCCAAACAAGGATCGTTCTCCGCATCg                                                                                                                 <  2:393288/228‑1 (MQ=255)                                                                                                                 tCTAACGTGCCACATCAATCGCAAGACGTTCTCGCTAAATTGTCTCCAGTGCAATTGAGTGCGTATAAACACAGCCTGGAATTATCGGTACAGACACCTTATATGCACAACGTGACCATTGTGTTACTGGACACTCCGCTGCACATGCTGACGAAGTACCCGCCATTGTAGAGTGTAAATCACAGTGTCGCCCAACACGCTCCGTTCTCCGCATCGCTTTCTcccccc                                                                                                     >  2:691942/1‑226 (MQ=255)                                                                                                                 tCTAACGTGCCACATCAATCGCAAGACGTTCTCGATAAAATGTCTCCAGAGCAATTGAGTGCGTATAAACACAGCCTGGTATTATCGGTACAGACACCTTATATGat                                                                                                                                                                                                                              >  1:41560/1‑107 (MQ=255)                                                                                                                 tCTAACGTGCCACATCAATCGCAAGACGTTCTCGATAAAATGTCGCCAGAGCAATTGAGTGCGTATAAACACAGCCTGGTATTATCGGTACAGACACCTTATATGATCATCGTGGCTATCGGGTTACTGGTCGCCCTGCTGATCATGCTGACGAAATTCACTGCATTGCAGAGTGATAATCACAGTGACGCCCAACAAGGATCGGTCTCCGCACCGCTTTCTCGcca                                                                                                      >  2:488709/1‑226 (MQ=255)                                                                                                                  cTAACGTGCCACATCAATCGCAAGACGTTCTCGATAAAATGTCTCCAGAGCAATTGAGTGCGTATAAACACAGCCTGGTATTATCGGTACAGACACCTTATATGatc                                                                                                                                                                                                                             <  2:41560/107‑1 (MQ=255)                                                                                                                  cTAACGGGCCACATCAATCGCAAGACGTTCTCGATAAAATGTCTCCAGAGCAATTGAGTGCGTATAAACACAGCCTGGTATTATCGGTACAGACACCTTATATGATCATCGTGGCTATCGTGTTACTGGTCGCCCTTCTGATCATGCTGACGAAATTCCCGGCATTGCaga                                                                                                                                                             <  1:633184/171‑1 (MQ=255)                                                                                                                             caTCAATCGCAAGACGTTCTCGATAAAATGTCTCCAGAGCAATTGAGTGCGTATAAACACAGCCTGGTATTATCGGTACAGACACCTTATATGATCATCGTGGCTATCGTTTTACTGGTCGCCCTGCTGTTCATGCTGACGAAATTCCCTGCATTGCAGATGGATATTCACAGTGACGCCAACCAAGGCCCGTTCTCCTCATCGCTTTCTCGCCTGGCGCGTATTCGc                                                                                         >  1:66371/1‑228 (MQ=255)                                                                                                                                      cAAGACGTTTCCGATAAAATGTCTCCAGAGCAATTAAGTGCGTATAAACACAGCCTGGTATTAGCGGTACAGACACCTTATATGATCATCGTGGCTATCGTGTTACTGGTCGCCCTGCTGATCATGCTGACGAAATTCCCGGCATTGCAGAGTGATAATCACAGTGACGCCAAACAAGGATCGTTCTCCGCATCGCTTTCTCGCCTGGCGCGTATTCGCCACTGgcgc                                                                                <  1:459077/228‑1 (MQ=255)                                                                                                                                           cGTTCTCGATAAAATGTCTCCAGAGAAATTGAGTGCGTTTATACCCAGCCTGGTAGTATCGGTACAGACACCTTATAGGATCATCGTGGCTACCGTGTTACTGGTCGCCCTGCTGAGCATGCTGACGAAATTCCCGGCATTGCAGAGTGATAATCACAGTGACGCCAAACAAGGATCGTTCTCCGCAGCGCTTTCTCGCCTGGCGCGTATTCGCCACTGGCGCTGGGc                                                                           <  1:248427/228‑1 (MQ=255)                                                                                                                                                         tttctcCAGAGCAATTGAGTGCTTTAAAACCCGGCCTGGGTTTATCGGTACCGACACCTTATTTGATAACCGTGGCTATCGTGTTACTGGTCGCCCTGGTGATCATGCTGACGAAACTCCCGGCATTGCAGAGTGATAATCACAGTGACGCCAAACAAGGATCGTTCTCCGCATCGCTTTCTCGCCTGGCGCGTATTCGCCACTGGCGCTGGGCGGTATTAGCGCAAt                                                             <  2:483403/226‑1 (MQ=255)                                                                                                                                                         tGTCTCCAGAGCAATTGAGTGCGTATAAACACAGCCTGGTATTATCGGTACAGACACCTTATATGATCATCGTGGCTATCGTGTTACTGGTCGCCCTGCTGATCATGCTGACGAAATTCCCGGCATTGCAGAGTGATAATCACAGTGACGCCAAACAAGGATCGTTCTCCGCATCGCTTTCTCGCCTGGCGCGTATTCGCCACTGGCGCTGGGCGGGATTAGCGCAAt                                                             >  1:500340/1‑228 (MQ=255)                                                                                                                                                                      ttttggTGGGTATAAACCCAGCTTGGTTTTATGGGTACAGACACCTTATATGATCATCGTGGCTATCGTGTTACTGGTCGCCCTGCTGATCATGCTGACGAAATTCCCGGCATTGCAGAGTGATAATCACAGTGACGCCAAACAAGGATCGTTCTCCGCATCGCTTTCTCGCCTGGCGCGTATTCGCCACTGGCGCTGGGCGGTATTAGCGCAATTCTGCTATGTCgg                                                <  2:364715/223‑1 (MQ=255)                                                                                                                                                                              cGTATAAACACAGCCTGGGATTATCGGTACAGACACCTTATATGATCATCGTGGCTATCGTGTTACTGGTCGCCCTGCTGATCATGCTGACGAAATTCCCGCCATTGCAGAGTGATAATCACAGCGACGCCAAACAAGGAGCGTTCTCCCCATACCTTGCTCGCCTGGCCGCAATTCGCCACTGGCGCCGTGCGGTCTTAGCGCCACTCTCCTATTTCGGCCCACaaa                                        >  2:88270/1‑228 (MQ=255)                                                                                                                                                                                           ccTGGTATTATCGGTACAGACACCTTATATGATCATCGTGGCTATCGTGTTACTGGTCGCCCTGCTGATCATGCTGACGAAATTCCCGGCATTGCAGAGTGATAATCACAGTGACGCCAAACAAGGATCGTTCTCCGCATCGCTTTCTCGCCTGGCGCGTATTCGCCACTGGCGCTGGGCGGTATTAGCGCAATTCTGCTTTGTCGGCGCCCAAACGGCCTGCTGGAg                           >  2:39559/1‑228 (MQ=255)                                                                                                                                                                                           ccTGGTATTATCGGTACAGACACCTTATATGATCATCGTGGCTATCGTGTTACTGGTCGCCCTGCTGATCATGCTGACGAAATTCCCGGCATTGCAGAGTGATAATCACAGTGACGCCAAACAAGGATCGTTCTCCGCATCGCTTTCTCGCCTGGCGCGTATTCGCCACTGGCGCTGGGCGGTATTAGCGCAATTCTGCTATGTCGGCGCACAAACGGCCTGCTGGAg                           >  1:443644/1‑228 (MQ=255)                                                                                                                                                                                           ccTGGTATTATCGGTACAGACACCTTATATGATCATCGTGGCTATCGTGTTACTGGTCGCCCTGCTGATCATGCTGACGAAATTCCCGGCATTGCAGAGTGATAATCACAGTGACGCCAAACAAGGATCGTTCTCCGCATCGCTTTCTCGCCTGGCGCGTATTCGCCACTGGCGCTGGGCGGCATTAGCGCAATTCTGCTATGTCGGCGCACAAACGGCCTGCTGGAg                           >  1:237384/1‑228 (MQ=255)                                                                                                                                                                                           ccTGGTATTATCGGTACAGACACCTTATATGATCATCGTGGCTATCGTGTTACTGGTCGCCCTGCTGATCATGCTGACGAAATTCCCGGCATTGCAGAGTGATAATCACAGTGACGCCAAACAAGGACCGTTCTCCGCATCGCTTTCTCGCCTGGTCCGTATTCGCCACTGGCGCTGGGCGGTATTAGCGCAACTCTGCTATGTCGCCGCACAAACGGCCGGCGGGAg                           >  2:39760/1‑228 (MQ=255)                                                                                                                                                                                                     tCGGTACAGACACCTTATATGATCATCGTGGCTATCGTGTTACTGGTCGCCCTGCTGATCATGCTGACGAAATTCCCGGCATTGCAGAGTGATAATCACAGTGACGCCAAACAAGGATCGTTCTCCGCATCGCTTTCTCGCCTGGCGCGTATTCGCCACTGGCGCTGGGCGGTATTAGCGCAATTCTGCTATGTCGGCGCACAAACGGCCTGCTGGAGCTATTTGAtt                 >  1:274610/1‑228 (MQ=255)                                                                                                                                                                                                         tACAGACACCTTATATGATCATCGTGGCTATCGTGTTACTGGTCGCCCTGCTGATCATGCTGACGAAATTCCCGGCATTGCAGAGTGATAATCACAGTGACGCCAAACAAGGATCGTTCTCCGCATCGCTTTCTCGCCTGGCGCGTATTCGCCACTGGCGCTggg                                                                            >  2:297566/1‑165 (MQ=255)                                                                                                                                                                                                         tACAGACACCTTATATGATCATCGTGGCTATCGTGTTACTGGTCGCCCTGCTGATCATGCTGACGAAATTCCCGGCATTGCAGAGTGATAATCACAGTGACGCCAAACAAGGATCGTTCTCCGCATCGCTTTCTCGCCTGGCGCGTATTCGCCACTGGCGCTGGGCGGTCTTAGCGCATTTCTGCTATGTCGGCGCACAAACGGCCTGCTGGAGCCATTTTATTCGCt             >  1:375837/1‑228 (MQ=255)                                                                                                                                                                                                          ccacaccccTTTTTTGACCATCGTGGCCATCGTGTTCCGGGGCCCCCGGCTGATCATGCTGACGAAATTCCCGGCATTGCAGAGGGATAATCACAGTGACGCCAAACAAGGATCGTTCTCCGCATCGCTTTCTCGCCTGGCGCGTATTCGCCACTGGCGCTGGGCGGTATTAGCGCAATTCTGCTATGTCGGCGCACAAACGGCCTGCTGGAGCTATTTGATTCGCTa            <  2:622662/221‑1 (MQ=255)                                                                                                                                                                                                          aCAGACTCCTTATATGATCATCGTGGCTATCGTGTTACTGGTCGCCCTGCTGATCATGCTGACGAAATTCCCGGCATTGCAGAGTGATAATCACAGTGACGCCAAACAAGGATCGTTCTCCGCATCGCTTTCTCGCCTGGCGCGTATTCGCCACTGGCGCTGGGc                                                                           <  1:297566/165‑1 (MQ=255)                                                                                                                                                                                                              acacCTTATATGATCATCGTGGCTATCGTGTTACTGGTCGCCCTGCTGATCATGCTGACGAAATTCCCGGCATTGCAGAGTGATAATCACAGTGACGCCAAACAAGGATCGTTCTCCGCATCGCTTTCTCGCCTGGCGCGTATTCGCCACTGGCGCTGGGCGGTATTAGCGCAATTCTGCTATGTCGGCGCACAAACGGCCTGCTGGAGCTATTTGATTCGCTACGCt        >  1:621853/1‑228 (MQ=255)                                                                                                                                                                                                                acCTTATATGATCATCGTGGCTATCGTGTTACTGGTCGCCCTGCTGATCATGCTGACGAAATTCCCGGCATTGCAGAGTGATAATCACAGTGACGCCAAACAAGGATCGTTCTCCGCATCGCTTTCTCGCCTGGCGCGTATTCGCCACTGGCGCTGGGCGGGATTAGCGCAATTCTGCTATGTCGGCGCACAAACGGCCTGCTGGCGCTATTTGCTTCGCTACGCTGt      >  1:175966/1‑228 (MQ=255)                                                                                                                                                                                                                   ttATATGATCATCGTGGCTATCGTGTTACTGGTCGCCCTGCTGATCATGCTGACGAAATTCCCGGCATTGCAGAGTTATAATCACAGTGACGCCAAACAAGGATCGTTCTCCGCATCGCTTTCTCGCCTGGCGCGAATTCGCCACTGGCGCTGGGCGGTATTAGCGCAATTCTGCTTTGTCGGCGCACAAACGGCCTGCTGGGGCTATTTGATTCGCTACCCTGTgga   >  1:687323/1‑228 (MQ=255)                                                                                                                                                                                                                   ttATATGATCATCGTGGCTATCGTGTTACTGGTCGCCCTGCTGATCATGCTGACGAAATTCCCGGCATTGCAGAGTGATAATCACAGTGACGCCAAACAAGGATCGTTCTCCGCATCGCTTTCTCGCCTGGCGCGTATTCGCCACTGGCGCTGGGCGGTATTAGCGCAATTCTGCTATGTCGGGGCACAAACGGCCTGCTGGAGCTATTTGATTCGCt             >  1:622662/1‑218 (MQ=255)                                                                                                                                                                                                                   ttATATGATCATCGTGGCTATCGTGTTACTGGTCGCCCTGCTGATCATGCTGACGAAATTCCCGGCATTGCAGAGAGATAATCCCAGTGACGAAAAACAAGGATCGTTCTCCGCATCGCTTTCTCGCCTGCCGCGTCTTCGCCACTGGCGCTGGGCGGCCTTAGCGCAATTCTGCTATGTCGGCGCACCAACGGCCTGCTGTAGCTGTTTGATTCGCTACACTGCaga   >  1:218678/1‑228 (MQ=255)                                                                                                                                                                                                                    tataTG‑TCATCGTGGCTATCGTGTTACTGGTCGCCCTGCTGATCATGCTGACGAAATTCCCGGCATTGCAGAGTGATAATCACAGTGACGCCAAACAAGGATCGTTCTCCGCATCGCTTTCTCGCCTGGCGCGTATTCGCCACTGGCGCTGGGCGGTATTAGCGCAATTCTGCTATGTCGGCGCACAAACGGCCTGCTGGAGCTATTTGATTCGCTACGCTGTagaa  <  1:172052/227‑1 (MQ=255)                                                                                                                                                                                                                          |                                                                                                                                                                                                                              TGTTACGGTATTAGGGCCGGAAAGTAGTGGTCACTTCCGCTTAAATCTTGCGCAAACATTTAACTCGTTTGGCGCAATTATCGCGGTTGTCTTTGGGCAAAGTCTTATTTTGTCTAACGTGCCACATCAATCGCAAGACGTTCTCGATAAAATGTCTCCAGAGCAATTGAGTGCGTATAAACACAGCCTGGTATTATCGGTACAGACACCTTATATGATCATCGTGGCTATCGTGTTACTGGTCGCCCTGCTGATCATGCTGACGAAATTCCCGGCATTGCAGAGTGATAATCACAGTGACGCCAAACAAGGATCGTTCTCCGCATCGCTTTCTCGCCTGGCGCGTATTCGCCACTGGCGCTGGGCGGTATTAGCGCAATTCTGCTATGTCGGCGCACAAACGGCCTGCTGGAGCTATTTGATTCGCTACGCTGTAGAA  >  NC_000913/2934657‑2935095

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