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‑‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATGCGCGG  >  NC_000913/1053411‑1053606
                                                                                                    |                                                                                                    
   cTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑‑ggg                                                                                                >  2:123876/1‑100 (MQ=255)
     tatcTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGATTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑‑ggggg                                                                                              <  1:443294/97‑1 (MQ=255)
     tAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑‑ggggg                                                                                              <  1:3298309/100‑1 (MQ=255)
     tAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑‑ggggg                                                                                              <  2:1630652/100‑1 (MQ=255)
     tAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑‑ggggg                                                                                              >  1:3405802/1‑100 (MQ=255)
      aaCTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑‑gggggg                                                                                             <  2:606133/100‑1 (MQ=255)
        cTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑‑gggggggc                                                                                           >  1:3135086/1‑98 (MQ=255)
         tGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑‑gggggggca                                                                                          <  1:1392398/100‑4 (MQ=255)
         tGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑‑gggggggc                                                                                           <  1:3041278/99‑3 (MQ=255)
          gtcTGAGATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑‑gggggggcag                                                                                         <  1:2672076/98‑5 (MQ=255)
               atatTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑‑GGGGGGCAGATAATg                                                                                   <  1:4121867/100‑1 (MQ=255)
                  ttGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAG‑‑GGGGGGCAGATAATGgc                                                                                 >  1:555522/1‑100 (MQ=255)
                  ttGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAG‑‑GGGGGGCAGATAATGgc                                                                                 <  2:877889/100‑1 (MQ=255)
                  ttGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑‑GGGGGGCAGATAATGgct                                                                                <  2:1547361/100‑1 (MQ=255)
                     aTTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAG‑‑GGGGGGCAGATAATGgctgc                                                                              >  1:2251234/1‑100 (MQ=255)
                     aTTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑‑GGGGGGCAGATAATGgctgct                                                                             <  2:2370545/100‑1 (MQ=255)
                      ttCAAATCGCGTTCAGGTCTTTCTTATGCAACTATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAG‑‑GGGGGGCAGATAATGgctgct                                                                             <  1:2493050/100‑1 (MQ=255)
                      ttCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAG‑‑GGGGGGCAGATAATGgctgct                                                                             >  1:4120637/1‑100 (MQ=255)
                        cAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑‑GGGGGGCAGATAATGgctgctgct                                                                          >  1:3897283/1‑100 (MQ=255)
                         aaaTCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAG‑‑GGGGGGCAGATAATGgctgctgct                                                                          >  2:2504458/1‑100 (MQ=255)
                           aTCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAG‑‑GGGGGGCAGATAATGGCTGCTGCTcc                                                                        >  1:4210464/1‑100 (MQ=255)
                           aTCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAG‑‑GGGGGGCAGATAATGGCTGCTGCTcc                                                                        <  2:4305761/100‑1 (MQ=255)
                           aTCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAG‑‑GGGGGGCAGATAATGGCTGCTGCTcc                                                                        >  1:1122114/1‑100 (MQ=255)
                             cgcgTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGt                                                                     <  1:514018/100‑1 (MQ=255)
                              gcgTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAG‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGt                                                                     <  2:711827/100‑1 (MQ=255)
                               cgTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAG‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGtt                                                                    >  1:1253527/1‑100 (MQ=255)
                                gTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑‑GGGGGGTCAGATAATGGCTGCTGCTCCAGtt                                                                    >  1:111110/1‑99 (MQ=255)
                                 ttCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAG‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTtgt                                                                  >  1:2196411/1‑100 (MQ=255)
                                   cAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAG‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTgc                                                                >  2:2733322/1‑100 (MQ=255)
                                    aGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTgcgc                                                              <  2:3551627/100‑1 (MQ=255)
                                     ggTCTTTCTTCTGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑‑GGGTGGTAGATAATGGCTGCTGCTCCAGTTGTGCGCa                                                             <  1:2892830/100‑1 (MQ=255)
                                     ggTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑‑GGGAGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCa                                                             >  1:1684816/1‑100 (MQ=255)
                                      gTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAG‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCa                                                             <  1:24322/100‑1 (MQ=255)
                                       tctttcttATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAG‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCAt                                                            <  2:890573/100‑1 (MQ=255)
                                        ctttcttATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAG‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATa                                                           <  1:840833/100‑1 (MQ=255)
                                        ctttcttATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAG‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATa                                                           >  2:2249981/1‑100 (MQ=255)
                                           tcttATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAgctg                                                       <  2:4102927/100‑1 (MQ=255)
                                            cttATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAG‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAgctg                                                       >  2:3909366/1‑100 (MQ=255)
                                            cttATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAG‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAgctg                                                       <  2:609465/100‑1 (MQ=255)
                                             ttATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAgctggc                                                     <  1:2675812/100‑1 (MQ=255)
                                              tATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAG‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAgctggc                                                     >  1:396575/1‑100 (MQ=255)
                                              tATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAG‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAgctggc                                                     <  2:2814175/100‑1 (MQ=255)
                                              tATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAG‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAgctggc                                                     <  2:3262451/100‑1 (MQ=255)
                                              tATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAG‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAgctggc                                                     >  2:4023264/1‑100 (MQ=255)
                                              tATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAG‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAgctggc                                                     >  2:4096492/1‑100 (MQ=255)
                                              tATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAG‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAgctggc                                                     >  1:567368/1‑100 (MQ=255)
                                              tATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAgctggct                                                    >  1:825862/1‑100 (MQ=255)
                                               aTGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAG‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAgctggct                                                    >  2:907441/1‑100 (MQ=255)
                                                tGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAG‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAgctggctg                                                   <  2:818435/100‑1 (MQ=255)
                                                  cAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAG‑‑GGGGGGCAGATAGTGGCTGCTGCTCCAGTTGTGCGCATAGCtggctggc                                                 <  2:1964776/100‑1 (MQ=255)
                                                   aaCCATGCTTTCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCg                                                <  1:998750/99‑1 (MQ=255)
                                                    aCCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGat                                              >  2:439756/1‑100 (MQ=255)
                                                      cATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGatat                                            <  1:4148256/100‑1 (MQ=255)
                                                       aTGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATg                                           <  2:238419/100‑1 (MQ=255)
                                                        tGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAG‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATg                                           <  2:3007003/100‑1 (MQ=255)
                                                          cTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAG‑‑GGGGGGCAGATAATGTCTGCTGCTCCAGTTGTGCGCATAGATGGCTGCAGatatccc                                         >  1:2583861/1‑97 (MQ=255)
                                                          cTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAG‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCa                                         >  2:1953982/1‑100 (MQ=255)
                                                            tCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAG‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCAtt                                       >  1:1693660/1‑100 (MQ=255)
                                                            tCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCAtt                                       >  1:2064026/1‑99 (MQ=255)
                                                               agagCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAG‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTccc                                    >  1:4071371/1‑100 (MQ=255)
                                                               agagCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCa                                   >  1:977851/1‑100 (MQ=255)
                                                               agagCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCa                                   <  1:2833923/100‑1 (MQ=255)
                                                                gagCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAg                                  <  1:4286231/100‑1 (MQ=255)
                                                                gagCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGGTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAg                                  <  2:1683886/100‑1 (MQ=255)
                                                                   cGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACt                               <  2:3160044/100‑1 (MQ=255)
                                                                    ggCACCACTGCGTTCAATTTCTTCGTTTGGAA‑‑‑GGTGGGGAAGATAATGGTTGCTGCTCCAGTTGTGCGCATAGCTGGCTGACGATATGCATTCCCAGACt                               <  1:3846486/100‑1 (MQ=255)
                                                                    ggCAACACTGCGTGTAATTTCTTCGTGTGGAAG‑‑GGGGGGCAGATAATGGCTGCTGCTCCCGTTGTGCGCCTAGCTGGCTGGCGCTATGCCTTCCCAGAc                                >  1:2998072/1‑99 (MQ=255)
                                                                    ggCAACACTGCGTGTAATTTCTTCGTGTGGAAG‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACt                               >  2:2250247/1‑100 (MQ=255)
                                                                      cAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACt                               >  1:1067443/1‑100 (MQ=255)
                                                                       aaCACTGCGTGTAATTTCTTCGTGTGGAAG‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGa                            <  2:66757/100‑1 (MQ=255)
                                                                        acacTGCGTGTAATTTCTTCGTGTGGAAG‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAg                           >  1:1372569/1‑100 (MQ=255)
                                                                        acacTGCGTGTAATTTCTTCGTGTGGAAG‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAg                           <  1:120573/100‑1 (MQ=255)
                                                                        acacTGCGTGTAATTTCTTCGTGTGGAAG‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAg                           <  1:759067/100‑1 (MQ=255)
                                                                         cacTGCGTGTAATTTCTTCGTGTGGAAG‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGc                          <  1:4096462/100‑1 (MQ=255)
                                                                         cacTGCGTGTAATTTCTTCGTGTGGAAG‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGc                          <  2:1384217/100‑1 (MQ=255)
                                                                          acTGCGTGTAATTTCTTCGTGTGGAAG‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCa                         <  2:1759508/100‑1 (MQ=255)
                                                                          acTGCGTGTAATTTCTTCGTGTGGAAG‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCa                         >  1:3475697/1‑100 (MQ=255)
                                                                           cTGCGTGTAATTTCTTCGTGTGGAA‑‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCTAGCAAc                       <  1:1126161/100‑1 (MQ=255)
                                                                             gCGTGTAATTTCTTCGTGTGGAAG‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACt                      >  2:2673697/1‑100 (MQ=255)
                                                                             gCGTGTAATTTCTTCGTGTGGAAG‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACt                      >  2:3431741/1‑100 (MQ=255)
                                                                                tttAATTTTTGCGTGTGGAA‑‑‑TGGGGGCAGATAATGTCTGCTGCTCCAGTTGTGCGCATAGCTGCCTGGCGATATGCATTCCCAGACTCGAGCAACTg                     <  1:4318031/95‑1 (MQ=255)
                                                                                tgtAATTTCTTCGTGTGGAAG‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCt                   >  1:689284/1‑100 (MQ=255)
                                                                                 gtAATTTCTTCGTGTGGAAG‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCtt                  >  1:3170163/1‑100 (MQ=255)
                                                                                  tAATTTCTTCGTGTGGAAG‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCttt                 <  1:559359/100‑1 (MQ=255)
                                                                                  tAATTTCTTCGTGTGGAAG‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCttt                 <  2:243432/100‑1 (MQ=255)
                                                                                   aaTTTCTTCGTGTGGAA‑‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCttttt               >  1:566321/1‑100 (MQ=255)
                                                                                   aaTTTCTTCGTGTGGAA‑‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCttttt               <  2:1856183/100‑1 (MQ=255)
                                                                                    aTTTCTTCGTGTGGAAG‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCttttn               <  2:201239/100‑2 (MQ=255)
                                                                                     cgatctcGTGTGGAAG‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGATGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCttt                 <  2:820978/92‑1 (MQ=255)
                                                                                      ttcttcGTGTGGAAG‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAg             <  2:4318177/100‑1 (MQ=255)
                                                                                      ttcttcGTGTGGAA‑‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGc            <  1:2033690/100‑1 (MQ=255)
                                                                                       tcttcGTGTGGAAG‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGc            <  2:2837242/100‑1 (MQ=255)
                                                                                       tcttcGTGTGGAAG‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGc            <  1:329628/100‑1 (MQ=255)
                                                                                        cttcGTGTGGAAG‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAggt           <  2:2979334/100‑3 (MQ=255)
                                                                                        cttcGTGTGGAAG‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCt           <  2:3163903/100‑1 (MQ=255)
                                                                                         ttcGTGTGGAAG‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTg          >  2:2581028/1‑100 (MQ=255)
                                                                                         ttcGTGTGGAAG‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTg          >  2:4050808/1‑100 (MQ=255)
                                                                                           cGTGTGGAAG‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGAt        <  2:2196421/100‑1 (MQ=255)
                                                                                             tgtgGAAG‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATgc      <  2:3552322/100‑1 (MQ=255)
                                                                                                ggAA‑‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATGCGCgg  >  2:4357022/1‑100 (MQ=255)
                                                                                                    |                                                                                                    
GTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAA‑‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATGCGCGG  >  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: