Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,171,857 A→T 16.0% 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,8570AT16.0% 450.0 / 66.7 179intergenic (+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 (83/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 = 7.00e-01

Only 100 of 195 total aligned reads displayed.

CGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATAC  >  NC_000913/3171758‑3171956
                                                                                                   |                                                                                                   
cGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTa                                                                                                     >  1:2056243/1‑100 (MQ=255)
cGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGAACTTCTCTGTCTATGAGAGCCGTTa                                                                                                     <  2:205991/99‑1 (MQ=255)
 gCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTaa                                                                                                    >  1:681948/1‑100 (MQ=255)
 gCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTATCTATGAGAGCCGTTaa                                                                                                    >  1:876704/1‑100 (MQ=255)
  ccAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTaaa                                                                                                   <  2:2804051/100‑1 (MQ=255)
    aaaTTGCATGGCATGAATGTTGAATTTGGGAATGCGGAAAAAGGTGGATTTGAGGAGAAGGTAAGATGGTAACTTGTTTGTCTATGAGAGACGTTaa                                                                                                    >  2:3581917/1‑97 (MQ=255)
       ttGCATGGCATGAATGTTGAATTTGGGAATGTGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTa                                                                                                     <  1:976567/93‑1 (MQ=255)
       ttGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTACGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACg                                                                                                <  2:2240771/98‑1 (MQ=255)
         gcatggcatgAATGTTTAATTTGGGAATGCGGATCAAGGTGGATTTGAGGCTAAGGTAAGCTGGTAACTTTTCTTTCTATGAGGTCCGGTAAAACg                                                                                                <  2:1574540/96‑1 (MQ=255)
           atggcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCtctc                                                                                          <  1:2863681/100‑1 (MQ=38)
             ggcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCAt                                                                                        >  1:3176574/1‑100 (MQ=255)
             ggcatgAATGTTGAATTTGGGAACGAGGACCAATGTGGATTTGAGGCGAAGGTAAGGTGCTATCTTCTCTTTCTCTGTGAGTCGTTTTAACGGCTCTCAt                                                                                        <  2:3174074/100‑1 (MQ=255)
              gcatgAATGTTGAATTTGGGGATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATa                                                                                       >  2:2391251/1‑100 (MQ=255)
              gcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTACTCTATCTATGAGAGCCGTTAAAACGACTCTCATa                                                                                       >  1:125495/1‑100 (MQ=255)
              gcatgAATGTTGAATTTGGGAATGCGAAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATa                                                                                       <  2:3734611/100‑1 (MQ=38)
               catgAATGTTGAATTTGGGAATGCGGAACAAGGTTGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATAg                                                                                      <  1:3926990/100‑1 (MQ=37)
               catgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATAg                                                                                      <  2:2857399/100‑1 (MQ=38)
                atgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATAGa                                                                                     <  1:1483656/100‑1 (MQ=38)
                atgAATGTTGAATTTGGGAATGCGGAAAAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATAg                                                                                      <  2:1060268/99‑1 (MQ=38)
                    aTGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGAtttt                                                                                 <  2:1415883/100‑1 (MQ=255)
                     tGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTa                                                                                >  2:3607088/1‑100 (MQ=255)
                      gTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATAGacagaga                                                                               <  1:2910794/100‑7 (MQ=25)
                      gTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTAc                                                                               <  1:2230561/100‑1 (MQ=255)
                      gTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTAc                                                                               >  2:567564/1‑100 (MQ=255)
                      gTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTAc                                                                               <  2:2622286/100‑1 (MQ=255)
                       ttGAATTTGGGTATGTGGAATAAGGTGGATTTGAGGAGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTAc                                                                               <  1:1821255/99‑1 (MQ=255)
                       ttGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACt                                                                              <  2:3501183/100‑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)
                             ttGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGc                                                                        <  1:2402226/100‑1 (MQ=255)
                                ggAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCaaa                                                                     >  1:593062/1‑100 (MQ=255)
                                 gAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCaaaa                                                                    >  2:2522908/1‑100 (MQ=255)
                                    tGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAatat                                                                 >  1:2971125/1‑100 (MQ=255)
                                    tGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAatat                                                                 >  2:735559/1‑100 (MQ=255)
                                      cGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATaa                                                               >  1:3377828/1‑100 (MQ=255)
                                      cGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATaa                                                               >  2:2014638/1‑100 (MQ=255)
                                       ggAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATaaa                                                              >  2:3236969/1‑100 (MQ=255)
                                        gAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAAc                                                             <  2:2037432/100‑1 (MQ=255)
                                        gAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAAc                                                             >  1:2348904/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)
                                              ggTGGATTTGAGGCGAAGGTAAGCTAGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAAATTTTACTAATAGCAAAATATAAAACGTc                                                         >  2:2659221/1‑98 (MQ=255)
                                               gTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTc                                                         >  2:894773/1‑97 (MQ=255)
                                                tGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCa                                                     <  2:1673142/100‑1 (MQ=255)
                                                    tttGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCaaaaa                                                 <  2:2218503/100‑1 (MQ=255)
                                                    tttGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCaaaaa                                                 >  2:3821458/1‑100 (MQ=255)
                                                      tGAGTCGAAGGTAAGCTGATAACTTCTCTGTTTATGAAAGCCGTTAAAACTACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAg                                               <  2:195181/100‑1 (MQ=255)
                                                       gAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGc                                              <  1:73835/100‑1 (MQ=255)
                                                         ggCGAAGGTAATCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGc                                              <  2:67136/98‑1 (MQ=255)
                                                         ggCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCGAAATATAAACCGTCCCCAAAAACGc                                              >  2:58192/1‑98 (MQ=255)
                                                         ggCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGcca                                            >  1:1471476/1‑100 (MQ=255)
                                                         ggCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGcca                                            >  2:45144/1‑100 (MQ=255)
                                                         ggCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGcca                                            >  2:96924/1‑100 (MQ=255)
                                                              aGGTAAGCTGGTAATTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCaccaac                                       <  1:1731376/100‑1 (MQ=255)
                                                               ggTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCaccaacc                                      <  2:2392605/100‑1 (MQ=255)
                                                               ggTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCaccaacc                                      <  2:3519953/100‑1 (MQ=255)
                                                               ggTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCaccaacc                                      >  2:3393276/1‑100 (MQ=255)
                                                               ggTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGAAAAAAAAAAAGCCaccaacc                                      <  2:3894423/100‑1 (MQ=255)
                                                                gTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCACAATATAAACCGTCCCCAAAAAAGCCaccaacca                                     <  2:2348882/100‑1 (MQ=255)
                                                                 tAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACcac                                    >  2:2908444/1‑100 (MQ=255)
                                                                 tAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACcac                                    >  1:2746212/1‑100 (MQ=255)
                                                                  aaGCTGGTCACTTCTCTGTCTATGAGAGCCGTTAAAACGACTATCATAGATTTTACTAATAGCCAAATATAAACAGTACACAAAAAATACACAAACcaca                                   >  2:622278/1‑100 (MQ=255)
                                                                    gCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACaaa                                 >  1:12278/1‑100 (MQ=255)
                                                                         tAACTTCTCAGCCTATGAGATCCGTTAAACCGACTCTAATACATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCATAAAcctg                             <  2:1044939/99‑3 (MQ=255)
                                                                          aaCTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGAc                           <  1:2071463/100‑1 (MQ=255)
                                                                          aaCTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGAc                           <  1:2418778/100‑1 (MQ=255)
                                                                           aCTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACa                          >  1:2582574/1‑100 (MQ=255)
                                                                            cTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACaa                         >  1:3586858/1‑100 (MQ=255)
                                                                            cTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACaa                         >  1:3222886/1‑100 (MQ=255)
                                                                             ttCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACaaa                        >  2:3247091/1‑100 (MQ=255)
                                                                               ctctGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATa                      >  2:1756251/1‑100 (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)
                                                                                     cTATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagc                >  1:3708302/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)
                                                                                       aTGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcga              >  1:3049579/1‑100 (MQ=38)
                                                                                        tGAGAGTCTTTTTAACGGCTCGCATATATTATACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagc                <  2:1970987/97‑1 (MQ=255)
                                                                                         gagaGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCgaga            >  2:2456577/1‑100 (MQ=255)
                                                                                           gagCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACAAACCACAAACCAGACAAATACCGAGAGAGAAt          >  2:2098904/1‑100 (MQ=255)
                                                                                             gCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATAc        <  1:2014656/100‑1 (MQ=255)
                                                                                              ccGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACa       <  1:96927/100‑1 (MQ=255)
                                                                                              ccGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACa       <  1:2179644/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)
                                                                                                   |                                                                                                   
CGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATAC  >  NC_000913/3171758‑3171956

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: