Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,171,863 A→G 18.4% intergenic (+30/+16) 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,8630AG18.4% 246.1 / 48.8 119intergenic (+30/+16)qseC/ygiZquorum sensing sensory histidine kinase in two‑component regulatory system with QseB/inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (42/55);  new base G (18/4);  total (60/59)
Fisher's exact test for biased strand distribution p-value = 1.68e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.70e-01

Only 100 of 119 total aligned reads displayed.

TTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAA  >  NC_000913/3171765‑3171961
                                                                                                  |                                                                                                  
ttGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGAc                                                                                                   <  1:2052091/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)
               gTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTAc                                                                                    >  2:691931/1‑100 (MQ=255)
                ttGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACt                                                                                   <  2:1680983/100‑1 (MQ=255)
                ttGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACt                                                                                   >  2:143486/1‑100 (MQ=255)
                  gAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTaa                                                                                 <  2:2056904/100‑1 (MQ=255)
                  gAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTaa                                                                                 <  2:1389907/100‑1 (MQ=255)
                     tttGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAg                                                                              <  2:2357602/100‑1 (MQ=255)
                     tttGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAg                                                                              <  2:1283439/100‑1 (MQ=255)
                     tttGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAg                                                                              <  2:740959/100‑1 (MQ=255)
                      ttGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGc                                                                             <  2:1098365/100‑1 (MQ=255)
                      ttGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGc                                                                             >  2:2564901/1‑100 (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)
                            aTGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAata                                                                       >  2:1646536/1‑100 (MQ=255)
                              cttccgatctgTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAAtata                                                                     <  2:2551942/90‑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)
                                        gTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATaagatcggaag                                                           >  1:2551923/1‑91 (MQ=255)
                                        gTGGATTTGAGGCGAAGGTAAGCTGGTAACTCCTCGGTCTATGAGAGCGGGTAAAACGATTCTCATAGTTTTTCCTAATAGCAAAACATCAACCTTcccc                                                           >  1:178588/1‑100 (MQ=255)
                                         tGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCa                                                          <  2:1241079/100‑1 (MQ=255)
                                         tGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCa                                                          <  2:1578498/100‑1 (MQ=255)
                                          ggATTTGAGGGCAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGAGTCTCATAGATTTTACTAATAGCAAAATATAAACCGTccc                                                            >  1:1175050/1‑97 (MQ=255)
                                              ttGAGGCGAAGGTAATCTGGTAACTTCGCTGTCTATGATAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCaaaaaa                                                     <  2:242003/100‑1 (MQ=255)
                                               tGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAg                                                    <  2:176807/100‑1 (MQ=255)
                                               tGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAg                                                    <  1:1346261/100‑1 (MQ=255)
                                                  ggCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGcc                                                  <  1:348638/99‑1 (MQ=255)
                                                  ggCGAAGGTAAGCGGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGcca                                                 <  1:54908/100‑1 (MQ=255)
                                                   gCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGccac                                                <  2:1994406/100‑1 (MQ=255)
                                                    cGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGcca                                                 >  2:1881320/1‑98 (MQ=255)
                                                       aGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCaccaac                                            >  2:2126918/1‑100 (MQ=255)
                                                        ggTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCaccaacc                                           <  1:688695/100‑1 (MQ=255)
                                                        ggTAAGCTGGTAACTTCGCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCaccaacc                                           <  2:2327952/100‑1 (MQ=255)
                                                          tAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCaccaacc                                           <  2:2122967/98‑1 (MQ=255)
                                                          tAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCAACCAACcac                                         >  1:2401589/1‑100 (MQ=255)
                                                             gCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACaaa                                      >  2:2036885/1‑100 (MQ=255)
                                                                ggTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCAACCAACCACAAAcaa                                   >  1:2592202/1‑98 (MQ=255)
                                                                  tAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGa                                 >  1:945469/1‑100 (MQ=255)
                                                                  tAACTTCTCTGTCTATGAGAGACGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACAGTAACCAAAAAAGCCACCAACCAAAAACCa                                   >  1:1609837/1‑98 (MQ=255)
                                                                   aaCTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGAc                                <  1:651854/100‑1 (MQ=255)
                                                                     cTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACaa                              >  2:266342/1‑100 (MQ=255)
                                                                     cTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACaa                              >  1:1972305/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)
                                                                             tCTATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACAAGACAAATACcga                       >  1:1407464/1‑99 (MQ=37)
                                                                              cTATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATAcc                         >  1:1620318/1‑96 (MQ=38)
                                                                              cTATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagc                     >  2:1309745/1‑100 (MQ=38)
                                                                              cTATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagc                     >  2:12123/1‑100 (MQ=38)
                                                                              cTATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTACCCAAAAAAGCCACCAACCACAAACCAGACAAAAACcgagc                     >  2:1358713/1‑100 (MQ=37)
                                                                              cTATGAGAGTCGGTTTAACGGCTCTTATAGATTTTACTAATATCAAAATATAAACCGTCGCCAAAACAGACACCAACGACAAACCAGTCAAATACcgagc                     >  1:1909313/1‑100 (MQ=255)
                                                                              cTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagc                     >  1:1556652/1‑100 (MQ=255)
                                                                              cTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagc                     >  2:878749/1‑100 (MQ=255)
                                                                                aTGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTCATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcga                   >  1:327997/1‑100 (MQ=38)
                                                                                aTGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcga                   >  2:1009432/1‑100 (MQ=38)
                                                                                aTGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcga                   >  1:947501/1‑100 (MQ=38)
                                                                                aTGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcga                   >  2:668225/1‑100 (MQ=38)
                                                                                aTGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcga                   >  2:2304010/1‑100 (MQ=38)
                                                                                aTGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCAACAAACCACAAACAAGACAAATACcgaaaga                   >  1:1641459/1‑96 (MQ=25)
                                                                                aTGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcga                   <  2:614402/100‑1 (MQ=255)
                                                                                aTGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcga                   <  2:1874940/100‑1 (MQ=255)
                                                                                 tGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcgag                  <  2:1372288/100‑1 (MQ=255)
                                                                                 tGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcgag                  <  1:39818/100‑1 (MQ=255)
                                                                                    gagCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAAt               <  2:2194899/100‑1 (MQ=255)
                                                                                    gagCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCAAACGAGGAt               >  1:192796/1‑100 (MQ=255)
                                                                                      gCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATAc             <  2:1035026/100‑1 (MQ=255)
                                                                                        cGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACat           <  1:160509/100‑1 (MQ=255)
                                                                                        cGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACat           >  2:843910/1‑100 (MQ=255)
                                                                                         gTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACata          <  1:2112205/100‑1 (MQ=255)
                                                                                         gTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACata          >  1:939349/1‑100 (MQ=255)
                                                                                         gTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACata          >  2:908191/1‑100 (MQ=255)
                                                                                            aaaaCGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATAc       >  2:1327516/1‑100 (MQ=255)
                                                                                             aaaCGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATAcc      <  2:1114880/100‑1 (MQ=255)
                                                                                              aaCGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATAccc     <  2:2271462/100‑1 (MQ=255)
                                                                                              aaCGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATAccc     <  1:11129/100‑1 (MQ=255)
                                                                                                 gACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCaaa  <  1:1646527/100‑1 (MQ=255)
                                                                                                  |                                                                                                  
TTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAA  >  NC_000913/3171765‑3171961

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: