Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,796,755	G→A	S259S (TCG→TCA)	waaC →	ADP‑heptose:LPS heptosyltransferase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,796,755	0	G	A	100.0%	118.4 / NA	34	S259S (TCG→TCA)	waaC	ADP‑heptose:LPS heptosyltransferase 1
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (16/18);  total (16/18)

TGCTATCGCACAGCATTTTCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCGGTGGATACGGGGTTAAGCCATTTAACGGCGGCACTGGATAGACCCAATATCACGGTTTATGGACCAACCGATCCGGGATTAATTGGTGGGTATGGGAAGAATCAGATGGTATGTAGGGCTCCAAGAGAAAATTT  >  NC_000913/3796467‑3796889                                                                                                                                                                                                                                                                                                 |                                                                                                                                       tGCTATTGCACAGCATTTCCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCAGTGGATACggg                                                                                                                             <  2:203574/300‑1 (MQ=255)     attgCACAGCATTTCCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACNGGCGAGAATTGATTGGTTTACTGGCTGATTCANGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCAGTGGATACGGGGTTa                                                                                                                         >  2:530226/4‑300 (MQ=255)        gCACAGCATTTCCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCAGTGGATACGGGGTTAAGc                                                                                                                      <  1:65380/300‑1 (MQ=255)             gCATTTCCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCAGTGGATACGGGGTTAAGCCAttt                                                                                                                 >  1:30288/1‑300 (MQ=255)             gCATTTCCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCAGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCAGTGGATACGGGGTTAAGCCAt                                                                                                                   >  2:154677/1‑298 (MQ=255)              cATTTCCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAANCACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCAGTGGATACGGGGTTAAGCCa                                                                                                                    <  1:517930/296‑1 (MQ=255)                     ttaCGAACCTGCCTACAGATGCTGGCGAATATGCCGTACTTCTTCATGCGACGACCCTTGATGATAAACACTGGCCGGAAGAACACCGTCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCTGCATGAGGAGGAACGGGCGAAACGACTGGCGGAGGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCATTGGATACGGGGTTAAGCCATTTAAcag                                                                                                            <  2:82169/295‑3 (MQ=255)                              tGCCTACAGATGCTGTCGAATATGCCGTATTTCTTCATGCGACGACCCTTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGAGTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCAGTGGATACGGGGTTAAGCCATTTAACAGCGGCACTGGAt                                                                                                 <  2:306323/299‑1 (MQ=255)                              tGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTACTGGCCGGGGCTAAATTTGTAGTGTCAGTGGATACGGGGTTAAGCCATTTAACAGCGGCACTGGAt                                                                                                 <  2:50924/299‑1 (MQ=255)                               gcccaCAGATGCTGGNCAATATGCCGTATTTCTTNATGCGACCACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCAGTGGATACGGGGTTAAGCCATTNAACAGCGGCACTGGAt                                                                                                 <  1:585857/294‑1 (MQ=255)                               gCCTACAGATGCTGGCNAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCAGTGGNTACGGGGTTAAGCCATTTNACAGCGGCACTGGATa                                                                                                >  1:582993/1‑299 (MQ=255)                               gCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTC‑‑GAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCAGTGGATACGGGGTTAAGCCATTTAACAGCGGCACTGGATAGAc                                                                                             >  1:42519/1‑300 (MQ=255)                                 cTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCAGTGGATACGGGGTTAAGCCATTTAACAGCGGCACTGGATAGAc                                                                                             >  1:173178/1‑300 (MQ=255)                                  tACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGGGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCAGTGGATACGGGGTTAAGCCATTTAACAGCGGCACTGGATAGAcc                                                                                            <  1:348600/300‑1 (MQ=255)                                  tACAGATGCTGGCGAATATGCCGTATGTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCAGTGGATACGGGGTTAAGCCATTTAACAGCGGCACTGGATAGAcc                                                                                            <  1:220446/300‑1 (MQ=255)                                       aTGCTGGCGAATATGCCGTAGTTCTTCCTGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGTCGAGAATTGATTGTTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCAGTGGATACGGGGTTAAGCCATTTAACAGCGGCACTGGATAGa                                                                                              <  2:40385/293‑1 (MQ=255)                                                     gCCGTATTTCTTCATGCGACGACCCGGGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGTGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCAGTGGATACGGGGTTAAGCCATTTAACAGCGGCACTGGATAGACCCAATATCACGGTTTATGGa                                                                         <  2:290429/300‑1 (MQ=255)                                                     gCCGTAGTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGTAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATCCATGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGTCTAAACGACTGGCGTAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGTGGCTAAATTTTTATTGTCAGTGGATACGGGGTTAAGCCATTTAACAGCGGCACTGGATAGACCCAATATCACGGTTTATGGa                                                                         <  2:290267/300‑1 (MQ=255)                                                      ccGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCAGTGGATACGGGGTTAAGCCATTTAACAGCGGCACTGGATAGACCCAATATCACGGTTTATGGAc                                                                        >  1:415561/1‑300 (MQ=255)                                                        gTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCAGTGGATACGGGGTTAAGCCATTTAACAGCGGCACTGGATAGACCCAATATCACGGTTTATGGAcct                                                                      >  1:686458/1‑299 (MQ=255)                                                        gTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCAGTGGATACGGGGTTAAGCCATTTAACAGCGGCACTGGATAGACCCAATATCAAGGTTTATGGAcct                                                                      >  1:686470/1‑299 (MQ=255)                                                                      gacgacCCTTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCAGTGGATACGGGGTTAAGCCATTTAACAGCGGCACTGGATAGACCCAATATCACGGTTTATGGACCTACCGATCCGGGAt                                                         <  1:162775/299‑1 (MQ=255)                                                                      gacgacCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACNGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCAGTGGATACGGGGTTAAGCCATTTAACAGCGGCACTGGATAGACCCAATATCACGGTTTATGGACCTACCGATCCGGGAt                                                         <  1:135315/299‑1 (MQ=255)                                                                          acCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCTGTGGGGCGCGCCGCATGAGGGAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTACTGGCCGGGGCTAAATTTGTAGTGTCAGTGGATACGGGGTTAAGCCATTTAACAGCGACACTGGATAGACCCAATATCACGGTTTATGGACCTACCGATCCGGGATTAAtt                                                    <  2:51849/300‑1 (MQ=255)                                                                                   gatAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCAGTGGATACGGGGTTAAGCCATTTAACAGCGGCACTGGATAGACCCAATATCACGGTTTATGGACCTACCGATCCGGGATTAATTGGTGGGTAt                                           >  1:330467/1‑300 (MQ=255)                                                                                   gatAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCAGTGGATACGGGGTTAAGCCATTTAACAGCGGCACTGGATAGACCCAATATCACGGTTTATGGACCTACCGATCCGGGATTAATTGGTGGGTAt                                           >  2:97857/1‑300 (MQ=255)                                                                                   gatAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCAGGGGCTAAATTTGTAGTGTCAGTGGATACGGGGTTAAGCCATTTAACAGCGGCACTGGATAGACCCAATATCACGGTTTATGGACCTACCGATCAGGGATTAGTTGGTGGGTAt                                           >  2:65401/1‑300 (MQ=255)                                                                                        acacTGGCCGTAAGAACACTGGCGAGAATTGATTGTTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGTCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCAGTGGATACGGGGTTAAGCCATTTAACAGCGGCACTGGATAGACCCAATATCACGGTTTATGGACCTACCGATCCGGGATTAATTGGTGGGTATGGgaa                                      <  2:586002/300‑1 (MQ=255)                                                                                        acacTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCAGGGGCTAAATTTGTAGTGTCAGTGGATACGGGGTTAAGCCATTTAACAGCGGCACTGGATAGACCCAATATCAAGGTTTATGGACCTACCGATCCGGGATTAATTGGTGGGTATGGgaa                                      >  2:444380/1‑300 (MQ=255)                                                                                                 ggAAGAACACTGTCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGAGTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCAGTGGATACGGGGTTAAGCCATTTAACAGCGGCACTGGATAGACCCAATATCACGGTTTATGGACCTACCGATCCGGGATTAATTGGTGGGTATGGGAAGAATCAGAt                             <  2:5656/300‑1 (MQ=255)                                                                                                  gaagaaCACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCAGTGGATACGGGGTTAAGCCATTTAACAGCGGCACTGGATAGACCCAATATCACGGTTTATGGACCTACCGATCCGGGATTAATTGGTGGGTATGGGAAGAATCAGAt                             >  2:694002/1‑299 (MQ=255)                                                                                                            ggCGAGGATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCAGTGGATACGGGGTTAAGCCATTTAACAGCGGCACTGGATAGACCCAATATCACGGTTTATGGACCTACCAATCCGGGATTAATTGGTGGGTATGGGAAGAATCAGATGGTTTGTa                     >  2:358408/1‑297 (MQ=255)                                                                                                                    ttgattgTTTTACTGGCTGATTCATGAAGACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGTGCTAAATTTGTAGTGTCAGTGGATACGGGGTTAAGCCATTTAACAGCGGCACTGGATAGACCCAATATCACGGTTTATGGACCTACCGATCCGGGATTAATTGGTGGGTATGGGAAGAATCAGATGGTTTGTAGGGCTccgggg          <  2:201325/300‑5 (MQ=255)                                                                                                                             ttACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTANTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCNTGCTGGGCGGGGCTAAATTTGTAGTGTCAGTGGATACGGGGTTAAGCCATTTAACAGCGGCACTGGATAGACCCAATATCACGGTTTATGGACCTACCGATCCGGGATTAATTGGTGGGTATGGGAAGAATCAGATGGTTTGTAGGGCTccggggaatgagtt  >  2:537383/1‑287 (MQ=255)                                                                                                                                                                                                                                                                                                 |                                                                                                                                       TGCTATCGCACAGCATTTTCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCGGTGGATACGGGGTTAAGCCATTTAACGGCGGCACTGGATAGACCCAATATCACGGTTTATGGACCAACCGATCCGGGATTAATTGGTGGGTATGGGAAGAATCAGATGGTATGTAGGGCTCCAAGAGAAAATTT  >  NC_000913/3796467‑3796889

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