Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 875,327 G→A *783* (TAG→TAA yliE → putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913875,3270GA100.0% 295.4 / NA 91*783* (TAG→TAAyliEputative membrane‑anchored cyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (44/47);  total (44/47)

TTGCGGAAGGCGTAGAGAGCGAAGACGTTGCGCGCGAATTAATCAAATTAGGATGTGTTCAGGCTCAGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGT  >  NC_000913/875191‑875459
                                                                                                                                        |                                                                                                                                    
ttGCGGAAGGCGTAGAGAGCGAAGACGTTGCGCGCGAATTAATCAAATTAGGATGTGTTCAGGCTCAGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAAtt                                                                                                                                    <  1:238821/139‑1 (MQ=255)
          cGTAGAGAGCGAAGACGTTGCGCGCGAATTAATCAAATTAGGATGTGTTCAGGCTCAGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTc                                                                                                                          >  1:494675/1‑139 (MQ=255)
            tAGAGAGCGAAGACGTTGCGCGCGAATTAATCAAATTAGGATGTGTTCAGGCTCAGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCa                                                                                                                        <  1:187229/139‑1 (MQ=255)
            tAGAGAGCGAAGACGTTGCGCGCGAATTAATCAAATTAGGATGTGTTCAGGCTCAGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCa                                                                                                                        <  2:494675/139‑1 (MQ=255)
                gagCGAAGACGTTGCGCGCGAATTAATCAAATTAGGATGTGTTCAGGCTCAGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAAt                                                                                                                    >  1:598728/1‑139 (MQ=255)
                   cGAAGACGTTGCGCGCGAATTAATCAAATTAGGATGTGTTCAGGCTCAGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCaa                                                                                                                 >  1:1712108/1‑139 (MQ=255)
                          gTTGCGCGCGAATTAATCAAATTAGGATGTGTTCAGGCTCAGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTc                                                                                                          <  1:1273682/139‑1 (MQ=255)
                            tGCGCGCGAATTAATCAAATTAGGATGTGTTCAGGCTCAGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGt                                                                                                        >  1:1119474/1‑139 (MQ=255)
                                cgcgAATTAATCAAATTAGGATGTGTTCAGGCTCAGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACtt                                                                                                    >  1:1656814/1‑139 (MQ=255)
                                     aTTAATCAAATTAGGATGTGTTCAGGCTCAGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTAcagt                                                                                               <  1:1828659/139‑1 (MQ=255)
                                     aTTAATCAAATTAGGATGTGTTCAGGCTCAGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTAcagt                                                                                               <  2:1656814/139‑1 (MQ=255)
                                        aaTCAAATTAGGATGTGTTCAGGCTCAGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTAcagtcag                                                                                            >  2:277797/1‑139 (MQ=255)
                                                tAGGATGTGTTCAGGCTCAGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGa                                                                                    <  2:1712108/139‑1 (MQ=255)
                                                tAGGATGTGTTCAGGCTCAGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATCTCCAGAATCAATAAGCTCGTACTTACAGTCAGTCTGCTGa                                                                                    <  2:1119474/139‑1 (MQ=255)
                                                   gATGTGTTCAGGCTCAGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGtt                                                                                                           >  2:43211/1‑113 (MQ=255)
                                                   gATGTGTTCAGGCTCAGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGtt                                                                                                           <  1:43211/113‑1 (MQ=255)
                                                       tgtTCAGGCTCAGGGGTATCTGTACCAGAACCCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAAtt                                                                                                                                    >  1:535257/1‑84 (MQ=255)
                                                       tgtTCAGGCTCAGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAAtt                                                                                                                                    <  2:535257/84‑1 (MQ=255)
                                                        gtTCAGGCTCAGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTa                                                                                                                              >  2:1670943/1‑89 (MQ=255)
                                                        gtTCAGGCTCAGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTa                                                                                                                              <  1:1670943/89‑1 (MQ=255)
                                                             ggCTCAGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTc                                                                                                                          >  1:692466/1‑88 (MQ=255)
                                                             ggCTCAGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTc                                                                                                                          <  2:692466/88‑1 (MQ=255)
                                                                 cAGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCg                                                                                                         >  1:1361797/1‑101 (MQ=255)
                                                                 cAGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCg                                                                                                         <  2:1361797/101‑1 (MQ=255)
                                                                 cAGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCAt                                                                          <  1:1619910/132‑1 (MQ=255)
                                                                 cAGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCAt                                                                          >  2:1619910/1‑132 (MQ=255)
                                                                   ggggTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCAt                                                                          <  1:1710367/130‑1 (MQ=255)
                                                                   ggggTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCAt                                                                          >  2:1710367/1‑130 (MQ=255)
                                                                   ggggTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTcagc                                                                 >  1:674338/1‑139 (MQ=255)
                                                                        aTCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTg                                                            >  1:1003832/1‑139 (MQ=255)
                                                                         tCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGt                                                                                                                               >  1:707585/1‑71 (MQ=255)
                                                                         tCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGt                                                                                                                               <  2:707585/71‑1 (MQ=255)
                                                                           tGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCt                                                         <  1:572961/139‑1 (MQ=255)
                                                                              aCCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCa                                                                                                                  >  2:1434257/1‑79 (MQ=255)
                                                                              aCCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCa                                                                                                                  <  1:1434257/79‑1 (MQ=255)
                                                                              aCCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTctgct                                                                                      <  2:342455/107‑1 (MQ=255)
                                                                              aCCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTctgct                                                                                      >  1:342455/1‑107 (MQ=255)
                                                                               ccAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCa                                                                                                                        >  1:1750239/1‑72 (MQ=255)
                                                                               ccAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCa                                                                                                                        <  2:1750239/72‑1 (MQ=255)
                                                                                  gAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAAt                                                                                                                >  1:1372561/1‑77 (MQ=255)
                                                                                  gAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAAt                                                                                                                <  2:1372561/77‑1 (MQ=255)
                                                                                  gAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACt                                                                                                     >  2:1368733/1‑88 (MQ=255)
                                                                                  gAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACt                                                                                                     <  1:1368733/88‑1 (MQ=255)
                                                                                  gAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGAt                                                  >  1:1391957/1‑139 (MQ=255)
                                                                                   aaaCCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGt                                                              >  1:1327943/1‑126 (MQ=255)
                                                                                   aaaCCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGt                                                              <  2:1327943/126‑1 (MQ=255)
                                                                                     aCCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTAcagt                                                                                               <  1:453454/91‑1 (MQ=255)
                                                                                     aCCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTAcagt                                                                                               >  2:453454/1‑91 (MQ=255)
                                                                                     aCCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTctgct                                                                                      <  1:1484596/100‑1 (MQ=255)
                                                                                     aCCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTctgct                                                                                      >  2:1484596/1‑100 (MQ=255)
                                                                                     aCCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTa                                               >  2:1811/1‑139 (MQ=255)
                                                                                       ccATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGc                                                       <  2:855634/129‑1 (MQ=255)
                                                                                       ccATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGc                                                       >  1:855634/1‑129 (MQ=255)
                                                                                       ccATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTaca                                             >  2:479303/1‑139 (MQ=255)
                                                                                       ccATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTaca                                             >  1:536569/1‑139 (MQ=255)
                                                                                         aTGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTctg                                                                                        <  1:546458/94‑1 (MQ=255)
                                                                                         aTGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACt                                           <  1:479303/139‑1 (MQ=255)
                                                                                         aTGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTACAGAATCAATAAGTTCGTACTTACAGTCAGTctg                                                                                        >  2:546458/1‑94 (MQ=255)
                                                                                           gCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCa                                         <  2:674338/139‑1 (MQ=255)
                                                                                              aTTCTCCGACTGTGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTTTCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAAt                                      <  1:74236/139‑1 (MQ=255)
                                                                                               ttCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATg                                                                                                                            <  2:918689/52‑1 (MQ=255)
                                                                                               ttCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATg                                                                                                                            >  1:918689/1‑52 (MQ=255)
                                                                                                 ctcCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGc                                                       <  2:343204/119‑1 (MQ=255)
                                                                                                 ctcCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGc                                                       >  1:343204/1‑119 (MQ=255)
                                                                                                 ctcCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGt                                   >  1:147390/1‑139 (MQ=255)
                                                                                                 ctcCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGt                                   >  1:301212/1‑139 (MQ=255)
                                                                                                    cGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCa                                                                                                                        >  2:800345/1‑51 (MQ=255)
                                                                                                    cGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCa                                                                                                                        <  1:800345/51‑1 (MQ=255)
                                                                                                    cGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCg                                                      <  1:1166639/117‑1 (MQ=255)
                                                                                                    cGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCg                                                      >  2:1166639/1‑117 (MQ=255)
                                                                                                    cGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTaaa                                <  2:1391957/139‑1 (MQ=255)
                                                                                                      ccTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAgg                              <  1:277797/139‑1 (MQ=255)
                                                                                                        tGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGaa                            <  2:536569/139‑1 (MQ=255)
                                                                                                         gggATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGc                                                       <  1:649761/111‑1 (MQ=255)
                                                                                                         gggATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGc                                                       >  2:649761/1‑111 (MQ=255)
                                                                                                         gggATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAAc                           >  1:1620819/1‑139 (MQ=255)
                                                                                                                aaGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGta                    <  1:1811/139‑1 (MQ=255)
                                                                                                                aaGTGGAAAATTAGTAAAAGAGTAATTTACGTATGGCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGta                    <  1:1584706/139‑1 (MQ=255)
                                                                                                                     gAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTcagca                                                                <  2:1238363/90‑1 (MQ=255)
                                                                                                                     gAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTcagca                                                                >  1:1238363/1‑90 (MQ=255)
                                                                                                                     gAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGcc               >  2:1025497/1‑139 (MQ=255)
                                                                                                                     gAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGcc               >  2:881344/1‑139 (MQ=255)
                                                                                                                     gAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGcc               >  2:263784/1‑139 (MQ=255)
                                                                                                                      aaaaTTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGc                                                       >  1:1449337/1‑98 (MQ=255)
                                                                                                                      aaaaTTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGc                                                       <  2:1449337/98‑1 (MQ=255)
                                                                                                                          ttAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGaaa          <  2:598728/139‑1 (MQ=255)
                                                                                                                              tAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGt      <  2:301212/139‑1 (MQ=255)
                                                                                                                              tAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGt      <  1:881344/139‑1 (MQ=255)
                                                                                                                                aaaGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCg                                                      >  1:1357526/1‑89 (MQ=255)
                                                                                                                                aaaGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCg                                                      <  2:1357526/89‑1 (MQ=255)
                                                                                                                                  agagTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGt  <  1:1025497/139‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                    
TTGCGGAAGGCGTAGAGAGCGAAGACGTTGCGCGCGAATTAATCAAATTAGGATGTGTTCAGGCTCAGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGT  >  NC_000913/875191‑875459

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