Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,332,868 G→T 16.2% R135R (CGC→CGA basS ← sensory histidine kinase in two‑component regulatory system with BasR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,332,8680GT16.2% 495.2 / 81.3 192R135R (CGC→CGAbasSsensory histidine kinase in two‑component regulatory system with BasR
Reads supporting (aligned to +/- strand):  ref base G (77/84);  new base T (17/14);  total (94/98)
Fisher's exact test for biased strand distribution p-value = 5.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.27e-01

Only 100 of 203 total aligned reads displayed.

CAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCC  >  NC_000913/4332769‑4332967
                                                                                                   |                                                                                                   
cagcagTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAgg                                                                                                     >  2:748312/1‑100 (MQ=255)
  gcagTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCATTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCg                                                                                                   <  2:2168009/100‑1 (MQ=255)
  gcagTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCg                                                                                                   >  1:11532/1‑100 (MQ=255)
      ttCCAGATGCAAACGCACCCCAGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTg                                                                                               <  2:349761/100‑1 (MQ=255)
        ccAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGAc                                                                                             >  1:1469809/1‑100 (MQ=255)
        ccAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGAc                                                                                             >  1:3518583/1‑100 (MQ=255)
         cAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGAcc                                                                                            >  1:2791720/1‑100 (MQ=255)
         cAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGAcc                                                                                            >  1:3066199/1‑100 (MQ=255)
         cAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGAcc                                                                                            >  1:2404642/1‑100 (MQ=255)
         cAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGAcc                                                                                            >  2:1554685/1‑100 (MQ=255)
          aGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCa                                                                                           <  2:2284193/100‑1 (MQ=255)
           gATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAg                                                                                          >  1:1459609/1‑100 (MQ=255)
            aTGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGTCGCCTGACCAgc                                                                                         <  1:877780/100‑2 (MQ=255)
             tGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGa                                                                                         >  1:411889/1‑99 (MQ=255)
             tGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGAt                                                                                        >  2:1853959/1‑100 (MQ=255)
             tGCAAACGCACCCCCGACAGTGGCGTTCGTAGTTCGTGCGCGACGTCCGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAgct                                                                                        >  2:3361385/1‑98 (MQ=255)
              gCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGTCGCCTGACCAgcac                                                                                       <  2:3247159/100‑4 (MQ=255)
              gCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGTCGCCTGACCAgcac                                                                                       <  2:2805919/100‑4 (MQ=255)
              gCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGTCGCCTGACCAgcac                                                                                       <  1:1254864/100‑4 (MQ=255)
               cAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGTCGCCTGACCAgcacg                                                                                      <  1:3645707/100‑5 (MQ=255)
               cAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGTCGCCTGACCAgcacg                                                                                      <  1:1644126/100‑5 (MQ=255)
               cAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCg                                                                                      >  2:2036139/1‑100 (MQ=255)
                aaaCGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGTCGCCTGACCAgcacgc                                                                                     <  2:899412/100‑6 (MQ=255)
                aaaCGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGt                                                                                     <  1:2119195/100‑1 (MQ=255)
                      aCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACg                                                                               <  2:2618634/100‑1 (MQ=255)
                       cccccGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGc                                                                              >  1:1570257/1‑100 (MQ=255)
                        ccccGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGcc                                                                             >  1:1805682/1‑100 (MQ=255)
                        ccccGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGcc                                                                             >  1:3521980/1‑100 (MQ=255)
                        ccccGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGcc                                                                             <  2:2608193/100‑1 (MQ=255)
                         cccGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCg                                                                            <  1:710664/100‑1 (MQ=255)
                           cGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGaa                                                                          <  2:1774716/100‑1 (MQ=255)
                            gCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCCCTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGGCGaaa                                                                         >  2:190301/1‑100 (MQ=255)
                              naGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCCTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAacc                                                                       >  2:567679/2‑100 (MQ=255)
                              cAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAacc                                                                       <  2:2024016/100‑1 (MQ=255)
                              cAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAacc                                                                       >  2:568024/1‑100 (MQ=255)
                               aGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAacca                                                                      >  2:2878019/1‑100 (MQ=255)
                               aGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAacca                                                                      >  2:846445/1‑100 (MQ=255)
                                 tGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAaccacc                                                                    <  2:2566960/100‑1 (MQ=255)
                                     gTTCGCAGTTCGTGCGAGACGTCAGCGGTAAACAACCTTTCGTTATATAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAGCGCCGAAACCACCGc                                                                  <  2:2944976/98‑1 (MQ=255)
                                       tCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCg                                                              <  1:3122343/100‑1 (MQ=255)
                                          cAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGAtt                                                           <  1:107204/100‑1 (MQ=255)
                                           aGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGAttt                                                          >  1:2383132/1‑100 (MQ=255)
                                               cGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAg                                                      <  2:2055162/100‑1 (MQ=255)
                                                gTGCGCGACGTCAGCGGTAAACATCCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAgg                                                     <  2:3333967/100‑1 (MQ=255)
                                                 tGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAggg                                                    <  2:444637/100‑1 (MQ=255)
                                                  gcgcGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGt                                                   >  1:768248/1‑100 (MQ=255)
                                                   cgcgACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTg                                                  >  1:2775530/1‑100 (MQ=255)
                                                       aCGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGgcgc                                              <  2:927596/100‑1 (MQ=255)
                                                        cGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAGCGCCGAAACCACCGCTTCGATTTCGAGGGTGgcgc                                              <  2:1078264/99‑1 (MQ=255)
                                                        cGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCt                                             >  1:1296005/1‑100 (MQ=255)
                                                         gTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCtg                                            >  2:1061252/1‑100 (MQ=255)
                                                           cAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCtgtg                                          <  2:7745/100‑1 (MQ=255)
                                                           cAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCtgtg                                          <  2:1570249/100‑1 (MQ=255)
                                                             gCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGaa                                        <  1:2044693/100‑1 (MQ=255)
                                                             gCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGaa                                        >  2:1479581/1‑100 (MQ=255)
                                                              cGGTAAACAACCTTTCGTTATCCAGCATGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAAt                                       >  1:2216694/1‑100 (MQ=255)
                                                               ggTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATg                                      <  2:3439765/100‑1 (MQ=255)
                                                               ggTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATg                                      <  1:2714363/100‑1 (MQ=255)
                                                                gTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATgg                                     <  1:3197722/100‑1 (MQ=255)
                                                                 tAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGc                                    <  1:3439590/100‑1 (MQ=255)
                                                                 tAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGc                                    <  1:2881281/100‑1 (MQ=255)
                                                                 tAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGc                                    >  2:1451184/1‑100 (MQ=255)
                                                                 tAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGc                                    >  2:535253/1‑100 (MQ=255)
                                                                 tAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGc                                    <  1:1233367/100‑1 (MQ=255)
                                                                  aaaCAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCa                                   <  1:956269/100‑1 (MQ=255)
                                                                   aacaacCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGGGTAACGCCGAAACCACCGCGTCGATGTCGAGGGTGGCGCTGTGAATGGCa                                   >  2:3039158/1‑99 (MQ=255)
                                                                       acCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATggg                              >  2:2087701/1‑100 (MQ=255)
                                                                        ccTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGACACCACCGCTTCGATTTCGAGGGTGGGGCTGTGAATGGCAATGGGc                             >  2:3026733/1‑100 (MQ=255)
                                                                        ccTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGc                             <  2:3624635/100‑1 (MQ=255)
                                                                            tCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGt                           <  2:1249323/98‑1 (MQ=255)
                                                                              gTTATCCTGCGTGCTGGTCAGGCGACTGTCTAGATTGTTGAGCGCCGAAACCGCCGCTTCGATTTCGTGGGTGGCGCTGTGAATGGCAATGGGCGTTaa                        <  1:1943092/99‑1 (MQ=255)
                                                                              gTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAg                       >  1:983513/1‑100 (MQ=255)
                                                                              gTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAg                       <  1:557959/100‑1 (MQ=255)
                                                                                tATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGtt                     >  1:3623594/1‑100 (MQ=255)
                                                                                 gttaaacgatcTGGTCAGTCGCCTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTg                    >  1:1883461/11‑100 (MQ=38)
                                                                                   ccAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTc                  >  2:74347/1‑100 (MQ=255)
                                                                                   ccAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTc                  <  1:2576046/100‑1 (MQ=255)
                                                                                    cAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTcg                 <  2:2685328/100‑1 (MQ=255)
                                                                                       cgatcTGGTCAGTCGCCTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTcggcg              >  2:1167296/5‑100 (MQ=255)
                                                                                       cgatcTGGTCAGTCGCCTGACCAGATCGTTTAACGCCGAAACCACAGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTcggcg              >  1:642375/5‑100 (MQ=255)
                                                                                        gatcTGGTCAGTCGCCTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTcggcgg             >  1:2859471/4‑100 (MQ=255)
                                                                                        gatcTGGTCAGTCGCCTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTcggcgg             >  1:293118/4‑100 (MQ=255)
                                                                                        gatcTGGTCAGTCGCCTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTcggcgg             >  1:309240/4‑100 (MQ=255)
                                                                                        gatcTGGTCAGTCGCCTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTcggcgg             >  1:46518/4‑100 (MQ=255)
                                                                                        gTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTcggcgg             <  2:3204314/100‑1 (MQ=255)
                                                                                         atcTGGTCAGTCGCCTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGt            >  1:1003657/3‑100 (MQ=255)
                                                                                          gCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTg           <  2:1219156/100‑1 (MQ=255)
                                                                                          gCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACAACCGCGTCGATTTCGAGGGTGGCGCTGTGAAGGGAAATGGGCGTTAAGTTATCGGCGGTg           >  2:1731131/1‑100 (MQ=255)
                                                                                           cTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTgc          <  1:2894814/100‑1 (MQ=255)
                                                                                           cTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTgc          >  1:2152977/1‑100 (MQ=255)
                                                                                            tGGTCAGGCGACTGACCAGATCGTTTAACGCCGAATTCTCCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGTAATGGGCGTTAAGTTGTCGGCGGTgc          <  1:2272021/99‑1 (MQ=255)
                                                                                            tGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTgcg         <  1:2341562/100‑1 (MQ=255)
                                                                                             ggTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGACTGGCAATGGGCGTTAAGTTGTCGGCGGTgcg         >  2:1683230/1‑99 (MQ=255)
                                                                                             ggTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTgcgc        >  2:1826561/1‑100 (MQ=255)
                                                                                             ggTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTgcgc        <  1:3417221/100‑1 (MQ=255)
                                                                                             ggTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTCGTCGGCGGTgcgc        >  1:1376882/1‑100 (MQ=255)
                                                                                               tCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTgcgcgc      <  1:2681310/100‑1 (MQ=255)
                                                                                                  ggCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTc   >  1:3484819/1‑100 (MQ=255)
                                                                                                   gCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTcc  <  2:1295149/100‑1 (MQ=255)
                                                                                                   gCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGCGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTcc  <  2:309412/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
CAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCC  >  NC_000913/4332769‑4332967

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