Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,171,857 A→T 15.5% 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,8570AT15.5% 572.8 / 86.6 219intergenic (+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 (86/99);  new base T (15/19);  total (101/118)
Fisher's exact test for biased strand distribution p-value = 8.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.16e-01

Only 100 of 234 total aligned reads displayed.

CGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATAC  >  NC_000913/3171758‑3171956
                                                                                                   |                                                                                                   
cGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTa                                                                                                     >  2:1643437/1‑100 (MQ=255)
  ccAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTaaa                                                                                                   >  2:1236518/1‑100 (MQ=255)
     aaTTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACg                                                                                                <  1:1250061/100‑1 (MQ=255)
      aTTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGa                                                                                               <  1:942282/100‑1 (MQ=255)
        tGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGAct                                                                                             >  1:1926708/1‑100 (MQ=255)
          catggcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGActct                                                                                           <  2:2142551/100‑1 (MQ=255)
             ggcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCAt                                                                                        <  1:1189426/100‑1 (MQ=38)
             ggcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCAt                                                                                        <  1:252828/100‑1 (MQ=255)
              gcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATa                                                                                       <  1:3679355/100‑1 (MQ=38)
              gcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATa                                                                                       <  1:3114103/100‑1 (MQ=38)
              gcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATa                                                                                       <  2:3097072/100‑1 (MQ=38)
               catgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAg                                                                                      >  1:467062/1‑100 (MQ=255)
               catgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAg                                                                                      >  1:1070402/1‑100 (MQ=255)
               catgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAg                                                                                      >  1:725041/1‑100 (MQ=255)
                atgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGa                                                                                     <  1:647401/100‑1 (MQ=255)
                 tgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATAGac                                                                                    <  1:1638566/100‑2 (MQ=37)
                  gagggTTGAAGTTGGGAATGCGGAACAAGGTGGATTTGGGGCGAAGGTAAGCGGGTAACTTCTGTGTCTATGAGAGCCGTTAAAACGACTCTCATAGAtt                                                                                   <  2:3124753/96‑1 (MQ=255)
                  gAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGAtt                                                                                   <  2:1432704/100‑1 (MQ=255)
                  gAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGAtt                                                                                   <  1:2640680/100‑1 (MQ=255)
                   aaTGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAATTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATAGacag                                                                                  <  1:338589/100‑4 (MQ=37)
                   aaTGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATAGacag                                                                                  <  2:1057138/100‑4 (MQ=37)
                   aaTGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGAtt                                                                                   <  2:698130/99‑1 (MQ=255)
                    aTGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGAtttt                                                                                 <  1:3993003/100‑1 (MQ=255)
                      gTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATAGacagat                                                                                >  2:1200812/1‑94 (MQ=37)
                      gTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTAc                                                                               <  2:1111635/100‑1 (MQ=255)
                      gTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTAc                                                                               >  2:1579091/1‑100 (MQ=255)
                      gTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTAc                                                                               <  2:1635347/100‑1 (MQ=255)
                      gTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTAc                                                                               <  1:3986155/100‑1 (MQ=255)
                      gTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTAc                                                                               >  1:3227341/1‑100 (MQ=255)
                      gTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTAc                                                                               >  1:379892/1‑100 (MQ=255)
                       ttGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACt                                                                              <  2:1543951/100‑1 (MQ=255)
                       ttGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACt                                                                              <  1:3511515/100‑1 (MQ=255)
                        tGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTa                                                                             >  1:893021/1‑100 (MQ=255)
                        tGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTa                                                                             >  2:2544589/1‑100 (MQ=255)
                         gAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTaa                                                                            >  1:55474/1‑100 (MQ=255)
                            tttGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAg                                                                         >  1:1058802/1‑100 (MQ=255)
                            tttGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAg                                                                         <  1:2203170/100‑1 (MQ=255)
                            tttGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAg                                                                         <  1:2372535/100‑1 (MQ=255)
                            nttGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTATGTCTATGAGCTCCGTTAAAAATACTCTCATAGATTTTACTAATAg                                                                         >  1:1010708/2‑100 (MQ=255)
                                ggAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCaaa                                                                     >  1:2253985/1‑100 (MQ=255)
                                 gAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCaaaa                                                                    <  2:2242263/100‑1 (MQ=255)
                                    tGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAatat                                                                 <  2:641511/100‑1 (MQ=255)
                                     gCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAAtata                                                                >  2:230831/1‑100 (MQ=255)
                                     gCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAAtata                                                                >  1:977676/1‑100 (MQ=255)
                                     gCGGAACAAGGTAGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAAtata                                                                >  1:2758693/1‑100 (MQ=255)
                                      cGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATaa                                                               >  1:3282350/1‑100 (MQ=255)
                                      cGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATaa                                                               >  1:2907295/1‑100 (MQ=255)
                                      cGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACAAATAGCAAAATATaa                                                               >  1:4133775/1‑100 (MQ=255)
                                        gAACAAGGGGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAAc                                                             <  1:1839503/100‑1 (MQ=255)
                                            aaGGTGGATTTGAGGCGAAGATAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTc                                                         >  2:2561478/1‑100 (MQ=255)
                                             aGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAACAGCAAAATATAAACCGTcc                                                        <  2:2877904/100‑1 (MQ=255)
                                              ggTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTccc                                                       >  1:3931642/1‑100 (MQ=255)
                                               gTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTccc                                                       >  2:754940/1‑99 (MQ=255)
                                                 ggATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACCCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCaa                                                    >  1:1160923/1‑100 (MQ=255)
                                                  gATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCaaa                                                   >  1:573225/1‑100 (MQ=255)
                                                  gATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCaaa                                                   <  2:3270458/100‑1 (MQ=255)
                                                  gATTTGAGGCGAAGGAAAGCCGGTAAGTTCTCTGTCTATGAGAGCCGTTACAGCGACGCTCATAGTTTTTACGAATTGCTAACTATAAACTGTGCCCaaa                                                   >  2:144606/1‑100 (MQ=255)
                                                      tGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAg                                               <  1:2263175/100‑1 (MQ=255)
                                                        aGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGcc                                             <  1:1456947/100‑1 (MQ=255)
                                                        aGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGcc                                             <  2:1997884/100‑1 (MQ=255)
                                                         ggCGAAGGTAAGCTGGTAACTTCTCTGTCTATGATAGCCGTTAAAACGACACTCATAGATTTTACTAATAGCAAAATATAAACCGTCACCAAAAAAGcca                                            <  2:3152951/100‑1 (MQ=255)
                                                         ggCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAGTATAAACCGTCCCCAAAAAAg                                               >  2:4051257/1‑97 (MQ=255)
                                                          gcaaagctaagctGGTACATTACATGTTTATGCGAGCCGTTAAAACCACTCTAATAGATTTTACTAATATTAAAATATAAACCTTTCCCATATAAgtc                                             <  2:2440953/95‑3 (MQ=255)
                                                            gAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGccacca                                         >  1:2861704/1‑100 (MQ=255)
                                                              aGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCaccaac                                       >  1:247120/1‑100 (MQ=255)
                                                                gTAAGCTTGTAACTTCTCTGTCTATGAGAGACGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCaccaac                                       <  2:3862573/98‑1 (MQ=255)
                                                                  aaGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACcaca                                   <  2:770030/100‑1 (MQ=255)
                                                                  aaGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACcaca                                   <  2:3000486/100‑1 (MQ=255)
                                                                  aaGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAAAAAAAAAAAAAAAAGCCACCAACcaca                                   <  2:1499211/100‑1 (MQ=255)
                                                                   aGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACaa                                  <  1:1532515/100‑1 (MQ=255)
                                                                    gCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACaaa                                 >  2:3988259/1‑100 (MQ=255)
                                                                    gCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACaaa                                 <  2:2729163/100‑1 (MQ=255)
                                                                     cTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAAc                                <  1:4094264/100‑1 (MQ=255)
                                                                     cTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAAc                                <  1:1643441/100‑1 (MQ=255)
                                                                       ggTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCa                              >  2:866591/1‑100 (MQ=255)
                                                                       ggTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCa                              <  2:2108091/100‑1 (MQ=255)
                                                                       ggTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCa                              <  2:1565005/100‑1 (MQ=255)
                                                                        gTAACTTCTCTATATATGAGAGCCGTTAAAACGACACTAAAAAATTATAATAATAAAAAAATAAAAAAAAAAAAAAAAAAAGCCACCAACCACAAAcc                               <  2:3773446/98‑1 (MQ=255)
                                                                          aaCTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGAc                           <  1:2688852/100‑1 (MQ=255)
                                                                            cTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACaa                         >  2:1224543/1‑100 (MQ=255)
                                                                            cTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACaa                         >  2:3481462/1‑100 (MQ=255)
                                                                               ctctGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATa                      >  2:1358788/1‑100 (MQ=255)
                                                                                 ctGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATAcc                    <  1:148728/100‑1 (MQ=255)
                                                                                 ctGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATAcc                    >  1:3848371/1‑100 (MQ=255)
                                                                                 ctGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATAcc                    >  1:373828/1‑100 (MQ=255)
                                                                                  tGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcg                   >  1:1793687/1‑100 (MQ=255)
                                                                                     cTATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagc                >  1:2448893/1‑100 (MQ=38)
                                                                                     cTATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagc                >  1:2112490/1‑100 (MQ=38)
                                                                                     cTATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagc                >  1:3168314/1‑100 (MQ=38)
                                                                                     cTATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACAAATAGCAAAATATAAACCGTCCTCAAAAAAGACACTAACCACCAACCAGACAAATCCcgaac                >  2:3390623/1‑98 (MQ=21)
                                                                                     cTATGAGAGTCGTTTTAACGGCTCTCATAGATTGTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCAACACACCACAAAACAGACAAATAc                     >  2:2775196/1‑95 (MQ=25)
                                                                                      tATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcga                  >  2:3612535/1‑97 (MQ=38)
                                                                                      tATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcg               >  1:593188/1‑100 (MQ=255)
                                                                                       aTGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcga              >  1:1734854/1‑100 (MQ=38)
                                                                                         gagaGCCGTTAAAACGACTCTCATAGATTTCACTGACACCAAAAGATAAACCGTCCCAAAGAAAGCCACCAACCACTAACCAGACAAATACCGAGCgaga            <  2:1866198/100‑1 (MQ=255)
                                                                                          agagCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGaa           >  1:1867560/1‑100 (MQ=255)
                                                                                             gCCGTTAAAACGCATATCATAGATTTTAATAATTGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACTTTCCAGACAAATACCGCGCGAGAATa         <  2:3557039/99‑1 (MQ=255)
                                                                                             gCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACAAAAAAAAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATAc        <  2:4149176/100‑1 (MQ=255)
                                                                                                gTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACata     <  2:2390491/100‑1 (MQ=255)
                                                                                                   aaaaCGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATAc  <  2:3077653/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 < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: