Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,011,430 A→G H263R (CAT→CGT)  metE → tetrahydropteroyltriglutamate methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,011,4300AG100.0% 843.8 / NA 258H263R (CAT→CGT) metEtetrahydropteroyltriglutamate methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (140/118);  total (140/118)

Only 100 of 267 total aligned reads displayed.

AGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCATGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTCTGATCAATG  >  NC_000913/4011340‑4011519
                                                                                          |                                                                                         
aGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTg                                                                                         >  1:5891/1‑93 (MQ=255)
 gggACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCgtgt                                                                                        >  2:3881267/1‑93 (MQ=255)
 gggACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCgtgt                                                                                        >  2:468465/1‑93 (MQ=255)
   gACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTg                                                                                      >  1:1879460/1‑93 (MQ=255)
    aCAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGa                                                                                     <  2:517744/93‑1 (MQ=255)
    aCAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGa                                                                                     >  2:1801697/1‑93 (MQ=255)
    aCAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCCCGACACGATAACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGa                                                                                     <  1:1685397/93‑1 (MQ=255)
      aGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGAcc                                                                                   <  1:845031/93‑1 (MQ=255)
      aGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGAcc                                                                                   <  2:6182204/93‑1 (MQ=255)
      aGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGAcc                                                                                   >  2:1532368/1‑93 (MQ=255)
       ggTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCt                                                                                  >  1:5221119/1‑93 (MQ=255)
         tGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCg                                                                                >  2:4329868/1‑93 (MQ=255)
           aagctgctgctgACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTa                                                                              <  2:4770912/90‑1 (MQ=255)
            aaCTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTAc                                                                             >  2:523416/1‑93 (MQ=255)
            aaCTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTAc                                                                             >  2:3180733/1‑93 (MQ=255)
            aaCTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTAc                                                                             >  2:1911344/1‑93 (MQ=255)
            aaCTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTAc                                                                             >  1:2730792/1‑93 (MQ=255)
             aCTGCTGCTGACCCCCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACa                                                                            <  1:3235248/93‑1 (MQ=255)
              ctgctgctgACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACAt                                                                           <  1:1367055/93‑1 (MQ=255)
               tgctgctgACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATg                                                                          <  2:5183018/93‑1 (MQ=255)
                 ctgctgACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGCTCTGCGTGTTGACCTCGTACATGGt                                                                        >  2:6481923/1‑93 (MQ=255)
                  tgctgACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTa                                                                       >  1:4342303/1‑93 (MQ=255)
                   gctgACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTaa                                                                      <  1:5475922/93‑1 (MQ=255)
                    ctgACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTaaa                                                                     <  1:3464833/93‑1 (MQ=255)
                      gACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGa                                                                   >  1:6303139/1‑93 (MQ=255)
                       accaccTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGAt                                                                  <  1:1229687/93‑1 (MQ=255)
                       accaccTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGAt                                                                  <  2:6431759/93‑1 (MQ=255)
                       accaccTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGAt                                                                  <  2:5873402/93‑1 (MQ=255)
                       accaccTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGAt                                                                  >  2:4232052/1‑93 (MQ=255)
                        ccaccTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATg                                                                 >  2:4634203/1‑93 (MQ=255)
                        ccaccTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATg                                                                 >  1:4742224/1‑93 (MQ=255)
                        ccaccTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATg                                                                 >  2:6395462/1‑93 (MQ=255)
                         caccTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGa                                                                >  1:3686615/1‑93 (MQ=255)
                         caccTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGa                                                                >  1:1529800/1‑93 (MQ=255)
                          cccTATTTTGAAGGCCTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGAc                                                               <  1:6030422/92‑1 (MQ=255)
                          accTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGAc                                                               <  1:6126008/93‑1 (MQ=255)
                          accTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTACGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGCACATGGTAAAGATGAc                                                               >  1:2857271/1‑93 (MQ=255)
                            cTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGt                                                             <  1:1297046/93‑1 (MQ=255)
                             tATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGtt                                                            >  2:2624755/1‑93 (MQ=255)
                              aTTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTg                                                           >  1:5275136/1‑93 (MQ=255)
                               ttttGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGc                                                          <  2:308868/93‑1 (MQ=255)
                               ttttGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGc                                                          >  2:3745780/1‑93 (MQ=255)
                                 ttGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTg                                                        <  1:1604582/93‑1 (MQ=255)
                                  tGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGa                                                       >  2:3809114/1‑93 (MQ=255)
                                    aaGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAAc                                                     <  2:5385462/93‑1 (MQ=255)
                                    aaGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAAc                                                     <  2:5148108/93‑1 (MQ=255)
                                     aGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACt                                                    <  1:1001944/93‑1 (MQ=255)
                                      ggCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTg                                                   <  2:2919755/93‑1 (MQ=255)
                                      ggCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTg                                                   >  2:5029492/1‑93 (MQ=255)
                                       gCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGc                                                  <  1:2072272/93‑1 (MQ=255)
                                       gCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGc                                                  <  2:6315979/93‑1 (MQ=255)
                                         gTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGcac                                                >  1:6265999/1‑93 (MQ=255)
                                         gTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGcac                                                >  1:4838501/1‑93 (MQ=255)
                                           aaCGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACaa                                              <  1:3057283/93‑1 (MQ=255)
                                           aaCGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACaa                                              >  1:570193/1‑93 (MQ=255)
                                           aaCGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACaa                                              <  2:1389728/93‑1 (MQ=255)
                                            aCGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAg                                             >  2:3301917/1‑93 (MQ=255)
                                             cGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAgc                                            >  2:187033/1‑93 (MQ=255)
                                             cGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAgc                                            <  2:5475557/93‑1 (MQ=255)
                                               ccAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAgcgc                                          >  2:2234663/1‑93 (MQ=255)
                                                cAACTCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCgcc                                         <  1:1415702/93‑1 (MQ=255)
                                                cAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCgcc                                         <  1:1582587/93‑1 (MQ=255)
                                                cAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCgcc                                         >  2:3385556/1‑93 (MQ=255)
                                                cAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCgcc                                         <  2:918207/93‑1 (MQ=255)
                                                cAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCgcc                                         >  2:3067888/1‑93 (MQ=255)
                                                 aaaTCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCgcct                                        <  1:3336264/93‑1 (MQ=255)
                                                  aaTCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCgcctg                                       >  2:1627926/1‑93 (MQ=255)
                                                   aTCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCgcctgc                                      >  2:1997787/1‑93 (MQ=255)
                                                   aTCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCgcctgc                                      >  2:6706787/1‑93 (MQ=255)
                                                     ctcGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCgcctgcct                                    <  1:1480507/93‑1 (MQ=255)
                                                       cGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTc                                  >  2:5677735/1‑93 (MQ=255)
                                                       cGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTc                                  <  1:3642063/93‑1 (MQ=255)
                                                         acacGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTg                                >  1:4031923/1‑93 (MQ=255)
                                                          cacGATTAGTGCGCTGCCTCTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGa                               >  1:6568163/1‑93 (MQ=255)
                                                          cacGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGa                               >  2:5511065/1‑93 (MQ=255)
                                                           acGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGAc                              <  2:209805/93‑1 (MQ=255)
                                                             gATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTg                            >  2:4048710/1‑93 (MQ=255)
                                                               ttACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGt                          >  2:3503727/1‑93 (MQ=255)
                                                               ttACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGt                          >  1:4956830/1‑93 (MQ=255)
                                                                tACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGtt                         <  1:6515584/93‑1 (MQ=255)
                                                                   tGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCt                      <  1:4583133/93‑1 (MQ=255)
                                                                      gcTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTc                   <  2:5727182/93‑1 (MQ=255)
                                                                       cTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCt                  >  2:1486823/1‑93 (MQ=255)
                                                                       cTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCt                  >  1:937185/1‑93 (MQ=255)
                                                                          ccTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCg               >  2:3703984/1‑93 (MQ=255)
                                                                            tGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCggg             <  1:3861484/93‑1 (MQ=255)
                                                                            tGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCggg             <  1:3186534/93‑1 (MQ=255)
                                                                            tGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCggg             <  1:2983414/93‑1 (MQ=255)
                                                                             gTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGt            >  1:103627/1‑93 (MQ=255)
                                                                             gTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGt            <  1:1264743/93‑1 (MQ=255)
                                                                              ttCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTc           >  2:2755533/1‑93 (MQ=255)
                                                                              ttCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTc           >  1:4275138/1‑93 (MQ=255)
                                                                              ttCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTc           >  1:3974306/1‑93 (MQ=255)
                                                                              ttCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTc           <  1:2627573/93‑1 (MQ=255)
                                                                                 aGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTCTGa        >  1:4728617/1‑93 (MQ=255)
                                                                                  gggTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTCTGAt       <  1:5385727/93‑1 (MQ=255)
                                                                                  gggTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTCTGAt       >  1:196357/1‑93 (MQ=255)
                                                                                     tCTGCGTGTTGACCTCGTACATGGGAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTCTGCTCaa    >  1:5189258/1‑93 (MQ=255)
                                                                                       tGCGTGTTGTGCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCGTGCCTTCTGACTGGTTGCTGTCTGCGGGTCTGATCAATg  <  1:751022/93‑1 (MQ=255)
                                                                                       tGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTCTGATCAATg  >  1:574349/1‑93 (MQ=255)
                                                                                          |                                                                                         
AGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCATGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTCTGATCAATG  >  NC_000913/4011340‑4011519

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