Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,171,852 C→T 16.5% 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.5% 430.0 / 53.5 189intergenic (+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 (90/63);  new base T (16/16);  total (109/80)
Fisher's exact test for biased strand distribution p-value = 4.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 195 total aligned reads displayed.

GAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACA  >  NC_000913/3171754‑3171951
                                                                                                  |                                                                                                   
gAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGcc                                                                                                    >  1:1323851/1‑100 (MQ=255)
gAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGcc                                                                                                    >  1:272490/1‑100 (MQ=255)
   tCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGTGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGtt                                                                                                 <  1:2979746/100‑1 (MQ=255)
     gtcAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGAAACAAGGTGGATTTGAGGCGAAGATAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTaa                                                                                               <  2:507790/98‑1 (MQ=255)
     gCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTaa                                                                                               >  2:2025445/1‑100 (MQ=255)
     gCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTaa                                                                                               >  1:681948/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)
           ttGCATGGCATGAATGTTGAATTTGGGAATGTGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTa                                                                                                <  1:976567/93‑1 (MQ=255)
           ttGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTACGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACg                                                                                           <  2:2240771/98‑1 (MQ=255)
           ttGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGAc                                                                                         <  1:291501/100‑1 (MQ=255)
            tGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGttttaacggct                                                                                        <  1:111639/100‑10 (MQ=255)
             gcatggnatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGActc                                                                                       >  2:3094079/1‑100 (MQ=255)
              catggcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGActct                                                                                      <  1:1543929/100‑1 (MQ=255)
              catggcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGActc                                                                                       <  2:2066830/99‑1 (MQ=255)
               atggcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACtctc                                                                                     <  2:1969325/100‑1 (MQ=255)
                 ggcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCAt                                                                                   >  1:3176574/1‑100 (MQ=255)
                  gcatgAATGTTGAATTTGGGGATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATa                                                                                  >  2:2391251/1‑100 (MQ=255)
                   catgAATGTTGAATTTGGGAATGCGGAACAAGGTTGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATAg                                                                                 <  1:3926990/100‑1 (MQ=37)
                   catgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATAg                                                                                 <  2:2857399/100‑1 (MQ=38)
                     tgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATAGac                                                                               <  2:383085/100‑2 (MQ=37)
                     tgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATAGac                                                                               <  1:1952507/100‑2 (MQ=37)
                     tgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGAt                                                                               >  2:3238982/1‑100 (MQ=255)
                      gAATGTTGAATTTGGGAATAAAGAACAAGATGGATTAGATGCGAAGGTAAGATGGTAACTTCTCTGTCTATGAGAGCAGTAAAAACGACTCCCATAg                                                                                 >  2:3669066/1‑97 (MQ=255)
                       aaTGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATAGacag                                                                             <  1:1470803/100‑4 (MQ=37)
                          gTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATAGacagaga                                                                          <  1:2153217/100‑7 (MQ=25)
                          gTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTAc                                                                          <  2:2622286/100‑1 (MQ=255)
                          gTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTAc                                                                          >  2:567564/1‑100 (MQ=255)
                          gTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTAc                                                                          <  1:2230561/100‑1 (MQ=255)
                            tGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTa                                                                        >  1:1233797/1‑100 (MQ=255)
                            tGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTa                                                                        >  1:1724712/1‑100 (MQ=255)
                             gAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTaa                                                                       >  2:122513/1‑100 (MQ=255)
                             gAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTaa                                                                       >  2:78761/1‑100 (MQ=255)
                               aTTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATa                                                                     <  1:78764/100‑1 (MQ=255)
                               aTTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATa                                                                     >  1:819852/1‑100 (MQ=255)
                               aTTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGAAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATa                                                                     >  2:2171281/1‑100 (MQ=255)
                                 ttGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGc                                                                   <  1:2402226/100‑1 (MQ=255)
                                    ggAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCaaa                                                                >  1:593062/1‑100 (MQ=255)
                                     gAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCaaaa                                                               >  2:2522908/1‑100 (MQ=255)
                                        tGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAatat                                                            >  2:735559/1‑100 (MQ=255)
                                         gCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTTTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAAtata                                                           <  1:1483177/100‑1 (MQ=255)
                                          cGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATaa                                                          >  1:3377828/1‑100 (MQ=255)
                                           ggAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATaaa                                                         <  2:3182246/100‑1 (MQ=255)
                                            gAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAAc                                                        >  1:2062214/1‑100 (MQ=255)
                                            gAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAAc                                                        >  1:2348904/1‑100 (MQ=255)
                                              aCAATGTGGATTTGAGGAGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCg                                                      <  1:971151/100‑1 (MQ=255)
                                              aCAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGACGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAAGCg                                                      >  2:1712373/1‑100 (MQ=255)
                                                aaGGTTGATTTGAGGCGAAAATAAGCTAGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATCTTACTAATAACAAAATATAAAcc                                                       >  2:1926642/1‑97 (MQ=255)
                                                  ggTGGATTTGAGGCGAAGGTAAGCTAGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAAATTTTACTAATAGCAAAATATAAAACGTc                                                    >  2:2659221/1‑98 (MQ=255)
                                                   gTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTc                                                    >  2:894773/1‑97 (MQ=255)
                                                    tGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCa                                                <  2:1673142/100‑1 (MQ=255)
                                                    tGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTATGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCa                                                >  2:1558025/1‑100 (MQ=255)
                                                    tGGATTTGAGACGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCa                                                >  2:2309602/1‑100 (MQ=255)
                                                      gATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCaaa                                              <  2:3128705/100‑1 (MQ=255)
                                                        tttGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCaaaaa                                            <  2:2218503/100‑1 (MQ=255)
                                                          tGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAg                                          >  2:1829293/1‑100 (MQ=255)
                                                           gAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGc                                         <  1:73835/100‑1 (MQ=255)
                                                             ggCGAAGGTAATCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGc                                         <  2:67136/98‑1 (MQ=255)
                                                             ggCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCGAAATATAAACCGTCCCCAAAAACGc                                         >  2:58192/1‑98 (MQ=255)
                                                             ggCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGcca                                       >  2:96924/1‑100 (MQ=255)
                                                                   ggTAAGCTGGTAACTTCTCTGTCTATGATAGTCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCaccaacc                                 <  1:1836685/100‑1 (MQ=255)
                                                                   ggTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCaccaacc                                 <  2:2392605/100‑1 (MQ=255)
                                                                   ggTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCaccaacc                                 <  2:133678/100‑1 (MQ=255)
                                                                   ggTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCaccaacc                                 <  2:3519953/100‑1 (MQ=255)
                                                                    gTAAGCTGGTATCTTCTCTGTCTATGAGAGTCGTTAAAACTACTCTCATAGATTTTACTACTAGCATCTTATTAACTGTTTCCAAAAAAGCCaccaacc                                 <  2:3434756/99‑1 (MQ=255)
                                                                    gTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCACAATATAAACCGTCCCCAAAAAAGCCaccaacca                                <  2:2348882/100‑1 (MQ=255)
                                                                    gTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCaccaacca                                >  2:1194940/1‑100 (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)
                                                                          tGGTAACTTCTCTGTCTATTCTAGCCGTTAAAACGATTCTCATAGATTTTAATAATAGTAAAATATAAACCGTCCCCAAAAAAGCCACCAACGACAAAcc                          <  2:479963/100‑1 (MQ=255)
                                                                          tGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAAcc                          >  2:3564579/1‑100 (MQ=255)
                                                                          tGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAAcc                          <  1:29083/100‑1 (MQ=255)
                                                                             tAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGa                       <  2:3358760/100‑1 (MQ=255)
                                                                             tAACTTCTCAGCCTATGAGATCCGTTAAACCGACTCTAATACATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCATAAAcctg                        <  2:1044939/99‑3 (MQ=255)
                                                                              aaCTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGAc                      <  1:2418778/100‑1 (MQ=255)
                                                                              aaCTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGAc                      >  2:1083302/1‑100 (MQ=255)
                                                                               aCTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACa                     >  1:2582574/1‑100 (MQ=255)
                                                                               aCTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACa                     >  2:1992156/1‑100 (MQ=255)
                                                                                cTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACaa                    >  1:3586858/1‑100 (MQ=255)
                                                                                cTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACaa                    >  1:3222886/1‑100 (MQ=255)
                                                                                 ttCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACaa                    <  2:3038141/99‑1 (MQ=255)
                                                                                   ctctGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACAAACCACAAACcaaacaaa                   >  2:2038274/1‑98 (MQ=255)
                                                                                      aatCTATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACAGTTCCCAAAAAAGCCAACAACCACAAAACAGAAAAATACcg              >  2:2738713/3‑100 (MQ=21)
                                                                                        tCTATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGACACCAACCACAAACCAGACAAATACcaag            >  2:3922453/1‑100 (MQ=37)
                                                                                         cTATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagc           >  1:3708302/1‑100 (MQ=38)
                                                                                         cTATGAGAGTCGTTTTAACGGCTCTCATAGAATTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAAACAAATACagaac           >  2:3877000/1‑98 (MQ=25)
                                                                                          tATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcg          >  2:836132/1‑100 (MQ=38)
                                                                                          tATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcg          >  1:2194089/1‑100 (MQ=38)
                                                                                          tATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgag            >  2:3818534/1‑98 (MQ=38)
                                                                                          tATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcg          >  1:17743/1‑100 (MQ=255)
                                                                                           aTGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcga         >  1:2787436/1‑100 (MQ=38)
                                                                                           aTGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcga         >  1:3049579/1‑100 (MQ=38)
                                                                                             gagaGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCAACAAAAAAGCAACCAACCACAAACCAGACAAATACAGAGCgaga       >  1:2029432/1‑100 (MQ=37)
                                                                                             gagaGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCgaga       >  2:2456577/1‑100 (MQ=255)
                                                                                             gagaGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCgaga       >  1:938799/1‑100 (MQ=255)
                                                                                              agagTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCAACAACCACAAACCAGACAAATACCGAGCGAGaa      >  2:3593193/1‑100 (MQ=38)
                                                                                               gagCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACAAACCACAAACCAGACAAATACCGAGAGAGAAt     >  2:2098904/1‑100 (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 < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: