Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,160,200	T→C	H316R (CAC→CGC)	ygiQ ←	radical SAM superfamily protein YgiQ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,160,200	0	T	C	100.0%	107.2 / NA	38	H316R (CAC→CGC)	ygiQ	radical SAM superfamily protein YgiQ
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (19/19);  total (19/19)

GCTCGGTGATAGAACAGAAAGAGCAGCCGCCAAAGCAGCCACGCATAATGTTGACCGAAAAACGGATCATTTCGTAAGCCGGAATACGGGCATTGCCATAGGCCGGATGTGGCACGCGCTTGTATGGCAGCGCAAAAACGCTGTCCATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGCAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGGCTGCACGGTTACGGCTTTGGCTTCCTGCTTTTTCGGTGCCACCGGTTTGTTATCCGCGCACGGCAAATCTTCACCATACGGATGCGGGATTGGGTCGATTTTTCCAGGGGTATCAAGACGGGTGGAATCC  >  NC_000913/3159955‑3160436                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                             gCTCGGTGATAGACCAGAAAGAGCAGCCGCCAAAGCAGCCACGCATAATGTGGCCCGAAAAACGGATCATTTCGTAAGCCGGAATACGGGTATTGCCATAGGCCGGATGTGGCACGCGCTTGTATGGCAGCGCAAAAACGCTGTCCATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGTTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCgcggtg                                                                                                                                                                                                                                          <  2:240680/250‑1 (MQ=255)     ggTGATAGAACAGAAAGAGCAGCCGCCAAAGCAGCCACGCATAATGTTGACCGAAAAACGGATCATTTCGTAAGCCGGAATACGGGCATTGCCATAGGCCGGATGTGGCACGCGCTTGTATGGCAGCGCAAAAACGCTGTCCATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGCGGGGCAGa                                                                                                                                                                                                                                      >  1:227502/1‑250 (MQ=255)           aGAACAGAAAGAGCAGCCGCCAAAGCAGCCACGCATAATGTTGACCGAAAAACGGATCATTTCGTAAGCCGGAATACGGGCATTGCCATAGGCCGGATGTGGCACGCGCTTGTATGGCAGCGCAAAAACGCTGTCCATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGCGGTGCAGAATAcgcg                                                                                                                                                                                                                               <  1:143661/251‑1 (MQ=255)               cAGAAAGAGCAGCCGCCAAAGCAGCCACGCATAATGTTGACCGAAAAACGGATCATTTCGTAAGCCGGAATACGGGCATTGCCATAGGCCGGATGTGGCACGCGCTTGTATGGCAGCGCAAAAACGCTGTCCATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGCGGTGCAGAATACGCGAAGc                                                                                                                                                                                                                           >  1:49987/1‑251 (MQ=255)                         agcCGCCAAAGCAGCCACGCATAATGTTGACCGAAAAACGGATCATTTCGTAAGCCGGAATACGGGCATTGCCATAGGCCGGATGTGGCACGCGCTTGTATGGCAGCGCAAAAACGCTGTCCATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGCGGGGCAGAATACGCGAAGCATGGGCGTAc                                                                                                                                                                                                                 >  2:84187/1‑251 (MQ=255)                             gccACAGCAGCAAAGCAGAATATTGACCGAAACACGGATCGGTTCGTAAGACGGAATGCGGGCCTTGACAGAGGCCGGATGTGGGGCGCGCGTGGGTGGCAGGTGAAAGACGGTGTCCATCTCGTCGGTAGAAAGCGGAATAGCAGTCAGGTTGATCCGCACATAGCGGTCGCCGTGTTTGGGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGCGGTGGAGAATACGCGAAGCAGGGGTGTACAGCa                                                                                                                                                                                                             <  2:188247/251‑1 (MQ=255)                              ccAAAGCAGCCACGCATAATGTTGACCGAAAAACGGATCATTTCGTAAGCCGGAATACGGGCATTGCCATAGGCCGGATGTGGCACGCGCTTGTATGGCAGCGCAAAAACGCTGTCCATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTGTCGCGGTGCAGAATACGCGAAGCATGGGCGTGCAGCAc                                                                                                                                                                                                            <  1:88753/251‑1 (MQ=255)                               cAAAGCAGCCACGCATAATGTTGACCGAAAAACGGATCATTTCGTAAGCCGGAATACGGGCATTGCCATAGGCCGGATGTGGCACGCGCTTGTATGGCAGCGCAAAAACGCTGTCCATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGCGGTGCAGAATACGCGAAGCATGGGCGTACAGCACt                                                                                                                                                                                                           >  1:246756/1‑251 (MQ=255)                                   gcagcCACGCATAATGTTGACCGAAAAACGGATCATTTCGTAAGCCGGAATACGGGCATTGCCATAGGCCGGATGTGGCACGCGCTTGTATGGCAGCGCAAAAACGCTGTCCATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGCGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTAt                                                                                                                                                                                                       <  1:84187/251‑1 (MQ=255)                                       cctcGCATAAGGTTGACCGAAAACAGGATCGTGTCGTTAGCCTGAATTCGGGCATGGCCATAGGCCGGATGTGGCACGCGCTGGTATGGCTGCGCAAAACCTCTGTCCATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATTCACACAGAGCGGTCTTCGTGTTTTTGTATCAATGCGGGTGCACAGCCGTGGTTGTTGTCGCGGTGCAGGATACGGTAAGTTTGGGCTTACAGCTCTTTATCGcc                                                                                                                                                                                                   <  2:159315/248‑1 (MQ=255)                                                aTGTTGACCGAAAAACGGATCATTTCGTAAGCCGGAATACGGGCATTGCCATAGGCCGGATGTGGCACGCGCTTGTATGGCAGCGCAAAAACGCTGTCCATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCCCAGCCTGGGTTGGTTTCGCGGTGCAGACTACGCGAAGCATGGGCGTACAGCCATTTATCGCCCTTCCCttt                                                                                                                                                                                          >  2:172095/1‑251 (MQ=255)                                                  gTTGACCGAAAAACGGATCATTTCGTAAGCCGGAATACGGGCATTGCCATAGGCCGGATGTGGCACGCGCTTGTATGGCAGCGCAAAAACGCTGTCCATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGCGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTtc                                                                                                                                                                                         <  1:222335/250‑1 (MQ=255)                                                       ccGAAAAACGGATCATTTCGTAAGCCGGAATACGGGCATTGCCATAGGCCGGATGTGGCACGCGCTTGTATGGCAGCGCAAAAACGCTGTCCATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGCGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGaaa                                                                                                                                                                                   <  1:109337/251‑1 (MQ=255)                                                                ggATCATTTCGTAAGCCGGAATACGGGCATTGCCATAGGCCGGATGTGGCACGCGCTTGTATGGCAGCGCAAAAACGCTGTCCATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGCGGTGCAGAATACGCGAAGCATGGGCGTACAGCACttt                                                                                                                                                                                                         >  1:202564/1‑220 (MQ=255)                                                                ggATCATTTCGTAAGCCGGAATACGGGCATTGCCATAGGCCGGATGTGGCACGCGCTTGTATGGCAGCGCAAAAACGCTGTCCATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGGTTTTGCATCAATGCTCGGGCACAGCCTGGGTTGGTTTCGCGGTGCAGAATACGCGAAGGATGGGCGTACAGCACttt                                                                                                                                                                                                         <  2:202564/220‑1 (MQ=255)                                                                ggATCATTTCGTAAGCCGGAATACGCGCATTGCCATAGGCCGGATATGGCACGCGCTTGTATGGCCGCGCAAAAACGCTGTCCATCTCTTCGGTAGAAAGCGGACTAGCAGGCGGGTTGATCCACACATACCGGCCGCCGTGTTTTTGCCTCAATGCGCGGGCACAACCTGGGTTGGTTTCGCGGTGCCGAATACGCGAAGCATGGGCGTACAGCACTTTTTCGCCCTTCCCTTTCTCGAAAGAAGgcccc                                                                                                                                                                          >  2:100217/1‑248 (MQ=255)                                                                     aTTTCGTAAGCCGGAATACGGGCATTGCCATAGGCCGGAGGTGGCACGCGCTTGTATGGCAGCGCAAAAACGCTGTCCATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGACTGGGTTGGTTTCGCGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTGTCTCGAAAGAAGGCAGCAacac                                                                                                                                                                     <  2:224110/251‑1 (MQ=255)                                                                           tAAGCCGGAAAACGGGCATTGCCATCGGCCGGATGTGGCCCGCGCTTGTATGGCAGCGCAAAAACGCTGACCATCTCTTCGGTAGAAAGAGGAATAGCAGGAGGGTTGATCCACAAATAGCGGTCGCCGTGTTTTTGCATCCATGAGCGGGCACAGCCTGGGTTGGTGTCGCGGCGCAGAATACGCGAAGCATGGGCGTACCGAACTTTATCCACCTTCACTTTCTCGAACGAAGGCAGCAACACGTCGGt                                                                                                                                                               >  2:40136/1‑251 (MQ=255)                                                                                 ggAATACGGGCATTGCCATAGGCCGGATGTGGCACGCGCTTGTATGGCAGCGCAAAAACGCTGTCCATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGCGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGCAGCAACACGGAGGTTTTTTc                                                                                                                                                         >  1:169219/1‑251 (MQ=255)                                                                                        gggCATTGCCATAGGCCGGATGTGGCACGCGCTTGTATGGAAGCGCAAAAACGCTGTCAATCTCTTCGGTAGAAAGCGGACTAGCAGGCGGGTTGATAAACACATAGCGGCCGCCGAGTTTTTGAATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGCGGCGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGCAGCAACACGTAGGTTTTTTCCCACGGt                                                                                                                                                  >  1:295201/1‑251 (MQ=255)                                                                                           cATTGCCCTAGGCCGGATGTGGCACGCGCTTGTATGGCAGCGCAAAACCGCTGTCCATCTCTTCGGTAGCAAGCGGAATAGCAGGCGGGTTGATCCCCAAATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGCGGTGCAGAAAAAGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAATGCAGCACCACGAAGGTTTTTTCCCACGGtttt                                                                                                                                               >  2:310385/1‑250 (MQ=255)                                                                                           cATTGCCATAGGCCGGATGTGGCACGCGCTTGTATGGCAGCGCAAAAACGCTGTCCATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGCGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGCAGCAACACGTAGGTTTTTTCCCACGGTTTc                                                                                                                                               <  2:246756/251‑1 (MQ=255)                                                                                                                                    gaaaaaaCGCTGGCCATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGCGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGCAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGGCTGCACGGTTACGGCTTTGGCTGCCTGCtt                                                                                                      <  2:158734/249‑1 (MQ=255)                                                                                                                                               gTCCATCTCTTCGGTAGACAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGCGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTAACTTTCTCGAAAGAAGGCAGCAACACGTCGGTTTTTTCCCACGGTTTCGGGCGCGGTGGCTGCACGGTTACGGCTTTGGCTTCCTGCTTTTTCGGTGcc                                                                                            >  2:51548/1‑250 (MQ=255)                                                                                                                                                    tctcTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGCGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGCAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGGCTGCACGGTTACGGCTTTGGCTTCCTGCTTTTTCGGTGCCACCGGt                                                                                      >  1:262790/1‑251 (MQ=255)                                                                                                                                                              aGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGGCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGCGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGCAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGGCTGCACGGTTACGGCTTTGGCTTCCTGCTTTTTCGGTGCCACCGGTTTGTTATCcg                                                                            <  2:169713/251‑1 (MQ=255)                                                                                                                                                                      gAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGCGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGCAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGGCTGCACGGTTACGGCTTTGGCTTCCTGCTTTTTCGGTGCCACCGGTTTGTTATCCGCGCACGGc                                                                    >  1:223922/1‑251 (MQ=255)                                                                                                                                                                       aaTAGCAGGCGGGTTGATCCACACATAGCGGTCGACGTGTCTTTGCATCAATGGGCGGGCACAGCCCGGGTTGGTTTCGCGGTGCAGAATACGCGCAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGCAGCAACACGTAGGCTTTTTCCCACGGTTTCGGGCGCGGTGGCTGCACGGTTACGGCTTTGGCTTCCTGCTTTTTCGGTGCCACCGGTTTGTTATCCGGGCACGGc                                                                    >  1:144358/1‑250 (MQ=255)                                                                                                                                                                               gCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGCGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGCAGCAACACGTAGGTTTTTTCCCTCGGTTTCGGGCGCGGTGGCTGCACGGTTACGGCTTTTGCTTCCTTCTTTTTCGTTGCCACCGGTTTGTTATCCGCGCACGGCAAATCTTCa                                                           <  1:285479/251‑1 (MQ=255)                                                                                                                                                                                      gATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGCGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGCAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGGCTGCACGGTTACGGCTTTGGCTTCCTGCTTTTTCGGTGCCACCGGTTTGTTATCCGCGCACGGCAAATCTTCCCCATACg                                                    >  1:57119/1‑251 (MQ=255)                                                                                                                                                                                                 aGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGCGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGCAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGGCTGCACGGTTACGGCTTTGGCTTCCTGCTTTTTCGGTGCCACCGGTTTGTTATCCGCGCACGGCAAATCTTCACCATACGGATGCGGGAt                                          <  1:112214/250‑1 (MQ=255)                                                                                                                                                                                                 aGCGGGCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGCGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGCAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGGCTGCACGGTTACGGCTTTGGCTTCCTGCTTTTTCGGTGCCACCGGTTTGTTATCCGCGCACGGCAAATCTTCACCATACGGATGCGGGAt                                          <  1:104594/250‑1 (MQ=255)                                                                                                                                                                                                       cGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGCGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGCAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGGCTGCACGGTTACGGCTTTGGCTTCCTGCTTTTTCGGTGCCACCGGTTTGTTATCCGCGCACGGCAAATCTTCACCATACGGATGCGGGATTGGGTCg                                   <  2:227502/251‑1 (MQ=255)                                                                                                                                                                                                               ttttGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGCGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGCAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGGCTGCACGGTTACGGCTTTGGCTTCCTGCTTTTTCGGTGCCACCGGTTTGTTATCCGCGCACGGCAAATCTTCACCATACGGATGCGGGATTGGGTCGATTTTTc                            <  2:57447/250‑1 (MQ=255)                                                                                                                                                                                                                tttGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGCGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGCAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGGCTGCACGGTTACGGCTTTGGCTTCCTGCTTTTTCGGTGCCACCGGTTTGTTATCCGCGCACGGCAAATCTTCACCATACGGATGCGGGATTGGGGCGATTTTTcc                           >  2:51049/1‑250 (MQ=255)                                                                                                                                                                                                                             gcgGGCACAGTCTGGGTTGGTTTCGCGGTGCAGAATACGCGAAGCAAGGGCGTACAGCACTTTATCGCCCATCAGTCTCTCGAAAGAAGGCAGCAACTCGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTTGCTGCTCGGTTACGGCTTTGGCGTCCTGCTATTTCGGTACCACCTGTTTGTTAGCCGCGCACGGCTAATCGTCGCGATTCGGATGCGGGATTGGGTCAATTTGTCCAGGGGTAGCAAGAc             <  1:252534/251‑1 (MQ=255)                                                                                                                                                                                                                                       ccTGTGTTGGTTTCGCGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGCAGCAACACGTAGGTTTTGTCCCACGGTTTCGGGCGCGGTGGCTGCACGGTTACGGCTGTGGCTGCCTGCTTTTTCGGTGCCTCCGGTTTTTTATCCGCGCACGGCAAATCTTCa                                                           <  2:303968/195‑1 (MQ=255)                                                                                                                                                                                                                                       ccTGGGTTGGTTTCGCGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGCAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGGCTGCACGGTTACGGCTTTGGCTTCCTGCTTTTTCGGTGCCACCGGTTTGTTATCCGCGCACGGCAAATCTTCa                                                           >  1:303968/1‑195 (MQ=255)                                                                                                                                                                                                                                         tGGGTTGGTTTCGCGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGCAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGGCTGCACGGTTACGGCTTTGGCTTCCTGCTTTTTCGGTGCCACCGGTTTGTTATCCGCGCACGGCAAATCTTCACCATACGGATGCGGGATTGGGTCGATTTTTCCAGGGGTATCAAGACGGGTGGAATcc  >  2:311891/1‑250 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                             GCTCGGTGATAGAACAGAAAGAGCAGCCGCCAAAGCAGCCACGCATAATGTTGACCGAAAAACGGATCATTTCGTAAGCCGGAATACGGGCATTGCCATAGGCCGGATGTGGCACGCGCTTGTATGGCAGCGCAAAAACGCTGTCCATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGCAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGGCTGCACGGTTACGGCTTTGGCTTCCTGCTTTTTCGGTGCCACCGGTTTGTTATCCGCGCACGGCAAATCTTCACCATACGGATGCGGGATTGGGTCGATTTTTCCAGGGGTATCAAGACGGGTGGAATCC  >  NC_000913/3159955‑3160436

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