Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	551,656	A→G	R313R (CGT→CGC)	purK ←	N5‑carboxyaminoimidazole ribonucleotide synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	551,656	0	A	G	100.0%	294.0 / NA	76	R313R (CGT→CGC)	purK	N5‑carboxyaminoimidazole ribonucleotide synthase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (33/43);  total (33/43)

GCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTACGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGATCATCACCGACGGATTATTCACCACTGGTTGCGGTAACGGCAGATCGGTAATCGCC  >  NC_000913/551509‑551803                                                                                                                                                    |                                                                                                                                                    gCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTgcgc                                                                                                                                                  >  1:537035/1‑151 (MQ=255) gCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTgcgc                                                                                                                                                  >  2:1143479/1‑151 (MQ=255)  cAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGcc                                                                                                                                                 >  1:948708/1‑151 (MQ=255)    aGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGcccc                                                                                                                                               <  1:334858/151‑1 (MQ=255)      aaTAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCcgg                                                                                                                                             >  1:433774/1‑151 (MQ=255)            gCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCggacggac                                                                                                                                       >  2:1011952/1‑151 (MQ=255)               ccAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGActtc                                                                                                                                    >  1:960563/1‑151 (MQ=255)               ccAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGActt                                                                                                                                     >  1:149258/1‑150 (MQ=255)                    tAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTg                                                                                                                                <  2:433777/150‑1 (MQ=255)                       ccGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGt                                                                                                                            >  2:421523/1‑151 (MQ=255)                         gAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTAc                                                                                                                          >  1:407472/1‑151 (MQ=255)                         gAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTAc                                                                                                                          <  1:163098/151‑1 (MQ=255)                          aaCTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTAcc                                                                                                                         >  2:311102/1‑151 (MQ=255)                           aCTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCa                                                                                                                        <  1:758625/151‑1 (MQ=255)                            cTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTAcc                                                                                                                         >  1:576066/1‑149 (MQ=255)                              tACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTg                                                                                                                     >  2:724452/1‑151 (MQ=255)                               aCTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTg                                                                                                                     <  1:703473/150‑1 (MQ=255)                                  ctGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCa                                                                                                                   <  1:1118599/149‑1 (MQ=255)                                    gcgcCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGAt                                                                                                                <  2:211680/150‑1 (MQ=255)                                      gcCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGc                                                                                                              >  1:1056605/1‑150 (MQ=255)                                        ccAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCAcc                                                                                                           <  1:842489/151‑1 (MQ=255)                                         cAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCAcc                                                                                                           >  2:212871/1‑150 (MQ=255)                                           aaTCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGc                                                                                                        <  1:506349/151‑1 (MQ=255)                                           aaTCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCCTTGTCGTACCAGTGCAGATGCACCAGc                                                                                                        <  2:983024/151‑1 (MQ=255)                                            aTCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCg                                                                                                       >  1:306040/1‑151 (MQ=255)                                            aTCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCg                                                                                                       <  1:147383/151‑1 (MQ=255)                                            aTCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCg                                                                                                       >  1:711087/1‑151 (MQ=255)                                             tCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCgg                                                                                                      >  1:116895/1‑151 (MQ=255)                                              cACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCgg                                                                                                      <  1:10127/150‑1 (MQ=255)                                              cACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGc                                                                                                     >  1:751191/1‑151 (MQ=255)                                               aCGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCa                                                                                                    <  2:604421/151‑1 (MQ=255)                                               aCGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCa                                                                                                    <  2:580366/151‑1 (MQ=255)                                                cgccgcTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAg                                                                                                   <  2:826596/151‑1 (MQ=255)                                                   cgcTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGttt                                                                                                <  2:537038/151‑1 (MQ=255)                                                    gcTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTc                                                                                               >  1:908569/1‑151 (MQ=255)                                                      tGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTc                                                                                               >  2:395523/1‑149 (MQ=255)                                                       ggCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCa                                                                                              <  2:1094239/149‑1 (MQ=255)                                                       ggCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAg                                                                                             <  1:49347/150‑1 (MQ=255)                                                         cATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCa                                                                                          <  1:992571/151‑1 (MQ=255)                                                            atTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGt                                                                                        <  1:32799/150‑1 (MQ=255)                                                            atTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTc                                                                                       >  2:588952/1‑151 (MQ=255)                                                              tCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCa                                                                                      <  1:892794/150‑1 (MQ=255)                                                              tCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCAt                                                                                     >  1:386199/1‑151 (MQ=255)                                                              tCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCCGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCa                                                                                      <  1:1004213/150‑1 (MQ=255)                                                               ccGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATa                                                                                    <  2:494035/151‑1 (MQ=255)                                                                cggcggCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATaa                                                                                   <  2:369005/151‑1 (MQ=255)                                                                 ggcggcAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAAt                                                                                  <  1:359484/151‑1 (MQ=255)                                                                  gcggcAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAAtt                                                                                 >  1:1027182/1‑151 (MQ=255)                                                                       agcagcGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACAt                                                                            >  1:900830/1‑151 (MQ=255)                                                                        gcagcGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATc                                                                           <  2:51831/151‑1 (MQ=255)                                                                            cGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGc                                                                         >  2:1107119/1‑149 (MQ=255)                                                                            cGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTa                                                                       <  2:834475/151‑1 (MQ=255)                                                                             gggATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTAc                                                                      <  2:130166/151‑1 (MQ=255)                                                                              ggATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTAcc                                                                     >  2:740570/1‑151 (MQ=255)                                                                               gATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCa                                                                    <  2:214351/151‑1 (MQ=255)                                                                                aTTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCaa                                                                   <  2:603513/151‑1 (MQ=255)                                                                                  tAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATc                                                                 >  1:486414/1‑151 (MQ=255)                                                                                      gCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATc                                                                 >  2:931064/1‑147 (MQ=255)                                                                                      gCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGAt                                                             <  1:588943/151‑1 (MQ=255)                                                                                         tCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGAtt                                                            <  1:951736/149‑1 (MQ=255)                                                                                          ccAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGa                                                          <  1:212868/150‑1 (MQ=255)                                                                                             gCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGatcat                                                      <  2:297396/151‑1 (MQ=255)                                                                                                   cAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGATCATCACCGa                                                <  2:285862/151‑1 (MQ=255)                                                                                                    aGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGATCATCACCGAc                                               <  2:790099/151‑1 (MQ=255)                                                                                                    aGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGATCATCACCGAc                                               >  1:855637/1‑151 (MQ=255)                                                                                                    aGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGATCATCACCGAc                                               <  2:867150/151‑1 (MQ=255)                                                                                                      tCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGATCATCACCGACgg                                             >  2:880244/1‑151 (MQ=255)                                                                                                           cgcgATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGATCATCACCGACGGattat                                        <  2:455963/151‑1 (MQ=255)                                                                                                                 gtgtCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGATCATCACCGACGGATTATTcacc                                   <  1:229164/150‑1 (MQ=255)                                                                                                                   gtCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGATCATCACCGACGGATTATTCACCACt                                <  1:699506/151‑1 (MQ=255)                                                                                                                       cTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGATCATCACCGACGGATTATTCACCACTGGtt                            >  2:726776/1‑151 (MQ=255)                                                                                                                            tgTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGATCATCACCGACGGATTATTCACCACTGGTTGCGGt                       >  1:966619/2‑151 (MQ=255)                                                                                                                             gTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCCCCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGATCATCACCGACGGATTATTCACCACTGGTTGCGGTa                      <  2:833403/151‑1 (MQ=255)                                                                                                                                  aTTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGATCATCACCGACGGATTATTCACCACTGGTTGCGGTAACGGc                 >  1:742150/1‑151 (MQ=255)                                                                                                                                           ccccACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGATCATCACCGACGGATTATTCACCACTGGTTGCGGTAACGGCAGATCGGTa        <  1:983268/151‑1 (MQ=255)                                                                                                                                              cACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGATCATCACCGACGGATTATTCACCACTGGTTGCGGTAACGGCAGATCGGTAATc     <  2:667043/151‑1 (MQ=255)                                                                                                                                                cTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGATCATCACCGACGGATTATTCACCACTGGTTGCGGTAACGGCAGATCGGTAATCg    >  2:879744/1‑150 (MQ=255)                                                                                                                                                 tttGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGATCATCACCGACGGATTATTCACCACTGGTTGCGGTAACGGCAGATCGGTAATCGcc  >  2:1148699/1‑151 (MQ=255)                                                                                                                                                    |                                                                                                                                                    GCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTACGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGATCATCACCGACGGATTATTCACCACTGGTTGCGGTAACGGCAGATCGGTAATCGCC  >  NC_000913/551509‑551803

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