Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,171,852 C→T 18.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,8520CT18.2% 264.8 / 41.8 121intergenic (+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 (43/54);  new base T (18/4);  total (63/58)
Fisher's exact test for biased strand distribution p-value = 1.81e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 122 total aligned reads displayed.

CGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATAC  >  NC_000913/3171753‑3171950
                                                                                                   |                                                                                                  
cGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGc                                                                                                    >  2:1531434/1‑100 (MQ=255)
  aaTCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCg                                                                                                  >  1:1756596/1‑100 (MQ=255)
   aTCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGt                                                                                                 <  1:2170955/100‑1 (MQ=255)
     cGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGt                                                                                                 <  1:996466/98‑1 (MQ=255)
      gCCAAATTGCATGGCATGAATGTTGAATTTGTGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTaa                                                                                              <  1:1631955/100‑1 (MQ=255)
      gCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTaa                                                                                              >  2:1280274/1‑100 (MQ=255)
      gCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTaa                                                                                              >  1:2549812/1‑100 (MQ=255)
      gCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTaa                                                                                              >  1:2259071/1‑100 (MQ=255)
      gCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTaa                                                                                              >  2:541045/1‑100 (MQ=255)
       ccAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTaaa                                                                                             >  1:989196/1‑100 (MQ=255)
        cAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTaaaa                                                                                            <  1:1039625/100‑1 (MQ=255)
         aaaTTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAc                                                                                           <  1:2256324/100‑1 (MQ=255)
             tGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGAct                                                                                       <  2:2270239/100‑1 (MQ=255)
             tGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGAGAAGGTAAGCTGGTAACTTCCCTGCCTATGACAGCCGTTAAACCGAc                                                                                        <  1:1677339/99‑1 (MQ=255)
              gcatggcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGActc                                                                                      <  2:1298475/100‑1 (MQ=255)
               natggcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCAAAGGTAAGCTGGGAACTTCACTGAATATGAGAGACGATAAAACGActct                                                                                     >  1:1931870/2‑100 (MQ=255)
                atggcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACtctc                                                                                    <  1:2098654/100‑1 (MQ=255)
                 tggcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCa                                                                                   <  1:601044/100‑1 (MQ=255)
                   gcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATa                                                                                 <  2:1152649/100‑1 (MQ=255)
                    catgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATAg                                                                                <  1:1556929/100‑1 (MQ=38)
                     atgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATAGa                                                                               <  1:64107/100‑1 (MQ=38)
                       gAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATAGaca                                                                             <  1:738466/100‑3 (MQ=37)
                       gAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGAtt                                                                             <  2:889770/100‑1 (MQ=255)
                       gAATGTTCAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGAtt                                                                             <  2:2559750/100‑1 (MQ=255)
                         aTGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGAtttt                                                                           <  1:1828661/100‑1 (MQ=255)
                          tGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATAGacagag                                                                          >  1:2299739/1‑95 (MQ=25)
                           gTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTAc                                                                         >  2:691931/1‑100 (MQ=255)
                            ttGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGTTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTAc                                                                         <  1:1436935/99‑1 (MQ=255)
                            ttGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACt                                                                        <  2:1680983/100‑1 (MQ=255)
                              gAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTaa                                                                      >  2:1714781/1‑100 (MQ=255)
                              gAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTaa                                                                      <  2:1389907/100‑1 (MQ=255)
                              gAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTaa                                                                      >  1:323057/1‑100 (MQ=255)
                              gAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTaa                                                                      <  2:2056904/100‑1 (MQ=255)
                                aTTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATa                                                                    >  2:761951/1‑100 (MQ=255)
                                 tttGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAg                                                                   <  2:1283439/100‑1 (MQ=255)
                                 tttGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAg                                                                   <  2:2357602/100‑1 (MQ=255)
                                  ttGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGc                                                                  >  2:2564901/1‑100 (MQ=255)
                                  ttGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGc                                                                  <  2:1098365/100‑1 (MQ=255)
                                  ttGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGc                                                                  <  1:1369287/100‑1 (MQ=255)
                                   tGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCa                                                                 >  2:11129/1‑100 (MQ=255)
                                     ggAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCaaa                                                               >  1:1578779/1‑100 (MQ=255)
                                     ggAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCaaa                                                               >  2:160510/1‑100 (MQ=255)
                                     ggAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTGCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCaaa                                                               >  1:902251/1‑100 (MQ=255)
                                        aTGCGGAACAAGGTGGATTTGAGGCGAATGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAata                                                            <  1:2031082/100‑1 (MQ=255)
                                            gtaaCAAGGTGTATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATaaa                                                        <  1:1224029/98‑1 (MQ=255)
                                                cAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGt                                                    >  2:651857/1‑100 (MQ=255)
                                                   ggTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTccc                                                 >  2:1403827/1‑100 (MQ=255)
                                                   ggTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTccc                                                 <  1:617367/100‑1 (MQ=255)
                                                   ggTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATCAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTccc                                                 <  1:1652560/100‑1 (MQ=255)
                                                    gTGGATTTGAGGCGAAGGTAAGCTGGTAACTCCTCGGTCTATGAGAGCGGGTAAAACGATTCTCATAGTTTTTCCTAATAGCAAAACATCAACCTTcccc                                                >  1:178588/1‑100 (MQ=255)
                                                     tGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCa                                               <  2:1578498/100‑1 (MQ=255)
                                                     tGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCa                                               <  2:1241079/100‑1 (MQ=255)
                                                      ggATTTGAGGGCAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGAGTCTCATAGATTTTACTAATAGCAAAATATAAACCGTccc                                                 >  1:1175050/1‑97 (MQ=255)
                                                      ggATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCaa                                              >  2:2403608/1‑100 (MQ=255)
                                                          ttGAGGCGAAGGTAATCTGGTAACTTCGCTGTCTATGATAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCaaaaaa                                          <  2:242003/100‑1 (MQ=255)
                                                           tGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAg                                         <  1:1346261/100‑1 (MQ=255)
                                                           tGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAg                                         <  2:176807/100‑1 (MQ=255)
                                                           tGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAg                                         >  2:1499004/1‑100 (MQ=255)
                                                              ggCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGcc                                       <  1:348638/99‑1 (MQ=255)
                                                               gCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGccac                                     <  2:1994406/100‑1 (MQ=255)
                                                                   aGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCaccaac                                 >  2:2126918/1‑100 (MQ=255)
                                                                    ggTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCaccaacc                                <  1:688695/100‑1 (MQ=255)
                                                                      tAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCAACCAACcac                              >  1:2401589/1‑100 (MQ=255)
                                                                         gCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACaaa                           >  2:2036885/1‑100 (MQ=255)
                                                                            ggTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCAACCAACCACAAAcaa                        >  1:2592202/1‑98 (MQ=255)
                                                                             gTAACTTCTCTGTCTATGAAAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAAc                          <  1:923895/97‑1 (MQ=255)
                                                                              tAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGa                      >  1:945469/1‑100 (MQ=255)
                                                                              tAACTTCTCTGTCTATGAGAGACGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACAGTAACCAAAAAAGCCACCAACCAAAAACCa                        >  1:1609837/1‑98 (MQ=255)
                                                                                 cTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACaa                   >  2:266342/1‑100 (MQ=255)
                                                                                 cTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACaa                   >  1:1972305/1‑100 (MQ=255)
                                                                                 cTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACaa                   >  2:120366/1‑100 (MQ=255)
                                                                                  agtaaaatCTATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACaaa                  >  2:2375839/8‑100 (MQ=25)
                                                                                    ctctGTCTACGAGAGCCGTTAAAACGACTCTCATAAATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAGCCACAAACCAGACAAATa                <  1:1801303/100‑1 (MQ=255)
                                                                                     tctGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATAc               >  1:2501950/1‑100 (MQ=255)
                                                                                      ctGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATAc               >  2:1202231/1‑99 (MQ=255)
                                                                                       aatCTATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcg             >  2:1968875/3‑100 (MQ=37)
                                                                                       aatCTATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcg             >  2:2355555/3‑100 (MQ=37)
                                                                                        gTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcga            <  2:945472/100‑1 (MQ=255)
                                                                                        atCTATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATTAACCGTACCCAAAAAAGCCACCAACCACAAAACAGACAAATACcga            >  1:2492622/2‑100 (MQ=25)
                                                                                         tCTATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgag           <  2:1641468/100‑1 (MQ=38)
                                                                                          cTATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagc          >  2:1309745/1‑100 (MQ=38)
                                                                                          cTATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTACCCAAAAAAGCCACCAACCACAAACCAGACAAAAACcgagc          >  2:1358713/1‑100 (MQ=37)
                                                                                          cTATGAGAGTCGGTTTAACGGCTCTTATAGATTTTACTAATATCAAAATATAAACCGTCGCCAAAACAGACACCAACGACAAACCAGTCAAATACcgagc          >  1:1909313/1‑100 (MQ=255)
                                                                                          cTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagc          >  2:878749/1‑100 (MQ=255)
                                                                                          cTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagc          >  1:1556652/1‑100 (MQ=255)
                                                                                           tATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcg         >  1:614399/1‑100 (MQ=255)
                                                                                            aTGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTCATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcga        >  1:327997/1‑100 (MQ=38)
                                                                                            aTGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcga        >  2:2304010/1‑100 (MQ=38)
                                                                                            aTGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcga        >  2:668225/1‑100 (MQ=38)
                                                                                            aTGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcga        >  2:1009432/1‑100 (MQ=38)
                                                                                            aTGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCAACAAACCACAAACAAGACAAATACcgaaaga        >  1:1641459/1‑96 (MQ=25)
                                                                                            aTGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcga        <  2:1874940/100‑1 (MQ=255)
                                                                                            aTGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcga        <  2:614402/100‑1 (MQ=255)
                                                                                             tGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcgag       >  2:33888/1‑100 (MQ=255)
                                                                                             tGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcgag       <  1:1948859/100‑1 (MQ=255)
                                                                                             tGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcgag       <  1:39818/100‑1 (MQ=255)
                                                                                             tGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcgag       >  2:2019749/1‑100 (MQ=255)
                                                                                             tGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcgag       <  2:1372288/100‑1 (MQ=255)
                                                                                                gagCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAAt    <  2:2194899/100‑1 (MQ=255)
                                                                                                  gCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATAc  <  2:1035026/100‑1 (MQ=255)
                                                                                                   |                                                                                                  
CGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATAC  >  NC_000913/3171753‑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: