Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,791,400 T→C Q61Q (CAA→CAG btuD ← ATP‑binding component of vitamin B12 transport system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,791,4000TC100.0% 600.5 / NA 188Q61Q (CAA→CAGbtuDATP‑binding component of vitamin B12 transport system
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (110/78);  total (110/78)

Only 100 of 191 total aligned reads displayed.

CCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGTTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCA  >  NC_000913/1791310‑1791488
                                                                                          |                                                                                        
  aGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGcc                                                                                      >  2:954330/1‑93 (MQ=255)
  aGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGcc                                                                                      >  1:2797680/1‑93 (MQ=255)
    aCCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGcccc                                                                                    <  1:1442428/93‑1 (MQ=255)
      cGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCcgc                                                                                  >  1:4379727/1‑93 (MQ=255)
      cGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCcgc                                                                                  >  2:6030201/1‑93 (MQ=255)
       ggCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCcgcg                                                                                 >  2:3992050/1‑93 (MQ=255)
       ggCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCcgcg                                                                                 >  2:5800520/1‑93 (MQ=255)
        gCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGa                                                                                >  2:5027756/1‑93 (MQ=255)
        gCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGa                                                                                >  2:4962332/1‑93 (MQ=255)
          gTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAAc                                                                              <  2:824312/93‑1 (MQ=255)
          gTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAAc                                                                              <  1:5742420/93‑1 (MQ=255)
          gTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTGTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAAc                                                                              >  2:6296536/1‑93 (MQ=255)
           ttGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACt                                                                             >  1:6225777/1‑93 (MQ=255)
             gCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGa                                                                           <  2:5620006/93‑1 (MQ=255)
             gCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGa                                                                           <  2:3261370/93‑1 (MQ=255)
              cAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGaa                                                                          <  2:458043/93‑1 (MQ=255)
              cAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGaa                                                                          >  1:1502929/1‑93 (MQ=255)
                aaCGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATg                                                                        <  2:4276500/93‑1 (MQ=255)
                 aCGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGc                                                                       <  2:1983128/93‑1 (MQ=255)
                   ggcggcGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCtt                                                                     >  2:5059963/1‑93 (MQ=255)
                   ggcggcGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCtt                                                                     <  2:1234789/93‑1 (MQ=255)
                   ggcggcGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCtt                                                                     <  2:5875992/93‑1 (MQ=255)
                      ggcgTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTccc                                                                  >  2:366293/1‑93 (MQ=255)
                      ggcgTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTccc                                                                  <  1:4639333/93‑1 (MQ=255)
                       gcgTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCt                                                                 >  2:6115900/1‑93 (MQ=255)
                       gcgTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCt                                                                 >  2:5163601/1‑93 (MQ=255)
                         gTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTa                                                               <  1:1800939/93‑1 (MQ=255)
                         gTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTa                                                               <  1:55647/93‑1 (MQ=255)
                         gTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTa                                                               >  2:5641004/1‑93 (MQ=255)
                         gTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTa                                                               <  1:420535/93‑1 (MQ=255)
                          tCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTTAATGCTTCCCTTAc                                                              >  2:1681388/1‑93 (MQ=255)
                          tCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTGCCTTAc                                                              >  2:4716410/1‑93 (MQ=255)
                          tCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTAc                                                              >  2:1618369/1‑93 (MQ=255)
                          tCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTAc                                                              >  2:4814823/1‑93 (MQ=255)
                          tCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTAc                                                              >  2:3568450/1‑93 (MQ=255)
                           ctgctgTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTAcc                                                             >  2:2475613/1‑93 (MQ=255)
                            tgctgTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCg                                                            <  1:1091118/93‑1 (MQ=255)
                             gctgTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGc                                                           >  2:2206928/1‑93 (MQ=255)
                             gctgTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGc                                                           >  1:3482365/1‑93 (MQ=255)
                             gctgTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGc                                                           <  2:3767536/93‑1 (MQ=255)
                              ctgTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCt                                                          >  1:3579717/1‑93 (MQ=255)
                              ctgTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCt                                                          >  1:5936301/1‑93 (MQ=255)
                               tgttgtGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTg                                                         <  1:6163029/93‑1 (MQ=255)
                               tgttgtGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTg                                                         <  1:3288971/93‑1 (MQ=255)
                                   gtgAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCa                                                     <  1:3918699/93‑1 (MQ=255)
                                    tgAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCAt                                                    >  1:198381/1‑93 (MQ=255)
                                     gAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTAACGCTGGTCAtt                                                   >  1:164435/1‑93 (MQ=255)
                                      aaaGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTc                                                  >  1:2689595/1‑93 (MQ=255)
                                      aaaGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTc                                                  >  2:4217142/1‑93 (MQ=255)
                                        aGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGCATGCTTCCCTTACCGCTGGTCATTCCg                                                >  1:385600/1‑93 (MQ=255)
                                        aGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCg                                                <  1:5670625/93‑1 (MQ=255)
                                        aGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCg                                                <  1:4820571/93‑1 (MQ=255)
                                         gATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCgg                                               <  2:902554/93‑1 (MQ=255)
                                          aTAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGc                                              >  2:176347/1‑93 (MQ=255)
                                            aGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCa                                            <  1:1062601/93‑1 (MQ=255)
                                             ggCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCAt                                           >  1:4025708/1‑93 (MQ=255)
                                               cgcgATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTc                                         <  2:1047734/93‑1 (MQ=255)
                                                gcgATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTATCGCTGGTCATTCCGGCCAttag                                        >  2:4512680/1‑91 (MQ=255)
                                                gcgATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTcg                                        <  1:2437672/93‑1 (MQ=255)
                                                gcgATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTcg                                        >  2:1139705/1‑93 (MQ=255)
                                                  gATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTcgcg                                      >  2:4628278/1‑93 (MQ=255)
                                                  gATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTcgcg                                      >  1:3575560/1‑93 (MQ=255)
                                                    tGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGcc                                    >  2:3620367/1‑93 (MQ=255)
                                                      cAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCCTTCCGGCCATTCGCGCCAg                                  >  1:6008042/1‑93 (MQ=255)
                                                       aGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGt                                 >  2:861939/1‑93 (MQ=255)
                                                        gcgcGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTa                                <  2:5379222/93‑1 (MQ=255)
                                                        gcgcGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTa                                >  2:4774359/1‑93 (MQ=255)
                                                        gcgcGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTa                                <  2:4898357/93‑1 (MQ=255)
                                                        gcgcGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTa                                >  2:3008085/1‑93 (MQ=255)
                                                         cgcgAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTaa                               <  1:881308/93‑1 (MQ=255)
                                                          gcgAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAg                              >  2:3355610/1‑93 (MQ=255)
                                                           cgAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAgg                             <  1:4982896/93‑1 (MQ=255)
                                                           cgAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAgg                             <  1:3160108/93‑1 (MQ=255)
                                                            gAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGt                            <  2:5182871/93‑1 (MQ=255)
                                                            gAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGt                            >  1:3252207/1‑93 (MQ=255)
                                                             aGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTa                           >  2:2067819/1‑93 (MQ=255)
                                                             aGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTa                           >  1:5377103/1‑93 (MQ=255)
                                                             aGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTa                           <  1:2986285/93‑1 (MQ=255)
                                                              gTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTAc                          <  1:3625519/93‑1 (MQ=255)
                                                               ttttGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTAct                         <  1:4897812/93‑1 (MQ=255)
                                                               ttttGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGGAct                         >  2:3384718/1‑93 (MQ=255)
                                                                   gttgCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTa                     <  1:4475723/93‑1 (MQ=255)
                                                                   gttgCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTa                     <  1:3539367/93‑1 (MQ=255)
                                                                        ggACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCcgc                <  1:105025/93‑1 (MQ=255)
                                                                           ccATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGcc             <  1:1224766/93‑1 (MQ=255)
                                                                            cATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCa            <  1:296168/93‑1 (MQ=255)
                                                                             aTGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCAt           >  2:3862771/1‑93 (MQ=255)
                                                                             aTGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCAt           >  2:1912478/1‑93 (MQ=255)
                                                                             aTGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCAt           <  1:3474011/93‑1 (MQ=255)
                                                                             aTGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCAt           <  2:2360460/93‑1 (MQ=255)
                                                                              tGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCAtt          <  1:5881876/93‑1 (MQ=255)
                                                                              tGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCAtt          >  2:2865179/1‑93 (MQ=255)
                                                                               gCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTc         <  2:2770525/93‑1 (MQ=255)
                                                                               gCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTc         >  2:3315014/1‑93 (MQ=255)
                                                                                cTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCg        <  1:426885/93‑1 (MQ=255)
                                                                                  tCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGc      <  1:3764440/93‑1 (MQ=255)
                                                                                    cAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGccc    >  1:3204433/1‑93 (MQ=255)
                                                                                     aGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGcccc   >  2:2706627/1‑93 (MQ=255)
                                                                                      gTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCcca  >  1:682017/1‑93 (MQ=255)
                                                                                      gTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCcca  >  1:6042758/1‑93 (MQ=255)
                                                                                          |                                                                                        
CCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGTTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCA  >  NC_000913/1791310‑1791488

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: