Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,171,857 A→T 17.6% intergenic (+24/+22) 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,8570AT17.6% 444.1 / 86.8 182intergenic (+24/+22)qseC/ygiZquorum sensing sensory histidine kinase in two‑component regulatory system with QseB/inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (76/73);  new base T (16/16);  total (93/89)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.75e-02

Only 100 of 190 total aligned reads displayed.

GCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATAC  >  NC_000913/3171759‑3171956
                                                                                                  |                                                                                                   
 ccAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTaaa                                                                                                   >  1:2776774/1‑100 (MQ=255)
 ccAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTaaa                                                                                                   >  2:978850/1‑100 (MQ=255)
   aaaTTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAc                                                                                                 <  2:3401814/100‑1 (MQ=255)
        gcatggcatgAATGTTGAATTTGGGAATGCGTACCACGGAGGATTTGAGGCGATGGTTAGCTGGTAACTTCTCTGTATCTGAGAGCGGTTaaa                                                                                                   >  1:1848840/1‑93 (MQ=255)
        gcatggcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGActc                                                                                            >  1:1453183/1‑100 (MQ=255)
        gcatggcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGActc                                                                                            >  2:2766396/1‑100 (MQ=255)
        gcatggcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGActc                                                                                            <  2:1832951/100‑1 (MQ=255)
        gcatggcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGAct                                                                                             >  1:3156445/1‑99 (MQ=255)
         catggcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGActct                                                                                           <  2:3653598/100‑1 (MQ=255)
            ggcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCan                                                                                        <  2:3552696/100‑2 (MQ=38)
              catgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGCGAGCCGTTAAAACGACTTTGATAg                                                                                      >  1:3532803/1‑100 (MQ=255)
              catgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATAg                                                                                      <  2:1080407/100‑1 (MQ=38)
              catgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATAg                                                                                      <  2:3744428/100‑1 (MQ=38)
               atgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATAGa                                                                                     <  2:2930622/100‑1 (MQ=38)
               atgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATAGa                                                                                     <  2:665599/100‑1 (MQ=38)
                  aaTGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGAt                                                                                    >  1:3749567/1‑98 (MQ=255)
                   aTGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATAGacaga                                                                                 <  2:59792/100‑5 (MQ=37)
                     gTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATAGacagaga                                                                               <  2:29859/100‑7 (MQ=25)
                     gTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTAc                                                                               <  1:1780966/100‑1 (MQ=255)
                      ttGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACt                                                                              <  2:2805705/100‑1 (MQ=255)
                        gAATTTGGGAATGCCGAACAAGATGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACt                                                                              <  2:561227/98‑1 (MQ=255)
                        gAATTCGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTaa                                                                            >  1:1214633/1‑100 (MQ=255)
                         aaTTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAAt                                                                           >  1:2283294/1‑100 (MQ=255)
                         aaTTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAAt                                                                           >  1:2764643/1‑100 (MQ=255)
                          aTTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATa                                                                          >  1:1505099/1‑100 (MQ=255)
                          aTTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATa                                                                          >  1:431242/1‑100 (MQ=255)
                          aTTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATa                                                                          >  2:3093565/1‑100 (MQ=255)
                          aTTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATa                                                                          <  1:238743/100‑1 (MQ=255)
                           tttGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAg                                                                         <  2:897493/100‑1 (MQ=255)
                              gggAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCaa                                                                      >  1:3551081/1‑100 (MQ=255)
                               ggAATGCGTAACAAGTGGGAATTGAGGCGAAGGTAAGCTGGTAACTCCTGTGTCTATTAGAGCAGCGAAAACGACTCTCATCGAATTTACTCATAGCa                                                                       >  1:2919013/1‑98 (MQ=255)
                               ggAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCaaa                                                                     >  1:1707752/1‑100 (MQ=255)
                                  aTGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTGCTATTAGTAAAatc                                                                  <  1:786506/100‑2 (MQ=255)
                                  aTGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAata                                                                  <  2:2589473/100‑1 (MQ=255)
                                  aTGCGGAACAAGGCGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAata                                                                  <  1:508555/100‑1 (MQ=255)
                                   tGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAatat                                                                 >  2:2011654/1‑100 (MQ=255)
                                   tGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAataa                                                                 >  2:2188178/1‑99 (MQ=255)
                                      ggAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATaaa                                                              <  1:2428842/100‑1 (MQ=255)
                                       naaCAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAAc                                                             >  1:2804196/2‑100 (MQ=255)
                                        aaCAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTGTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAAcc                                                            <  1:2947141/100‑1 (MQ=255)
                                            aGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTcn                                                        <  1:71777/100‑2 (MQ=255)
                                            aGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTc                                                         >  1:144599/1‑99 (MQ=255)
                                             ggTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTccc                                                       >  2:1698261/1‑100 (MQ=255)
                                                     tGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAg                                               <  1:606347/100‑1 (MQ=255)
                                                     tGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAg                                               <  2:3105669/100‑1 (MQ=255)
                                                       aGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAg                                               >  1:3695801/1‑98 (MQ=255)
                                                        ggCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGcca                                            >  2:2664445/1‑100 (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:2701229/100‑1 (MQ=255)
                                                              ggTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCaccaacc                                      <  2:2839486/100‑1 (MQ=255)
                                                                tAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACcac                                    <  2:3345193/100‑1 (MQ=255)
                                                                tAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACcac                                    >  1:1110341/1‑100 (MQ=255)
                                                                tAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACcac                                    >  2:2608584/1‑100 (MQ=255)
                                                                tAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACcac                                    >  1:334150/1‑100 (MQ=255)
                                                                    cTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAAc                                <  2:3436187/100‑1 (MQ=255)
                                                                     tGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAAcc                               >  1:753775/1‑100 (MQ=255)
                                                                     tGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAAcc                               <  1:1194736/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                             <  1:2573870/100‑1 (MQ=255)
                                                                        tAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGa                            <  2:3137281/100‑1 (MQ=255)
                                                                         aaCTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGAc                           <  2:203590/100‑1 (MQ=255)
                                                                         aaCTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGAc                           <  2:1646683/100‑1 (MQ=255)
                                                                         aaCTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGAc                           <  2:1025121/100‑1 (MQ=255)
                                                                               tctGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATAc                     <  1:397130/100‑1 (MQ=255)
                                                                               aaaatCTATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATGAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACaaa                        >  1:1834877/5‑97 (MQ=25)
                                                                               aaaatCTATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATAc                     >  2:1135770/5‑100 (MQ=25)
                                                                                ctGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATAcc                    <  2:2383839/100‑1 (MQ=255)
                                                                                 aatCTATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCTCCAAAAAAGCCACCAACCACAAACCAGACAAATACcg                   >  1:2810430/3‑100 (MQ=37)
                                                                                 aatCTATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcg                   >  2:2071797/3‑100 (MQ=37)
                                                                                  atCTATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcga                  >  2:500266/2‑100 (MQ=37)
                                                                                   tCTATGAGAGCCGTTAAAACGCCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgag                 <  2:2515727/100‑1 (MQ=255)
                                                                                   ncTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgag                 >  1:1568185/2‑100 (MQ=255)
                                                                                    cTATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagc                >  2:1308879/1‑100 (MQ=38)
                                                                                    cTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagc                >  2:1530375/1‑100 (MQ=255)
                                                                                    cTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagc                <  2:3419320/100‑1 (MQ=255)
                                                                                     tATGAGAGCCGTTAAAACGCCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCACAAAACAGACACCACCCACAAACCAGACAAATACcgagcg               >  1:160196/1‑100 (MQ=255)
                                                                                     tATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcg               >  2:2151701/1‑100 (MQ=255)
                                                                                     tATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagc                <  1:1881565/99‑1 (MQ=255)
                                                                                     tATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagag               <  2:2625375/100‑3 (MQ=255)
                                                                                     tATCAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATAcc                    <  1:19442/95‑1 (MQ=255)
                                                                                      aTGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCGAAATATAAACCGTCCCCAAAAAAGCCACCAACAACAAACCAGACAAATACcgagcga              >  1:2507804/1‑100 (MQ=37)
                                                                                      aTGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCAACCAACAACAAACAAGACAAATACagag                 >  1:2828273/1‑97 (MQ=25)
                                                                                       tGAGAGCCGTTCAAACGAATCTCATAGATTTTACTAATTGTAAAATATAAACTGTCACCGAAAAAGCTATCAACGACAAACCAGACAAATATcgagcgag             >  1:1062952/1‑100 (MQ=255)
                                                                                       tGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCTTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcgag             <  1:2148485/100‑1 (MQ=255)
                                                                                       tGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAAGAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcgag             >  1:2072792/1‑100 (MQ=255)
                                                                                           agCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATa         <  1:2006220/100‑1 (MQ=255)
                                                                                           agCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATa         <  1:3304173/100‑1 (MQ=255)
                                                                                               gTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACata     >  1:1274060/1‑100 (MQ=255)
                                                                                               gTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCACAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACata     >  1:1074800/1‑100 (MQ=255)
                                                                                                 tAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACAtata   >  2:2955815/1‑100 (MQ=255)
                                                                                                 tAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACAtata   >  1:2862299/1‑100 (MQ=255)
                                                                                                  aaaaCGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACAtata   >  1:785794/1‑99 (MQ=255)
                                                                                                  aaaaCGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATAc  <  1:1331590/100‑1 (MQ=255)
                                                                                                  |                                                                                                   
GCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATAC  >  NC_000913/3171759‑3171956

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: