Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,171,863 A→G 16.2% 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,8630AG16.2% 397.6 / 63.0 186intergenic (+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 (91/62);  new base G (16/14);  total (107/79)
Fisher's exact test for biased strand distribution p-value = 5.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.30e-01

Only 100 of 188 total aligned reads displayed.

TTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAA  >  NC_000913/3171765‑3171962
                                                                                                  |                                                                                                   
  gcatggnatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGActc                                                                                                  >  2:3094079/1‑100 (MQ=255)
   catggcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGActct                                                                                                 <  1:1543929/100‑1 (MQ=255)
   catggcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGActct                                                                                                 >  2:2369164/1‑100 (MQ=255)
   catggcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGActc                                                                                                  <  2:2066830/99‑1 (MQ=255)
    atggcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCtctc                                                                                                <  1:2863681/100‑1 (MQ=38)
    atggcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACtctc                                                                                                <  2:1969325/100‑1 (MQ=255)
       gcatgAATGTTGAATTTGGGGATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATa                                                                                             >  2:2391252/1‑100 (MQ=255)
       gcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTACTCTATCTATGAGAGCCGTTAAAACGACTCTCATa                                                                                             >  1:125495/1‑100 (MQ=255)
         atgAATGTTGAATTTGGGAATGCGGAAAAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATAg                                                                                            <  2:1060268/99‑1 (MQ=38)
          tgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATAGac                                                                                          <  1:1952507/100‑2 (MQ=37)
          tgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGAt                                                                                          >  2:3238982/1‑100 (MQ=255)
           gAATGTTGAATTTGGGAATAAAGAACAAGATGGATTAGATGCGAAGGTAAGATGGTAACTTCTCTGTCTATGAGAGCAGTAAAAACGACTCCCATAg                                                                                            >  2:3669066/1‑97 (MQ=255)
            aaTGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATAGacag                                                                                        <  1:1470803/100‑4 (MQ=37)
             aTGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGAtttt                                                                                       <  2:1415883/100‑1 (MQ=255)
               gTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATAGacagaga                                                                                     <  1:2910794/100‑7 (MQ=25)
               gTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTAc                                                                                     <  1:2230561/100‑1 (MQ=255)
               gTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTAc                                                                                     <  2:2627317/100‑1 (MQ=255)
                ttGAATTTGGGTATGTGGAATAAGGTGGATTTGAGGAGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTAc                                                                                     <  1:1821255/99‑1 (MQ=255)
                 tGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTa                                                                                   >  1:1724712/1‑100 (MQ=255)
                  gAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTaa                                                                                  >  2:78761/1‑100 (MQ=255)
                    aTTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATa                                                                                >  2:495130/1‑100 (MQ=255)
                    aTTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGAAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATa                                                                                >  2:2171281/1‑100 (MQ=255)
                     tttGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAg                                                                               >  1:787444/1‑100 (MQ=255)
                      ttGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGc                                                                              <  2:2834474/100‑1 (MQ=255)
                         ggAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCATTAAAACGACTCTCATAGATTTTACTAATAGCaaa                                                                           >  1:260979/1‑100 (MQ=255)
                          gAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCaaaa                                                                          >  2:2522908/1‑100 (MQ=255)
                             tGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAatat                                                                       >  1:2971125/1‑100 (MQ=255)
                             tGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAatat                                                                       >  2:735559/1‑100 (MQ=255)
                             tGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAatat                                                                       <  1:666567/100‑1 (MQ=255)
                              gCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTTTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAAtata                                                                      <  1:1483177/100‑1 (MQ=255)
                               cGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATaa                                                                     >  1:3377828/1‑100 (MQ=255)
                                ggAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATaaa                                                                    >  2:3236969/1‑100 (MQ=255)
                                ggAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATaaa                                                                    <  2:3182246/100‑1 (MQ=255)
                                 gAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAAc                                                                   <  2:2037432/100‑1 (MQ=255)
                                 gAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAAc                                                                   >  1:2348904/1‑100 (MQ=255)
                                 gAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAACCGTTAAAACGAATATCATAGATTTTACTAATAGCAAAATATAAAc                                                                   >  1:1230201/1‑100 (MQ=255)
                                   aCAATGTGGATTTGAGGAGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCg                                                                 <  1:971151/100‑1 (MQ=255)
                                   aCAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGACGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAAGCg                                                                 >  2:1712373/1‑100 (MQ=255)
                                   aCAAGGTGGATTTGAAGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGACGTTAAAACGACTCTCATACATTTTAATCATAGCAAAATAAAAccc                                                                  >  2:1259268/1‑99 (MQ=255)
                                    cAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGt                                                                >  2:3849459/1‑100 (MQ=255)
                                     aaGGTTGATTTGAGGCGAAAATAAGCTAGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATCTTACTAATAACAAAATATAAAcc                                                                  >  2:1926642/1‑97 (MQ=255)
                                       ggTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTacc                                                             >  2:2324590/1‑100 (MQ=255)
                                        gTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTc                                                               >  2:894773/1‑97 (MQ=255)
                                         tGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCa                                                           <  2:1673142/100‑1 (MQ=255)
                                         tGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTATGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCa                                                           >  2:1558025/1‑100 (MQ=255)
                                           gATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCaaa                                                         <  2:3128705/100‑1 (MQ=255)
                                             tttGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGTAAAATATAAACCGTCCCCaaaaa                                                       >  1:554718/1‑100 (MQ=255)
                                             tttGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCaaaaa                                                       <  2:2218503/100‑1 (MQ=255)
                                             tttGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCaaaaa                                                       >  2:3821458/1‑100 (MQ=255)
                                               tgtggCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTAAAAACGACTCTCCTAGATTTTACTAATAGCAAAATATAAACCGTCCCCaaaaaa                                                      >  2:562847/4‑99 (MQ=255)
                                               tGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAg                                                     >  2:1829293/1‑100 (MQ=255)
                                                  ggCGAAGGTAATCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGc                                                    <  2:67136/98‑1 (MQ=255)
                                                  ggCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCGAAATATAAACCGTCCCCAAAAACGc                                                    >  2:58192/1‑98 (MQ=255)
                                                  ggCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAgaga                                                  >  1:685750/1‑97 (MQ=255)
                                                  ggCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGcca                                                  >  2:45144/1‑100 (MQ=255)
                                                  ggCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGcca                                                  >  2:96924/1‑100 (MQ=255)
                                                  ggCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTACACAAAAAAGcaa                                                  >  1:827398/1‑98 (MQ=255)
                                                     gAAGGTAAGCTGTTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTACTCGCACAATATAAACCGTCACCAAAAAAGccacc                                                <  2:221977/99‑1 (MQ=255)
                                                        ggTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCaccaacc                                            <  2:3519953/100‑1 (MQ=255)
                                                        ggTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGAAAAAAAAAAAGCCaccaacc                                            <  2:3894423/100‑1 (MQ=255)
                                                         gTAAGCTGGTATCTTCTCTGTCTATGAGAGTCGTTAAAACTACTCTCATAGATTTTACTACTAGCATCTTATTAACTGTTTCCAAAAAAGCCaccaacc                                            <  2:3434756/99‑1 (MQ=255)
                                                             gCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACaaa                                       >  1:12278/1‑100 (MQ=255)
                                                               tGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAAcc                                     >  2:3564579/1‑100 (MQ=255)
                                                                ggTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCa                                    >  1:3137797/1‑100 (MQ=255)
                                                                 gTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCa                                    <  1:793502/99‑1 (MQ=255)
                                                                  tAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGa                                  <  2:3358760/100‑1 (MQ=255)
                                                                  tAACTTCTCAGCCTATGAGATCCGTTAAACCGACTCTAATACATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCATAAAcctg                                   <  2:1044939/99‑3 (MQ=255)
                                                                   aaCTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGAc                                 <  1:2071463/100‑1 (MQ=255)
                                                                    aCTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACa                                >  2:1992156/1‑100 (MQ=255)
                                                                    aCTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACa                                <  1:2369186/100‑1 (MQ=255)
                                                                    aCTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACa                                >  2:5785/1‑100 (MQ=255)
                                                                      ttCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACaa                               <  2:3038141/99‑1 (MQ=255)
                                                                       tctcTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAAt                             >  2:2903228/1‑100 (MQ=255)
                                                                       gtaaaatCTATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAAt                             >  1:3181956/7‑100 (MQ=25)
                                                                        ctctGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATa                            >  2:1756251/1‑100 (MQ=255)
                                                                           aatCTATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACAGTTCCCAAAAAAGCCAACAACCACAAAACAGAAAAATACcg                         >  2:2738713/3‑100 (MQ=21)
                                                                            gTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATAcc                          >  2:1205431/1‑98 (MQ=255)
                                                                             tCTATAAGAGACGTTAAGACGATTCTCACAGATTTTACTAATAGCAAAATATAAAACTTCCCCAAAAAAGCCACCTTCCACAAACCAGACAAATACcgag                       <  2:108966/100‑1 (MQ=255)
                                                                              cTATGCGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAAAATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAAAAcc                          >  2:2713088/1‑96 (MQ=37)
                                                                              cTATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACAACAAACAAAAAAAATACcga                        >  1:3536713/1‑98 (MQ=25)
                                                                              cTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagc                      >  1:3696285/1‑100 (MQ=255)
                                                                               tATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcg                     >  2:836132/1‑100 (MQ=38)
                                                                               tATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcg                     >  1:17743/1‑100 (MQ=255)
                                                                               tATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcg                     <  1:2495563/100‑1 (MQ=255)
                                                                                aTGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcga                    >  1:3049579/1‑100 (MQ=38)
                                                                                  gagaGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCAACAAAAAAGCAACCAACCACAAACCAGACAAATACAGAGCgaga                  >  1:2029432/1‑100 (MQ=37)
                                                                                  gagaGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCgaga                  >  2:2456577/1‑100 (MQ=255)
                                                                                  gagaGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCgaga                  >  1:2776263/1‑100 (MQ=255)
                                                                                  gagaGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCgaga                  >  1:938799/1‑100 (MQ=255)
                                                                                   agagTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCAACAACCACAAACCAGACAAATACCGAGCGAGaa                 >  2:3593193/1‑100 (MQ=38)
                                                                                    gagCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAAt                >  1:2405166/1‑100 (MQ=255)
                                                                                     agCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAAtn               <  1:3819811/100‑2 (MQ=255)
                                                                                       ccGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACa             <  1:2179644/100‑1 (MQ=255)
                                                                                        cGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCATCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAACGAAAATACat            >  1:1451447/1‑100 (MQ=255)
                                                                                         gTTAAAACGACTCTCATAGATTTTACTATTAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAAACAAATAACGAGCGTGAATACata           <  2:2977694/100‑1 (MQ=255)
                                                                                          ttttaaCGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACatat          <  1:892808/96‑1 (MQ=255)
                                                                                              aaCGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATAccc      >  2:3810414/1‑100 (MQ=255)
                                                                                              aaCGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATAc        >  1:1883799/1‑98 (MQ=255)
                                                                                              aaCGACTCTCATAGATTTTACAAATAGCAAAATATAAACCGTCCCCAAAAAAAACACCAACCACAAACCAGAGAAATCCCAAGCGAGAATACATATAccc      >  2:591469/1‑100 (MQ=255)
                                                                                                  aCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCaaa   <  2:3773785/99‑1 (MQ=255)
                                                                                                  |                                                                                                   
TTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAA  >  NC_000913/3171765‑3171962

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