Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,140,222 T→G I33L (ATT→CTT)  hybF ← protein involved with the maturation of hydrogenases 1 and 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,140,2220TG100.0% 226.1 / NA 61I33L (ATT→CTT) hybFprotein involved with the maturation of hydrogenases 1 and 2
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (33/28);  total (33/28)

CCACCACCTGGCTGCAATCCCAGCACCAAGCCTGGGCGGGTTTATAGACGATATGTAAATCGCACCCTTGCGCCACCGTTCCGTGGCAGACAATTTCAAAACTAAAACGGACGGCGCTCTCCTCAACGCAGGAGAGCGCGCCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGT  >  NC_000913/3140078‑3140365
                                                                                                                                                |                                                                                                                                               
ccaccaCCTGGCTGCAATCCCAGCACCAAGCCTGGGCGGGTTTATAGACGATATGTAAATCGCACCCTTGCGCCACCGTTCCGTGGCAGACAATTTCAAAACTAAAACGGACGGCGCTCTCCTCAACGCAGGAGAGCGCGCCAAGTTc                                                                                                                                              <  2:411270/148‑1 (MQ=255)
   ccaccTGGCTGCAATCCCAGCACCAAGCCTGGGCGGGTTTATAGACGATATGTAAATCGCACCCTTGCGCCACCGTTCCGTGGCAGACAATTTCAAAACTAAAACGGACGGCGCTCTCCTCAACGCAGGAGAGCGCGCCAAGTTCCAGcc                                                                                                                                         <  1:673130/150‑1 (MQ=255)
   ccaccTGGCTGCAATCCCAGCACCAAGCCTGGGCGGGTTTATAGACGATATGTAAATCGCACCCTTGCGCCACCGTTCCGTGGCAGACAATTTCAAAACTAAAACGGACGGCGCTCTCCTCAACGCAGGAGAGCGCGCCAAGTTCCAGCca                                                                                                                                        <  2:711133/151‑1 (MQ=255)
    caccTGGCTGCAATCCCAGCACCAAGCCTGGGCGGGTTTATAGACGATATGTAAATCGCACCCTTGCGCCACCGTTCCGTGGCAGACAATTTCAAAACTAAAACGGACGGCGCTCTCCTCAACGCAGGAGAGCGCGCCAAGTTCCAGCcac                                                                                                                                       >  1:432947/1‑151 (MQ=255)
        tGGCTGCAATCCCAGCACCAAGCCTGGGCGGGTTTATAGACGATATGTAAATCGCACCCTTGCGCCACCGTTCCGTGGCAGACAATTTCAAAACTAAAACGGACGGCGCTCTCCTCAACGCAGGAGAGCGCGCCAAGTTCCAGCCACAcgg                                                                                                                                   >  2:407157/1‑151 (MQ=255)
          gCTGCAATCCCAGCACCAAGCCTGGGCGGGTTTATAGACGATATGTAAATCGCACCCTTGCGCCACCGTTCCGTGGCAGACAATTTCAAAACTAAAACGGACGGCGCTCTCCTCAACGCAGGAGAGCGCGCCAAGTTCCAGCCACAcggcg                                                                                                                                 <  1:324370/151‑1 (MQ=255)
             gCAATCCCAGCACCAAGCCTGGGCGGGTTTATAGACGATATGTAAATCGCACCCTTGCGCCACCGTTCCGTGGCAGACAATTTCAAAACTAAAACGGACGGCGCTCTCCTCAACGCAGGAGAGCGCGCCAAGTTCCAGCCACACGGCGGTg                                                                                                                              <  1:827292/151‑1 (MQ=255)
                  cccAGCACCAAGCCTGGGCGGGTTTATAGACGATATGTAAATCGCACCCTTGCGCCACCGTTCCGTGGCAGACAATTTCAAAACTAAAACGGACGGCGCTCTCCTCAACGCAGGAGAGCGCGCCAAGTTCCAGCCACACGGCGGTGAcgcg                                                                                                                         <  2:517563/151‑1 (MQ=255)
                       cACCAAGCCTGGGCGGGTTTATAGACGATATGTAAATCGCACCCTTGCGCCACCGTTCCGTGGCAGACAATTTCAAAACTAAAACGGACGGCGCTCTCCTCAACGCAGGAGAGCGCGCCAAGTTCCAGCCACACGGCGGTGACGCGCTTa                                                                                                                     >  1:371594/1‑150 (MQ=255)
                        aCCAAGCCTGGGCGGGTTTATAGACGATATGTAAATCGCACCCTTGCGCCACCGTTCCGTGGCAGACAATTTCAAAACTAAAACGGACGGCGCTCTCCTCAACGCAGGAGAGCGCGCCAAGTTCCAGCCACACGGCGGTGACGCGCTTaa                                                                                                                    >  2:534349/1‑150 (MQ=255)
                         ccAAGCCTGGGCGGGTTTATAGACGATATGTAAATCGCACCCTTGCGCCACCGTTCCGTGGCAGACAATTTCAAAACTAAAACGGACGGCGCTCTCCTCAACGCAGGAGAGCGCGCCAAGTTCCAGCCACACGGCGGTGACGCGCTTAAc                                                                                                                   >  2:615893/1‑150 (MQ=255)
                         ccAAGCCTGGGCGGGTTTATAGACGATATGTAAATCGCACCCTTGCGCCACCGTTCCGTGGCAGACAATTTCAAAACTAAAACGGACGGCGCTCTCCTCAACGCAGGAGAGCGCGCCAAGTTCCAGCCACACGGCGGTGACGCGCTTAACa                                                                                                                  <  2:454123/151‑1 (MQ=255)
                             gCCTGGGCGGGTTTATAGACGATATGTAAATCGCACCCTTGCGCCACCGTTCCGTGGCAGACAATTTCAAAACTAAAACGGACGGCGCTCTCCTCAACGCAGGAGAGCGCGCCAAGTTCCAGCCACACGGCGGTGACGCGCTTAACATCg                                                                                                               >  1:433792/1‑150 (MQ=255)
                              ccTGGGCGGGTTTATAGACGATATGTAAATCGCACCCTTGCGCCACCGTTCCGTGGCAGACAATTTCAAAACTAAAACGGACGGCGCTCTCCTCAACGCAGGAGAGCGCGCCAAGTTCCAGCCACACGGCGGTGACGCGCTTAACATCGtg                                                                                                             <  1:703803/151‑1 (MQ=255)
                                tGGGCGGGTTTATAGACGATATGTAAATCGCACCCTTGCGCCACCGTTCCGTGGCAGACAATTTCAAAACTAAAACGGACGGCGCTCTCCTCAACGCAGGAGAGCGCGCCAAGTTCCAGCCACACGGCGGTGACGCGCTTAACATCg                                                                                                               >  1:353721/1‑147 (MQ=255)
                                tGGGCGGGTTTATAGACGATATGTAAATCGCACCCTTGCGCCACCGTTCCGTGGCAGACAATTTCAAAACTAAAACGGACGGCGCTCTCCTCAACGCAGGAGAGCGCGCCAAGTTCCAGCCACACGGCGGTGACGCGCTTAACATCGtgc                                                                                                            >  2:421541/1‑150 (MQ=255)
                                         ttATAGACGATATGTAAATCGCACCCTTGCGCCACCGTTCCGTGGCAGACAATTTCAAAACTAAAACGGACGGCGCTCTCCTCAACGCAGGAGAGCGCGCCAAGTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCt                                                                                                  >  2:665169/1‑151 (MQ=255)
                                           ataGACGATATGTAAATCGCACCCTTGCGCCCCCGTTCCGTGGCAGACAATTTCAAAACTAAAACGGACGGCGCTCGCCTCAACGCAGGAGAGCGCGCCAAGTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTg                                                                                                 <  1:269914/150‑1 (MQ=255)
                                           ataGACGATATGTAAATCGCACCCTTGCGCCACCGTTCCGTGGCAGACAATTTCAAAACTAAAACGGACGGCGCTCTCCTCAACGCAGGAGAGCGCGCCAAGTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGc                                                                                                >  2:119808/1‑151 (MQ=255)
                                            taGACGATATGTAAATCGCACCCTTGCGCCACCGTTCCGTGGCAGACAATTTCAAAACTAAAACGGACGGCGCTCTCCTCAACGCAGGAGAGCGCGCCAAGTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGcc                                                                                               >  2:367461/1‑151 (MQ=255)
                                              gACGATATGTAAATCGCACCCTTGCGCCACCGTTCCGTGGCAGACAATTTCAAAACTAAAACGGACGGCGCTCTCCTCAACGCAGGAGAGCGCGCCAAGTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGt                                                                                             >  1:232929/1‑151 (MQ=255)
                                                    atGTAAATCGCACCCTTGCGCCACCGTTCCGTGGCAGACAATTTCAAAACTAAAACGGACGGCGCTCTCCTCAACGCAGGAGAGCGCGCCAAGTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGAt                                                                                        >  1:750872/1‑150 (MQ=255)
                                                    atGTAAATCGCACCCTTGCGCCACCGTTCCGTGGCAGACAATTTCAAAACTAAAACGGACGGCGCTCTCCTCAACGCAGGAGAGCGCGCCAAGTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGAt                                                                                        >  1:436361/1‑150 (MQ=255)
                                                    atGTAAATCGCACCCTTGCGCCACCGTTCCGTGGCAGACAATTTCAAAACTAAAACGGACGGCGCTCTCCTCAACGCAGGAGAGCGCGCCAAGTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATa                                                                                       <  1:275768/151‑1 (MQ=255)
                                                      gTAAATCGCACCCTTGCGCCACCGTTCCGTGGCAGACAATTTCAAAACTAAAACGGACGGCGCTCTCCTCAACGCAGGAGAGCGCGCCAAGTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAAt                                                                                     >  1:718946/1‑151 (MQ=255)
                                                        aaaTCGCACCCTTGCGCCACCGTTCCGTGGCAGACAATTTCAAAACTAAAACGGACGGCGCTCTCCTCAACGCAGGAGAGCGCGCCAAGTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAAt                                                                                     <  2:62222/149‑1 (MQ=255)
                                                           tCGCACCCTTGCGCCACCGTTCCGTGGCAGACAATTTCAAAACTAAAACGGACGGCGCTCTCCTCAACGCAGGAGAGCGCGCCAAGTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATACTTTCa                                                                                 >  2:68711/1‑150 (MQ=255)
                                                            cGCACCCTTGCGCCACCGTTCCGTGGCAGACAATTTCAAAACTAAAACGGACGGCGCTCTCCTCAACGCAGGAGAGCGCGCCAAGTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCaa                                                                                <  2:291185/150‑1 (MQ=255)
                                                                 ccTTGCGCCACCGTTCCGTGGCAGACAATTTCAAAACTAAAACGGACGGCGCTCTCCTCAACGCAGGAGAGCGCGCCAAGTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGgcgc                                                                          <  2:803342/151‑1 (MQ=255)
                                                                   ttGCGCCACCGTTCCGTGGCAGACAATTTCAAAACTAAAACGGACGGCGCTCTCCTCAACGCAGGAGAGCGCGCCAAGTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGct                                                                         >  2:122236/1‑150 (MQ=255)
                                                                        ccACCGTTCCGTGGCAGACAATTTCAAAACTAAAACGGACGGCGCTCTCCTCAACGCAGGAGAGCGCGCCAAGTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGc                                                                    >  2:599397/1‑150 (MQ=255)
                                                                           ccGTTCCGTGGCAGACAATTTCAAAACTAAAACGGACGGCGCTCTCCTCAACGCAGGAGAGCGCGCCAAGTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCaaa                                                                 <  1:407142/150‑1 (MQ=255)
                                                                                ccGTGGCAGACAATTTCAAAATTAAAACGGACGGCGCTCTCCTCAACGCAGGAGAGCGCGCCAAGTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGAc                                                            >  2:850327/1‑150 (MQ=255)
                                                                                ccGTGGCAGACAATTTCAAAACTAAAACGGACGGCGCTCTCCTCAACGCAGGAGAGCGCGCCAAGTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGAc                                                            >  1:134554/1‑150 (MQ=255)
                                                                                 cGTGGCAGACAATTTCAAAACTAAAACGGACGGCGCTCTCCTCAACGCAGGAGAGCGCGCCAAGTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACaa                                                          <  1:465172/151‑1 (MQ=255)
                                                                                 cGTGGCAGACAATTTAAAAACTAAAACGGGCGGCGCTCGCCTCAACGCAGGAGAGCGCGCCAAGTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACaa                                                          <  1:57218/151‑1 (MQ=255)
                                                                                  gTGGCAGACAATTTCAAAACTAAAACGGACGGCGCTCTCCTCAACGCAGGAGAGCGCGCCAAGTTCCAGCCACACGGCGGTGACGCGCTTAACAACGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACaa                                                          >  1:651645/1‑150 (MQ=255)
                                                                                    ggCAGACAATTTCAAAACTAAAACGGACGGCGCTCTCCTCAACGCAGGAGAGCGCGCCAAGTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTc                                                       >  1:629931/1‑151 (MQ=255)
                                                                                     gCAGACAATTTCAAAACTAAAACGGACGGCGCTCTCCTCAACGCAGGAGAGCGCGCCAAGTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTc                                                       >  2:173266/1‑150 (MQ=255)
                                                                                        gACAATTTCAAAACTAAAACGGACGGCGCTCTCCTCAACGCAGGAGAGCGCGCCAAGTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCCCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGCGACAACTCATg                                                    >  1:485160/1‑150 (MQ=255)
                                                                                           aaTTTCAAAACTAAAACGGACGGCGCTCTCCTCAACGCAGGAGAGCGCGCCAAGTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCAtt                                                >  2:849831/1‑151 (MQ=255)
                                                                                              ttCAAAACTAAAACGGACGGCGCTCTCCTCAACGCAGGAGAGCGCGCCAAGTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTc                                             >  1:522501/1‑151 (MQ=255)
                                                                                                cAAAACTAAAACGGACGGCGCTCTCCTCAACGCAGGAGAGCGCGCCAAGTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCg                                            <  2:529592/150‑1 (MQ=255)
                                                                                                 aaaaCTAAAGCGGACGGCGCTCTCCTCAACGCAGGAGAGCGCGCCAAGTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCg                                            <  1:538370/149‑1 (MQ=255)
                                                                                                     cTAAAACGGACGGCGCTCTCCTCAACGCAGGAGAGCGCGCCAAGTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTc                                      >  2:815209/1‑151 (MQ=255)
                                                                                                      tAAAACGGACGGCGCTCTCCTCAACGCAGGAGAGCGCGCCAAGTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTcc                                     >  2:143712/1‑151 (MQ=255)
                                                                                                       aaaaCGGACGGCGCTCTCCTCAACGCAGGAGAGCGCGCCAAGTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTc                                      <  1:272641/149‑1 (MQ=255)
                                                                                                       aaaaCGGACGGCGCTCTCCTCAACGCAGGAGAGCGCGCCAAGTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCt                                    <  1:408478/151‑1 (MQ=255)
                                                                                                         agcggacggCGCTCTCCTCAACGCAGGAGAGCGCGCCAAGTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGc                                  <  2:276627/149‑1 (MQ=255)
                                                                                                         aaCGGACGGCGCTCTCCTCAACGCAGGAGAGCGCGCCAAGTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTcc                                     >  1:814499/1‑148 (MQ=255)
                                                                                                             gacggCGCTCTCCTCAACGCAGGAGAGCGCGCCAAGTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGAcc                              <  2:789359/151‑1 (MQ=255)
                                                                                                               cgccgcTCTCCTCAACGCAGGAGAGCGCGCCAAGTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCaa                            <  2:369129/148‑1 (MQ=255)
                                                                                                                     tctcCTCAACGCAGGAGAGCGCGCCAAGTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAAcagc                       >  2:184149/1‑150 (MQ=255)
                                                                                                                        cctcAACGCAGGAGAGCGCGCCAAGTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAgg                    <  2:163105/150‑1 (MQ=255)
                                                                                                                            aaCGCAGGAGAGCGCGCCAAGTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGgcgc                <  2:444986/150‑1 (MQ=255)
                                                                                                                              cGCAGGAGAGCGCGCCAAGTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGAc             <  2:43422/151‑1 (MQ=255)
                                                                                                                               gCAGGAGAGCGCGCCAAGTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACg            >  1:36034/1‑151 (MQ=255)
                                                                                                                                     gagCGCGCCAAGTTCCAGCCTCACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGc      >  1:762053/1‑151 (MQ=255)
                                                                                                                                     gagCGCGCCAAGTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTg       <  2:555447/150‑1 (MQ=255)
                                                                                                                                         gcgcGAAGTTCCAGCCACGCGGGGGTGACGCGCTTAACGTTGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCgt    <  1:585722/149‑1 (MQ=255)
                                                                                                                                         gcgcCAAGTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCgtgt  >  1:582191/1‑151 (MQ=255)
                                                                                                                                                |                                                                                                                                               
CCACCACCTGGCTGCAATCCCAGCACCAAGCCTGGGCGGGTTTATAGACGATATGTAAATCGCACCCTTGCGCCACCGTTCCGTGGCAGACAATTTCAAAACTAAAACGGACGGCGCTCTCCTCAACGCAGGAGAGCGCGCCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGT  >  NC_000913/3140078‑3140365

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