Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,639,286 A→G V180A (GTG→GCG)  gcpE ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,639,2860AG99.7% 949.7 / NA 300V180A (GTG→GCG) gcpEhypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (155/144);  total (156/144)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 312 total aligned reads displayed.

GCTTCGGTGATCCCCAGATGCAACGGCTGATCGATCTGTTTTGCCAGCAAACGATAAGACTCAACAGCGAGGAAGACGTCAGACGCTTTCACGCTGACTTTGAACTGATCGAAGTTCAGGCGATCGAGATGATCAACATGACGCATGGCAGATTCCAGCAACGCCTGCGGCGTCGGTTCGC  >  NC_000913/2639196‑2639376
                                                                                          |                                                                                          
gCTTCGGTGATCCCCAGATGCAACGGCTGATCGATCTGTTTTGCCAGCAAACGATAAGACTCAACAGCGAGGAAGACGTCAGACGCTTTcgcg                                                                                          >  2:1959532/1‑93 (MQ=255)
 cTTCGGTGATCCCCAGATGCAACGGCTGATCGATCTGTTTTGCCAGCAAACGATAAGACTCAACAGCGAGGAAGACGTCAGACGCTTTCgcgc                                                                                         >  1:7536918/1‑93 (MQ=255)
  ttCGGTGATCCCCAGATGCAACGGCTGATCGATCTGTTTTGCCAGCAAACGATAAGACTCAACAGCGAGGAAGACGTCAGACGCTTTCGCGCt                                                                                        <  1:4879964/93‑1 (MQ=255)
   tCGGTGATCCCCAGATGCAACGGCTGATCGATCTGTTTTGCCAGCAAACGATAAGACTCAACAGCGAGGAAGACGTCAGACGCTTTCGCGCTg                                                                                       >  2:2888654/1‑93 (MQ=255)
      gTGATCCCCAGATGCAACGGCTGATCGATCTGTTTTGCCAGCAAACGATAAGACTCAACAGCGAGGAAGACGTCAGACGCTTTCGCGCTGACt                                                                                    >  2:4827179/1‑93 (MQ=255)
       tGATCCCCAGATGCAACGGCTGATCGATCTGTTTTGCCAGCAAACGATAAGACTCAACAGCGAGGAAGACGTCAGACGCTTTCGCGCTGACtt                                                                                   >  2:6734061/1‑93 (MQ=255)
        gATCCCCAGATGCAACGGCTGATCGATCTGTTTTGCCAGCAAACGATAAGACTCAACAGCGAGGAAGACGTCAGACGCTTTCGCGCTGACttt                                                                                  <  1:6551428/93‑1 (MQ=255)
          tCCCCAGATGCAACGGCTGATCGATCTGTTTTGCCAGCAAACGATAAGACTCAACAGCGAGGAAGACGTCAGACGCTTTCGCGCTGACTTTGa                                                                                >  2:543227/1‑93 (MQ=255)
            cccAGATGCAACGGCTGATCGATCTGTTTTGCCAGCAAACGATAAGACTCAACAGCGAGGAAGACGTCAGACGCTTTCGCGCTGACTTTGAAc                                                                              >  1:5179388/1‑93 (MQ=255)
             ccAGATGCAACGGCTGATCGATCTGTTTTGCCAGCAAACGATAAGACTCAACAGCGAGGAAGACGTCAGACGCTTTCGCGCTGACTTTGAACt                                                                             >  1:1456927/1‑93 (MQ=255)
             ccAGATGCAACGGCTGATCGATCTGTTTTGCCAGCAAACGATAAGACTCAACAGCGAGGAAGACGTCAGACGCTTTCGCGCTGACTTTGAACt                                                                             <  2:313351/93‑1 (MQ=255)
              cAGATGCAACGGCTGATCGATCTGTTTTGCCAGCAAACGATAAGACTCAACAGCGAGGAAGACGTCAGACGCTTTCGCGCTGACTTTGAACTg                                                                            >  2:85208/1‑93 (MQ=255)
                gATGCAACGGCTGATCGATCTGTTTTGCCAGCAAACGATAAGACTCAACAGCGAGGAAGACGTCAGACGCTTTCGCGCTGACTTTGAACTGAt                                                                          <  1:4461769/93‑1 (MQ=255)
                gATGCAACGGCTGATCGATCTGTTTTGCCAGCAAACGATAAGACTCAACAGCGAGGAAGACGTCAGACGCTTTCGCGCTGACTTTGAACTGAt                                                                          >  1:7042402/1‑93 (MQ=255)
                gATGCAACGGCTGATCGATCTGTTTTGCCAGCAAACGATAAGACTCAACAGCGAGGAAGACGTCAGACGCTTTCGCGCTGACTTTGAACTGAt                                                                          <  2:7648489/93‑1 (MQ=255)
                   gCAACGGCTGATCGATCTGTTTTGCCAGCAAACGATAAGACTCAACAGCGAGGAAGACGTCAGACGCTTTCGCGCTGACTTTTGACTGATCGa                                                                       >  1:1045719/1‑93 (MQ=255)
                   gCAACGGCTGATCGATCTGTTTTGCCAGCAAACGATAAGACTCAACAGCGAGGAAGACGTCAGACGCTTTCGCGCTGACTTTGAACTGATCGa                                                                       <  2:914689/93‑1 (MQ=255)
                   gCAACGGCTGATCGATCTGTTTTGCCAGCAAACGATAAGACTCAACAGCGAGGAAGACGTCAGACGCTTTCGCGCTGACTTTGAACTGATCGa                                                                       >  1:1735695/1‑93 (MQ=255)
                    cAACGGCTGATCGATCTGTTTTGCCAGCAAACGATAAGACTCAACAGCGAGGAAGACGTCAGACGCTTTCGCGCTGACTTTGAACTGATCGaa                                                                      <  2:1230024/93‑1 (MQ=255)
                     aaCGGCTGATCGATCTGTTTTGCCAGCAAACGATAAGACTCAACAGCGAGGAAGACGTCAGACGCTTTCGCGCTGACTTTGAACTGATCGAAg                                                                     >  2:3354676/1‑93 (MQ=255)
                     aaCGGCTGATCGATCTGTTTTGCCAGCAAACGATAAGACTCAACAGCGAGGAAGACGTCAGACGCTTTCGCGCTGACTTTGAACTGATCGAAg                                                                     >  2:5093693/1‑93 (MQ=255)
                     aaCGGCTGATCGATCTGTTTTGCCAGCAAACGATAAGACTCAACAGCGAGGAAGACGTCAGACGCTTTCGCGCTGACTTTGAACTGATCGAAg                                                                     >  2:7498436/1‑93 (MQ=255)
                      agcgCTGATCGATCTGTTATGCCAGCAAACGATAAGACTCAACCGCGAGCAGGACGTCAGGCGCTTTCGCGCTGACTTTGAACTGATCGAAGt                                                                    <  2:5196140/90‑1 (MQ=255)
                      aCGGCTGATCGATCTGTTTTGCCAGCAAACGATAAGACTCAACAGCGAGGAAGACGTCAGACGCTTTCGCGCTGACTTTGAACTGATCGAAGt                                                                    >  1:5370010/1‑93 (MQ=255)
                       cGGCTGATCGATCTGTTTTGCCAGCAAACGATAAGACTCAACAGCGAGGAAGACGTCAGACGCTTTCGCGCTGACTTTGAACTGATCGAAGtt                                                                   <  1:5470757/93‑1 (MQ=255)
                        ggCTGATCGATCTGTTTTGCCAGCAAACGATAAGACTCAACAGCGAGGAAGACGTCAGACGCTTTCGCGCTGACTTTGAACTGATCGAAGTTc                                                                  >  2:3690054/1‑93 (MQ=255)
                        ggCTGATCGATCTGTTTTGCCAGCAAACGATAAGACTCAACAGCGAGGAAGACGTCAGACGCTTTCGCGCTGACTTTGAACTGATCGAAGTTc                                                                  <  2:3203454/93‑1 (MQ=255)
                           tGATCGATCTGTTTTGCCAGCAAACGATAAGACTCAACAGCGAGGAAGACGTCAGACGCTTTCGCGCTGACTTTGAACTGATCGAAGTTCAgg                                                               <  1:3995535/93‑1 (MQ=255)
                            gatcgatcTGTTTTGCCAGCAAACGATAAGACTCAACAGCGAGGAAGACGTCAGACGCTTTCGCGCTGACTTTGAACTGATCGAAGTTCAggg                                                              >  2:1884192/1‑92 (MQ=255)
                            gatcgatcTGTTTTGCCAGCAAACGATAAGACTCAACAGCGAGGAAGACGTCAGACGCTTTCGCGCTGACTTTGAACTGATCGAAGTTCAGGc                                                              <  2:5825746/93‑1 (MQ=255)
                            gatcgatcTGTTTTGCCAGCAAACGATAAGACTCAACAGCGAGGAAGACGTCAGACGCTTTCGCGCTGACTTTGAACTGATCGAAGTTCAGGc                                                              <  2:4736611/93‑1 (MQ=255)
                               cgatcTGTTTTGCCAGCAAACGATAAGACTCAACAGCGAGGAAGACGTCAGACGCTTTCGCGCTGACTTTGAACTGATCGAAGTTCAGGCGAt                                                           <  1:5070099/93‑1 (MQ=255)
                                gatcTGTTTTGCCAGCAAACGATAAGACTCAACAGCGAGGAAGACGTCAGACGCTTTCGCGCTGACTTTGAACTGATCGAAGTTCAGGCGCTc                                                          >  1:6718078/1‑93 (MQ=255)
                                  tcTGTTTTGCCAGCAAACGATAAGACTCAACAGCGAGGAAGACGTCAGACGCTTTCGCGCTGACTTTGAACTGATCGAAGTTCAGGCGATCga                                                        >  2:1478225/1‑93 (MQ=255)
                                  tcTGTTTTGCCAGCAAACGATAAGACTCAACAGCGAGGAAGACGTCAGACGCTTTCGCGCTGACTTTGAACTGATCGAAGTTCAGGCGATCga                                                        >  1:7931860/1‑93 (MQ=255)
                                  tcTGTTTTGCCAGCAAACGATAAGACTCAACAGCGAGGAAGACGTCAGACGCTTTCGCGCTGACTTTGAACTGATCGAAGTTCAGGCGATCga                                                        >  2:172002/1‑93 (MQ=255)
                                  tcTGTTTTGCCAGCAAACGATAAGACTCAACAGCGAGGAAGACGTCAGACGCTTTCGCGCTGACTTTGAACTGATCGAAGTTCAGGCGATCga                                                        >  1:2566626/1‑93 (MQ=255)
                                    tGTTTTGCCAGCAAACGATAAGACTCAACAGCGAGGAAGACGTCAGACGCTTTCGCGCTGACTTTGAACTGATCGAAGTTCAGGCGATCgaga                                                      >  2:1149799/1‑93 (MQ=255)
                                     gTTTTGCCAGCAAACGATAAGACTCAACAGCGAGGAAGACGTCAGACGCTTTCGCGCTGACTTTGAACTGATCGAAGTTCAGGCGATCGAgat                                                     >  1:617058/1‑93 (MQ=255)
                                       tttGCCAGCAAACGATAAGACTCAACAGCGAGGAAGACGTCAGACGCTTTCGCGCTGACTTTGAACTGATCGAAGTTCAGGCGATCGAgatga                                                   >  1:6420343/1‑93 (MQ=255)
                                       tttGCCAGCAAACGATAAGACTCAACAGCGAGGAAGACGTCAGACGCTTTCGCGCTGACTTTGAACTGATCGAAGTTCAGGCGATCGAgatga                                                   <  1:4560232/93‑1 (MQ=255)
                                       tttGCCAGCAAACGATAAGACTCAACAGCGAGGAAGACGTCAGACGCTTTCGCGCTGACTTTGAACTGATCGAAGTTCAGGCGATCGAgatga                                                   <  1:2250845/93‑1 (MQ=255)
                                         tGCCAGCAAACGATAAGACTCAACAGCGAGGAAGACGTCAGACGCTTTCGCGCTGACTTTGAACTGATCGAAGTTCAGGCGATCGAGATGATc                                                 <  2:7175640/93‑1 (MQ=255)
                                            cAGCAAACGATAAGACTCAACAGCGAGGAAGACGTCAGACGCTTTCGCGCTGACTTTGAACTGATCGAAGTTCAGGCGATCGAGATGATCAAc                                              >  2:7824050/1‑93 (MQ=255)
                                             aGCAAACGATAAGACTCAACAGCGAGGAAGACGTCAGACGCTTTCGCGCTGACTTTGAACTGATCGAAGTTCAGGCGATCGAGATGATCAACa                                             <  2:4361118/93‑1 (MQ=255)
                                              gCAAACGATAAGACTCAACAGCGAGGAAGACGTCAGACGCTTTCGCGCTGACTTTGAACTGATCGAAGTTCAGGCGATCGAGATGATCAACAt                                            >  2:4543698/1‑93 (MQ=255)
                                              gCAAACGATAAGACTCAACAGCGAGGAAGACGTCAGACGCTTTCGCGCTGACTTTGAACTGATCGAAGTTCAGGCGATCGAGATGATCAACAt                                            >  1:4059433/1‑93 (MQ=255)
                                              gCAAACGATAAGACTCAACAGCGAGGAAGACGTCAGACGCTTTCGCGCTGACTTTGAACTGATCGAAGTTCAGGCGATCGAGATGATCAACAt                                            <  1:2000685/93‑1 (MQ=255)
                                                aaaCGATAAGACTCAACAGCGAGGAAGACGTCAGACGCTTTCGCGCTGACTTTGAACTGATCGAAGTTCAGGCGATCGAGATGATCAACATGa                                          <  1:2752572/93‑1 (MQ=255)
                                                 aaCGATAAGACTCAACAGCGAGGAAGACGTCAGACGCTTTCGCGCTGACTTTGAACTGATCGAAGTTCAGGCGATCGAGATGATCAACATGAc                                         <  2:1337512/93‑1 (MQ=255)
                                                 aaCGATAAGACTCAACAGCGAGGAAGACGTCAGACGCTTTCGCGCTGACTTTGAACTGATCGAAGTTCAGGCGATCGAGATGATCAACATGAc                                         >  2:3575422/1‑93 (MQ=255)
                                                 aaCGATAAGACTCAACAGCGAGGAAGACGTCAGACGCTTTCGCGCTGACTTTGAACTGATCGAAGTTCAGGCGATCGAGATGATCAACATGAc                                         >  1:2203140/1‑93 (MQ=255)
                                                  aCGATAAGACTCAACAGCGAGGAAGACGTCAGACGCTTTCGCGCTGACTTTGAACTGATCGAAGTTCAGGCGATCGAGATGATCAACATGACg                                        >  1:3345419/1‑93 (MQ=255)
                                                  aCGATAAGACTCAACAGCGAGGAAGACGTCAGACGCTTTCGCGCTGACTTTGAACTGATCGAAGTTCAGGCGATCGAGATGATCAACATGACg                                        >  1:5348682/1‑93 (MQ=255)
                                                    gATAAGACTCAACAGCGAGGAAGACGTCAGACGCTTTCGCGCTGACTTTGAACTGATCGAAGTTCAGGCGATCGAGATGATCAACATGACGCa                                      <  2:6377215/93‑1 (MQ=255)
                                                    gATAAGACTCAACAGCGAGGAAGACGTCAGACGCTTTCGCGCTGACTTTGAACTGATCGAAGTTCAGGCGATCGAGATGATCAACATGACGCa                                      <  2:5524380/93‑1 (MQ=255)
                                                       aaGACTCAACAGCGAGGAAGACGTCAGACGCTTTCGCGCTGACTTTGAACTGATCGAAGTTCAGGCGATCGAGATGATCAACATGACGCATgg                                   <  2:3548097/93‑1 (MQ=255)
                                                        aGACTCAACAGCGAGGAAGACGGCAGACGCTTTCGCGCTGACCTTGAACTGATCGAAGGTCAGGCGATCGAGATGATCAACATGAAGCATTGc                                  >  2:7082200/1‑93 (MQ=255)
                                                         gACTCAACAGCGAGGAAGACGTCAGACGCTTTCGCGCTGACTTTGAACTGATCGAAGTTCAGGCGATCGAGATGATCAACATGACGCATGGCa                                 <  2:5321854/93‑1 (MQ=255)
                                                         gACTCAACAGCGAGGAAGACGTCAGACGCTTTCGCGCTGACTTTGAACTGATCGAAGTTCAGGCGATCGAGATGATCAACATGACGCATGGCa                                 <  2:6974115/93‑1 (MQ=255)
                                                          aCTCAACAGCGAGGAAGACGTCAGACGCTTTCGCGCTGACTTTGAACTGATCGAAGTTCAGGCGATCGAGATGATCAACATGACGCATGGCAg                                >  2:4486667/1‑93 (MQ=255)
                                                          aCTCAACAGCGAGGAAGACGTCAGACGCTTTCGCGCTGACTTTGAACTGATCGAAGTTCAGGCGATCGAGATGATCAACATGACGCATGGCAg                                >  2:6619134/1‑93 (MQ=255)
                                                           cTCAACAGCGAGGAAGACGTCAGACGCTTTCGCGCTGACTTTGAACTGATCGAAGTTCAGGCGATCGAGATGATCAACATGACGCATGGCAGa                               >  1:5655821/1‑93 (MQ=255)
                                                             cAACAGCGAGGAAGACGTCAGACGCTTTCGCGCTGACTTTGAACTGATCGAAGTTCAGGCGATCGAGATGATCAACATGACGCATGGCAGAtt                             >  1:3658131/1‑93 (MQ=255)
                                                              aaCAGCGAGGAAGACGTCAGACGCTTTCGCGCTGACTTTGAACTGATCGAAGTTCAGGCGATCGAGATGATCAACATGACGCATGGCAGATTc                            >  1:354625/1‑93 (MQ=255)
                                                              aaCAGCGAGGAAGACGTCAGACGCTTTCGCGCTGACTTTGAACTGATCGAAGTTCAGGCGATCGAGATGATCAACATGACGCATGGCAGATTc                            <  2:3986380/93‑1 (MQ=255)
                                                               aCAGCGAGGAAGACGTCAGACGCTTTCGCGCTGACTTTGAACTGATCGAAGTTCAGGCGATCGAGATGATCAACATGACGCATGGCAGATTcc                           >  2:6173379/1‑93 (MQ=255)
                                                                cAGCGAGGAAGACGTCAGACGCTTTCGCGCTGACTTTGAACTGATCGAAGTTCAGGCGATCGAGATGATCAACATGACGCATGGCAGATTCCa                          <  2:1024739/93‑1 (MQ=255)
                                                                cAGCGAGGAAGACGTCAGACGCTTTCGCGCTGACTTTGAACTGATCGAAGTTCAGGCGATCGAGATGATCAACATGACGCATGGCAGATTCCa                          <  2:4851088/93‑1 (MQ=255)
                                                                   cGAGGAAGACGTCAGACGCTTTCGCGCTGACTTTGAACTGATCGAAGTTCAGGCGATCGAGATGATCAACATGACGCATGGCAGATTCCAGCa                       >  1:1571418/1‑93 (MQ=255)
                                                                     aGGAAGACGTCAGACGCTTTCGCGCTGACTTTGAACTGATCGAAGTTCAGGCGATCGAGATGATCAACATGACGCATGGCAGATTCCAGCAAc                     >  2:6312847/1‑93 (MQ=255)
                                                                     aGGAAGACGTCAGACGCTTTCGCGCTGACTTTGAACTGATCGAAGTTCAGGCGATCGAGATGATCAACATGACGCATGGCAGATTCCAGCAAc                     >  2:4126614/1‑93 (MQ=255)
                                                                     aGGAAGACGTCAGACGCTTTCGCGCTGACTTTGAACTGATCGAAGTTCAGGCGATCGAGATGATCAACATGACGCATGGCAGATTCCAGCAAc                     <  1:1231868/93‑1 (MQ=255)
                                                                     aGGAAGACGTCAGACGCTTTCGCGCTGACTTTGAACTGATCGAAGTTCAGGCGATCGAGATGATCAACATGACGCATGGCAGATTCCAGCAAc                     >  1:6067283/1‑93 (MQ=255)
                                                                      ggAAGACGTCAGACGCTTTCGCGCTGACTTTGAACTGATCGAAGTTCAGGCGATCGAGATGATCAACATGACGCATGGCAGATTCCAGCAACg                    >  1:7391913/1‑93 (MQ=255)
                                                                      ggAAGACGTCAGACGCTTTCGCGCTGACTTTGAACTGATCGAAGTTCAGGCGATCGAGATGATCAACATGACGCATGGCAGATTCCAGCAACg                    >  2:3784155/1‑93 (MQ=255)
                                                                      ggAAGACGTCAGACGCTTTCGCGCTGACTTTGAACTGATCGAAGTTCAGGCGATCGAGATGATCAACATGACGCATGGCAGATTCCAGCAACg                    >  1:4596499/1‑93 (MQ=255)
                                                                      ggAAGACGTCAGACGCTTTCGCGCTGACTTTGAACTGATCGAAGTTCAGGCGATCGAGATGATCAACATGACGCATGGCAGATTCCAGCAACg                    >  2:28626/1‑93 (MQ=255)
                                                                      ggAAGACGTCAGACGCTTTCGCGCTGACTTTCAACTGATCGAACTTCATGCGATCGAGACGATCAACATGACCCATGGCAGATCCCAGCAACg                    >  1:4698871/1‑93 (MQ=255)
                                                                        aaGACGTCAGACGCTTTCGCGCTGACTTTGAACTGATCGAAGTTCAGGCGATCGAGATGATCAACATGACGCATGGCAGATTCCAGCAACGcc                  <  2:27168/93‑1 (MQ=255)
                                                                          gACGTCAGACGCTTTCGCGCTGACTTTGAACTGATCGAAGTTCAGGCGATCGAGATGATCAACATGACGCATGGCAGATTCCAGCAACGCCTg                <  2:4258085/93‑1 (MQ=255)
                                                                          gACGTCAGACGCTTTCGCGCTGACTTTGAACTGATCGAAGTTCAGGCGATCGAGATGATCAACATGACGCATGGCAGATTCCAGCAACGCCTg                >  2:2066371/1‑93 (MQ=255)
                                                                            cGTCAGACGCTTTCGCGCTGACTTTGAACTGATCGAAGTTCAGGCGATCGAGATGATCAACATGACGCATGGCAGATTCCAGCAACGCCTgcg              >  1:5551013/1‑93 (MQ=255)
                                                                            cGTCAGACGCTTTCGCGCTGACTTTGAACTGATCGAAGTTCAGGCGATCGAGATGATCAACATGACGCATGGCAGATTCCAGCAACGCCTgcg              <  1:765067/93‑1 (MQ=255)
                                                                             gTCAGACGCTTTCGCGCTGACTTTGAACTGATCGAAGTTCAGGCGATCGAGATGATCAACATGACGCATGGCAGATTCCAGCAACGCCTgcgg             >  1:5722997/1‑93 (MQ=255)
                                                                             gTCAGACGCTTTCGCGCTGACTTTGAACTGATCGAAGTTCAGGCGATCGAGATGATCAACATGACGCATGGCAGATTCCAGCAACGCCTgcgg             >  2:726314/1‑93 (MQ=255)
                                                                             gTCAGACGCTTTCGCGCTGACTTTGAACTGATCGAAGTTCAGGCGATCGAGATGATCAACATGACGCATGGCAGATTCCAGCAACGCCTgcgg             >  2:4581055/1‑93 (MQ=255)
                                                                              tCAGACGCTTTCGCGCTGACTTTGAACTGATCGAAGTTCAGGCGATCGAGATGATCAACATGACGCATGGCAGATTCCAGCAACGCCTgcggc            >  2:3233355/1‑93 (MQ=255)
                                                                              tCAGACGCTTTCGCGCTGACTTTGAACTGATCGAAGTTCAGGCGATCGAGATGATCAACATGACGCATGGCAGATTCCAGCAACGCCTgcggc            >  1:858385/1‑93 (MQ=255)
                                                                               cAGACGCTTTCGCGCTGACTTTGAACTGATCGAAGTTCAGGCGATCGAGATGATCAACATGACGCATGGCAGATTCCAGCAACGCCTgcggcg           >  2:7925147/1‑93 (MQ=255)
                                                                               cAGACGCTTTCGCGCTGACTTTGAACTGATCAAAGTTCAGGCGATCGAGATGATCAACATGACGCATGGCAGATTCCAGCAACGCCTgcggcg           >  1:3131344/1‑93 (MQ=255)
                                                                                aGACGCTTTCGCGCTGACTTTGAACTGATCGAAGTTCAGGCGATCGAGATGATCAACATGACGCATGGCAGATTCCAGCAACGCCTGCGGCGt          >  2:5792877/1‑93 (MQ=255)
                                                                                 gACGCTTTCGCGCTGACTTTGAACTGATCGAAGTTCAGGCGATCGAGATGATCAACATGACGCATGGCAGATTCCAGCAACGCCTGCGGCGTc         >  1:7222035/1‑93 (MQ=255)
                                                                                 gACGCTTTCGCGCTGACTTTGAACTGATCGAAGTTCAGGCGATCGAGATGATCAACATGACGCATGGCAGATTCCAGCAACGCCTGCGGCGTc         <  2:2248919/93‑1 (MQ=255)
                                                                                 gACGCTTTCGCGCTGACTTTGAACTGATCGAAGTTCAGGCGATCGAGATGATCAACATGACGCATGGCAGATTCCAGCAACGCCTGCGGCGTc         <  1:580873/93‑1 (MQ=255)
                                                                                  aCGCTTTCGCGCTGACTTTGAACTGATCGAAGTTCAGGCGATCGAGATGATCAACATGACGCATGGCAGATTCCAGCAACGCCTGCGGCGTCg        >  1:6855363/1‑93 (MQ=255)
                                                                                      tttCGCGCTGACTTTGAACTGATCGAAGTTCAGGCGATCGAGATGATCAACATGACGCATGGCAGATTCCAGCAACGCCTGCGGCGTCGGTTc    >  1:4371261/1‑93 (MQ=255)
                                                                                       ttCGCGCTGACTTTGAACTGATCGAAGTTCAGGCGATCGAGATGATCAACATGACGCATGGCAGATTCCAGCAACGCCTGCGGCGTCGGTTCg   >  1:3870165/1‑93 (MQ=255)
                                                                                       ttCGCGCTGACTTTGAACTGATCGAAGTTCAGGCGATCGAGATGATCAACATGACGCATGGCAGATTCCAGCAACGCCTGCGGCGTCGGTTCg   >  2:739906/1‑93 (MQ=255)
                                                                                        tCACGCTGACTTTGAACTGATCGAAGTTCCGGCGATCGAGATGATCAACATGACTCATGGCAGATTCCAGCAACTCCTGCGGCGTCGGTTTGc  >  2:2611660/1‑93 (MQ=255)
                                                                                          |                                                                                          
GCTTCGGTGATCCCCAGATGCAACGGCTGATCGATCTGTTTTGCCAGCAAACGATAAGACTCAACAGCGAGGAAGACGTCAGACGCTTTCACGCTGACTTTGAACTGATCGAAGTTCAGGCGATCGAGATGATCAACATGACGCATGGCAGATTCCAGCAACGCCTGCGGCGTCGGTTCGC  >  NC_000913/2639196‑2639376

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