Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,064,506	T→C	E239E (GAA→GAG)	cobS ←	cobalamin synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,064,506	0	T	C	100.0%	242.6 / NA	61	E239E (GAA→GAG)	cobS	cobalamin synthase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (36/25);  total (36/25)

CCCAGGCTACCAACAGGTTTGAGTAACCCGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGTTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCGAGGATGAAAATCGCCACCATCGTCACCACCATAGCAGCCACACCATGCATACCCGGCAATAATACGGCAGC  >  NC_000913/2064358‑2064649                                                                                                                                                     |                                                                                                                                                cccAGGCTACCAACAGGTTTGAGTAACCCGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTc                                                                                                                                               >  1:173451/1‑151 (MQ=255) cccAGGCTACCAACAGGTTTGAGTAACCCGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTc                                                                                                                                               >  1:421928/1‑151 (MQ=255)  ccAGGCTACCAACAGGTTTGAGTAACCCGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTCa                                                                                                                                              <  1:669565/151‑1 (MQ=255)   cAGGCTACCAACAGGTTTGAGTAACCCGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTCAc                                                                                                                                             >  2:991696/1‑151 (MQ=255)    aGGCTACCAACAGGTTTGAGTAACCCGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTCAcc                                                                                                                                            >  2:1106434/1‑151 (MQ=255)      gCTACCAACAGGTTTGAGTAACCCGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTCACCaa                                                                                                                                          <  1:235390/151‑1 (MQ=255)        tACCAACAGGTTTGAGTAACCCGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTCACCAAGc                                                                                                                                        >  1:406625/1‑151 (MQ=255)         aCCAACAGGTTTGAGTAACCCGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTCACCAAg                                                                                                                                         <  1:415654/149‑1 (MQ=255)              cAGGTTTGAGTAACCCGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATg                                                                                                                                  <  1:753907/151‑1 (MQ=255)              cAGGTTTGAGTAACCCGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATg                                                                                                                                  >  2:384186/1‑151 (MQ=255)                     gAGTAACCCGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTc                                                                                                                           <  1:869872/151‑1 (MQ=255)                      aGTAACCCGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTcc                                                                                                                          <  2:803250/151‑1 (MQ=255)                             cGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTa                                                                                                                   <  1:836418/151‑1 (MQ=255)                                 aaTATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATcac                                                                                                               <  1:993348/151‑3 (MQ=255)                                                  gcgACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGc                                                                                              >  1:214600/1‑151 (MQ=255)                                                   cgACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCg                                                                                             >  2:255336/1‑151 (MQ=255)                                                         atcGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGtt                                                                                       >  1:187635/3‑151 (MQ=255)                                                         aGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGt                                                                                        <  2:842159/150‑1 (MQ=255)                                                             gTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTgag                                                                                   <  2:441950/151‑1 (MQ=255)                                                               aGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTgaga                                                                                  >  1:592314/1‑150 (MQ=255)                                                               aGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGaa                                                                                 <  2:399816/151‑1 (MQ=255)                                                                   tCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTg                                                                             >  2:1164195/1‑151 (MQ=255)                                                                        ggCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCg                                                                         <  2:380195/150‑1 (MQ=255)                                                                             ggATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCGAGGATg                                                                   >  2:1041237/1‑151 (MQ=255)                                                                                  gTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCGAGGATGAAAAt                                                              >  1:609379/1‑151 (MQ=255)                                                                                  gTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCGAGGATGAAAAt                                                              >  1:144056/1‑151 (MQ=255)                                                                                    aTTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCGAGGATGAAAATCg                                                            >  2:823866/1‑151 (MQ=255)                                                                                     ttCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGA‑CAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCGAGGATGAAAATCGcc                                                          >  2:149356/1‑151 (MQ=255)                                                                                      tCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCGAGGATGAAAATCGcc                                                          >  1:858270/1‑151 (MQ=255)                                                                                        aGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCGAGGATGAAAATCGccac                                                        >  2:435756/1‑151 (MQ=255)                                                                                          tAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCGAGGATGAAAATCGccacc                                                       >  1:506986/1‑150 (MQ=255)                                                                                           aaaTCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGA‑CAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCGAGGATGAAAATCGCCACCATc                                                    <  1:149354/151‑1 (MQ=255)                                                                                            aaTCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCGAGGATGAAAATCGCCACCATc                                                    <  2:552335/151‑1 (MQ=255)                                                                                              tCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCGAGGATGAAAATCGCCACCATCGt                                                  >  1:127456/1‑151 (MQ=255)                                                                                               cGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCGAGGATGAAAATCGCCACCATCGTc                                                 >  2:663468/1‑151 (MQ=255)                                                                                                ggCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCGAGGATGAAAATCGCCACCATCGTca                                                <  2:406628/151‑1 (MQ=255)                                                                                                      tATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCGAGGATGAAAATCGCCACCATCGTCaccacc                                           >  1:720631/1‑150 (MQ=255)                                                                                                       aTTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCGAGGATGAAAATCGCCACCATCGTCACCACCAt                                         >  1:858834/1‑151 (MQ=255)                                                                                                       aTTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTCACCAAGCGCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCGAGGATGAAAATCGCCACCATCGTCACCACCAt                                         <  1:375242/151‑1 (MQ=255)                                                                                                         ttGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCGAGGATGAAAATCGCCACCATCGTCACCACCATa                                        <  2:909706/150‑1 (MQ=255)                                                                                                          tGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCGAGGATGAAAATCGCCACCATCGTCACCACCATag                                       <  1:677859/150‑1 (MQ=255)                                                                                                           gCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCACGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCGAGGATGAAAATCGCCACCATCGTCACCACCATagc                                      >  1:710476/1‑150 (MQ=255)                                                                                                            cATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCGAGGATGAAAATCGCCACCATCGTCACCACCATagc                                      >  1:256561/1‑149 (MQ=255)                                                                                                               ctcACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCGAGGATGAAAATCGCCACCATCGTCACCACCATAGCAGcc                                  >  2:750086/1‑150 (MQ=255)                                                                                                                   cGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTCTGGCTGCTCCCAGCGTACCGCCCGTTTGCCCGCCCCGCGTACGTTTGAGAAGCTGACCGAGGATGAAAATCGCCACCATCGTCACCACCATAGAAGCCACAcc                             >  1:289037/1‑151 (MQ=255)                                                                                                                      ctctTTATAGCAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCGAGGATGAAAATCGCCACCATCGTCACCACCATAGCAGCCACACCATg                          >  1:1183294/1‑151 (MQ=255)                                                                                                                          ttATAGCAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCGAGGATGAAAATCGCCACCATCGTCACCACCATAGCAGCCACACCATGCAt                       >  1:1084687/1‑150 (MQ=255)                                                                                                                          ttATAGCAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCGAGGATGAAAATCGCCACCATCGTCACCACCATAGCAGCCACACCATGCAt                       >  2:705654/1‑150 (MQ=255)                                                                                                                            ataGCAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCGAGGATGAAAATCGCCACCATCGTCACCACCATAGCAGCCACACCATGCATAcc                    >  2:114461/1‑151 (MQ=255)                                                                                                                            ataGCAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCGAGGATGAAAATCGCCACCATCGTCACCACCATAGCAGCCACACCATGCATAcc                    <  2:386767/151‑1 (MQ=255)                                                                                                                              aGCAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCGAGGATGAAAATCGCCACCATCGTCACCACCATAGCAGCCACACCATGCATACCCg                  >  2:541606/1‑151 (MQ=255)                                                                                                                                cAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCGAGGATGAAAATCGCCACCATCGTCACCACCATAGCAGCCACACCATGCATACCCGGc                <  1:498294/151‑1 (MQ=255)                                                                                                                                   tGCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCGAGGATGAAAATCGCCACCATCGTCACCACCATAGCAGCCACACCATGCATACCCGGCaat             >  1:919467/1‑151 (MQ=255)                                                                                                                                    gCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCGAGGATGAAAATCGCCACCATCGTCACCACCATAGCAGCCACACCATGCATACCCGGCaata            <  1:723299/151‑1 (MQ=255)                                                                                                                                      cagcagGAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCGAGGATGAAAATCGCCACCATCGTCACCACCATAGCAGCCACACCATGCATACCCGGCaataat          >  2:117514/1‑151 (MQ=255)                                                                                                                                         cagGAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCGAGGATGAAAATCGCCACCATCGTCACCACCATAGCAGCCACACCATGCATACCCGGCAATAATACg       <  2:706564/151‑1 (MQ=255)                                                                                                                                           ggAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCGAGGATGAAAATCGCCACCATCGTCACCACCATAGCAGCCACACCATGCATACCCGGCAATAATACGgc     <  1:797619/151‑1 (MQ=255)                                                                                                                                            gAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCGAGGATGAAAATCGCCACCATCGTCACCACCATAGCAGCCACACCATGCATACCCGGCAATAATACGgca    <  2:653118/151‑1 (MQ=255)                                                                                                                                              aGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCGAGGATGAAAATCGCCACCATCGTCACCACCATAGCAGCCACACCATGCATACCCGGCAATAATACGgcagc  <  2:542530/151‑1 (MQ=255)                                                                                                                                               gACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCGAGGATGAAAATCGCCACCATCGTCACCACCATAGCAGCCACACCATGCATACCCGGCAATAATACGgcagc  >  1:408640/1‑150 (MQ=255)                                                                                                                                               gACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGGATCGCCCGTTTGCCCGCCCAGTGTACGTTTGAGAAGCTGACCGAGGATGAAAATCGCCACCATCGTCACCACCATAGCAGCCACACCATGCCTGCACCCCAATAATACgaaagc  >  1:850300/1‑145 (MQ=255)                                                                                                                                                     |                                                                                                                                                CCCAGGCTACCAACAGGTTTGAGTAACCCGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGTTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCGAGGATGAAAATCGCCACCATCGTCACCACCATAGCAGCCACACCATGCATACCCGGCAATAATACGGCAGC  >  NC_000913/2064358‑2064649

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