Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,171,858 A→T 16.2% intergenic (+25/+21) qseC → / ← ygiZ quorum sensing sensory histidine kinase in two‑component regulatory system with QseB/inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,171,8580AT16.2% 324.5 / 46.8 128intergenic (+25/+21)qseC/ygiZquorum sensing sensory histidine kinase in two‑component regulatory system with QseB/inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (65/42);  new base T (11/10);  total (76/52)
Fisher's exact test for biased strand distribution p-value = 4.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 134 total aligned reads displayed.

GCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACC  >  NC_000913/3171759‑3171957
                                                                                                   |                                                                                                   
gCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTaa                                                                                                     >  1:1002275/1‑100 (MQ=255)
gCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTaa                                                                                                     <  2:83770/100‑1 (MQ=255)
gCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTaa                                                                                                     >  1:160202/1‑100 (MQ=255)
gCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTaa                                                                                                     >  2:2307487/1‑100 (MQ=255)
gCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTaa                                                                                                     >  1:1431050/1‑100 (MQ=255)
    aaTTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAc                                                                                                  <  2:2143705/99‑1 (MQ=255)
         catggcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGActct                                                                                            <  2:1479530/100‑1 (MQ=255)
         catggcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGActct                                                                                            >  1:2401495/1‑100 (MQ=255)
         catggcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGActct                                                                                            >  1:1554228/1‑100 (MQ=255)
         catggcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGActct                                                                                            >  1:1554134/1‑100 (MQ=255)
         catggcatgAATGTTGAATTTGGGAATGCGGAACAAGCTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGActct                                                                                            >  2:1364635/1‑100 (MQ=255)
           tggcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCa                                                                                          <  1:82456/100‑1 (MQ=38)
           tggcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCa                                                                                          <  2:1938317/100‑1 (MQ=255)
            ggcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCAt                                                                                         <  2:2124903/100‑1 (MQ=38)
            ggcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCAt                                                                                         <  1:354125/100‑1 (MQ=255)
             gcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATa                                                                                        <  1:645397/100‑1 (MQ=38)
              catgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATAg                                                                                       <  1:788369/100‑1 (MQ=38)
              catgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAg                                                                                       >  1:2175068/1‑100 (MQ=255)
              catgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAg                                                                                       >  2:2382755/1‑100 (MQ=255)
               atgAATGTTGACTTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATAGa                                                                                      <  2:1213408/100‑1 (MQ=37)
               atgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATAGa                                                                                      <  1:491072/100‑1 (MQ=38)
                   aTGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGATTCTCATAGAtttt                                                                                  <  1:1115190/100‑1 (MQ=255)
                   aTGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGAtttt                                                                                  <  2:1964742/100‑1 (MQ=255)
                     gTTGAATTTGGGAATGCGGAACAAGGTTGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTAc                                                                                <  2:1835722/100‑1 (MQ=255)
                     gTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATAGacagaga                                                                                <  1:481005/100‑7 (MQ=25)
                     gTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTAc                                                                                >  1:2134829/1‑100 (MQ=255)
                      ttGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACt                                                                               <  1:2048746/100‑1 (MQ=255)
                      ttGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACt                                                                               >  1:759345/1‑100 (MQ=255)
                      ttGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACt                                                                               <  1:1856543/100‑1 (MQ=255)
                       tGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTa                                                                              >  2:2519568/1‑100 (MQ=255)
                          aTTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATa                                                                           >  1:888071/1‑100 (MQ=255)
                           tttGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTGACGTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATa                                                                           >  2:919501/1‑99 (MQ=255)
                             tGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCa                                                                        >  1:2297625/1‑100 (MQ=255)
                                gAATGCGGAACAAGGTGGATTTGAGGCGAATGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCaaaa                                                                     <  2:752094/100‑1 (MQ=255)
                                gAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCaaaa                                                                     <  2:2160593/100‑1 (MQ=255)
                                 aaTGCGGAACAAGGTGGGTTTGAGGCGAAGATAAGCTGGTAACTTCTGTGTCTATGAGAGCCGTTAAAACGACGCTCATAGATTTTACTAATAGCaaaa                                                                     <  2:2386098/99‑1 (MQ=255)
                                  aTGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAata                                                                   >  1:116184/1‑100 (MQ=255)
                                    gCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAAtata                                                                 >  1:836810/1‑100 (MQ=255)
                                     cGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATaa                                                                >  2:1740211/1‑100 (MQ=255)
                                     cGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATaa                                                                >  2:1493627/1‑100 (MQ=255)
                                       gAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAAc                                                              >  2:2366316/1‑100 (MQ=255)
                                         aCAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCg                                                            >  1:1518230/1‑100 (MQ=255)
                                         aCAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCg                                                            >  1:2361758/1‑100 (MQ=255)
                                            aGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCTAAATATAAACCGTcc                                                         >  1:2518619/1‑100 (MQ=255)
                                             ggTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTccc                                                        >  1:1454917/1‑100 (MQ=255)
                                                ggATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCaa                                                     >  2:1770015/1‑100 (MQ=255)
                                                    ttGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCaaaaaa                                                 >  1:1782932/1‑100 (MQ=255)
                                                      gAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAg                                                <  1:2447413/99‑1 (MQ=255)
                                                        ggCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGcca                                             >  1:957126/1‑100 (MQ=255)
                                                        ggCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGcca                                             >  2:1132880/1‑100 (MQ=255)
                                                             aGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCaccaac                                        >  1:1044946/1‑100 (MQ=255)
                                                              ggTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCaccaacc                                       <  1:2441758/100‑1 (MQ=255)
                                                               gTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCaccaacca                                      >  1:566647/1‑100 (MQ=255)
                                                               gTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCaccaacca                                      >  1:711463/1‑100 (MQ=255)
                                                                tAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACcac                                     >  1:291953/1‑100 (MQ=255)
                                                                tAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACcac                                     >  2:81980/1‑100 (MQ=255)
                                                                tAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACcac                                     >  2:1311771/1‑100 (MQ=255)
                                                                 aaGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACcaca                                    <  2:227427/100‑1 (MQ=255)
                                                                   gCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACaaa                                  >  2:2561039/1‑100 (MQ=255)
                                                                   gCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACaaa                                  >  1:2075727/1‑100 (MQ=255)
                                                                   gCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACAGTCACCAAACAAGCCACCACCCACaaa                                  <  2:161525/100‑1 (MQ=255)
                                                                    cTGGTAACTTCTCTGTCTATGAGAGCCTTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAAc                                 <  2:1214388/100‑1 (MQ=255)
                                                                         aaCTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGAc                            <  1:2175301/100‑1 (MQ=255)
                                                                         aaCTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGAc                            <  2:1804486/100‑1 (MQ=255)
                                                                         aaCTTCTCTGTCTATGAGAGCCGTCAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAg                              <  2:1764700/98‑1 (MQ=255)
                                                                          aCTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACa                           <  2:1014338/100‑1 (MQ=255)
                                                                          aCTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACa                           <  1:979821/100‑1 (MQ=255)
                                                                             tctcTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAAt                        <  1:350093/100‑1 (MQ=255)
                                                                              tcttGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATa                       <  1:1965474/97‑1 (MQ=255)
                                                                               tctGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATAc                      >  1:800519/1‑100 (MQ=255)
                                                                                 tGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATaaga                    >  2:1965460/1‑97 (MQ=255)
                                                                                 tGTCTATGAGAACCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATaaga                    >  2:1967623/1‑97 (MQ=255)
                                                                                   tCTATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgag                  >  1:1033964/1‑100 (MQ=38)
                                                                                   tCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgag                  <  1:1386550/100‑1 (MQ=255)
                                                                                   tCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAACAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgag                  >  1:1097754/1‑100 (MQ=255)
                                                                                    cTATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagc                 >  1:1691144/1‑100 (MQ=38)
                                                                                    cTATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagc                 >  1:1824572/1‑100 (MQ=38)
                                                                                    cTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagc                 >  2:1296925/1‑100 (MQ=255)
                                                                                    cTAGGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagc                 <  2:551306/100‑1 (MQ=255)
                                                                                     tATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGCCCCCAAAAAAGCCACCAACCACAAACCACACAAAAACctagc                 >  2:719675/1‑99 (MQ=37)
                                                                                     tATGAGAGTCGTTTTAACGGCTCTCATAGATCTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcga                   >  2:1956655/1‑97 (MQ=38)
                                                                                     tATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcg                <  1:1917757/100‑1 (MQ=255)
                                                                                     tATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACAAACCACAAACCAGACAAATACcgagcg                >  2:1471913/1‑100 (MQ=255)
                                                                                     tATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGACACCAACCACAAACCAGACAAATACcgagcg                >  2:1735692/1‑100 (MQ=255)
                                                                                      aTGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcga               >  1:1341330/1‑100 (MQ=38)
                                                                                      aTGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcga               >  2:268571/1‑100 (MQ=38)
                                                                                      aTGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgag                  >  2:1397442/1‑97 (MQ=38)
                                                                                      aTGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGAACAATATAAACCGTCCCCAAAAAAGCCACAAACCACAAAACAAACAACTACcgag                  >  2:1470638/1‑97 (MQ=25)
                                                                                        gaggGCCGTTAAAACGACTCTGATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCgaga             <  1:2412755/100‑1 (MQ=255)
                                                                                         agagCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGaa            >  1:684287/1‑100 (MQ=255)
                                                                                           agCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATa          <  2:1504043/100‑1 (MQ=255)
                                                                                           agCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATa          <  2:571766/100‑1 (MQ=255)
                                                                                            gCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATAc         >  1:1499743/1‑100 (MQ=255)
                                                                                             ccGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACa        <  1:2307501/100‑1 (MQ=255)
                                                                                             ccGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACa        <  2:917029/100‑1 (MQ=255)
                                                                                             ccGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACa        <  1:2057142/100‑1 (MQ=255)
                                                                                               gTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACata      >  1:924072/1‑100 (MQ=255)
                                                                                               gTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACata      >  2:2239988/1‑100 (MQ=255)
                                                                                               gTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACaag      >  2:165496/1‑98 (MQ=255)
                                                                                                   aaaCGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATAcc  >  1:2539683/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
GCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACC  >  NC_000913/3171759‑3171957

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