Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,141,722 A→G L153P (CTG→CCG)  hybB ← putative cytochrome Ni/Fe component of hydrogenase‑2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,141,7220AG97.5% 371.3 / 0.0 121L153P (CTG→CCG) hybBputative cytochrome Ni/Fe component of hydrogenase‑2
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base G (51/67);  total (53/68)
Fisher's exact test for biased strand distribution p-value = 5.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.16e-01

Only 100 of 122 total aligned reads displayed.

ATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGTTCACGTTGAAGTGA  >  NC_000913/3141632‑3141811
                                                                                          |                                                                                         
aTGGTCGGCAGCAGCGCACCGAGGGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGt                                                                                         <  2:2625436/93‑1 (MQ=255)
 tGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTg                                                                                        <  1:1179185/93‑1 (MQ=255)
 tGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTg                                                                                        <  2:3098880/93‑1 (MQ=255)
  ggTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGc                                                                                       >  1:3176755/1‑93 (MQ=255)
  ggTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACAGTGc                                                                                       >  1:4734546/1‑93 (MQ=255)
   gTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGcc                                                                                      >  1:1722725/1‑93 (MQ=255)
   gTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGCGcc                                                                                      <  1:3510940/93‑1 (MQ=255)
    gcGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCg                                                                                     <  1:3330331/92‑1 (MQ=255)
      ggCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGa                                                                                   <  1:2495458/93‑1 (MQ=255)
          gcagcGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCaa                                                                               <  2:1320313/93‑1 (MQ=255)
          gcagcGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCaa                                                                               <  2:1470389/93‑1 (MQ=255)
           cagcGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCaaa                                                                              <  1:1258365/93‑1 (MQ=255)
           cagcGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGAGCaaa                                                                              <  2:1687597/93‑1 (MQ=255)
           cagcGCACCGAGCGCGATGATGAAGAACATCACCTTGCTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCaaa                                                                              <  1:3086644/93‑1 (MQ=255)
            agcGCACCGAGCGCGATGATGAAGAACATCCCCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAAc                                                                             >  1:1623805/1‑93 (MQ=255)
             gcgcACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAAcc                                                                            >  1:233158/1‑92 (MQ=255)
             gcgcACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACt                                                                            >  2:4144400/1‑93 (MQ=255)
             gcgcACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACt                                                                            <  2:168972/93‑1 (MQ=255)
             gcgcACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACt                                                                            >  1:376432/1‑93 (MQ=255)
                 cccGAGCGCCATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAg                                                                        <  1:3633892/92‑1 (MQ=255)
                 aCCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAg                                                                        <  1:4495089/93‑1 (MQ=255)
                 aCCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAg                                                                        >  2:3411485/1‑93 (MQ=255)
                 aCCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAg                                                                        >  2:2551197/1‑93 (MQ=255)
                  ccGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGt                                                                       <  2:1397276/93‑1 (MQ=255)
                  ccGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGt                                                                       <  2:1523668/93‑1 (MQ=255)
                   cGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTg                                                                      <  2:2385148/93‑1 (MQ=255)
                     aGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGcc                                                                    <  2:1116716/93‑1 (MQ=255)
                     aGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGcc                                                                    >  1:2694486/1‑93 (MQ=255)
                     aGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGcc                                                                    >  1:1880013/1‑93 (MQ=255)
                      gcgcGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCa                                                                   <  2:64797/93‑1 (MQ=255)
                         cgATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATcc                                                                >  1:1061783/1‑92 (MQ=255)
                           atgatgAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACg                                                              <  1:2834633/93‑1 (MQ=255)
                            tgatgaAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGc                                                             <  2:3132456/93‑1 (MQ=255)
                            tgatgaAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGc                                                             >  1:2244239/1‑93 (MQ=255)
                             gatgaAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGcc                                                            <  2:4369152/93‑1 (MQ=255)
                              atgaAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCg                                                           <  1:1821796/93‑1 (MQ=255)
                               tgaAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGa                                                          <  1:4381199/93‑1 (MQ=255)
                               tgaAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGa                                                          <  2:4418571/93‑1 (MQ=255)
                                gaagaaCATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGat                                                         >  1:3020165/1‑93 (MQ=255)
                                gaagaaCATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGat                                                         <  2:704552/93‑1 (MQ=255)
                                 aagaaCATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGata                                                        >  2:1827814/1‑93 (MQ=255)
                                 aagaaCATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGata                                                        >  1:2559985/1‑93 (MQ=255)
                                 aagaaCATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGata                                                        <  2:3923341/93‑1 (MQ=255)
                                  agaaCATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGatat                                                       >  1:1413031/1‑93 (MQ=255)
                                   gaaCATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGAtata                                                      <  1:4529866/93‑1 (MQ=255)
                                     aCATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGa                                                    <  1:1664360/93‑1 (MQ=255)
                                      cATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGAt                                                   >  2:2222957/1‑93 (MQ=255)
                                       aTCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATg                                                  <  2:944995/93‑1 (MQ=255)
                                       aTCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATg                                                  <  2:1569600/93‑1 (MQ=255)
                                       aTCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATg                                                  <  2:63732/93‑1 (MQ=255)
                                       aTCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATg                                                  >  1:1989094/1‑93 (MQ=255)
                                       aTCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATg                                                  <  1:4037159/93‑1 (MQ=255)
                                        tCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATgg                                                 >  2:2837753/1‑93 (MQ=255)
                                         cACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGt                                                >  1:2023270/1‑93 (MQ=255)
                                         cACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGt                                                <  2:2572515/93‑1 (MQ=255)
                                            cTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCAt                                             >  2:203027/1‑93 (MQ=255)
                                                 ttAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGa                                        <  1:826641/93‑1 (MQ=255)
                                                   aGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGAcc                                      >  2:1772549/1‑93 (MQ=255)
                                                    gTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGAccg                                     >  2:1782195/1‑93 (MQ=255)
                                                     tCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGAccgc                                    >  2:3671462/1‑93 (MQ=255)
                                                       gCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGAccgccg                                  <  1:457304/93‑1 (MQ=255)
                                                        cTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGt                                 >  1:2217630/1‑93 (MQ=255)
                                                         tGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGtc                                <  1:4195305/93‑1 (MQ=255)
                                                         tGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGtc                                <  1:463889/93‑1 (MQ=255)
                                                         tGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGtc                                >  1:466845/1‑93 (MQ=255)
                                                         tGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGtc                                <  2:1487606/93‑1 (MQ=255)
                                                         tGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGtc                                >  2:905732/1‑93 (MQ=255)
                                                          gTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGtct                               <  2:2209473/93‑1 (MQ=255)
                                                          gTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGtct                               <  2:142579/93‑1 (MQ=255)
                                                          gTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGtct                               <  1:669981/93‑1 (MQ=255)
                                                            aGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCg                             <  2:1944851/93‑1 (MQ=255)
                                                            aGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCg                             >  1:677058/1‑93 (MQ=255)
                                                            aGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCg                             <  1:3030219/93‑1 (MQ=255)
                                                             gCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGa                            <  2:748119/93‑1 (MQ=255)
                                                               gACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAAc                          <  2:3815952/93‑1 (MQ=255)
                                                                acacCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACa                         >  2:3021106/1‑93 (MQ=255)
                                                                 cacCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAg                        >  1:4719134/1‑93 (MQ=255)
                                                                  acCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGt                       <  1:728613/93‑1 (MQ=255)
                                                                    cTTCCACCCCAGACGTTCAAACGGTGCCGGAGCGAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTAc                     <  1:4593013/93‑1 (MQ=255)
                                                                      tCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCg                   <  2:484282/93‑1 (MQ=255)
                                                                       ccACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCTATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGa                  <  1:1751248/93‑1 (MQ=255)
                                                                        cACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAg                 <  1:363921/93‑1 (MQ=255)
                                                                         aCCCCAGACGTTCAAACGGTGCCGGAGCCCACTCCAGTGCCATCCCGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGt                >  1:2340875/1‑93 (MQ=255)
                                                                         aCCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGt                >  1:2278182/1‑93 (MQ=255)
                                                                          ccccAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGtt               <  2:3152265/93‑1 (MQ=255)
                                                                            ccAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGTTCa             <  1:1238352/93‑1 (MQ=255)
                                                                             cAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCNGNNTCGAACAGTACCNAGTTCan            >  1:1590123/1‑92 (MQ=255)
                                                                             cAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGTTCAc            <  2:837730/93‑1 (MQ=255)
                                                                              tgACGATTAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATTGATGGTCATACAGACCGCCGTCTCGAACAGCACCGAGTTCACg           >  2:4220001/2‑93 (MQ=255)
                                                                               gACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGTTCACGt          >  2:2946608/1‑93 (MQ=255)
                                                                               gACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGTTCACGt          >  1:980371/1‑93 (MQ=255)
                                                                               gACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCAGCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGTTCACGt          <  2:3043771/93‑1 (MQ=255)
                                                                                  gTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGTTCACGTTGa       >  1:515997/1‑93 (MQ=255)
                                                                                   ttCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGTTCACGTTGaa      >  1:3928036/1‑93 (MQ=255)
                                                                                   ttCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGTTCACGTTGaa      <  2:2571281/93‑1 (MQ=255)
                                                                                   ttCAAACAGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGTTCACGTTGaa      >  2:3112322/1‑93 (MQ=255)
                                                                                     cAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGTTCACGTTGAAGt    <  2:4632848/93‑1 (MQ=255)
                                                                                     cAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGTTCACGTTGAAGt    >  1:474827/1‑93 (MQ=255)
                                                                                      aaaCGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGTTCACGTTGAAGTg   <  2:3800774/93‑1 (MQ=255)
                                                                                       aaCGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGTTCACGTTGAAGTGa  <  1:4253713/93‑1 (MQ=255)
                                                                                          |                                                                                         
ATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGTTCACGTTGAAGTGA  >  NC_000913/3141632‑3141811

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