Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,171,858 A→T 16.1% intergenic (+25/+21) 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,8580AT16.1% 442.1 / 70.3 179intergenic (+25/+21)qseC/ygiZquorum sensing sensory histidine kinase in two‑component regulatory system with QseB/inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (84/66);  new base T (16/13);  total (100/79)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.08e-01

Only 100 of 194 total aligned reads displayed.

GCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACC  >  NC_000913/3171759‑3171957
                                                                                                   |                                                                                                   
gcgaaaTTGCATGACATGAATGTAAAATTTGGGAATGCCGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATAAGAGCCGTTaa                                                                                                     >  2:2416965/4‑100 (MQ=255)
gCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTaa                                                                                                     >  2:2025445/1‑100 (MQ=255)
gCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTATCTATGAGAGCCGTTaa                                                                                                     >  1:876704/1‑100 (MQ=255)
 ccAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTaaa                                                                                                    <  2:2804051/100‑1 (MQ=255)
  cAAATTGCATTGCATGAATTTTTAATTTGTCAATGCCGAACAAGGTTGATTTGAGGCTAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTaaaa                                                                                                   <  2:3512800/100‑1 (MQ=255)
  cAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTaaaa                                                                                                   >  1:3256261/1‑100 (MQ=255)
   aaaTTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAc                                                                                                  <  2:2541865/100‑1 (MQ=255)
   aaaTTGCATGGCATGAATGTTGAATTTGGGAATGCGGAAAAAGGTGGATTTGAGGAGAAGGTAAGATGGTAACTTGTTTGTCTATGAGAGACGTTaa                                                                                                     >  2:3581917/1‑97 (MQ=255)
      ttGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTACGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACg                                                                                                 <  2:2240771/98‑1 (MQ=255)
      ttGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGAc                                                                                               <  1:291501/100‑1 (MQ=255)
        gcatggcatgAATGTTTAATTTGGGAATGCGGATCAAGGTGGATTTGAGGCTAAGGTAAGCTGGTAACTTTTCTTTCTATGAGGTCCGGTAAAACg                                                                                                 <  2:1574540/96‑1 (MQ=255)
         catggcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGActct                                                                                            >  2:2369164/1‑100 (MQ=255)
          atggcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACtctc                                                                                           <  2:1969325/100‑1 (MQ=255)
            ggcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCa                                                                                          <  2:1021427/99‑1 (MQ=255)
            ggcatgAATGTTGAATTTGGGAACGAGGACCAATGTGGATTTGAGGCGAAGGTAAGGTGCTATCTTCTCTTTCTCTGTGAGTCGTTTTAACGGCTCTCAt                                                                                         <  2:3174074/100‑1 (MQ=255)
             gcatgAATGTTGAATTTGGGGATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATa                                                                                        >  2:2391252/1‑100 (MQ=255)
             gcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTACTCTATCTATGAGAGCCGTTAAAACGACTCTCATa                                                                                        >  1:125495/1‑100 (MQ=255)
             gcatgAATGTTGAATTTGGGAATGCGAAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATa                                                                                        <  2:3734611/100‑1 (MQ=38)
              catgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATAg                                                                                       <  2:2857399/100‑1 (MQ=38)
               atgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATAGa                                                                                      <  1:1483656/100‑1 (MQ=38)
               atgAATGTTGAATTTGGGAATGCGGAAAAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATAg                                                                                       <  2:1060268/99‑1 (MQ=38)
                tgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGAt                                                                                     >  2:3238982/1‑100 (MQ=255)
                 gAATGTTGAATTTGGGAATAAAGAACAAGATGGATTAGATGCGAAGGTAAGATGGTAACTTCTCTGTCTATGAGAGCAGTAAAAACGACTCCCATAg                                                                                       >  2:3669066/1‑97 (MQ=255)
                   aTGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGAtttt                                                                                  <  2:1415883/100‑1 (MQ=255)
                     gTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATAGacagaga                                                                                <  1:2153217/100‑7 (MQ=25)
                     gTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTAc                                                                                >  2:567564/1‑100 (MQ=255)
                     gTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTAc                                                                                <  2:2627317/100‑1 (MQ=255)
                      ttGAATTTGGGTATGTGGAATAAGGTGGATTTGAGGAGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTAc                                                                                <  1:1821255/99‑1 (MQ=255)
                      ttGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACt                                                                               >  1:1947268/1‑100 (MQ=255)
                        gAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTaa                                                                             >  2:78761/1‑100 (MQ=255)
                        gAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACATCTCTGACTATGAGAGCAGTTAAAACGACTCTCATAGATTTTACTaa                                                                             >  2:1865304/1‑100 (MQ=255)
                         aaTTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAAt                                                                            <  2:3083145/100‑1 (MQ=255)
                          aTTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATa                                                                           <  1:78764/100‑1 (MQ=255)
                          aTTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATa                                                                           >  1:819852/1‑100 (MQ=255)
                          aTTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATa                                                                           >  2:495130/1‑100 (MQ=255)
                          aTTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGAAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATa                                                                           >  2:2171281/1‑100 (MQ=255)
                           tttGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAg                                                                          >  1:787444/1‑100 (MQ=255)
                                   tGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAatat                                                                  <  1:666567/100‑1 (MQ=255)
                                   tGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAatat                                                                  >  2:735559/1‑100 (MQ=255)
                                       gAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAAc                                                              >  1:2062214/1‑100 (MQ=255)
                                       gAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAAc                                                              <  1:914558/100‑1 (MQ=255)
                                       gAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAACCGTTAAAACGAATATCATAGATTTTACTAATAGCAAAATATAAAc                                                              >  1:1230201/1‑100 (MQ=255)
                                        aaCAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGCGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAAcc                                                             <  2:3526996/100‑1 (MQ=255)
                                         aCAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGACGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAAGCg                                                            >  2:1712373/1‑100 (MQ=255)
                                         aCAAGGTGGATTTGAAGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGACGTTAAAACGACTCTCATACATTTTAATCATAGCAAAATAAAAccc                                                             >  2:1259268/1‑99 (MQ=255)
                                          cAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGt                                                           >  2:3849459/1‑100 (MQ=255)
                                             ggTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTacc                                                        >  2:2324590/1‑100 (MQ=255)
                                              gTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTc                                                          >  2:894773/1‑97 (MQ=255)
                                               tGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCa                                                      <  2:1673142/100‑1 (MQ=255)
                                                   tttGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGTAAAATATAAACCGTCCCCaaaaa                                                  >  1:554718/1‑100 (MQ=255)
                                                   tttGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCaaaaa                                                  <  2:2218503/100‑1 (MQ=255)
                                                   tttGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCaaaaa                                                  >  2:3821458/1‑100 (MQ=255)
                                                     tGAGTCGAAGGTAAGCTGATAACTTCTCTGTTTATGAAAGCCGTTAAAACTACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAg                                                <  2:195181/100‑1 (MQ=255)
                                                        ggCGAAGGTAATCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGc                                               <  2:67136/98‑1 (MQ=255)
                                                        ggCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCGAAATATAAACCGTCCCCAAAAACGc                                               >  2:58192/1‑98 (MQ=255)
                                                        ggCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAgaga                                             >  1:685750/1‑97 (MQ=255)
                                                        ggCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGcca                                             >  2:96924/1‑100 (MQ=255)
                                                        ggCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGcca                                             >  2:45144/1‑100 (MQ=255)
                                                        ggCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGcca                                             >  1:1471476/1‑100 (MQ=255)
                                                           gAAGGTGAGCTGTTGATTTCTCGGTCTATGAGAGCCGTTAAATCGCCTCTCATAGATGTTATTAATCGCATCATCTAAACTGTCCCCATAAAAGccac                                            <  2:2516232/98‑1 (MQ=255)
                                                            aaGGTAAGCTGGTAACTTCTCTGTCTATGATAGCCGTTAAAACGATTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGccacca                                          <  1:3014774/99‑1 (MQ=255)
                                                             aGGTAAGCTGGTAATTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCaccaac                                        <  1:1731376/100‑1 (MQ=255)
                                                              ggTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGAAAAAAAAAAAGCCaccaacc                                       <  2:3894423/100‑1 (MQ=255)
                                                               gTAAGCTGGTATCTTCTCTGTCTATGAGAGTCGTTAAAACTACTCTCATAGATTTTACTACTAGCATCTTATTAACTGTTTCCAAAAAAGCCaccaacc                                       <  2:3434756/99‑1 (MQ=255)
                                                                tAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACcac                                     >  1:2746212/1‑100 (MQ=255)
                                                                tAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACcac                                     >  2:2908444/1‑100 (MQ=255)
                                                                 aaGCTGGTCACTTCTCTGTCTATGAGAGCCGTTAAAACGACTATCATAGATTTTACTAATAGCCAAATATAAACAGTACACAAAAAATACACAAACcaca                                    >  2:622278/1‑100 (MQ=255)
                                                                   gCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACaaa                                  >  1:12278/1‑100 (MQ=255)
                                                                     tGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAAcc                                <  1:29083/100‑1 (MQ=255)
                                                                      ggTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCa                               >  1:3137797/1‑100 (MQ=255)
                                                                        tAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGa                             <  2:3358760/100‑1 (MQ=255)
                                                                         aaCTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGAc                            <  1:2418778/100‑1 (MQ=255)
                                                                          aCTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACa                           >  1:2582574/1‑100 (MQ=255)
                                                                          aCTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACa                           <  1:2369186/100‑1 (MQ=255)
                                                                          aCTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACa                           >  2:1992156/1‑100 (MQ=255)
                                                                             gtaaaatCTATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAAt                        >  1:3181956/7‑100 (MQ=25)
                                                                              ctctGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATa                       >  2:1756251/1‑100 (MQ=255)
                                                                              ctctGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACAAACCACAAACcaaacaaa                         >  2:2038274/1‑98 (MQ=255)
                                                                                cagTCTATGAGAGCCGTTAAAACGCCTCTCATAGATAAAAAAAAAAACAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATAcc                     <  2:3014148/98‑1 (MQ=255)
                                                                                 aatCTATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACAGTTCCCAAAAAAGCCAACAACCACAAAACAGAAAAATACcg                    >  2:2738713/3‑100 (MQ=21)
                                                                                  gTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATAcc                     >  2:1205431/1‑98 (MQ=255)
                                                                                   tCTATAAGAGACGTTAAGACGATTCTCACAGATTTTACTAATAGCAAAATATAAAACTTCCCCAAAAAAGCCACCTTCCACAAACCAGACAAATACcgag                  <  2:108966/100‑1 (MQ=255)
                                                                                    cTATGCGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAAAATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAAAAcc                     >  2:2713088/1‑96 (MQ=37)
                                                                                    cTATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACAACAAACAAAAAAAATACcga                   >  1:3536713/1‑98 (MQ=25)
                                                                                    cTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagc                 >  1:3696285/1‑100 (MQ=255)
                                                                                     tATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcg                >  1:2194089/1‑100 (MQ=38)
                                                                                     tATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcg                >  2:836132/1‑100 (MQ=38)
                                                                                     tATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCAAAAACCAGACAAATACcgagcg                >  1:1598045/1‑100 (MQ=255)
                                                                                      aTGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcga               >  1:2787436/1‑100 (MQ=38)
                                                                                      aTGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCATCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcga               >  1:2242066/1‑100 (MQ=38)
                                                                                         agagTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCAACAACCACAAACCAGACAAATACCGAGCGAGaa            >  2:3593193/1‑100 (MQ=38)
                                                                                          gagCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAAt           >  1:2405166/1‑100 (MQ=255)
                                                                                          gagCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACAAACCACAAACCAGACAAATACCGAGAGAGAAt           >  2:2098904/1‑100 (MQ=255)
                                                                                           agCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAAtn          <  1:3819811/100‑2 (MQ=255)
                                                                                            gCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATAc         <  1:2014656/100‑1 (MQ=255)
                                                                                             ccGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACa        <  1:96927/100‑1 (MQ=255)
                                                                                             ccGTCAAAACGACTCGCATAGAGTTGAATAATAGCAAAATATAAACCGTCCNCAAGNAAGCCACNAANCACAAACCAGACANATACCGAGCGAGAATa          <  2:3098665/98‑1 (MQ=255)
                                                                                              cGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCATCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAACGAAAATACat       >  1:1451447/1‑100 (MQ=255)
                                                                                               gTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACata      <  2:51369/100‑1 (MQ=255)
                                                                                                   aaaCGATTCTCATAGATTTTACTAGTAGCAAAATATAAACAGTCCACAAAAAAGCCACCAACCACAAACCAGACAATTACCGAGCGAGAATACATATAcc  <  2:319679/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
GCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACC  >  NC_000913/3171759‑3171957

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