Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,088,010	2 bp→TT	intergenic (+150/‑273)	metK → / → galP	methionine adenosyltransferase/galactose:H(+) symporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,088,010	0	C	T	100.0%	133.0 / NA	42	intergenic (+150/‑274)	metK/galP	methionine adenosyltransferase/galactose:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (14/28);  total (14/28)
*	NC_000913	3,088,011	0	C	T	100.0%	141.8 / NA	42	intergenic (+151/‑273)	metK/galP	methionine adenosyltransferase/galactose:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (14/28);  total (14/28)

TGGTCGTGAACATTTCCCGTGGGAAAAAACCGACAAAGCGCAGCTGCTGCGC‑GATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAACCGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTC  >  NC_000913/3087788‑3088252                                                                                                                                                                                                                                ||                                                                                                                                                                                                                                                  tggccgtaaacgttgccctggtgAAAAAACCGCCAAAGGGCAGCTGCTG‑GCTGTTTCGTCCGGTCTGAAGTAATCTTTCTTCACCGGCGTACAAAGGCCCGCTCCGCGCGGCCCTTTGTCTTTTGGATAGGCGTTCACTCCGCATCCGTCAAAAAACCCGCCCGCTCAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAATTGATTACACCAACCACAACAGACAAAg                                                                                                                                                                                                                         <  2:387842/228‑1 (MQ=255)  ggTCGTGAACAGTTCCCGGGGGAAAAACCCGACAAAGCGCACCTGCTGCGC‑GATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAATTGATTACACCAACCACAACAGACAAAGAt                                                                                                                                                                                                                       <  2:508221/251‑1 (MQ=255)     cgttaaCATTTCCCGTGGGAAAAAACCGACAAATCCCAGCTGCCGCGC‑GATGCGCCCGGGCTGAAGTAATCTGTCTTCCCCTGCGTTCAAAGGCCAGCCTCCCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAATTGATTACACCAACCACAACAGACAAAGATTTg                                                                                                                                                                                                                    <  1:220513/247‑1 (MQ=255)     cGTGAACATTTCCCGTGGGAAAAAACCGACAAAGCGCAGCTGCTGCGC‑GATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAATTGATTACACCAACcaca                                                                                                                                                                                                                                   >  1:8690/1‑236 (MQ=255)     cGTGAACATTTCCCGTGGGAAAAAACCGACAAAGCGCAGCTGCTGCGC‑GATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAATTGATTACACCAACCACAACAGACAAAGAttt                                                                                                                                                                                                                     <  1:631500/250‑1 (MQ=255)           cATTCCCCGTGGGAAAAAACCGACAAAGCGCAGCTGCTGCGC‑GATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAATTGATTACACCAACCACAACAGACAAAGATTTGTAatat                                                                                                                                                                                                              <  1:164550/251‑1 (MQ=255)             tttCCCGTGGGAAAAAACCGACAAAGCGCAGCTGCTGCGC‑GTTGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAACAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAATTGATTACACCAACCACAACAGCCAAAGATTTGTAATAttt                                                                                                                                                                                                            >  2:181447/1‑251 (MQ=255)                    tGGGAAAAAACCGACAAAGCGCAGCTGCTGCGC‑GATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAATTGCTTACACCAACCACACCAGACAAAGATTTGTAATATTTTCAtatt                                                                                                                                                                                                     >  2:270712/1‑251 (MQ=255)                     gggAAAAAACCGACAAAGCGCAGCTGCTGCGC‑GATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAATTGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCAtatta                                                                                                                                                                                                    <  1:261651/251‑1 (MQ=255)                        aaaaaaCCGACAAAGCGAAGATGCTGTGC‑GATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGTCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCTCAGATTATCATCACAACTGAAATTGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCAtattatta                                                                                                                                                                                                 <  1:394752/251‑1 (MQ=255)                           aaaCCGACAAAGCGCAGCTGCTGCGC‑GATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAATTGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTc                                                                                                                                                                                              <  1:99180/251‑1 (MQ=255)                              ccGACAAAGCGCAGCTGCTGCGC‑GATGCTGCCGGTCTGAAGCAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTATTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACAGCCCTCCCAATAACATCATTCTTCCTGATCACGTTTCCCCGCAGATTATCATCACAACTGAAATTGCTTACACCAACCACAACAGACAAAGATTTGTAATATTTTCCTATTATTATTCGGt                                                                                                                                                                                           >  2:353514/1‑251 (MQ=255)                                                tgcGC‑GGTGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAATTGATTACACCAACCACAACAGACAAAGCTTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACAtt                                                                                                                                                                         >  2:152020/1‑251 (MQ=255)                                                 gcgc‑GATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAATTGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACAttt                                                                                                                                                                        >  2:220180/1‑251 (MQ=255)                                                   gc‑gATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAATTGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCt                                                                                                                                                                      <  1:244225/251‑1 (MQ=255)                                                        ggctgccgttatgaataaaTCTTTCTTCCCCTGCGTTCAAAGGCCAGCCCGGCGCTGGCCTTTTTCTTTTGGATAGGCGTCCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACAGCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAATTGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCttgtc                                                                                                                                                                  <  2:99462/234‑4 (MQ=255)                                                              cGGTTTGAAGTAATCTTTCTTCACTGGCGTTCAAAGGCCAGCCTCGCGTTGGCCTTTTTCTTTTTGAAAGGCTTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAATTGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTaaa                                                                                                                                                            <  2:583366/251‑1 (MQ=255)                                                                gTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAATTGATTACACCAACCACAACAGACAACGATTTGTAATATTTTCATATTATTATTCGGttt                                                                                                                                                                                         >  2:419255/1‑220 (MQ=255)                                                                gTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGCCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAATTGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGttt                                                                                                                                                                                         <  1:419255/220‑1 (MQ=255)                                                                 ccTGAAGTAATCTTTCTTCCCCGGGGTTAAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCTCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAATTGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTc                                                                                                                                                         <  2:10287/250‑1 (MQ=255)                                                                                                               gtgcttttttCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATGCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAATTGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAatt                                                                                                           <  1:224521/246‑1 (MQ=255)                                                                                                                    tttttCTTGTGGATAGGCGTTCACGCCGCATCCTGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAATTGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATTATATAACTCAattatt                                                                                                        <  2:707815/249‑1 (MQ=255)                                                                                                                     ttttCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAATTGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTAtttt                                                                                                      <  1:408297/250‑1 (MQ=255)                                                                                                                            ttggtaaggatttCACGCCCCACCCGGCATAAAAACCTCCCGCAAAACAACCTCATTCTTCCTGATCACGTTTCACCGGAGAGTATCATCGCAACTGAAATTGATTACACCAACCACAACAGACAAAGATGTGTAATATTTTCATATTATTATTCGGTTTTCACATTTGTTACATGTCTTTTCAGTAAAGTCTTAATTGTAGATAACAGCGTTTAATCTATGATGATATAACTCAAATATGTTCATGctc                                                                                               <  2:694632/239‑3 (MQ=255)                                                                                                                                           cgccgcATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTGACCGCAGATTATCATCCCAACTGAAATTGATTACACCAACCACAACAGACAAAGATTTGTaa                                                                                                                                                                                                                 >  2:402580/1‑121 (MQ=255)                                                                                                                                           cgccgcAGCCGGCAAAAAAACCGCCCGCACAATAACATCATCCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAATTGATTACACCAACCACAACAGACAAAGATTTGTaa                                                                                                                                                                                                                 <  1:402580/121‑1 (MQ=255)                                                                                                                                                  tCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAATTGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCt                                                                                                                              <  2:239459/197‑1 (MQ=255)                                                                                                                                                  tCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAATTGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCt                                                                                                                              >  1:239459/1‑197 (MQ=255)                                                                                                                                                        aaaaaaaCCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAATTGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTGCATGCACTTAAATCATAACTAAGATAAATGTTAgtg                                                                  <  2:692180/251‑1 (MQ=255)                                                                                                                                                           aaaaCCGCACTCACAATAACACCATTCTTCCATATCACGTTTCACCGCAGTTTAACATCACAACTAAAATTGATTACACCAACCACAATAGACAAAGACTTGTAATATTTTCATATTATTATTCGGTTTTCACAGGTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACCGTGACCATTACTAAGATAAATGTTAGTGTat                                                               <  2:682529/251‑2 (MQ=255)                                                                                                                                                             aaCCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCCCCGCAGATTATCATCACAACTGAAATTGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCGTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGc                                                             <  2:331590/251‑1 (MQ=255)                                                                                                                                                                    cgcACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAATTGATTACACCAACCACAACAGTCAAAGAGTTGTAATATTTTCATATTATTATTCGGtttt                                                                                                                                                                                        <  2:101869/121‑1 (MQ=255)                                                                                                                                                                    cgcACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAATTGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGtttt                                                                                                                                                                                        >  1:101869/1‑121 (MQ=255)                                                                                                                                                                                      ttcCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAATTGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTc                                     >  2:167401/1‑250 (MQ=255)                                                                                                                                                                                       tcCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAATTGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTca                                    <  1:57278/250‑1 (MQ=255)                                                                                                                                                                                            aTCACGTTTCACCGCAGATTATCATCACAACTGAAATTGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCa                                                                                                                                                >  2:489361/1‑137 (MQ=255)                                                                                                                                                                                            aTCACGTTTCACCGCAGATTATCATCACAACTGAAATTGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCa                                                                                                                                                <  1:489361/137‑1 (MQ=255)                                                                                                                                                                                              cACGTTTCACCGCAGATTATCATCACAACTGAAATTGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTCACATCttt                             <  1:486328/250‑1 (MQ=255)                                                                                                                                                                                                      aCCGCAGATTATCATCACAACTGAAATTGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTAcgtcgt                    <  1:374081/251‑1 (MQ=255)                                                                                                                                                                                                         gCAGATTATCATCACAACTGAAATTGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACt                 >  1:482757/1‑251 (MQ=255)                                                                                                                                                                                                           aGATTATCATCACAACTGAAATTGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTc                >  1:659215/1‑250 (MQ=255)                                                                                                                                                                                                                        aaCTGAAATTGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTc  <  1:589319/251‑1 (MQ=255)                                                                                                                                                                                                                                ||                                                                                                                                                                                                                                                  TGGTCGTGAACATTTCCCGTGGGAAAAAACCGACAAAGCGCAGCTGCTGCGC‑GATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAACCGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTC  >  NC_000913/3087788‑3088252

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