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%	231.3 / NA	62	R313R (CGT→CGC)	purK	N5‑carboxyaminoimidazole ribonucleotide synthase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (27/35);  total (27/35)

ATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTACGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGATCATCACCGACGGATTATTCACCACTGGTTGCGGTAACGGCAGATC  >  NC_000913/551515‑551793                                                                                                                                              |                                                                                                                                          aTAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCcgg                                                                                                                                   >  1:441973/1‑150 (MQ=255)      gCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCggacggac                                                                                                                             >  1:74870/1‑151 (MQ=255)       cACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGAct                                                                                                                            <  1:1013793/151‑1 (MQ=255)        aCCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGActt                                                                                                                           >  1:671722/1‑151 (MQ=255)             ttAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGActtctt                                                                                                                        >  1:501745/1‑149 (MQ=255)              tAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTg                                                                                                                      <  1:1010301/150‑1 (MQ=255)               acccGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTACTTTGTc                                                                                                                    <  2:759448/149‑1 (MQ=255)                 ccGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCg                                                                                                                   <  2:194618/150‑1 (MQ=255)                        tACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTg                                                                                                           >  2:62351/1‑151 (MQ=255)                              gcgcCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAg                                                                                                        <  1:142992/148‑1 (MQ=255)                              gcgcCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGa                                                                                                       <  1:464363/149‑1 (MQ=255)                              gcgcCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATg                                                                                                     >  2:20550/1‑151 (MQ=255)                                    aaaTCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCAcc                                                                                                 >  1:890980/1‑149 (MQ=255)                                        cACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCgg                                                                                            <  1:991681/150‑1 (MQ=255)                                         aCGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCa                                                                                          <  1:1056088/151‑1 (MQ=255)                                          cgccgcTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAg                                                                                         <  2:500049/151‑1 (MQ=255)                                             cgcTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGttt                                                                                      <  1:172381/151‑1 (MQ=255)                                                 ggCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCa                                                                                    <  2:750831/149‑1 (MQ=255)                                                    atatTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAg                                                                               <  2:934397/151‑1 (MQ=255)                                                      atTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGt                                                                              <  1:470515/150‑1 (MQ=255)                                                      atTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTc                                                                             <  2:69889/151‑1 (MQ=255)                                                          cggcggCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTGCAGCCAGTCATaa                                                                         <  2:594597/151‑1 (MQ=255)                                                           ggcggcAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAAt                                                                        <  1:382343/151‑1 (MQ=255)                                                            gcggcAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAAtt                                                                       <  1:480284/151‑1 (MQ=255)                                                             cggcAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTc                                                                      >  1:837856/1‑151 (MQ=255)                                                                cagcagCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTcaca                                                                   <  1:806903/151‑1 (MQ=255)                                                                cagcagCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTcac                                                                    >  1:60022/1‑150 (MQ=255)                                                                  gcagcGGGATTAAGGCTTCCAGCGTCGCAGTCAGACTCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACTGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACAt                                                                  <  2:249632/150‑1 (MQ=255)                                                                        ggATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTAcc                                                           >  1:848145/1‑151 (MQ=255)                                                                        ggATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTAcc                                                           >  2:31782/1‑151 (MQ=255)                                                                          aTTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGCCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCa                                                          >  1:208890/1‑150 (MQ=255)                                                                           ttAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAAt                                                        >  1:616817/1‑151 (MQ=255)                                                                               ggCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATc                                                       >  2:115154/1‑148 (MQ=255)                                                                                   tCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTg                                                 >  1:370318/1‑150 (MQ=255)                                                                                    ccAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGat                                               >  2:583214/1‑151 (MQ=255)                                                                                         gTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGatcatc                                           >  2:228023/1‑150 (MQ=255)                                                                                         gTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGatcatc                                           >  1:309125/1‑150 (MQ=255)                                                                                              aGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGATCATCACCGAc                                     <  1:142796/151‑1 (MQ=255)                                                                                              aGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGATCATCACCGAc                                     >  1:470817/1‑151 (MQ=255)                                                                                                tCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGATCATCACCGACgg                                   >  2:970124/1‑151 (MQ=255)                                                                                                 cAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGATCATCACCGACGGa                                  <  2:715853/151‑1 (MQ=255)                                                                                                  aGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGATCATCACCGACGGat                                 >  1:139326/1‑151 (MQ=255)                                                                                                  aGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGATCATCACCGACGGa                                  >  2:661019/1‑150 (MQ=255)                                                                                                   gACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGATCATCACCGACGGatt                                <  2:806903/151‑1 (MQ=255)                                                                                                     cgcgATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGATCATCACCGACGGattat                              <  2:994683/151‑1 (MQ=255)                                                                                                         aTGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGATCATCACCGACGGATTATTcac                          <  2:999596/151‑1 (MQ=255)                                                                                                            tgtCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGATCATCACCGACGGATTATTcaccac                       >  2:640186/1‑151 (MQ=255)                                                                                                             gtCGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGATCATCACCGACGGATTATTCACCACt                      <  2:837295/151‑1 (MQ=255)                                                                                                               cGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGATCATCACCGACGGATTATTCACCACTgg                    <  2:537784/151‑1 (MQ=255)                                                                                                               cGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGATCATCACCGACGGATTATTCACCACTgg                    >  1:359474/1‑151 (MQ=255)                                                                                                               cGCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGATCATCACCGACGGATTATTCACCACTg                     <  2:293642/150‑1 (MQ=255)                                                                                                                gCTGTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCCACCAATCAGATTGATCATCACCGACGGATTATTCACCACTGGt                   <  1:217851/151‑1 (MQ=255)                                                                                                                   gTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGATCATCACCGACGGATTATTCACCACTGGtt                  <  1:1015934/149‑1 (MQ=255)                                                                                                                   gTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGATCATCACCGACGGATTATTCACCACTGGtt                  <  1:516188/149‑1 (MQ=255)                                                                                                                   gTCGGTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGATCATCACCGACGGATTATTCACCACTGGTTGc                <  1:583246/151‑1 (MQ=255)                                                                                                                      ggTCAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGATCATCACCGACGGATTATTCACCACTGGTTGCGGt             <  1:1015931/151‑1 (MQ=255)                                                                                                                         cAAATTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGATCATCACCGACGGATTATTCACCACTGGTTGCGGTAAc          >  1:989453/1‑151 (MQ=255)                                                                                                                          aaaTTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGATCATCACCGACGGATTATTCACCACTGGTTGCGGTAACg         <  2:60021/151‑1 (MQ=255)                                                                                                                            aTTCAGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGATCATCACCGACGGATTATTCACCACTGGTTGCGGTAACGGc       <  2:747269/151‑1 (MQ=255)                                                                                                                                aGATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGATCATCACCGACGGATTATTCACCACTGGTTGCGGTAACGGc       >  1:521660/1‑147 (MQ=255)                                                                                                                                 gATGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGATCATCACCGACGGATTATTCACCACTGGTTGCGGTAACGGCAGATc  <  2:941176/151‑1 (MQ=255)                                                                                                                                   tGCCCCACTTTGCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGATCATCACCGACGGATTATTCACCACTGGTTGCGGTAACGGc       >  1:89564/1‑144 (MQ=255)                                                                                                                                              |                                                                                                                                          ATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTACGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGATCATCACCGACGGATTATTCACCACTGGTTGCGGTAACGGCAGATC  >  NC_000913/551515‑551793

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