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,7220AG99.7% 912.7 / NA 284L153P (CTG→CCG) hybBputative cytochrome Ni/Fe component of hydrogenase‑2
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (123/160);  minor base C (1/0);  total (124/160)
Fisher's exact test for biased strand distribution p-value = 4.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 300 total aligned reads displayed.

ATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGTTCACGTTGAAGTGACCC  >  NC_000913/3141632‑3141814
                                                                                          |                                                                                            
 tGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTg                                                                                           >  2:156687/1‑93 (MQ=255)
 tGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTg                                                                                           >  2:2308304/1‑93 (MQ=255)
  ggTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGc                                                                                          <  2:7826055/93‑1 (MQ=255)
      ggCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGa                                                                                      >  1:637361/1‑93 (MQ=255)
       gcagcaGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAg                                                                                     <  2:2694826/93‑1 (MQ=255)
        cagcagCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGc                                                                                    <  1:8806183/93‑1 (MQ=255)
         agcagcGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCa                                                                                   >  1:7181361/1‑93 (MQ=255)
          gcagcGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCaa                                                                                  >  1:1435912/1‑93 (MQ=255)
          gcagcGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCaa                                                                                  >  1:1808213/1‑93 (MQ=255)
          gcagcGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCaa                                                                                  >  1:4828128/1‑93 (MQ=255)
           cagcGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCaaa                                                                                 >  2:317092/1‑93 (MQ=255)
               gcACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTcc                                                                             >  1:312193/1‑93 (MQ=255)
                  ccGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGt                                                                          >  1:8377793/1‑93 (MQ=255)
                  ccGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGt                                                                          >  2:8274579/1‑93 (MQ=255)
                  ccGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGt                                                                          <  2:6246269/93‑1 (MQ=255)
                  ccGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGt                                                                          <  1:9337781/93‑1 (MQ=255)
                  ccGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACCGTGCCGGAGCACACTCCAGt                                                                          >  2:8826405/1‑93 (MQ=255)
                   cGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTg                                                                         <  1:5768607/93‑1 (MQ=255)
                      acgcgATTATGAAAAACATCCCCTTTTTTAGTCCCCGTAGCGACACCTTCGACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCa                                                                      <  1:7228416/92‑1 (MQ=255)
                       cgcgATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCAt                                                                     <  1:9964224/93‑1 (MQ=255)
                        gcgATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATc                                                                    >  1:2526880/1‑93 (MQ=255)
                        gcgATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATc                                                                    <  2:1892826/93‑1 (MQ=255)
                         cgATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCa                                                                   <  2:7266985/93‑1 (MQ=255)
                         caatgatgAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCa                                                                   <  1:2618693/91‑1 (MQ=255)
                          gatgatGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCAc                                                                  <  2:3373680/93‑1 (MQ=255)
                          gatgatGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCAc                                                                  <  2:5215545/93‑1 (MQ=255)
                          catgatgAAGAACATCTCCTTGTTTAGTCGCTGTAGCGACCCCTTCCACCCCAGACGTTCAAACGGTTCCGGAGCAAACTCCAGTGCCATCAc                                                                  <  2:5134533/92‑1 (MQ=255)
                           atgatgAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACg                                                                 >  1:55219/1‑93 (MQ=255)
                           atgatgAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACg                                                                 >  2:1692334/1‑93 (MQ=255)
                              atgaAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCg                                                              >  2:9684284/1‑93 (MQ=255)
                               tgaAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCCCCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGa                                                             <  1:2795918/93‑1 (MQ=255)
                               tgaAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGa                                                             >  2:305346/1‑93 (MQ=255)
                               tgaAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGa                                                             <  2:8078673/93‑1 (MQ=255)
                               tgaAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGa                                                             <  2:8422745/93‑1 (MQ=255)
                                gaagaaCATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGat                                                            <  1:1212045/93‑1 (MQ=255)
                                 aagaaCATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGata                                                           <  1:694805/93‑1 (MQ=255)
                                     aCATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGa                                                       <  1:1964283/93‑1 (MQ=255)
                                       aTCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATg                                                     <  1:6778437/93‑1 (MQ=255)
                                        tCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATgg                                                    >  1:7159626/1‑93 (MQ=255)
                                        tCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATgg                                                    >  2:9548328/1‑93 (MQ=255)
                                        tCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATgg                                                    <  2:8740524/93‑1 (MQ=255)
                                         cACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGt                                                   <  1:3685250/93‑1 (MQ=255)
                                          aCCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTc                                                  >  1:6042343/1‑93 (MQ=255)
                                           ccTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCa                                                 <  2:6689018/93‑1 (MQ=255)
                                            cTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCAt                                                >  1:967204/1‑93 (MQ=255)
                                             ttGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATa                                               <  1:5397384/93‑1 (MQ=255)
                                             ttGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATa                                               >  1:8902554/1‑93 (MQ=255)
                                              tGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATAc                                              <  2:6573702/93‑1 (MQ=255)
                                              tGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATAc                                              >  1:3166644/1‑93 (MQ=255)
                                               gTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACa                                             >  2:7585059/1‑93 (MQ=255)
                                               gTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACa                                             >  2:2781135/1‑93 (MQ=255)
                                                  tAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGAc                                          <  1:10038249/93‑1 (MQ=255)
                                                  tAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGAc                                          <  2:387982/93‑1 (MQ=255)
                                                  tAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGAc                                          <  1:4216353/93‑1 (MQ=255)
                                                    gTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGAccg                                        >  1:740911/1‑93 (MQ=255)
                                                     tCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGAccgc                                       >  1:6542451/1‑93 (MQ=255)
                                                      cGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGAccgcc                                      >  2:5797779/1‑93 (MQ=255)
                                                      cGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGAccgcc                                      <  1:2992887/93‑1 (MQ=255)
                                                       gCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGAccgccg                                     >  2:3247786/1‑93 (MQ=255)
                                                        cTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGt                                    <  2:8890190/93‑1 (MQ=255)
                                                         tGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGtc                                   <  1:4502330/93‑1 (MQ=255)
                                                         tGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGtc                                   <  2:7299203/93‑1 (MQ=255)
                                                         tGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGtc                                   <  1:2196619/93‑1 (MQ=255)
                                                           tAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGtctc                                 >  1:3490172/1‑93 (MQ=255)
                                                           tAGCGACAAATTCCACCCCAGACGTTCAAACGGTGCCGGCGCAAACTCCAGTGCCATCACGCCGATATAGAAGGTCATACCGACCGCCGtctc                                 <  2:8971460/93‑1 (MQ=255)
                                                              cGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCCATATAGATGGTCATACAGACCGCCCTCTCGac                              >  1:793048/1‑92 (MQ=255)
                                                               gACACCTCCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAAc                             <  1:5642266/93‑1 (MQ=255)
                                                                cccccTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACa                            <  2:5855674/90‑1 (MQ=255)
                                                                acacCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAATCTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACa                            <  2:4870091/93‑1 (MQ=255)
                                                                acacCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACa                            <  2:8526833/93‑1 (MQ=255)
                                                                acacCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACa                            >  1:7660436/1‑93 (MQ=255)
                                                                   ccTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTa                         <  1:7892899/93‑1 (MQ=255)
                                                                   ccTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTa                         >  2:595073/1‑93 (MQ=255)
                                                                   ccTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTa                         >  2:5343284/1‑93 (MQ=255)
                                                                    cTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTAc                        >  2:4732825/1‑93 (MQ=255)
                                                                     ttCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTGcc                       <  2:5609056/93‑1 (MQ=255)
                                                                     ttCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTAcc                       >  2:2494382/1‑93 (MQ=255)
                                                                        cACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAg                    <  1:2869376/93‑1 (MQ=255)
                                                                         cccccAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGt                   <  2:7439435/92‑1 (MQ=255)
                                                                         aCCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGt                   <  2:4401209/93‑1 (MQ=255)
                                                                          ccccAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGtt                  >  1:2493248/1‑93 (MQ=255)
                                                                           cccAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGAACCGacatc                 >  2:6035095/1‑89 (MQ=255)
                                                                            ccAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGTTCa                >  1:9790271/1‑93 (MQ=255)
                                                                            ccAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGTTCa                <  1:3935882/93‑1 (MQ=255)
                                                                               gACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGTTCACGt             <  1:4079167/93‑1 (MQ=255)
                                                                               gACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGTTCACGt             >  2:2014245/1‑93 (MQ=255)
                                                                                gcGTTCAACCGGTGCCGTAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGTTCACGtt            <  2:3589328/92‑1 (MQ=255)
                                                                                aCGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGTTCACGtt            <  2:4631270/93‑1 (MQ=255)
                                                                                aCGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGTTCACGtt            <  1:2011932/93‑1 (MQ=255)
                                                                                 cGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGTTCACGTTg           >  1:5035732/1‑93 (MQ=255)
                                                                                 cGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGTTCACGTTg           <  2:7139041/93‑1 (MQ=255)
                                                                                   ttCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGTTCACGTTGaa         >  1:5298824/1‑93 (MQ=255)
                                                                                   ttCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGTTCACGTTGaa         <  2:9241767/93‑1 (MQ=255)
                                                                                   ttCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGTTCACGTTGaa         >  1:566928/1‑93 (MQ=255)
                                                                                    tCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGTTCACGTTGAAg        <  2:1107470/93‑1 (MQ=255)
                                                                                     taaaCGGTGCCGGAGCAATCTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGTTCACGTTGAAGt       <  1:6226030/92‑1 (MQ=255)
                                                                                     cAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATAGAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGTTCACGTTGAAGt       <  1:8169416/93‑1 (MQ=255)
                                                                                      aaaCGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGTTCACGTTGAAGTg      >  2:6402192/1‑93 (MQ=255)
                                                                                       aaCGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGTTCACGTTGAAGTGa     >  2:4497124/1‑93 (MQ=255)
                                                                                          aGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGTTCACGTTGAAGTGAccc  <  2:4312265/93‑1 (MQ=255)
                                                                                          |                                                                                            
ATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGTTCACGTTGAAGTGACCC  >  NC_000913/3141632‑3141814

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