Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,053,516:1 (G)6→7 100% coding (2663/2745 nt) torS ← hybrid sensory histidine kinase in two‑component regulatory system with TorR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,053,5101.G68.8% 336.6 / 167.3 138coding (2669/2745 nt)torShybrid sensory histidine kinase in two‑component regulatory system with TorR
Reads supporting (aligned to +/- strand):  ref base . (17/26);  new base G (52/43);  total (69/69)
Fisher's exact test for biased strand distribution p-value = 1.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 153 total aligned reads displayed.

GTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑‑GG‑GGGG‑CAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATGCGCGG  >  NC_000913/1053411‑1053606
                                                                                                    |                                                                                                    
   cTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑‑gg‑g                                                                                                >  2:123876/1‑100 (MQ=255)
     tatcTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGATTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑‑gg‑ggg                                                                                              <  1:443294/97‑1 (MQ=255)
     tAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑‑gg‑ggg                                                                                              <  1:3298309/100‑1 (MQ=255)
     tAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑‑gg‑ggg                                                                                              <  2:1630652/100‑1 (MQ=255)
     tAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑‑gg‑ggg                                                                                              >  1:3405802/1‑100 (MQ=255)
      aaCTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑‑gg‑gggg                                                                                             <  2:606133/100‑1 (MQ=255)
        cTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑‑gg‑gggggc                                                                                           >  1:3135086/1‑98 (MQ=255)
         tGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑‑gg‑gggggca                                                                                          <  1:1392398/100‑4 (MQ=255)
         tGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑‑gg‑gggggc                                                                                           <  1:3041278/99‑3 (MQ=255)
          gtcTGAGATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑‑gg‑gggggcag                                                                                         <  1:2672076/98‑5 (MQ=255)
               atatTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑‑GG‑GGGG‑CAGATAATg                                                                                   <  1:4121867/100‑1 (MQ=255)
                  ttGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAG‑‑GG‑GGGG‑CAGATAATGgc                                                                                 >  1:555522/1‑100 (MQ=255)
                  ttGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAG‑‑GG‑GGGG‑CAGATAATGgc                                                                                 <  2:877889/100‑1 (MQ=255)
                  ttGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑‑GG‑GGGG‑CAGATAATGgct                                                                                <  2:1547361/100‑1 (MQ=255)
                     aTTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAG‑‑GG‑GGGG‑CAGATAATGgctgc                                                                              >  1:2251234/1‑100 (MQ=255)
                     aTTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑‑GG‑GGGG‑CAGATAATGgctgct                                                                             <  2:2370545/100‑1 (MQ=255)
                      ttCAAATCGCGTTCAGGTCTTTCTTATGCAACTATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAG‑‑GG‑GGGG‑CAGATAATGgctgct                                                                             <  1:2493050/100‑1 (MQ=255)
                      ttCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAG‑‑GG‑GGGG‑CAGATAATGgctgct                                                                             >  1:4120637/1‑100 (MQ=255)
                        cAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑‑GG‑GGGG‑CAGATAATGgctgctgct                                                                          >  1:3897283/1‑100 (MQ=255)
                         aaaTCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAG‑‑GG‑GGGG‑CAGATAATGgctgctgct                                                                          >  2:2504458/1‑100 (MQ=255)
                           aTCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAG‑‑GG‑GGGG‑CAGATAATGGCTGCTGCTcc                                                                        >  1:4210464/1‑100 (MQ=255)
                           aTCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAG‑‑GG‑GGGG‑CAGATAATGGCTGCTGCTcc                                                                        <  2:4305761/100‑1 (MQ=255)
                           aTCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAG‑‑GG‑GGGG‑CAGATAATGGCTGCTGCTcc                                                                        >  1:1122114/1‑100 (MQ=255)
                             cgcgTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑‑GG‑GGGG‑CAGATAATGGCTGCTGCTCCAGt                                                                     <  1:514018/100‑1 (MQ=255)
                              gcgTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAG‑‑GG‑GGGG‑CAGATAATGGCTGCTGCTCCAGt                                                                     <  2:711827/100‑1 (MQ=255)
                               cgTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAG‑‑GG‑GGGG‑CAGATAATGGCTGCTGCTCCAGtt                                                                    >  1:1253527/1‑100 (MQ=255)
                                gTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑‑GG‑GGGGTCAGATAATGGCTGCTGCTCCAGtt                                                                    >  1:111110/1‑99 (MQ=255)
                                 ttCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAG‑‑GG‑GGGG‑CAGATAATGGCTGCTGCTCCAGTtgt                                                                  >  1:2196411/1‑100 (MQ=255)
                                   cAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAG‑‑GG‑GGGG‑CAGATAATGGCTGCTGCTCCAGTTGTgc                                                                >  2:2733322/1‑100 (MQ=255)
                                    aGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑‑GG‑GGGG‑CAGATAATGGCTGCTGCTCCAGTTGTgcgc                                                              <  2:3551627/100‑1 (MQ=255)
                                     ggTCTTTCTTCTGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑‑GG‑GTGG‑TAGATAATGGCTGCTGCTCCAGTTGTGCGCa                                                             <  1:2892830/100‑1 (MQ=255)
                                     ggTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑‑GG‑GAGG‑CAGATAATGGCTGCTGCTCCAGTTGTGCGCa                                                             >  1:1684816/1‑100 (MQ=255)
                                      gTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAG‑‑GG‑GGGG‑CAGATAATGGCTGCTGCTCCAGTTGTGCGCa                                                             <  1:24322/100‑1 (MQ=255)
                                       tctttcttATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAG‑‑GG‑GGGG‑CAGATAATGGCTGCTGCTCCAGTTGTGCGCAt                                                            <  2:890573/100‑1 (MQ=255)
                                        ctttcttATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAG‑‑GG‑GGGG‑CAGATAATGGCTGCTGCTCCAGTTGTGCGCATa                                                           <  1:840833/100‑1 (MQ=255)
                                        ctttcttATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAG‑‑GG‑GGGG‑CAGATAATGGCTGCTGCTCCAGTTGTGCGCATa                                                           >  2:2249981/1‑100 (MQ=255)
                                           tcttATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑‑GG‑GGGG‑CAGATAATGGCTGCTGCTCCAGTTGTGCGCATAgctg                                                       <  2:4102927/100‑1 (MQ=255)
                                            cttATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAG‑‑GG‑GGGG‑CAGATAATGGCTGCTGCTCCAGTTGTGCGCATAgctg                                                       >  2:3909366/1‑100 (MQ=255)
                                            cttATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAG‑‑GG‑GGGG‑CAGATAATGGCTGCTGCTCCAGTTGTGCGCATAgctg                                                       <  2:609465/100‑1 (MQ=255)
                                             ttATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑‑GG‑GGGG‑CAGATAATGGCTGCTGCTCCAGTTGTGCGCATAgctggc                                                     <  1:2675812/100‑1 (MQ=255)
                                              tATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAG‑‑GG‑GGGG‑CAGATAATGGCTGCTGCTCCAGTTGTGCGCATAgctggc                                                     >  1:396575/1‑100 (MQ=255)
                                              tATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAG‑‑GG‑GGGG‑CAGATAATGGCTGCTGCTCCAGTTGTGCGCATAgctggc                                                     <  2:2814175/100‑1 (MQ=255)
                                              tATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAG‑‑GG‑GGGG‑CAGATAATGGCTGCTGCTCCAGTTGTGCGCATAgctggc                                                     <  2:3262451/100‑1 (MQ=255)
                                              tATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAG‑‑GG‑GGGG‑CAGATAATGGCTGCTGCTCCAGTTGTGCGCATAgctggc                                                     >  2:4023264/1‑100 (MQ=255)
                                              tATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAG‑‑GG‑GGGG‑CAGATAATGGCTGCTGCTCCAGTTGTGCGCATAgctggc                                                     >  2:4096492/1‑100 (MQ=255)
                                              tATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAG‑‑GG‑GGGG‑CAGATAATGGCTGCTGCTCCAGTTGTGCGCATAgctggc                                                     >  1:567368/1‑100 (MQ=255)
                                              tATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑‑GG‑GGGG‑CAGATAATGGCTGCTGCTCCAGTTGTGCGCATAgctggct                                                    >  1:825862/1‑100 (MQ=255)
                                               aTGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAG‑‑GG‑GGGG‑CAGATAATGGCTGCTGCTCCAGTTGTGCGCATAgctggct                                                    >  2:907441/1‑100 (MQ=255)
                                                tGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAG‑‑GG‑GGGG‑CAGATAATGGCTGCTGCTCCAGTTGTGCGCATAgctggctg                                                   <  2:818435/100‑1 (MQ=255)
                                                  cAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAG‑‑GG‑GGGG‑CAGATAGTGGCTGCTGCTCCAGTTGTGCGCATAGCtggctggc                                                 <  2:1964776/100‑1 (MQ=255)
                                                   aaCCATGCTTTCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑‑GG‑GGGG‑CAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCg                                                <  1:998750/99‑1 (MQ=255)
                                                    aCCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑‑GG‑GGGG‑CAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGat                                              >  2:439756/1‑100 (MQ=255)
                                                      cATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑‑GG‑GGGG‑CAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGatat                                            <  1:4148256/100‑1 (MQ=255)
                                                       aTGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑‑GG‑GGGG‑CAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATg                                           <  2:238419/100‑1 (MQ=255)
                                                        tGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAG‑‑GG‑GGGG‑CAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATg                                           <  2:3007003/100‑1 (MQ=255)
                                                          cTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAG‑‑GG‑GGGG‑CAGATAATGTCTGCTGCTCCAGTTGTGCGCATAGATGGCTGCAGatatccc                                         >  1:2583861/1‑97 (MQ=255)
                                                          cTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAG‑‑GG‑GGGG‑CAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCa                                         >  2:1953982/1‑100 (MQ=255)
                                                            tCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAG‑‑GG‑GGGG‑CAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCAtt                                       >  1:1693660/1‑100 (MQ=255)
                                                            tCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑‑GG‑GGGG‑CAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCAtt                                       >  1:2064026/1‑99 (MQ=255)
                                                               agagCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAG‑‑GG‑GGGG‑CAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTccc                                    >  1:4071371/1‑100 (MQ=255)
                                                               agagCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑‑GG‑GGGG‑CAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCa                                   >  1:977851/1‑100 (MQ=255)
                                                               agagCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑‑GG‑GGGG‑CAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCa                                   <  1:2833923/100‑1 (MQ=255)
                                                                gagCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑‑GG‑GGGG‑CAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAg                                  <  1:4286231/100‑1 (MQ=255)
                                                                gagCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑‑GG‑GGGG‑CAGATAATGGCTGCTGCTCCAGGTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAg                                  <  2:1683886/100‑1 (MQ=255)
                                                                   cGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑‑GG‑GGGG‑CAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACt                               <  2:3160044/100‑1 (MQ=255)
                                                                    ggCACCACTGCGTTCAATTTCTTCGTTTGGAA‑‑‑GGTGGGG‑AAGATAATGGTTGCTGCTCCAGTTGTGCGCATAGCTGGCTGACGATATGCATTCCCAGACt                               <  1:3846486/100‑1 (MQ=255)
                                                                    ggCAACACTGCGTGTAATTTCTTCGTGTGGAAG‑‑GG‑GGGG‑CAGATAATGGCTGCTGCTCCCGTTGTGCGCCTAGCTGGCTGGCGCTATGCCTTCCCAGAc                                >  1:2998072/1‑99 (MQ=255)
                                                                    ggCAACACTGCGTGTAATTTCTTCGTGTGGAAG‑‑GG‑GGGG‑CAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACt                               >  2:2250247/1‑100 (MQ=255)
                                                                      cAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGG‑GGGG‑CAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACt                               >  1:1067443/1‑100 (MQ=255)
                                                                       aaCACTGCGTGTAATTTCTTCGTGTGGAAG‑‑GG‑GGGG‑CAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGa                            <  2:66757/100‑1 (MQ=255)
                                                                        acacTGCGTGTAATTTCTTCGTGTGGAAG‑‑GG‑GGGG‑CAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAg                           >  1:1372569/1‑100 (MQ=255)
                                                                        acacTGCGTGTAATTTCTTCGTGTGGAAG‑‑GG‑GGGG‑CAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAg                           <  1:120573/100‑1 (MQ=255)
                                                                        acacTGCGTGTAATTTCTTCGTGTGGAAG‑‑GG‑GGGG‑CAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAg                           <  1:759067/100‑1 (MQ=255)
                                                                         cacTGCGTGTAATTTCTTCGTGTGGAAG‑‑GG‑GGGG‑CAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGc                          <  1:4096462/100‑1 (MQ=255)
                                                                         cacTGCGTGTAATTTCTTCGTGTGGAAG‑‑GG‑GGGG‑CAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGc                          <  2:1384217/100‑1 (MQ=255)
                                                                          acTGCGTGTAATTTCTTCGTGTGGAAG‑‑GG‑GGGG‑CAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCa                         <  2:1759508/100‑1 (MQ=255)
                                                                          acTGCGTGTAATTTCTTCGTGTGGAAG‑‑GG‑GGGG‑CAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCa                         >  1:3475697/1‑100 (MQ=255)
                                                                           cTGCGTGTAATTTCTTCGTGTGGAA‑‑‑GG‑GGGG‑CAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCTAGCAAc                       <  1:1126161/100‑1 (MQ=255)
                                                                             gCGTGTAATTTCTTCGTGTGGAAG‑‑GG‑GGGG‑CAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACt                      >  2:2673697/1‑100 (MQ=255)
                                                                             gCGTGTAATTTCTTCGTGTGGAAG‑‑GG‑GGGG‑CAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACt                      >  2:3431741/1‑100 (MQ=255)
                                                                                tttAATTTTTGCGTGTGGAA‑‑‑TG‑GGGG‑CAGATAATGTCTGCTGCTCCAGTTGTGCGCATAGCTGCCTGGCGATATGCATTCCCAGACTCGAGCAACTg                     <  1:4318031/95‑1 (MQ=255)
                                                                                tgtAATTTCTTCGTGTGGAAG‑‑GG‑GGGG‑CAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCt                   >  1:689284/1‑100 (MQ=255)
                                                                                 gtAATTTCTTCGTGTGGAAG‑‑GG‑GGGG‑CAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCtt                  >  1:3170163/1‑100 (MQ=255)
                                                                                  tAATTTCTTCGTGTGGAAG‑‑GG‑GGGG‑CAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCttt                 <  1:559359/100‑1 (MQ=255)
                                                                                  tAATTTCTTCGTGTGGAAG‑‑GG‑GGGG‑CAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCttt                 <  2:243432/100‑1 (MQ=255)
                                                                                   aaTTTCTTCGTGTGGAA‑‑‑GG‑GGGG‑CAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCttttt               >  1:566321/1‑100 (MQ=255)
                                                                                   aaTTTCTTCGTGTGGAA‑‑‑GG‑GGGG‑CAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCttttt               <  2:1856183/100‑1 (MQ=255)
                                                                                    aTTTCTTCGTGTGGAAG‑‑GG‑GGGG‑CAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCttttn               <  2:201239/100‑2 (MQ=255)
                                                                                     cgatctcGTGTGGAAG‑‑GG‑GGGG‑CAGATAATGGCTGCTGCTCCAGATGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCttt                 <  2:820978/92‑1 (MQ=255)
                                                                                      ttcttcGTGTGGAAG‑‑GG‑GGGG‑CAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAg             <  2:4318177/100‑1 (MQ=255)
                                                                                      ttcttcGTGTGGAA‑‑‑GG‑GGGG‑CAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGc            <  1:2033690/100‑1 (MQ=255)
                                                                                       tcttcGTGTGGAAG‑‑GG‑GGGG‑CAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGc            <  2:2837242/100‑1 (MQ=255)
                                                                                       tcttcGTGTGGAAG‑‑GG‑GGGG‑CAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGc            <  1:329628/100‑1 (MQ=255)
                                                                                        cttcGTGTGGAAG‑‑GG‑GGGG‑CAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAggt           <  2:2979334/100‑3 (MQ=255)
                                                                                        cttcGTGTGGAAG‑‑GG‑GGGG‑CAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCt           <  2:3163903/100‑1 (MQ=255)
                                                                                         ttcGTGTGGAAG‑‑GG‑GGGG‑CAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTg          >  2:2581028/1‑100 (MQ=255)
                                                                                         ttcGTGTGGAAG‑‑GG‑GGGG‑CAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTg          >  2:4050808/1‑100 (MQ=255)
                                                                                           cGTGTGGAAG‑‑GG‑GGGG‑CAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGAt        <  2:2196421/100‑1 (MQ=255)
                                                                                             tgtgGAAG‑‑GG‑GGGG‑CAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATgc      <  2:3552322/100‑1 (MQ=255)
                                                                                                ggAA‑‑‑GG‑GGGG‑CAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATGCGCgg  >  2:4357022/1‑100 (MQ=255)
                                                                                                    |                                                                                                    
GTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑‑GG‑GGGG‑CAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATGCGCGG  >  NC_000913/1053411‑1053606

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