Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,064,756	A→G	L156S (TTG→TCG)	cobS ←	cobalamin 5'‑phosphate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,064,756	0	A	G	100.0%	148.2 / NA	42	L156S (TTG→TCG)	cobS	cobalamin 5'‑phosphate synthase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (26/16);  total (26/16)

GCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCGAGGATGAAAATCGCCACCATCGTCACCACCATAGCAGCCACACCATGCATACCCGGCAATAATACGGCAGCAAAAATAGCGGCTAAGCCGAGTGTGACGCAGGTTTGTCGCCCATCAATTTTGC‑CAATGAACACATTGCCAAGCCCTTCTTCCCGCGCATAACGATGACGATACATCAATAATGCGGCAGTCCCGCGGCTGACCGCACATGCTGCCGCCAGCGATGCAAGGATCGATTCGCCACGCAGAGCCAACTCGCTTAACACCAGAATCTTTGCCAGTACCACAAAAATCAATGCCAGACCGCCGTGGGTGCCTAAACGGCTATCACGCATGATTTCCAACATGCGATCGCGGCTACGTGCAGAAAATACGCCATCACAGGTATCAGCGAGGCCGTCCAGATGGAATCCCCCGGTC  >  NC_000913/2064514‑2064998                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                   gCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCGAGGATGAAAATCGCCACCATCGTCACCACCATAGCAGCCACACCATGCATACCCGGCAATAATACGGCAGCAAAAATAGCGGCTAAGCCGAGTGTGACGCAGGTTTGTCGCCCATCAATTTTGC‑CAATGAACACATTGCCAAGCCCTTCTTCCCGCGCATAACGATGGCGATACATCGATAATGCg                                                                                                                                                                                                                                            <  1:328581/251‑1 (MQ=255)      aTGGCTGCTCCCAGCGTATCGCCAGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCGAGGATGAAAATCGCCACCATCGTCACCACCATAGCAGCCACACCATGCATACCCGGCAATAATACGGCAGCAAAAATAGCGGCTAAGCCGAGTGTGACGCAGGTTTGTCGCCCATCAATTTTGC‑CAATGAACACATTGCCAAGCCCTTCTTCCCGCGCATAACGATGACGATACATCGATAATGCGGCAGt                                                                                                                                                                                                                                       >  1:281109/1‑251 (MQ=255)         gctgctCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTCCGAACCTGACCGAGGCTGAAAATCGCCCCCATCGTCACCACCATAGCAGCCACACCATGCATACCCGGCAATAATACGGCAGCAAAAATAGCGGCTAAGCCGAGTGTGACGCAGGTTTGTCGCCCATCAATTTTGC‑CAATGAACACATTGCCAAGCCCTTCTTCCCGCGCATAACGATGACGATACATCGATAATGCGGCAGTccc                                                                                                                                                                                                                                    >  1:188805/1‑251 (MQ=255)               cccAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCGAGGATGAAAATCGCCACCATCGTCACCACCATAGCAGCCACACCATGCATACCCGGCAATAATACGGCAGCAAAAATAGCGGCTAAGCCGAGTGTGACGCAGGTTTGTCGCCCATCAATTTTGC‑CAATGAACACATTGCCAAGCCCTTCTTCCCGCGCATAACGATGACGATACATCGATAATGCGGCAGTCCCGCGGCt                                                                                                                                                                                                                              <  2:60379/251‑1 (MQ=255)                 cAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCGAGGATGAAAATCGCCACCATCGTCACCACCATAGCAGACACACCATGCATACCCGGCAATAATACGGCAGCAAAAATAGCGGCTAAGCCGAGTGTGACGCAGGTTTGTCGCCCATCAATTTTGC‑CAATGAACACCTTGCCAAGCCCTTCTTCCCGCGCATAACGATGACGATACATCGATAATGCGGCAGGCCcgctgatg                                                                                                                                                                                                                             >  1:232851/1‑245 (MQ=255)                    cGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCGAGGATGAAAATCGCCACCATCGTCACCACCATAGCAGCCACACCATGCATACCCGGCAATAATACGGCAGCAAAAATAGCGGCTAAGCCGAGTGTGACGCAGGTTTGTCGCCCATCAATTTTGC‑CAATGAACACATTGCCAAGCCCTTCTTCCCGCGCATAACGATGACGATACATCGATAATGCGGCAGTCCCGCGGCTGACCg                                                                                                                                                                                                                         >  2:37906/1‑251 (MQ=255)                     gTATCGCTCGTTTGCCAGCGCAGCGTACGTTCGAGCAGCCGACCGAGGATGAAAATAGCCACCATCGTCACCACCATAGCAGCCACACCATGCATACCCGGCAATAATCCGGCAGCACACATAGCGGCTACACCGAGTGCGACGAAGGTTTGTCGCCCATCAATTTTGC‑CAATGAACACATTGCCAAGCCCTTCTTCCCGCGCATAACGATGACGATACAACGATAATGCGGCAGCCCCGAGGCTGACCGc                                                                                                                                                                                                                        >  2:44536/1‑251 (MQ=255)                                tggcccgcccAGCGTACGTTTGAGAAGATGACCGAGGATGAAAATCGCCACCATCGTCACCACCATAGCAGCCACACCATGCATACCCGGCAATAATACGGCAGCAAAAAGAGCGGCTAAGCCGAGTGTGACGCAGGTTTGTCGCCCATCAAGTTTGC‑CAATGAACACATTGCCAAGCCCTTCTTCCCGCGCATAACGATGACGATACATCGAGAATGCGGCAGTCACGCGGCTGACCGCACATGCTgcc                                                                                                                                                                                                              <  2:224310/248‑1 (MQ=255)                                  gcccgcccAGCGTACGTTTGAGAAGCTGACCGAGGATGAAAATCGCCACCATCGTCACCACCATAGCAGCCACACCATGCATACCCGGCAATAATACGGCAGCAAAAATAGCGGCTAAGCCGAGTGTGACGCAGGTTTGTCGCCCATCAATTTTGC‑CAATGAACACATTGCCAAGCCCTTCTTCCCGCGCATAACGATGACGATACATCGATAATGCGGCAGTCCCGCGGCTGACCGCACATGCTgccgcc                                                                                                                                                                                                           >  1:294924/1‑251 (MQ=255)                                        ccAGCGTACGTTTGAGAAGCTGACCGAGGATGAAAATCGCCACCATCGTCACCACCATAGCAGCCACACCATGCATACCCGGCAATAATACGGCAGCAAAAATAGCGGCTAAGCCGAGTGTGACGCAGGTTTGTCGCCCATCAATTTTGC‑CAATGAACACATTGCCAAGCCCTTCTTCCCGCGCATAACGATGACGATACATCGATAATGCGGCAGTCCCGCGGCTGACCGCACATGCTGCCGCCAGCGAt                                                                                                                                                                                                     >  1:214225/1‑251 (MQ=255)                                        ccAGCGTACGTTTGAGAAGCTGACCGAGGATGAAAATCGCCACCATCGTCACCACCATAGCAGCCACACCATGCATACCCGGCAATAATACGGCAGCAAAAATAGCGGCTAAGCCGAGTGTGACGCAGGTTTGTCGCCCATCAATTTTGC‑CAATGAACACATTGCCAAGCCCTTCTTCCCGCGCATAACGATGACGATACATCGATAATGCGGCAGTCCCGCGGCTGACCGCACATGCTGCCGCCAGCGAt                                                                                                                                                                                                     >  1:291899/1‑251 (MQ=255)                                          aGCGTACGTTTGAGAAGCTGACCGAGGATGAAAATCGCCACCATCGTCCCCACCATCGCAGCCACACCATGCATACCCGGCAATAATACGGCCGCAAAAATAGCGGCTAAGCCGAGTGTGACGCAGGTTTGTCGCCCATCAATTTTGC‑CAATGAACACATTGCCAAGCCCTTCTTCCCGCGCATAACGATGACGATACATCGATAATGCGGCAGTCCCGCGGCTGACCGCACATGCTGCCGCCAGCGATGc                                                                                                                                                                                                   >  1:354625/1‑251 (MQ=255)                                           gCGTACGTTTGAGAAGCTGACCGAGGATGAAAATCGCCACCATCGTCACCACCATAGCAGCCACACCATGCATACCCGGCAATAATACGGCAGCAAAAATAGCGGCTAAGCCGAGTGTGACGCAGGTTTGTCGCCCATCAATTTTGC‑CAATGAACACATTGCCAAGCCCTTCTTCCCGCGCATAACGATGACGATACATCGATAATGCGGCAGTCCCGCGGCTGACCGCACATGCTGCCGCCAGCGATGCa                                                                                                                                                                                                  <  1:323769/251‑1 (MQ=255)                                                cGTTTGAGAAGCTGACCGAGGATGAAAATCGCCACCATCGTCACCACCATAGCAGCCACACCATGCATACCCGGCAATAATACGGCAGCAAAAATAGCGGCTAAGCCGAGTGTGACGCAGGTTTGTCGCCCATCAATTTTGC‑CAATGAACACATTGCCAAGCCCTTCTTCCCGCGCATAACGATGACGATACATCGATAATGCGGCAGTCCCGCGGCTGACCGCACATGCTGCCGCCAGCGATGCAAGGAt                                                                                                                                                                                             <  1:354022/251‑1 (MQ=255)                                                 ggttGAGAAGCTGACCGAGGATGAAAATCGCCACCATCGTCACCACCATAGCAGCCACACCATGCATACCCGGCAATAATACGGCAGCAAAAATAGCGGCTAAGCCGAGTGTGACGCAGGTTTGTCGCCCATCAATTTTGC‑CAATGAACACATTGCCAAGCCCTTCTTCCCGCGCATAACGATGACGATACATCGATAATGCGGCAGTCCCGCGGCTGACCGCACATGCTGCCGCCAGCGATGCAAGGATc                                                                                                                                                                                            <  2:248649/249‑1 (MQ=255)                                                                                 cccaTCGTCACCACCATAGCAGCCACACCATGCATACCCGGCAATAATACGGCAGCAAAAATAGCGGCTAAGCCGAGTGTGACGCAGGTTTGTCGCCCATCAATTTTGC‑CAATGAACACATTGCCAAGCCCTTCTTCCCGCGCATAACGATGACGATACATCGATAATGCGGCAGTCCCGCGGCTGACCGCACATGCTGCCGCCAGCGATGCAAGGATCGATTCGCCACGCAGAGCCAACTCGCTTAacac                                                                                                                                                            <  1:56654/250‑1 (MQ=255)                                                                                       gTCACCACCATAGCAGCCACACCATGCATACCCGGCAATAATACGGCAGCAAAAATAGCGGCTAAGCCGAGTGTGACGCAGGTTTGTCGCCCATCAATTTTGC‑CAATGAACACATTGCCAAGCCCTTCTTCCCGCGCATAACGATGACGATACATCGATAATGCGGCAGTCCCGCGGCTGACCGCACATGCTGCCGCCAGCGATGCAAGGATCGATTCGCCACGCAGAGCCAACTCGCTTAACACCAGAAt                                                                                                                                                      >  1:189185/1‑251 (MQ=255)                                                                                                  agcagcCACACCATGCATACCCGGCAATAATACGGCAGCAAAAATAGCGGCTAAGCCGAGTGTGACGCAGGTTTGTCGCCCATCAATTTTGC‑CAATGAACACATTGCCAAGCCCTTCTTCCCGCGCATAACGATGACGATACATCGATAATGCGGCAGTCCCGCGGCTGACCGCACATGCTGCCGCCAGCGATGCAAGGATCGATTCGCCACGCAGAGCCAACTCGCTTAACACCAGAATCTTTGCCAGt                                                                                                                                            >  1:159232/1‑250 (MQ=255)                                                                                                     agcCACACCATGCATACCCGGCAATAATACGGCAGCAAAAATAGCGGCTAAGCCGAGTGGGACGCAGGTTTGTCGCCCATCAATTTTGA‑CAATGAACACATTGCCAAGCCCTTCTTCCCGCGCATAACGATGACGATACATCGATAATGCGGCAGTCCCGCGGCTGACCGCACATGCTGCCGCCAGCGATGCAAGGATCGATTCGCCACGCAGAGCCAACTCGCTTAACACCAGAATCTTTGCCAGtccc                                                                                                                                         >  1:31471/1‑247 (MQ=255)                                                                                                      gcCACACCATGCATACCCGGCAATAATACGGCAGCAAAAATAGCGGCTAAGCCGAGTGTGACGCAGGTTTGTCGCCCATCAATTTTGC‑CAATGAACACATTGCCAAGCCCTTCTTCCCGCGCATAACGATGACGATCCATCGATAATGCGGCAGTCCCGCGGCTGACCGCACATGCTGCCGCCAGCGATGCAAGGATCGATTCGCCACGCAGAGCCAACTCGCTTAACACCAGAATCTTTGCCAGTAcccc                                                                                                                                       <  2:283599/251‑3 (MQ=255)                                                                                                                 cATACCCGGCAATAATACGGCAGCAAAAATAGCGGCTAAGCCGAGTGTGACGCAGGTTTGTCGCCCATCAATTTTGC‑CAATGAACACATTGCCAAGCCCTTCTTCCCGCGCATAACGATGACGATACATCGATAATGCGGCAGTCCCGCGGCTGACCGCACATGCTGCCGCCAGCGATGCAAGGATCGATTCGCCACGCAGAGCCAACTCGCTTAACACCAGAATCTTTGCCAGTACCACAAAAATCAATg                                                                                                                            >  1:7803/1‑251 (MQ=255)                                                                                                                  aTACCCGGCAATAATACGGCAGCAAAAATAGCGGCTAAGCCGAGTGTGACGCAGGTTTGTCGCCCATCAATTTTGC‑CAATGAACACATTGCCAAGCCCTTCTTCCCGCGCATAACGATGACGATACATCGATAATGCGGCAGTCCCGCGGCTGACCGCACATGCTGCCGCCAGCGATGCAAGGATCGATTCGCCACGCAGAGCCAACTCGCTTAACACCAGAATCTTTGCCAGTACCACAAAAATCAATGc                                                                                                                           >  1:346181/1‑251 (MQ=255)                                                                                                                                   ggCAGCAAAAATAGCGGCTAAGCCGAGTGTGACGCAGGTTTGTCGCCCATCAATTTTGC‑CAATGAACACATTGCCAAGCCCTTCTTCCCGCGCATAACGATGACGATACATCGATAATGCGGCAGTCCCGCGGCTGACCGCACATGCTGCCGCCAGCGATGCAAGGATCGATTCGCCACGCAGAGCCAACTCGCTTAACACCAGAATCTTTGCCAGTACCACAAAAATCAATGCCAGACCGCCGTGGGTGc                                                                                                          >  2:262684/1‑251 (MQ=255)                                                                                                                                        caAAAATAGCGGCTAAGCCGAGACTGACGCAGGTCTGTCGCCCATAAATATTGC‑CAATGAACACAATGCCAAGACCTTCCTCAAGCGAATAACGATGAAGATAAATAGATAATGCGGCAGTCCCGCGGCTGACCGCACATGCTGCCGCCAGCGATGCAAGGATCGATTCGCCACGAAGAGACAACTCGCTTAACACCAGAATCTTTGCCAGTACCACAAAAATCAATGCCAGACCGCCGTGGGTGCATaa                                                                                                      >  1:349238/1‑250 (MQ=255)                                                                                                                                             aTAGCGGCTAAGCCGAGTGTGACGCAGGTTTGTCGCCCATCAATTTTGC‑CAATGAACACATTGCCAAGCCCTTCTTCCCGCGCATAACGATGACGATACATCGATAATGCGGCAGTCCCGCGGCTGACCGCACATGCTGCCGCCAGCGATGCAAGGATCGATTCGCCCCGCAGAGCCAACTCGCTTAACACCAGAATCTTTGCCAGTACCACAAAAATCAATGCCAGACCGCCGTGGGTGCCTAAACGGCt                                                                                                >  1:286532/1‑251 (MQ=255)                                                                                                                                                     tAAGCCGAGTGTGACGCAGGTTTGTCGCCCATCAATTTTGC‑CAATGAACACATTGCCAAGCCCTTCTTCCCGCGCATAACGATGACGATACATCGATAATGCGGCAGTCCCGCGGCTGACCGCACATGCTGCCGCCAGCGATGCAAGGATCGATTCGCCACGCAGAGCCAACTCGCTTAACACCAGAATCTTTGCCAGTACCACAAAAATCAATGCCAGACCGCCGTGGGTGCCTAAACGGCTATCACGc                                                                                         >  1:77822/1‑250 (MQ=255)                                                                                                                                                        gCCGAGTGTGACGCAGGTTTGTCGCCCATCAATTTTGC‑CAATGAACACATTGCCAAGCCCTTCTTCCCGCGCATAACGATGACGATACATCGATAATGCGGCAGTCCCGCGGCTGACCGCACATGCTGCCGCCAGCGATGCAAGGATCGATTCGCCACGCAGAGCCAACTCGCTTAACACCAGAATCTTTGCCAGTACCACAAAAATCAATGCCAGACCGCCGTGGGTGCCTAAACGGCTATCACGCATg                                                                                      >  1:353389/1‑250 (MQ=255)                                                                                                                                                            aGTGTGACGCAGGTTTGTCGCCCATCAATTTTGC‑CAATGAACACATTGCCAAGCCCTTCTTCCCGCGCATAACGATGACGATACATCGATAATGCGGCAGTCCCGCGGCTGACCGCACATGCTGCCGCCAGCGATGCAAGGATCGATTCGCCACGCAGAGCCAACTCGCTTAACACCAGAATCTTTGCCAGTACCACAAAAATCAATGCCAGACCGCCGTGGGTGCCTAAACGGCTATCACGCATGATTTc                                                                                 >  1:174825/1‑251 (MQ=255)                                                                                                                                                                    gcacgggtGTCGCCCATCAATTT‑GCTCAATGCACACATTGACAAGCCCTTCTTCCCGGGTATTGCGATTACGATATATCGATAATGCGGCATTCCCGCGGCTGACCGCACATGCTGCCGCCAGCGATGCACGGATCGATTCGCCACGCTGAGCCAACTCGCTTAACACCAGAATCTTTGCCAGTCCCTCAAAAATAATTTCCAGGCCGCCGTGGGTGCCTAAACGGCTATCACGCATGATTTCCAACATGc                                                                         <  2:361288/244‑1 (MQ=255)                                                                                                                                                                        gTTTGTCGCCCATCAATTTTGC‑CAATGAACACATTGCCAAGCCCTTCTTCCCGCGCATAACGATGACGATACATCGATAATGCGGCAGTCCCGCGGCTGACCGCACATGCTGCCGCCAGCGATGCAAGGATCGATTCGCCACGCAGAGCCAACTCGCTTAACACCAGAATCTTTGCCAGTACCACAAAAATCAATGCCAGACCGCCGTGGGTGCCTAAACGGCTATCACGCATGATTTCCAACATGCGATc                                                                     <  1:349909/251‑1 (MQ=255)                                                                                                                                                                        gTTTGTCGCCCATCAATTTTGC‑CAATGAACACATTGCCAAGCCCTTCTTCCCGCGCATAACGATGACGATACATCGATAATGCGGCAGTCCCGCGGCTGACCGCACATGCTGCCGCCAGCGATGCAAGGATCGATTCGCCACGCAGAGCCAACTCGCTTAACACCAGAATCTTTGCCAGTACCACAAAAATCAATGCCAGACCGCCGTGGGGGCCTAAACGGCTATCACGCATGATTTCCAACATGCGATc                                                                     <  1:220998/251‑1 (MQ=255)                                                                                                                                                                            gTCGCCCATCAATTTTGC‑CAATGAACACATTGCCAAGCCCTTCTTCCCGCGCATAACGATGACGATACATCGATAATGCGGCAGTCCCGCGGCTGACCGCACATGCTGCCGCCAGCGATGCAAGGATCGATTCGCCACGCAGAGCCAACTCGCTTAACACCAGAATCTTTGCCAGTACCACAAAAATCAATGCCAGACCGCCGTGGGTGCCTAAACGGCTATCACGCATGATTTCCAACATGCGATcgcg                                                                  <  2:291899/250‑1 (MQ=255)                                                                                                                                                                              cGCCCATCAATTTTGC‑CAATGAACACATTGCCAAGCCCTTCTTCCCGCGCATAACGATGACGATACATCGATAATGCGGCAGTCCCGCGGCTGACCGCACATGCTGCCGCCAGCGATGCAAGGATCGATTCGCCACGCAGAGCCAACTCGCTTAACACCAGAATCTTTGCCAGTACCACAAAAATCAATGCCAGACCGCCGTGGGTGCCTAAACGGCTATCACGCATGATTTCCAACATGCGATCGCGGCt                                                               <  1:230740/251‑1 (MQ=255)                                                                                                                                                                                 ccATCAATTTTGC‑CAATGAACACATTGCCAAGCCCTTCTTCCCGCGCATAACGATGACGATACATCGATAATGCGGCAGTCCCGCGGCTGACCGCACATGCTGCCGCCAGCGATGCAAGGATCGATTCGCCACGCAGAGCCAACTCGCTTAACACCAGAATCTTTGCCAGTACCACAAAAATCAATGCCAGACCGCCGTGGGTGCCTAAACGGCTATCACGCATGATTTCCAACATGCGATCGCGGCTACg                                                            >  1:127862/1‑251 (MQ=255)                                                                                                                                                                                  cATCAATTTTGC‑CAATGAACACATTGCCAAGCCCTTCTTCCCGCGCATAACGATGACGATACATCGATAATGCGGCAGTCCCGCGGCTGACCGCACATGCTGCCGCCAGCGATGCAAGGATCGATTCGCCACGCAGAGCCAACTCGCTTAACACCAGAATCTTTGCCAGTACCACAAAAATCAATGCCAGACCGCCGTGGGTGCCTAAACGGCTATCACGCATGATTTCCAACATGCGATCGCGGCTACGt                                                           <  1:37906/251‑1 (MQ=255)                                                                                                                                                                                   aTCAATTTTGC‑CAATGAACACATTGCCAAGCCCTTCTTCCCGCGCATAACGATGACGATACATCGATAATGCGGCAGTCCCGCGGCTGACCGCACATGCTGCCGCCAGCGATGCAAGGATCGATTCGCCCCGCAGAGCCAACTCGCTTAACACCAGAATCTTTGCCAGTACCACACAAATCAATGCCAGACCGCCGTGGGTGCCTAAACGGCTATCACGCATGATTTCCAACATGCGATCGCGGCTACGTg                                                          >  1:48502/1‑251 (MQ=255)                                                                                                                                                                                   aTCAATTTTGC‑CAATGAACACATTGCCAAGCCCTTCTTCCCGCGCATAACGATGACGATACATCGATAATGCGGCAGTCCCGCGGCTGACCGCACATGCTGCCGCCAGCGATGCAAGGATCGATTCGCCACGCAGAGCCAACTCGCTTAACACCAGAATCTTTGCCAGTACCACAAAAATCAATGCCAGACCGCCGTGGGTGCCTAAACGGCTATCACGCATGATTTCCAACATGCGATCGCGGCTACGTg                                                          <  1:234420/251‑1 (MQ=255)                                                                                                                                                                                                    aaCACATTGCCAAGCCCTTCTTCCCGCGCATAACGATGACGATACATCGATAATGCGGCAGTCCCGCGGCTGACCGCACATGCTGCCGCCAGCGATGCAAGGATCGATTCGCCACGCAGAGCCAACTCGCTTAACACCAGAATCTTTGCCAGTACCACAAAAATCAATGCCAGACCGCCGTGGGTGCCTAAACGGCTATCACGCATGATTTCCAACATGCGATCGCGGCTACGTGCAGAAAATAAGCCATc                                          >  2:219049/1‑251 (MQ=255)                                                                                                                                                                                                          ttGCCAAGCCCTTCTTCCCGCGCATAACGATGACGATACATCGATAATGCGGCAGTCCCGCGGCTGACCGCACATGCTGCCGCCAGCGATGCAAGGATCGATTCGCCACGCAGAGCCAACTCGCTTAACACCAGAATCTTTGCCAGTACCACAAAAATCAATGCCAGACCGCCGTGGCTGCCTAAACGGCTATCACGCATGATTTCCAACATGCGATCGCGGCTACGTGCAGAAAATACGCCATCACAGGt                                    >  1:71412/1‑251 (MQ=255)                                                                                                                                                                                                                          cccGCGCATAACGATGACGATACATCGATAATGCGGCAGTCCCGCGGCTGACCGCACATGCTGCCGCCAGCGATGCAAGGATCGATTCGCCACGCAGAGCCAACTCGCTTAACACCAGAATCTTTGCCAGTACCACAAAAATCAATGCCAGACCGCCGTGGGTGCCTAAACGGCTATCACGCATGATTTCCAACATGCGATCGCGGCTACGTGCAGAAAATACGCCATCACAGGTATCAGCGAGGCCGTcc                    <  1:360681/251‑1 (MQ=255)                                                                                                                                                                                                                                      gATGACGATACATCGATAATGCGGCAGTCCCGCGGCTGACCGCACATGCTGCCGCCAGCGATGCAAGGATCGATTCGCCACGCAGAGCCAACTCGCTTAACACCAGAATCTTTGCCAGTACCACAAAAATCAATGCCAGACCGCCGTGGGTGCCTAAACGGCTATCACGCATGATTTCCAACATGCGATCGCGGCTACGTGCAGCAAATACGCCATCACAGGTATCAGCGAGGCCGTCCAGATGGAATccc        >  1:169848/1‑251 (MQ=255)                                                                                                                                                                                                                                            gATACATCGATAATGCGGCAGTCCCGCGGCTGACCGCACATGCTGCCGCCAGCGATGCAAGGATCGATTCGCCACGCAGAGCCAACTCGCTTAACACCAGAATCTTTGCCAGTACCACAAAAATCAATGCCAGACCGCCGTGGGTGCCTAAACGGCTATCACGCATGATTTCCAACATGCGATCGCGGCTACGTGCAGAAAATACGCCATCACAGGTATCAGCGAGGCCGTCCAGATGGAATCCCCCGGtc  >  2:290652/1‑251 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                   GCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCGAGGATGAAAATCGCCACCATCGTCACCACCATAGCAGCCACACCATGCATACCCGGCAATAATACGGCAGCAAAAATAGCGGCTAAGCCGAGTGTGACGCAGGTTTGTCGCCCATCAATTTTGC‑CAATGAACACATTGCCAAGCCCTTCTTCCCGCGCATAACGATGACGATACATCAATAATGCGGCAGTCCCGCGGCTGACCGCACATGCTGCCGCCAGCGATGCAAGGATCGATTCGCCACGCAGAGCCAACTCGCTTAACACCAGAATCTTTGCCAGTACCACAAAAATCAATGCCAGACCGCCGTGGGTGCCTAAACGGCTATCACGCATGATTTCCAACATGCGATCGCGGCTACGTGCAGAAAATACGCCATCACAGGTATCAGCGAGGCCGTCCAGATGGAATCCCCCGGTC  >  NC_000913/2064514‑2064998

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