Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,171,852 C→T 16.0% intergenic (+19/+27) 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,8520CT16.0% 439.2 / 61.0 186intergenic (+19/+27)qseC/ygiZquorum sensing sensory histidine kinase in two‑component regulatory system with QseB/inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (81/74);  new base T (17/13);  total (98/88)
Fisher's exact test for biased strand distribution p-value = 6.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 189 total aligned reads displayed.

GAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACA  >  NC_000913/3171754‑3171951
                                                                                                  |                                                                                                   
  aTCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGt                                                                                                  <  2:3610232/100‑1 (MQ=255)
   tCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGtt                                                                                                 <  1:2444877/100‑1 (MQ=255)
    cGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGttt                                                                                                <  2:1981935/100‑2 (MQ=255)
    cGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTa                                                                                                >  2:2144121/1‑100 (MQ=255)
     gCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTaa                                                                                               >  1:806117/1‑100 (MQ=255)
       cAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTaaaa                                                                                             >  2:2487999/1‑100 (MQ=255)
         aaTTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAc                                                                                            <  2:357777/99‑1 (MQ=255)
         aaTTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACg                                                                                           <  1:2232593/100‑1 (MQ=255)
             gcatggcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGActc                                                                                       >  1:1068454/1‑100 (MQ=255)
             gcatggcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGActc                                                                                       >  1:1405966/1‑100 (MQ=255)
              catggcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGActct                                                                                      >  2:2537522/1‑100 (MQ=255)
              catggcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGActct                                                                                      >  1:3665220/1‑100 (MQ=255)
              catggcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGAcaga                                                                                      >  1:2346571/1‑97 (MQ=255)
                tggcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCa                                                                                    >  1:1633753/1‑100 (MQ=255)
                tggcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCa                                                                                    >  2:1689093/1‑100 (MQ=255)
                  gcatgAATGTTGAGTTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGGAAGCGGGTAACGGCACTGTGTGGGAGAGACCTTCAGAAGACCCTCAt                                                                                   <  2:3074256/99‑1 (MQ=255)
                  gcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATa                                                                                  >  2:3713453/1‑100 (MQ=255)
                   catgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATAg                                                                                 <  2:2811941/100‑1 (MQ=38)
                   catgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATAg                                                                                 <  1:519154/100‑1 (MQ=38)
                   catgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATAg                                                                                 <  2:1002717/100‑1 (MQ=38)
                   catgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAg                                                                                 <  1:825224/100‑1 (MQ=255)
                   cataAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATa                                                                                  >  2:1165584/1‑99 (MQ=255)
                     tgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGAt                                                                               >  2:2748772/1‑100 (MQ=255)
                     tgAATGTTGAATTTGGGAATGCGGAACAAGGTAGATTTGCGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGAt                                                                               <  2:3239825/100‑1 (MQ=255)
                      gAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATAGaca                                                                              <  1:3381302/100‑3 (MQ=37)
                          gTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTAc                                                                          <  1:1256602/100‑1 (MQ=255)
                          gTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTAc                                                                          <  1:549266/100‑1 (MQ=255)
                          gTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTAc                                                                          <  1:631522/100‑1 (MQ=255)
                          gTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTAc                                                                          <  2:560990/100‑1 (MQ=255)
                          gTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTAc                                                                          >  1:1580186/1‑100 (MQ=255)
                          gTTGAATTTGGGAATGCGCAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCGGTTAAAACGACTCTCATAGATTTTAc                                                                          >  2:455932/1‑100 (MQ=255)
                           ttGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACt                                                                         <  2:1789163/100‑1 (MQ=255)
                             gAATTTGGGAATGCGGAACAAAGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTaa                                                                       <  2:232190/100‑1 (MQ=255)
                                tttGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAg                                                                    <  1:767907/100‑1 (MQ=255)
                                 ttGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGc                                                                   >  1:3419307/1‑100 (MQ=255)
                                    ggAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCaaa                                                                >  2:3147238/1‑100 (MQ=255)
                                     gAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCaaaa                                                               <  2:167553/100‑1 (MQ=255)
                                     gAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCaaaa                                                               <  1:630938/100‑1 (MQ=255)
                                     gAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCaaaa                                                               >  1:216072/1‑100 (MQ=255)
                                     gAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCaaaa                                                               <  1:2388875/100‑1 (MQ=255)
                                        tGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTTATAGCATAatat                                                            >  2:2952361/1‑100 (MQ=255)
                                          cGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATaa                                                          >  2:686711/1‑100 (MQ=255)
                                           ggAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATaaa                                                         >  1:1224040/1‑100 (MQ=255)
                                            gAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAAc                                                        <  2:3665193/100‑1 (MQ=255)
                                            gAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAAc                                                        >  1:543165/1‑100 (MQ=255)
                                             aaCAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAAcc                                                       >  1:3600896/1‑100 (MQ=255)
                                                aaGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTc                                                    <  2:150992/100‑1 (MQ=255)
                                                  ggTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATTAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCg                                                      >  1:1344049/1‑96 (MQ=255)
                                                  ggTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTccc                                                  >  2:338214/1‑100 (MQ=255)
                                                  ggTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTccc                                                  <  2:2738878/100‑1 (MQ=255)
                                                      gATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCaaa                                              >  2:967127/1‑100 (MQ=255)
                                                      gATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCaaa                                              <  1:282993/100‑1 (MQ=255)
                                                        tttGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCaaaaa                                            >  1:1457786/1‑100 (MQ=255)
                                                          tGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAg                                          <  1:3169871/100‑1 (MQ=255)
                                                          tGAGGCGAAGGTAAGCTGGTAACTTATATGTCTATGAGAGTCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAg                                          <  2:3804226/100‑1 (MQ=255)
                                                             ggCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGcca                                       <  2:285019/100‑1 (MQ=255)
                                                             ggCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGcca                                       >  2:1409196/1‑100 (MQ=255)
                                                             ggCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGcc                                        >  2:2365838/1‑99 (MQ=255)
                                                               cGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGccacc                                     >  1:375777/1‑100 (MQ=255)
                                                                gAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGccacca                                    <  2:2572162/100‑1 (MQ=255)
                                                                  aGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCaccaac                                  <  1:1824734/100‑1 (MQ=255)
                                                                   ggTAAGCTGGTAACTGCTCTGTCGATGAGAGTCGTTAAAACGACTCTCATAGATTTTACATATAGCAAAATATCAACCGTACCCTAAAAAGCCaccatcc                                 <  2:3292766/100‑1 (MQ=255)
                                                                    gTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCaccaacca                                <  2:2230655/100‑1 (MQ=255)
                                                                       aGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACaa                             <  2:1628753/100‑1 (MQ=255)
                                                                       aGCTGGTAACTTCTCTGTCTATGAGAGCCATTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACaa                             <  1:3808762/100‑1 (MQ=255)
                                                                        gCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACaaa                            >  2:1052290/1‑100 (MQ=255)
                                                                        gCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACaaa                            >  2:1808079/1‑100 (MQ=255)
                                                                           gttAACTTCTCTGTCTATGAGAGCCGTTAAAACAACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACTACAAAcc                          <  2:2000680/97‑1 (MQ=255)
                                                                           ggTAACTGCTCTGTCTATTGGAGCCGTTAAATCGACACTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAATAAAGCCACCAACCACAAACCa                         <  2:886633/100‑1 (MQ=255)
                                                                            gTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAg                        <  1:2099566/100‑1 (MQ=255)
                                                                            gTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAg                        <  1:1359863/100‑1 (MQ=255)
                                                                             tAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAg                        <  2:822691/99‑1 (MQ=255)
                                                                              aaCTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGAc                      >  1:1960011/1‑100 (MQ=255)
                                                                              aaCTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGAc                      <  1:1030097/100‑1 (MQ=255)
                                                                               acttccatctCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGAc                      <  1:2621405/90‑1 (MQ=255)
                                                                                cTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACaa                    >  1:1910375/1‑100 (MQ=255)
                                                                                   taaaatCTATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATa                 >  2:3494911/6‑100 (MQ=25)
                                                                                     ctGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATAcc               >  1:2498445/1‑100 (MQ=255)
                                                                                      aatCTATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcg              >  1:2709069/3‑100 (MQ=37)
                                                                                       gTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcga             >  1:2673355/1‑100 (MQ=255)
                                                                                       gTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcga             <  1:1889829/100‑1 (MQ=255)
                                                                                       gTCTATGAAAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcga             <  2:3210149/100‑1 (MQ=255)
                                                                                       atCTATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAAtcacga             >  2:1529484/2‑95 (MQ=37)
                                                                                        tCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgag            <  1:2054233/100‑1 (MQ=255)
                                                                                         cTATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagc           >  1:1601782/1‑100 (MQ=38)
                                                                                          tATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcg          >  1:3258528/1‑100 (MQ=255)
                                                                                          tATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcg          >  2:1329759/1‑100 (MQ=255)
                                                                                           aTGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcga         >  1:2968255/1‑100 (MQ=38)
                                                                                           aTGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcga         >  1:1992720/1‑100 (MQ=38)
                                                                                           aTGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcg          >  1:2116328/1‑99 (MQ=38)
                                                                                           aTGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgaacga         >  2:872424/1‑100 (MQ=38)
                                                                                            tGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcgag        >  1:278724/1‑100 (MQ=38)
                                                                                             gagaGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCgaga       >  2:165364/1‑100 (MQ=38)
                                                                                             gagaGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcga         >  1:1691589/1‑98 (MQ=255)
                                                                                             gagaGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCgaga       >  1:3004291/1‑100 (MQ=255)
                                                                                             gagaGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCgaga       >  1:2892710/1‑100 (MQ=255)
                                                                                               gagCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAAt     <  2:521435/100‑1 (MQ=255)
                                                                                               gagCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCACAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAAt     >  2:3217508/1‑100 (MQ=255)
                                                                                                agCCGTTAAAACGACTCTCATAGATTTTACTAATAACAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATa    <  2:1344040/100‑1 (MQ=255)
                                                                                                  ccGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACa  <  2:1289499/100‑1 (MQ=255)
                                                                                                  |                                                                                                   
GAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACA  >  NC_000913/3171754‑3171951

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: