Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,934,872	Δ1 bp	8.3%	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,872	0	A	.	8.3%	133.3 / 5.4	48	coding (640/1317 nt)	fucP	L‑fucose transporter
Reads supporting (aligned to +/- strand):  ref base A (23/19);  new base . (2/2);  total (27/21)
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.61e-01

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

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: ‑