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_0009132,064,5060TC100.0% 223.6 / NA 61E239E (GAA→GAGcobScobalamin synthase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (34/27);  total (34/27)

GCTACCAACAGGTTTGAGTAACCCGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGTTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCGAGGATGAAAATCGCCACCATCGTCACCACCATAGCAGCCACACCATGCATACCCGGCAATAATACGGCAGC  >  NC_000913/2064363‑2064649
                                                                                                                                               |                                                                                                                                               
gCTACCAACAGGTTTGAGTAACCCGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTCACCaa                                                                                                                                          <  2:612225/151‑1 (MQ=255)
gCTACCAACAGGTTTGAGTAACCCGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTCACCa                                                                                                                                           >  1:756399/1‑150 (MQ=255)
 cTACCAACAGGTTTGAGTAACCCGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTCACCa                                                                                                                                           >  1:1024719/1‑149 (MQ=255)
    ccAACAGGTTTGAGTAACCCGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTCACCAAGctc                                                                                                                                      <  1:973779/151‑1 (MQ=255)
     cAACAGGTTTGAGTAACCCGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTCACCAAGctct                                                                                                                                     >  2:23111/1‑151 (MQ=255)
        cAGGTTTGAGTAACCCGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATg                                                                                                                                  >  2:983366/1‑151 (MQ=255)
        cAGGTTTGAGTAACCCGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATg                                                                                                                                  >  1:421169/1‑151 (MQ=255)
         aGGTTTGAGTAACCCGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATgg                                                                                                                                 <  1:393751/151‑1 (MQ=255)
                aGTAACCCGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCCGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTcc                                                                                                                          <  2:743117/151‑1 (MQ=255)
                aGTAACCCGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTcc                                                                                                                          >  2:792026/1‑151 (MQ=255)
                   aaCCCGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTccc                                                                                                                         >  1:516512/1‑149 (MQ=255)
                   aaCCCGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTccc                                                                                                                         >  1:110851/1‑149 (MQ=255)
                      ccGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGt                                                                                                                    <  1:311200/151‑1 (MQ=255)
                       cGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGTAAGAATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACAAGCTCACCAAGCTCTATGGCTGCTCCCAGCg                                                                                                                     >  2:677746/1‑149 (MQ=255)
                        gTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTAt                                                                                                                  <  2:496506/151‑1 (MQ=255)
                              atGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGccc                                                                                                             <  2:255492/150‑1 (MQ=255)
                               tGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGt                                                                                                           >  1:1055347/1‑151 (MQ=255)
                                gCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGtt                                                                                                          <  2:291049/151‑1 (MQ=255)
                                  cGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCANCAGGAAGACCAGCTCACCAAGCTCTATNGNTGCNNNCNGCGTAt                                                                                                                  <  2:256617/141‑1 (MQ=255)
                                  cGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCAAGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTg                                                                                                        >  2:256445/1‑151 (MQ=255)
                                           cgcgACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTgcccgccc                                                                                                 >  1:181951/1‑149 (MQ=255)
                                           cgcgACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTgcccgccc                                                                                                 >  2:583131/1‑149 (MQ=255)
                                            gcgACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGc                                                                                              <  2:653961/151‑1 (MQ=255)
                                                           aaTCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGc                                                                               >  1:751042/1‑151 (MQ=255)
                                                                aTGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGAcc                                                                          >  2:30276/1‑151 (MQ=255)
                                                                 tGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCg                                                                         >  1:838036/1‑151 (MQ=255)
                                                                      gggATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCGAgg                                                                      <  2:280117/149‑1 (MQ=255)
                                                                        gATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCGAGGATGa                                                                  <  2:769557/151‑1 (MQ=255)
                                                                              aTTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCGAGGATGAAAATCg                                                            >  1:136784/1‑151 (MQ=255)
                                                                                  aGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAg                                                                                <  1:65531/127‑1 (MQ=255)
                                                                                  aGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAg                                                                                >  2:65530/1‑127 (MQ=255)
                                                                                   gTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCGAGGATGAAAATCGccacc                                                       <  2:747371/151‑1 (MQ=255)
                                                                                    tAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCGAGGATGAAAATCGccacc                                                       >  1:967215/1‑150 (MQ=255)
                                                                                     aaaTCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCGAGGATGAAAATCGccacca                                                      <  2:421121/150‑1 (MQ=255)
                                                                                       aTCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCCGCGTACGTTTGAGAAGCTGACCGAGGATGAAAATCGCCACCATCg                                                   >  1:962035/1‑151 (MQ=255)
                                                                                       aTCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCGAGGATGAAAATCGCCACCATCg                                                   <  2:679956/151‑1 (MQ=255)
                                                                                          ggCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCGAGGATGAAAATCGCCACCATCGTca                                                <  2:488326/151‑1 (MQ=255)
                                                                                           gCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCGAGGATGAAAATCGCCACCATCGTcac                                               <  1:258714/151‑1 (MQ=255)
                                                                                            cAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCGAGGATGAAAATCGCCACCATCGTcacc                                              <  2:21845/151‑1 (MQ=255)
                                                                                                 aTTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCGAGGATGAAAATCGCCACCATCGTCACCACCAt                                         <  1:156335/151‑1 (MQ=255)
                                                                                                   tttcATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATAGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCGAGGATGAAAATCGCCACCATCGTCACCACCAt                                         >  2:619681/4‑149 (MQ=255)
                                                                                                    tGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCGAGGATGAAAATCGCCACCATCGTCACCACCATagc                                      >  1:429425/1‑151 (MQ=255)
                                                                                                          tcACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCGAGGATGAAAATCGCCACCATCGTCACCACCATAGCAGCcac                                >  1:1051371/1‑151 (MQ=255)
                                                                                                           cACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCGAGGATGAAAATCGCCACCATCGTCACCACCATAGCAGCcac                                >  1:199830/1‑150 (MQ=255)
                                                                                                               tctcTTTATAGCAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCGAGGATGAAAATCGCCACCATCGTCACCACCATAGCAGCCACACCa                            >  2:697295/1‑150 (MQ=255)
                                                                                                                    ttATAGCAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCGAGGATGAAAATCGCCACCATCGTCACCACCATAGCAGCCACACCATGCAt                       <  1:707830/150‑1 (MQ=255)
                                                                                                                     tataGCAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCGAGGATGAAAATCGCCACCATCGTCACCACCATAGCAGCCACACCATGCATAc                     >  2:191086/1‑151 (MQ=255)
                                                                                                                     tataGCAGAGCCAGCAGGAAGCCCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCGCACCGAGGATGAAAATCGCCACCATCGTCACCACCATAGCAGCCACACCATGCAt                       <  2:191448/149‑1 (MQ=255)
                                                                                                                       taGCAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCGAGGATGAAAATCGCCACCATCGTCACCACCATAGCAGCCACACCATGCATAccc                   >  1:619795/1‑151 (MQ=255)
                                                                                                                        aGCAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCGAGGATGAAAATCGCCACCATCGTCACCACCATAGCAGCCACACCATGCATAccc                   <  1:356886/150‑1 (MQ=255)
                                                                                                                        aGCAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCGAGGATGAAAATCGCCACCATCGTCACCACCATAGCAGCCACACCATGCATACCCg                  <  1:242611/151‑1 (MQ=255)
                                                                                                                         gCAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCGAGGATGAAAATCGCCACCATCGTCACCACCATAGCAGCCACACCATGCATACCCgg                 >  1:12932/1‑151 (MQ=255)
                                                                                                                          cAGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCGAGGATGAAAATCGCCACCATCGTCACCACCATAGCAGCCACACCATGCATACCCGGc                >  2:699719/1‑151 (MQ=255)
                                                                                                                           aGTGGCAACAGGGAAAACAGCTCACCAAGCTCTATTGCTGCTCCCCGCGGATCGCCCGGTTGCCCCCCCCGCGGACGTTTTAGAAACTGACCGAGGATTAAAAACGCCACCATCGTCACCACCATAGCAGCCCCACCCAGCATAacccg                 >  2:418723/1‑147 (MQ=255)
                                                                                                                           aGTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCGAGGATGAAAATCGCCACCATCGTCACCACCATAGCAGCCACACCATGCATACCCGGc                >  2:694984/1‑150 (MQ=255)
                                                                                                                            gTGCCAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCGAGGATGAAAATCGCCACCATCGTCACCACCATAGCAGCCACACCATGCATACCCGGCa               <  2:384670/150‑1 (MQ=255)
                                                                                                                               ccAGCAGGAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCGAGGATGAAAATCGCCACCATCGTCACCACCATAGCAGCCACACCATGCATACCCGGCaata            >  2:328076/1‑150 (MQ=255)
                                                                                                                                cagcagGAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCGAGGATGAAAATCGCCACCATCGTCACCACCATAGCAGCCACACCATGCATACCCGGCaataa           >  2:763676/1‑150 (MQ=255)
                                                                                                                                 agcagGAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCGAGGATGAAAATCGCCACCATCGTCACCACCATAGCAGCCACACCATGCATACCCGGCAataata         <  2:296989/151‑1 (MQ=255)
                                                                                                                                      gAAGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCGAGGATGAAAATAGCCACCATCGTCACCACCATAGCAGCCACACCCTGCATACCCGGCAATAATAc        >  2:163356/1‑147 (MQ=255)
                                                                                                                                        aGACCAGCTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCGAGGATGAAAATCGCCACCATCGTCACCACCATAGCAGCCACACCATGCATACCCGGCAATAATACGgcagc  <  2:806815/151‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                               
GCTACCAACAGGTTTGAGTAACCCGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGTTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCGAGGATGAAAATCGCCACCATCGTCACCACCATAGCAGCCACACCATGCATACCCGGCAATAATACGGCAGC  >  NC_000913/2064363‑2064649

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: