Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,171,863 A→G 17.2% intergenic (+30/+16) qseC → / ← ygiZ quorum sensing sensory histidine kinase in two‑component regulatory system with QseB/inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,171,8630AG17.2% 414.6 / 77.6 192intergenic (+30/+16)qseC/ygiZquorum sensing sensory histidine kinase in two‑component regulatory system with QseB/inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (81/76);  new base G (16/17);  total (98/94)
Fisher's exact test for biased strand distribution p-value = 8.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.62e-01

Only 100 of 193 total aligned reads displayed.

GCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAA  >  NC_000913/3171767‑3171962
                                                                                                |                                                                                                   
gcatggcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGActc                                                                                                  <  2:1832951/100‑1 (MQ=255)
 catggcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGActct                                                                                                 <  2:272972/100‑1 (MQ=255)
   tggcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCa                                                                                               <  1:1808417/100‑1 (MQ=38)
     gcatgAATGTTGAATTTGGGAATGCGGCACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATa                                                                                             <  2:216613/100‑1 (MQ=255)
      catgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGCGAGCCGTTAAAACGACTTTGATAg                                                                                            >  1:3532803/1‑100 (MQ=255)
      catgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATAg                                                                                            <  2:3744428/100‑1 (MQ=38)
       atgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATAGa                                                                                           <  2:2894994/100‑1 (MQ=38)
       atgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATAGa                                                                                           <  2:665599/100‑1 (MQ=38)
         gAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATAGaca                                                                                         <  2:3571328/100‑3 (MQ=37)
          aaTGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATAGacag                                                                                        <  1:2833408/100‑4 (MQ=37)
          aaTGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGAt                                                                                          >  1:3749567/1‑98 (MQ=255)
           aTGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATAGacaga                                                                                       <  2:59792/100‑5 (MQ=37)
            tGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTa                                                                                      >  2:2809894/1‑100 (MQ=255)
             gTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGAGTAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGAAGTTAAAACGACTCTAATAGATTTTAc                                                                                     >  1:3650821/1‑100 (MQ=255)
              ttGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACt                                                                                    >  1:1256364/1‑100 (MQ=255)
               tGATTTTGGGAATGTGGTACAAGGTGGATTTGAGGCGAAAGTAAGCTGGTATCTTCTCTGTCTATGAAAGCCGTTAAAACGACTCTCATAGATTTTACTa                                                                                   <  1:239408/100‑1 (MQ=255)
               tGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTa                                                                                   <  2:2070735/100‑1 (MQ=255)
                gAATTTGGGAATGCCGAACAAGATGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACt                                                                                    <  2:561227/98‑1 (MQ=255)
                   tttGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAg                                                                               <  2:897493/100‑1 (MQ=255)
                       ggAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCaaa                                                                           >  1:1707752/1‑100 (MQ=255)
                          aTGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAata                                                                        <  2:2589473/100‑1 (MQ=255)
                          aTGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTGTACTAATAGCAAAata                                                                        >  2:966172/1‑100 (MQ=255)
                              ggAACAAGTTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTATATCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATaaa                                                                    <  1:1890339/100‑1 (MQ=255)
                              ggAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATaaa                                                                    <  2:104091/100‑1 (MQ=255)
                              ggAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATaaa                                                                    <  2:2783702/100‑1 (MQ=255)
                               naaCAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAAc                                                                   >  1:2804196/2‑100 (MQ=255)
                               gAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAAc                                                                   >  1:277521/1‑100 (MQ=255)
                               gAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAAc                                                                   <  2:2077863/100‑1 (MQ=255)
                                    aGGTGGGTTTGAGGCGAAGGTAAGCGGGTAACTTCTCTGTCTATGAGAGCCGTTAAACCGACTCGCTTAGATTTTACTAATAGCAAAATATATACCGTcc                                                              <  1:1876406/100‑1 (MQ=255)
                                     ggTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTccc                                                             <  2:1214641/100‑1 (MQ=255)
                                     ggTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTccc                                                             >  2:1698261/1‑100 (MQ=255)
                                     ggTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTccc                                                             <  1:1324347/100‑1 (MQ=255)
                                          aTTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCaaaa                                                        >  2:3525447/1‑100 (MQ=255)
                                             tGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAg                                                     >  1:1842858/1‑100 (MQ=255)
                                             tGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAg                                                     <  2:3105669/100‑1 (MQ=255)
                                               aGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAg                                                     >  1:3695801/1‑98 (MQ=255)
                                                ggCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGcca                                                  >  2:2664445/1‑100 (MQ=255)
                                                ggCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGACGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGcca                                                  >  1:3015696/1‑100 (MQ=255)
                                                 gCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGccac                                                 <  2:2396183/100‑1 (MQ=255)
                                                 gCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGccac                                                 <  2:2001820/100‑1 (MQ=255)
                                                   gAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGccacca                                               >  2:3248117/1‑100 (MQ=255)
                                                   gAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGccacca                                               >  2:3776104/1‑100 (MQ=255)
                                                    aaGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCaccaa                                              <  1:1585815/100‑1 (MQ=255)
                                                     aGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCaccaac                                             >  2:3187285/1‑100 (MQ=255)
                                                      ggTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCaccaacc                                            <  2:2018525/100‑1 (MQ=255)
                                                        tAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACcac                                          >  1:334150/1‑100 (MQ=255)
                                                        tAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACcac                                          >  2:2608584/1‑100 (MQ=255)
                                                        tAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACcac                                          >  1:1110341/1‑100 (MQ=255)
                                                        tAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACcac                                          <  2:3345193/100‑1 (MQ=255)
                                                            cTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAAc                                      <  2:3436187/100‑1 (MQ=255)
                                                             tGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAAcc                                     <  1:1194736/100‑1 (MQ=255)
                                                             tGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCAAAAAAAAAAAAGCCACCAACCACAAAcc                                     <  1:1449484/100‑1 (MQ=255)
                                                             tGGTAACTTCGCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCACCAAAAAAGCCACCAACCACAAAcc                                     <  1:3425126/100‑1 (MQ=255)
                                                              ggTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCa                                    >  2:755597/1‑100 (MQ=255)
                                                              ggTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCa                                    >  2:2463284/1‑100 (MQ=255)
                                                              ggTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCa                                    >  1:2701522/1‑100 (MQ=255)
                                                               gTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAg                                   <  2:3079678/100‑1 (MQ=255)
                                                               gTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAg                                   <  1:2573870/100‑1 (MQ=255)
                                                               gTAACTTCTCTGTCTATCAGACCCGTTAAAACGACACTCACAGAGGTTACGAATAGCAACATCTAAACGGTCCCCAAAAAAGCCACCACCCACAAACCAg                                   <  1:1373250/100‑1 (MQ=255)
                                                                tAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGa                                  >  1:460275/1‑100 (MQ=255)
                                                                tAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGa                                  <  2:3137281/100‑1 (MQ=255)
                                                                 aaCTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGAc                                 <  1:642904/100‑1 (MQ=255)
                                                                 aaCTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGAc                                 <  2:1646683/100‑1 (MQ=255)
                                                                 aaCTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGACTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGAc                                 <  1:2716709/100‑1 (MQ=255)
                                                                  aCTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACa                                >  1:426093/1‑100 (MQ=255)
                                                                  aCTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACa                                >  2:2954801/1‑100 (MQ=255)
                                                                        ctGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATAcc                          <  2:233957/100‑1 (MQ=255)
                                                                        ctGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATAcc                          <  1:2206089/100‑1 (MQ=255)
                                                                          atCTATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcga                        >  2:500266/2‑100 (MQ=37)
                                                                           tCTATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgag                       >  1:2492448/1‑100 (MQ=38)
                                                                           tCTATGAGAGCCGTTAAAACGCCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgag                       <  2:2515727/100‑1 (MQ=255)
                                                                            cTATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATAccaa                        >  1:584498/1‑96 (MQ=37)
                                                                            cTATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgag                       >  1:2939356/1‑99 (MQ=38)
                                                                            cTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagc                      <  2:3419320/100‑1 (MQ=255)
                                                                             tATTAGATCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGACCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcg                     <  1:388326/100‑1 (MQ=255)
                                                                             tATGAGAGCCGTTAAAACGCCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCACAAAACAGACACCACCCACAAACCAGACAAATACcgagcg                     >  1:160196/1‑100 (MQ=255)
                                                                             tATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcg                     >  2:1559842/1‑100 (MQ=255)
                                                                             tATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagag                     <  2:2625375/100‑3 (MQ=255)
                                                                             tATCAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATAcc                          <  1:19442/95‑1 (MQ=255)
                                                                              aTGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcga                    >  2:3346590/1‑100 (MQ=38)
                                                                              aTGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcga                    <  1:733355/100‑1 (MQ=255)
                                                                               tGAGAGCCGTTCAAACGAATCTCATAGATTTTACTAATTGTAAAATATAAACTGTCACCGAAAAAGCTATCAACGACAAACCAGACAAATATcgagcgag                   >  1:1062952/1‑100 (MQ=255)
                                                                               tGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcgag                   >  2:2975145/1‑100 (MQ=255)
                                                                               tGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcgag                   <  2:1110349/100‑1 (MQ=255)
                                                                               tGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAAGAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcgag                   >  1:2072792/1‑100 (MQ=255)
                                                                                gagaGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCgaga                  >  1:955111/1‑100 (MQ=255)
                                                                                gagaGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCgaga                  >  1:2909430/1‑100 (MQ=255)
                                                                                gagaGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAACCACCAACCACAAACCAGACAAATACcgagcgag                   >  1:2300942/1‑99 (MQ=255)
                                                                                   agCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATa               <  1:3304173/100‑1 (MQ=255)
                                                                                   agCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATa               <  1:2006220/100‑1 (MQ=255)
                                                                                   agCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATa               >  1:2820412/1‑100 (MQ=255)
                                                                                    gCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATAc              >  1:1562714/1‑100 (MQ=255)
                                                                                      cGTTAAAACGACTCTCATAGTTTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATAc              <  2:1503470/98‑1 (MQ=255)
                                                                                       gTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCACAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACata           >  1:1074800/1‑100 (MQ=255)
                                                                                          aaaaCGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACAtata         >  1:785794/1‑99 (MQ=255)
                                                                                            aaCGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATAccc      >  1:3159234/1‑100 (MQ=255)
                                                                                            aaCGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATAccc      <  1:1455953/100‑1 (MQ=255)
                                                                                             aCGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCa     <  1:1759124/100‑1 (MQ=255)
                                                                                             aCGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCa     <  1:1753779/100‑1 (MQ=255)
                                                                                               gACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCaaa   <  1:1500618/100‑1 (MQ=255)
                                                                                                |                                                                                                   
GCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAA  >  NC_000913/3171767‑3171962

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: