Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,171,852 C→T 17.2% 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,8520CT17.2% 424.0 / 69.3 186intergenic (+19/+27)qseC/ygiZquorum sensing sensory histidine kinase in two‑component regulatory system with QseB/inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (78/74);  new base T (16/16);  total (96/90)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.87e-01

Only 100 of 187 total aligned reads displayed.

GAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATAC  >  NC_000913/3171754‑3171950
                                                                                                  |                                                                                                  
     gCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTaa                                                                                              >  1:1210365/1‑100 (MQ=255)
     gCCAAATTGCATGGCATGAATGGTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGtt                                                                                                <  1:574146/98‑1 (MQ=255)
     gCCAAATGGCATGGTATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGGTAAGGTAAGCTGGTAACTCCTCTGTCTATGAGAGCCGtt                                                                                                <  1:2394652/98‑1 (MQ=255)
      ccAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTaaa                                                                                             >  1:2776774/1‑100 (MQ=255)
      ccAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTaaa                                                                                             >  2:978850/1‑100 (MQ=255)
        aaaTTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAc                                                                                           <  1:778365/100‑1 (MQ=255)
        aaaTTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAc                                                                                           <  2:3401814/100‑1 (MQ=255)
             gcatggcatgAATGTTGAATTTGGGAATGCGTACCACGGAGGATTTGAGGCGATGGTTAGCTGGTAACTTCTCTGTATCTGAGAGCGGTTaaa                                                                                             >  1:1848840/1‑93 (MQ=255)
             gcatggcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGActc                                                                                      <  2:1832951/100‑1 (MQ=255)
             gcatggcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGActc                                                                                      >  2:2766396/1‑100 (MQ=255)
              catggcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGActct                                                                                     <  2:3653598/100‑1 (MQ=255)
               atggcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACtctc                                                                                    <  1:1722086/100‑1 (MQ=255)
                 ggcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCan                                                                                  <  2:3552696/100‑2 (MQ=38)
                   catgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATAg                                                                                <  2:1080407/100‑1 (MQ=38)
                    atgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATAGa                                                                               <  2:665599/100‑1 (MQ=38)
                    atgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATAGa                                                                               <  2:2894994/100‑1 (MQ=38)
                       aaTGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGAt                                                                              >  1:3749567/1‑98 (MQ=255)
                        aTGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATAGacaga                                                                           <  2:59792/100‑5 (MQ=37)
                         tGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTa                                                                          >  2:3750235/1‑100 (MQ=255)
                          gTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATAGacagaga                                                                         <  2:29859/100‑7 (MQ=25)
                          gTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTAc                                                                         <  1:1780966/100‑1 (MQ=255)
                          gTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTAc                                                                         >  1:3337879/1‑100 (MQ=255)
                           ttGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACt                                                                        >  1:1256364/1‑100 (MQ=255)
                           ttGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACt                                                                        <  2:2805705/100‑1 (MQ=255)
                             gAATTTGGGAATGCCGAACAAGATGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACt                                                                        <  2:561227/98‑1 (MQ=255)
                              aaTTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAAt                                                                     >  1:2764643/1‑100 (MQ=255)
                              aaTTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAAt                                                                     >  1:2283294/1‑100 (MQ=255)
                               aTTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATa                                                                    <  1:238743/100‑1 (MQ=255)
                               aTTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATa                                                                    >  2:3093565/1‑100 (MQ=255)
                               aTTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATa                                                                    >  1:1505099/1‑100 (MQ=255)
                               aTTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATa                                                                    >  1:431242/1‑100 (MQ=255)
                                tttGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAg                                                                   <  2:3174964/100‑1 (MQ=255)
                                tttGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAg                                                                   <  2:897493/100‑1 (MQ=255)
                                  tGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTGTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCa                                                                 >  1:1590438/1‑100 (MQ=255)
                                   gggAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCaa                                                                >  1:3551081/1‑100 (MQ=255)
                                       aTGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAata                                                            <  2:2589473/100‑1 (MQ=255)
                                       aTGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAat                                                             >  2:3234381/1‑99 (MQ=255)
                                       aTGCGGAACAAGGCGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAata                                                            <  1:508555/100‑1 (MQ=255)
                                        tGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAataa                                                           >  2:2188178/1‑99 (MQ=255)
                                           ggAACAAGTTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTATATCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATaaa                                                        <  1:1890339/100‑1 (MQ=255)
                                           ggAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATaaa                                                        <  1:1109437/100‑1 (MQ=255)
                                           ggAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATaaa                                                        <  2:104091/100‑1 (MQ=255)
                                           ggAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATaaa                                                        <  1:2428842/100‑1 (MQ=255)
                                            naaCAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAAc                                                       >  1:2804196/2‑100 (MQ=255)
                                            gAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAAc                                                       >  2:397133/1‑100 (MQ=255)
                                                 aGGTGGGTTTGAGGCGAAGGTAAGCGGGTAACTTCTCTGTCTATGAGAGCCGTTAAACCGACTCGCTTAGATTTTACTAATAGCAAAATATATACCGTcc                                                  <  1:1876406/100‑1 (MQ=255)
                                                 aGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTcn                                                  <  1:71777/100‑2 (MQ=255)
                                                  ggTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTccc                                                 >  2:1698261/1‑100 (MQ=255)
                                                          tGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAg                                         <  1:606347/100‑1 (MQ=255)
                                                            aGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAg                                         >  1:3695801/1‑98 (MQ=255)
                                                             ggCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGcca                                      >  2:2664445/1‑100 (MQ=255)
                                                             ggCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGACGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGcca                                      >  1:3015696/1‑100 (MQ=255)
                                                              gCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGccac                                     <  2:2001820/100‑1 (MQ=255)
                                                                   ggTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCaccaacc                                <  2:2018525/100‑1 (MQ=255)
                                                                   ggTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCaccaacc                                <  2:2839486/100‑1 (MQ=255)
                                                                   ggTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCaccaacc                                <  2:2701229/100‑1 (MQ=255)
                                                                     tAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAGTATAAACCGTCCCCAAAAAAGCCACCAACcac                              <  1:3468423/100‑1 (MQ=255)
                                                                     tAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACcac                              <  2:3345193/100‑1 (MQ=255)
                                                                      aaGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACcaca                             >  1:1151411/1‑100 (MQ=255)
                                                                         cTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAAc                          <  2:3436187/100‑1 (MQ=255)
                                                                          tGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAAcc                         >  1:753775/1‑100 (MQ=255)
                                                                          tGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCAAAAAAAAAAAAGCCACCAACCACAAAcc                         <  1:1449484/100‑1 (MQ=255)
                                                                          tGGTAACTTCGCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCACCAAAAAAGCCACCAACCACAAAcc                         <  1:3425126/100‑1 (MQ=255)
                                                                           ggTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCa                        >  1:2701522/1‑100 (MQ=255)
                                                                           ggTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCa                        >  2:2463284/1‑100 (MQ=255)
                                                                           ggTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCa                        >  2:755597/1‑100 (MQ=255)
                                                                            gTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAg                       <  1:2573870/100‑1 (MQ=255)
                                                                            gTAACTTCTCTGTCTATCAGACCCGTTAAAACGACACTCACAGAGGTTACGAATAGCAACATCTAAACGGTCCCCAAAAAAGCCACCACCCACAAACCAg                       <  1:1373250/100‑1 (MQ=255)
                                                                              aaCTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGAc                     <  2:203590/100‑1 (MQ=255)
                                                                               aCTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACa                    >  2:2954801/1‑100 (MQ=255)
                                                                               aCTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACa                    <  1:1703270/100‑1 (MQ=255)
                                                                                    tctGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATAc               <  1:397130/100‑1 (MQ=255)
                                                                                    aaaatCTATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATGAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACaaa                  >  1:1834877/5‑97 (MQ=25)
                                                                                    aaaatCTATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATAc               >  2:1135770/5‑100 (MQ=25)
                                                                                     ctGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATAcc              <  2:2383839/100‑1 (MQ=255)
                                                                                      tGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATAcc              >  1:1975194/1‑99 (MQ=255)
                                                                                      aatCTATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCTCCAAAAAAGCCACCAACCACAAACCAGACAAATACcg             >  1:2810430/3‑100 (MQ=37)
                                                                                      aatCTATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATa                >  1:3180334/3‑97 (MQ=37)
                                                                                      aatCTATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcg             >  2:2071797/3‑100 (MQ=37)
                                                                                       atCTATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcga            >  2:500266/2‑100 (MQ=37)
                                                                                        tCTATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgag           >  2:3264940/1‑100 (MQ=38)
                                                                                        tCTATGAGAGCCGTTAAAACGCCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgag           <  2:2515727/100‑1 (MQ=255)
                                                                                        tCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgag           >  2:312632/1‑100 (MQ=255)
                                                                                        ncTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgag           >  1:1568185/2‑100 (MQ=255)
                                                                                         cTATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagc          >  2:1308879/1‑100 (MQ=38)
                                                                                         cTATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgag           >  1:2939356/1‑99 (MQ=38)
                                                                                         cTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagc          <  2:3419320/100‑1 (MQ=255)
                                                                                          tATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcg         >  1:1483743/1‑100 (MQ=38)
                                                                                          tATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcg         >  2:1559842/1‑100 (MQ=255)
                                                                                          tATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcg         >  2:2151701/1‑100 (MQ=255)
                                                                                          tATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagag         <  2:2625375/100‑3 (MQ=255)
                                                                                           aTGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCGAAATATAAACCGTCCCCAAAAAAGCCACCAACAACAAACCAGACAAATACcgagcga        >  1:2507804/1‑100 (MQ=37)
                                                                                           aTGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcg         >  2:2606995/1‑99 (MQ=38)
                                                                                           aTGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCAACCAACAACAAACAAGACAAATACagag           >  1:2828273/1‑97 (MQ=25)
                                                                                           aTGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcga        <  1:200695/100‑1 (MQ=255)
                                                                                            tGAGAGCCGTTCAAACGAATCTCATAGATTTTACTAATTGTAAAATATAAACTGTCACCGAAAAAGCTATCAACGACAAACCAGACAAATATcgagcgag       >  1:1062952/1‑100 (MQ=255)
                                                                                            tGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCTTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcgag       <  1:2148485/100‑1 (MQ=255)
                                                                                             gagaGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCgaga      >  1:2909430/1‑100 (MQ=255)
                                                                                             gagaGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCgaga      >  1:955111/1‑100 (MQ=255)
                                                                                                 gCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATAc  >  1:3423590/1‑100 (MQ=255)
                                                                                                  |                                                                                                  
GAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATAC  >  NC_000913/3171754‑3171950

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: