Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,843,138 G→A P453S (CCG→TCG)  hycE ← large subunit of hydrogenase 3

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,843,1380GA99.7% 1231.5 / NA 341P453S (CCG→TCG) hycElarge subunit of hydrogenase 3
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (129/211);  total (129/212)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 346 total aligned reads displayed.

GGATATCATCGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGGCAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACGCACTTTCAGACGGGAAATAACGTCGCAGCCCT  >  NC_000913/2843048‑2843227
                                                                                          |                                                                                         
ggATATCATCGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACa                                                                                         >  1:5574240/1‑93 (MQ=255)
 gATATCATCGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAg                                                                                        >  1:2340181/1‑93 (MQ=255)
   tatCATCGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGt                                                                                      >  1:808537/1‑93 (MQ=255)
   tatCATCGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGt                                                                                      >  2:5809302/1‑93 (MQ=255)
   tatCATCGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGt                                                                                      >  2:5999750/1‑93 (MQ=255)
    atcatcGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGtt                                                                                     <  1:4282656/93‑1 (MQ=255)
     tcatcGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTa                                                                                    <  1:604702/93‑1 (MQ=255)
     tcatcGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTa                                                                                    <  1:6842338/93‑1 (MQ=255)
     tcatcGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTa                                                                                    <  2:4068268/93‑1 (MQ=255)
     tcatcGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTa                                                                                    <  1:119190/93‑1 (MQ=255)
     tcatcGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTa                                                                                    <  2:7659128/93‑1 (MQ=255)
      catcGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCAGCAGTGGGCCACCCGACAGGTTAt                                                                                   >  2:4285988/1‑93 (MQ=255)
        tcGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATcc                                                                                 <  1:530875/93‑1 (MQ=255)
        tcGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATcc                                                                                 <  1:2606420/93‑1 (MQ=255)
        tcGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATcc                                                                                 >  2:5570677/1‑93 (MQ=255)
        tcGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATcc                                                                                 <  2:2303645/93‑1 (MQ=255)
         cgccgcGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCa                                                                                <  1:4733659/93‑1 (MQ=255)
         cgccgcGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCa                                                                                >  2:748576/1‑93 (MQ=255)
           ccgcGGCGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATTTAGGGAAAGCCTTCCCCCATCAGTGGGCCACCCGACAGGTTATCCAGa                                                                              <  1:1526117/93‑1 (MQ=255)
            cgcgcgGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCTACCATCAGTGGGCCACCCGACAGGTTATCCAGAc                                                                             >  2:6993598/1‑93 (MQ=255)
            cgcgcgGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGAc                                                                             <  1:6242339/93‑1 (MQ=255)
              cgcgGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCg                                                                           <  2:253001/93‑1 (MQ=255)
               gcggcgCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGt                                                                          <  1:746978/93‑1 (MQ=255)
               gcggcgCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGt                                                                          >  1:7203310/1‑93 (MQ=255)
               gcggcgCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGt                                                                          <  1:6600273/93‑1 (MQ=255)
                cggcgCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTa                                                                         <  1:2759451/93‑1 (MQ=255)
                  gcgcTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGt                                                                       <  1:4050359/93‑1 (MQ=255)
                       tCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATc                                                                  <  1:3166579/93‑1 (MQ=255)
                       tCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATc                                                                  >  1:7173108/1‑93 (MQ=255)
                         ggCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGCCCGTAGTCGATCAt                                                                <  2:4890786/93‑1 (MQ=255)
                         ggCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCAt                                                                <  1:5773724/93‑1 (MQ=255)
                          gCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATg                                                               >  2:627222/1‑93 (MQ=255)
                          gCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATg                                                               <  1:98525/93‑1 (MQ=255)
                          gCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATg                                                               >  2:6818602/1‑93 (MQ=255)
                            caaGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGtt                                                             <  1:2495342/92‑1 (MQ=255)
                            aaaGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGtt                                                             <  2:3379358/93‑1 (MQ=255)
                               gCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAg                                                          >  1:7083855/1‑93 (MQ=255)
                                cccAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCCCCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAgc                                                         <  2:2147129/93‑1 (MQ=255)
                                 ccAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAgcg                                                        >  2:2074710/1‑93 (MQ=255)
                                  cAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAgcgc                                                       >  2:2915709/1‑93 (MQ=255)
                                    gcgcAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCgg                                                     >  2:287910/1‑93 (MQ=255)
                                      gcAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGta                                                   <  1:3731197/93‑1 (MQ=255)
                                      gcAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGta                                                   >  2:5027608/1‑93 (MQ=255)
                                      gcAAAGCGGTGCGGAATGTAAGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGta                                                   <  1:6096255/93‑1 (MQ=255)
                                      gcAAAGCGGTACGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCGCCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGta                                                   <  1:2688038/93‑1 (MQ=255)
                                       cAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGtat                                                  <  2:4742973/93‑1 (MQ=255)
                                       cAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGtat                                                  <  2:4348795/93‑1 (MQ=255)
                                        aaaGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGtata                                                 <  2:3699090/93‑1 (MQ=255)
                                         aaGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAg                                                <  2:2665462/93‑1 (MQ=255)
                                          aGCGTTGCGGAATGTAGGTAAAGCCTTCCCCCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGa                                               <  2:1513474/93‑1 (MQ=255)
                                          aGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGa                                               <  1:75644/93‑1 (MQ=255)
                                          aGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGa                                               <  2:2337573/93‑1 (MQ=255)
                                           gCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGAc                                              <  2:7768348/93‑1 (MQ=255)
                                             ggTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATATACtt                                            <  1:7642722/93‑1 (MQ=255)
                                              gTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTc                                           <  1:6709561/93‑1 (MQ=255)
                                               tGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCg                                          >  1:2158123/1‑93 (MQ=255)
                                                 gggAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGtt                                        <  1:3921945/92‑1 (MQ=255)
                                                 cGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGtt                                        <  2:5757831/93‑1 (MQ=255)
                                                 cGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGtt                                        <  1:5036173/93‑1 (MQ=255)
                                                  ggAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTg                                       >  1:1116560/1‑93 (MQ=255)
                                                   gAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGa                                      >  1:6394789/1‑93 (MQ=255)
                                                     aTGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATa                                    <  1:5908426/93‑1 (MQ=255)
                                                     aTGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATa                                    <  1:5421505/93‑1 (MQ=255)
                                                       gTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACg                                  >  2:6242876/1‑93 (MQ=255)
                                                        tAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACGc                                 >  2:4207432/1‑93 (MQ=255)
                                                          ggTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACGCAc                               <  1:5265445/93‑1 (MQ=255)
                                                             aaaGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACGCACttt                            <  1:7843959/93‑1 (MQ=255)
                                                             aaaGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACGCACttt                            >  1:3419842/1‑93 (MQ=255)
                                                              aaGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACGCACTTTc                           <  1:1395972/93‑1 (MQ=255)
                                                              aaGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACGCACTTTc                           >  1:2230611/1‑93 (MQ=255)
                                                               aGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACGCACTTTCa                          <  2:792771/93‑1 (MQ=255)
                                                               aGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACGCACTTTCa                          <  1:1248213/93‑1 (MQ=255)
                                                               aGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACGCACTTTCa                          >  2:3909701/1‑93 (MQ=255)
                                                               aGCCTCCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACGCACTTTCa                          <  2:1439760/93‑1 (MQ=255)
                                                                 ccTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACGCACTTTCAGa                        <  1:1953219/93‑1 (MQ=255)
                                                                 ccTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACGCACTTTCAGa                        <  2:6588613/93‑1 (MQ=255)
                                                                  cTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACGCACTTTCAGAc                       <  1:1016535/93‑1 (MQ=255)
                                                                   ttCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACGCACTTTCAGACg                      >  2:1653827/1‑93 (MQ=255)
                                                                   ttCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACGCACTTTCAGACg                      >  1:4881367/1‑93 (MQ=255)
                                                                       accaTCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACGCACTTTCAGACGGGaa                  >  2:5927646/1‑93 (MQ=255)
                                                                       accaTCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACGCACTTTCAGACGGGaa                  >  1:3454811/1‑93 (MQ=255)
                                                                        ccaTCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACGCACTTTCAGACGGGaaa                 >  1:1499702/1‑93 (MQ=255)
                                                                         caTCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACGCACTTTCAGACGGGAAAt                <  1:655017/93‑1 (MQ=255)
                                                                         caTCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACGCACTTTCAGACGGGAAAt                <  1:1803101/93‑1 (MQ=255)
                                                                          aTCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACGCACTTTCAGACGGGAAATa               <  2:4182098/93‑1 (MQ=255)
                                                                          aTCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACGCACTTTCAGACGGGAAATa               <  2:4843171/93‑1 (MQ=255)
                                                                          aTCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACGCACTTTCAGACGGGAAATa               <  2:6128471/93‑1 (MQ=255)
                                                                          aTCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACGCACTTTCAGACGGGAAATa               >  2:3151559/1‑93 (MQ=255)
                                                                          aTCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACGCACTTTCAGACGGGAAATa               <  2:4443267/93‑1 (MQ=255)
                                                                          aTCAGTGGACCACCCGACAGGTTATCCAGTCCGTAGTCGATCATGTTCAGCGCGGTATAGCCTTCGTTGATACGCGCTTTCAGACGGGAAATa               <  2:7184213/93‑1 (MQ=255)
                                                                            cAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACGCACTTTCAGACGGGAAATAAc             <  2:5382730/93‑1 (MQ=255)
                                                                             aGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACGCACTTTCAGACGGGAAATAACg            <  2:1840990/93‑1 (MQ=255)
                                                                               tGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACGCACTTTCAGACGGGAAATAACGTc          >  2:6715673/1‑93 (MQ=255)
                                                                                    cACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACGCACTTTCAGACGGGAAATAACGTCGCAGc     >  1:2972488/1‑93 (MQ=255)
                                                                                     aCCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACGCACTTTCAGACGGGAAATAACGTCGCAGcc    >  2:3335535/1‑93 (MQ=255)
                                                                                      cccGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACGCACTTTCAGACGGGAAATAACGTCGCAGccc   <  2:6690705/93‑1 (MQ=255)
                                                                                      cccGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACGCACTTTCAGACGGGAAATAACGTCGCAGccc   <  1:5459217/93‑1 (MQ=255)
                                                                                      cccGACAGGTTAGCCAGTCCGTAGTCGATCATGTTCGGCTCGTTATAGACTTCGTTGATACGCACTTTCAGACGGGAAATAACGTCGCAGccc   <  2:3602746/93‑1 (MQ=255)
                                                                                       ccGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACGCACTTTCAGACGGGAAATAACGTCGCAGCCCt  <  1:3706295/93‑1 (MQ=255)
                                                                                       ccGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACGCACTTTCAGACGGGAAATAACGTCGCAGCCCt  <  1:4449927/93‑1 (MQ=255)
                                                                                          |                                                                                         
GGATATCATCGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGGCAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACGCACTTTCAGACGGGAAATAACGTCGCAGCCCT  >  NC_000913/2843048‑2843227

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