Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,815,292	G→C	R159R (CGC→CGG)	pyrE ←	orotate phosphoribosyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,815,292	0	G	C	100.0%	244.1 / NA	69	R159R (CGC→CGG)	pyrE	orotate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (35/34);  total (35/34)

TTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGCAATCAGGTCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAA‑TCTCGCCGCGCCCGCGTTCCTGACGATCGAGCGAAATCAACACGCCAGCAAGCGTCGCGCCATTGGCCTGAATAATCTCCATCGACTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCTTCTTTGCGGTTAAAGCAGTACGGCAGGTCCAGGTCGTGATGCTCCGCCAGTGCCACAGCGGTTGTGGTGGC  >  NC_000913/3815050‑3815537                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                      ttgccGCCCGCGGGCCGTTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTAAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAA‑TCTCGCCGCGCCCCCGTTCCt                                                                                                                                                                                                                                                <  2:335408/246‑1 (MQ=255)    aCGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAA‑TCTCGCCGCGCCCCCGTTCCTGACg                                                                                                                                                                                                                                            <  1:334445/251‑1 (MQ=255)     cGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAA‑TCTCGCCGCGCCCCCGTTCCTGACg                                                                                                                                                                                                                                            >  1:43009/1‑250 (MQ=255)          gCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAA‑TCTCGCCGCGCCCCCGTTCCTGACGATcaag                                                                                                                                                                                                                                      <  1:352322/251‑4 (MQ=255)           gtgccctctattggCGGCACAGTGGAAGTAAAATAACGCCGGAGAACTTTCCATCCGGCGAGTTTCTTTAATCGCAAAACTCTTCGCGTGAGTCCTTAACCGCCTCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGTTAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAA‑TCTCGCCGCGCCCCCGTTCCTGACGATcaag                                                                                                                                                                                                                                      <  2:194984/237‑4 (MQ=255)               ccGCTAATTGCGGCACAGTTGCATTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAA‑TCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGaaa                                                                                                                                                                                                                                 <  2:345093/251‑1 (MQ=255)               ccGCTAATTGCGGCACAGTTGAAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAA‑TCTCGACGCGCCCCCGTTCCTGACGATCAAGCGaa                                                                                                                                                                                                                                  >  2:423298/1‑250 (MQ=255)                          ggCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAA‑TCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGcc                                                                                                                                                                                                                       >  1:432589/1‑250 (MQ=255)                                        attaggACGCCGGATGACTTTCCAGCCGGCGAGTTTCTTTAAACGCCAATCTCTTCGGGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATATGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCGGAATCGCCGAAA‑TCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAGGCGTCGCGcc                                                                                                                                                                                                        <  2:422798/246‑1 (MQ=255)                                        aaTATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAA‑TCTCGCCGCGCCCCCGTTCCTGCCGATCAAGCGAAATCAATAAGCCAGCAAGCGTCGCGcc                                                                                                                                                                                                        >  2:151919/1‑251 (MQ=255)                                             gACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAA‑TCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCAtt                                                                                                                                                                                                     >  2:285109/1‑249 (MQ=255)                                                    ggtGACTTTTCATCCGGGGAGTTTCTTTAAANCCCAAACTCTTCGCGATCGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAA‑TCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCTCCATTAGCCTGaat                                                                                                                                                                                            <  2:212996/249‑1 (MQ=255)                                                     aTGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAA‑TCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGaat                                                                                                                                                                                            >  1:346154/1‑250 (MQ=255)                                                       gACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAA‑TCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGaataat                                                                                                                                                                                         <  1:285109/251‑1 (MQ=255)                                                            ttCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAA‑TCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTcc                                                                                                                                                                                     <  1:380916/250‑1 (MQ=255)                                                                 ccGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAA‑TCTCGCCGCGCCCCCGTTCCTGACGATCACGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGat                                                                                                                                                                               >  2:51858/1‑250 (MQ=255)                                                                   ggCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAA‑TCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGattc                                                                                                                                                                             >  1:26589/1‑248 (MQ=255)                                                                     cGAGTTTCTTTAACCGCCAATCTCTTCGCGATAGGCCTTACCCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTGCCAGATAGGCAATCAGATCTTTCAGGGTGACGATAGAGATCACTTTGCAGTTTTAATCACGCTCAACTTCCTGAATCGCCGAAA‑TCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTcgc                                                                                                                                                                           <  2:423434/251‑1 (MQ=255)                                                                              ttAAACGCCAACCTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAA‑TCTCGCCGCGCCCCCGTTCCTGACTATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCAGCGATTCGCGAATCg                                                                                                                                                                     <  1:353073/248‑1 (MQ=255)                                                                                aaaCGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAA‑TCTCGCCGCGCCCCCGTTCCTGACg                                                                                                                                                                                                                                            >  2:253456/1‑175 (MQ=255)                                                                                aaaCGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAA‑TCTCGCCGCGCCCCCGTTCCTGACg                                                                                                                                                                                                                                            <  1:253456/175‑1 (MQ=255)                                                                                    gCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAA‑TCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCcgg                                                                                                                                                            >  1:298604/1‑251 (MQ=255)                                                                                       aaaCTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAA‑TCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCcggcgg                                                                                                                                                         >  2:322680/1‑251 (MQ=255)                                                                                        aaCTCTTCGCGATAGGCCTTAACCGCCGCCCGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAA‑TCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCcggcggg                                                                                                                                                        >  1:138404/1‑250 (MQ=255)                                                                                                     aGGCTTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAA‑TCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCt                                                                                                                                           <  2:346322/251‑1 (MQ=255)                                                                                                        ccTTAACCGCCGCCAGATGTTCCGACATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAAGGAGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGACGAAA‑TCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGttctaaactcatcttac                                                                                                                                        >  1:196220/1‑235 (MQ=255)                                                                                                            aaCCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAA‑TCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGc                                                                                                                                     >  1:46368/1‑250 (MQ=255)                                                                                                             aCCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTAAACTTCCTGAATCGCCGAAA‑TCTCGCCGCGCCCCCGTTCCTGTCGATCAAGCGAAATCAATACGCCCGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACCTCTTCTACCGGCAt                                                                                                                                   >  1:86476/1‑251 (MQ=255)                                                                                                                gccgccAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAA‑TCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTAc                                                                                                                                <  1:51858/251‑1 (MQ=255)                                                                                                                  cgccAGATGTTCCGCCATGTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAA‑TCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTGCCGGCGGTGAGCACATCATCTACCAGCATTAcgc                                                                                                                              <  2:180019/251‑1 (MQ=255)                                                                                                                      agctGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCTGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAA‑TCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTcc                                                                                                                          <  2:243854/248‑1 (MQ=255)                                                                                                                            tCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAA‑TCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTa                                                                                                                     <  2:72601/250‑1 (MQ=255)                                                                                                                               gCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAA‑TCTCGCCGCGCACCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCCTCTACCAGCATTCCGCGTCCTTGTAAcgc                                                                                                                 >  2:291383/1‑251 (MQ=255)                                                                                                                                 cATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAA‑TCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTTCGCGTCCTTGTAACgcgc                                                                                                               >  2:46016/1‑251 (MQ=255)                                                                                                                                  aTTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAA‑TCTCCCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGAACACATCATCTACCCGCATTACGCGTCCTTGTAACGCGCt                                                                                                              >  1:296541/1‑251 (MQ=255)                                                                                                                                      ccGGCTTCTCTTCCAGATAGGCAATCAGCTCTTTCAGGGTGATGATAGAGATCACTTCGCAGTTGTAATCACGCTCAACTTCCTGACTCGTCGAAA‑TCTCGCCGCGCCCCCGTTCCTGACGAACAAGCTAAATCAAGACGCCAGCAAGCGGCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAAcgctgtacca                                                                                                          <  2:46884/251‑8 (MQ=255)                                                                                                                                            tctcTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAA‑TCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCCTCTACCAGCATTTCGCGTCCTTGTAACGCGCTacccacc                                                                                                       >  2:282872/1‑248 (MQ=255)                                                                                                                                            tctcTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAA‑TCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTaccaacc                                                                                                       >  1:199728/1‑248 (MQ=255)                                                                                                                                                    agataggCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAA‑TCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTgccgcc                                                                                            <  1:111541/245‑1 (MQ=255)                                                                                                                                                     gataggcaatcagagttttgaGGGGGATGATAGGGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAA‑TCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTGCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCt                                                                                           <  2:197223/231‑1 (MQ=255)                                                                                                                                                      ataggCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAA‑TCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCtt                                                                                          <  2:346154/247‑1 (MQ=255)                                                                                                                                                          gCAATCAGATCTTGCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAA‑TCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCAcc                                                                                      <  1:384498/251‑1 (MQ=255)                                                                                                                                                               cagatCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAA‑TCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGCCCTTGTAACGCGCCACCAACCAGATTTCCGCCTTCCCCGTggg                                                                                 >  2:240550/5‑250 (MQ=255)                                                                                                                                                                agatCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAA‑TCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCCTCTACCAGCCTTACGCGTCCTTGTAACGCGCTACCCACCAGATTGCCGCCTTCACCGTGGTc                                                                                >  2:243618/4‑251 (MQ=255)                                                                                                                                                                  atCTTTCAGGGTGATGATAGAGATCCCTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAA‑TCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCt                                                                                                                         <  2:241270/207‑1 (MQ=255)                                                                                                                                                                  atCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAA‑TCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCt                                                                                                                         >  1:241270/2‑208 (MQ=255)                                                                                                                                                                       tCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAA‑TCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATTCCATCTCCTAACAGCATTACGTGTCCTTTTCACCCGCTACCACNCAGACTGCCACCCTCACCGTGGTCttttccc                                                                         >  2:223825/1‑248 (MQ=255)                                                                                                                                                                        cAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAA‑TCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCCTCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGccc                                                                        >  2:291027/1‑249 (MQ=255)                                                                                                                                                                         aGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAA‑TCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGcctc                                                                       <  1:38870/251‑4 (MQ=255)                                                                                                                                                                           ggTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAA‑TCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGcctctt                                                                     >  1:146882/1‑251 (MQ=255)                                                                                                                                                                           ggTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAA‑TCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGcctctt                                                                     <  1:329988/251‑1 (MQ=255)                                                                                                                                                                             tgatgaTAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAA‑TCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAAtct                                                                                                                                                                                       <  2:254776/135‑1 (MQ=255)                                                                                                                                                                             tgatgaTAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAA‑TCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAAtct                                                                                                                                                                                       >  1:254776/1‑135 (MQ=255)                                                                                                                                                                               atgatAGAGATCACTTTGCAGTTGTAATCACGCTCACCTTCCTGAATCGCCGACA‑TCTCGCCGCGCCCCAGTTCCTGACGATCAAACGAAATCAATACGCCAGCAAGCGTCGCGCCATCAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGAGGGGATCAAATCATCTACCAGCATTACGCGTCCCTTTAACGACCTCCCAACCAGATCGACCCCGTCACCGCGGACTTTTGcctatttccg                                                                 >  2:266748/1‑242 (MQ=255)                                                                                                                                                                               atgatAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAC‑TCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCCCATCATCTCCAAGCATTAAGCGTCCTTGTAACGCGCTACCCACCAGATTGCCGCCTTCACCCGGGTCTTTTGCCTCTTTGCg                                                                 >  2:88163/1‑251 (MQ=255)                                                                                                                                                                                   tAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGACA‑TCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGACTCGCCGTTCCGGCGGTGATTACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTCCCCACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGttc                                                             >  2:88290/1‑250 (MQ=255)                                                                                                                                                                                        aTCACTTTGCATGTGTAATCACGCTCAACGTCCTTAATCGCCGAAA‑TCTCGCCTCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCTCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCCTTCCGGCGGTGATCACATCATCTACCTGCAGTACGCGTCCTTGTACCGCGCTGCCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCGCTTTTCGGTTGAAGCa                                                        <  2:242623/251‑1 (MQ=255)                                                                                                                                                                                           aCTTTGCAGTTGTAATCACGCTCATCTTCCTGAATCGCCGAAAGTC‑CGCCGCGCCCCGGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAATCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGGATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTGTTGCCTCGTTGCGGTTAAAGCAGTa                                                     <  2:84347/251‑1 (MQ=255)                                                                                                                                                                                               tGCAGTTGTAATCACGCTCAACTCCCTGAATCGCCGAAA‑TCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCTGTTAAAGCAGTACgg                                                  <  2:265619/250‑1 (MQ=255)                                                                                                                                                                                                       tAATCACGCTCAACTTCATGAATCGCCGAAA‑TCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTTCCGCCTTCACCGTGGTCTTTTTCCCCTTTGCGGTTAAAGCAGTACGGcagggcc                                           >  2:367969/1‑246 (MQ=255)                                                                                                                                                                                                         atcccGCTTAACTTCCTTAATCGCCGAAA‑TCTCGCCGCGCCCCCGGTCCTGACTATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGTCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCTGTTAAAGCAGTACGGcaggtcaaggt                                       <  1:35227/247‑1 (MQ=255)                                                                                                                                                                                                              ccTCAACTTCCTGAACCGCCGAAA‑TCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGa                                  <  1:423237/250‑1 (MQ=255)                                                                                                                                                                                                                      tCCTGAATCGCCGAAA‑TCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGGCAAGGGCATGATGCTcccc                          >  1:422325/1‑249 (MQ=255)                                                                                                                                                                                                                          gAATCGCCGAAA‑TCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAgc                      <  2:146882/251‑2 (MQ=255)                                                                                                                                                                                                                           aaTCGCCGAAA‑TCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGAGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGcagggccagggtatggtggcccgcccccg                     >  2:375105/1‑232 (MQ=255)                                                                                                                                                                                                                            aTCGCCGAAA‑TCTCCCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAgcgc                    <  2:191190/251‑4 (MQ=255)                                                                                                                                                                                                                                     a‑TCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTGCTGGTGCTCCggccgcgccaccggggg           >  1:382223/1‑235 (MQ=255)                                                                                                                                                                                                                                        ctcGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGGCAAGATGGTCCGCCAGCGccccggcggggg         >  2:327145/1‑241 (MQ=255)                                                                                                                                                                                                                                               cgcCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAgtggtcgc  <  1:9585/251‑4 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                      TTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGCAATCAGGTCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAA‑TCTCGCCGCGCCCGCGTTCCTGACGATCGAGCGAAATCAACACGCCAGCAAGCGTCGCGCCATTGGCCTGAATAATCTCCATCGACTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCTTCTTTGCGGTTAAAGCAGTACGGCAGGTCCAGGTCGTGATGCTCCGCCAGTGCCACAGCGGTTGTGGTGGC  >  NC_000913/3815050‑3815537

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