Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 657,533 A→T 15.5% intergenic (+32/+22) cspE → / ← flc constitutive cold shock family transcription antitermination protein; negative regulator of cspA transcription; RNA melting protein; ssDNA‑binding protein/fluoride efflux channel, dual topology membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913657,5330AT15.5% 346.4 / 53.6 135intergenic (+32/+22)cspE/flcconstitutive cold shock family transcription antitermination protein; negative regulator of cspA transcription; RNA melting protein; ssDNA‑binding protein/fluoride efflux channel, dual topology membrane protein
Reads supporting (aligned to +/- strand):  ref base A (54/60);  new base T (5/16);  total (59/76)
Fisher's exact test for biased strand distribution p-value = 5.63e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.00e-01

Only 100 of 140 total aligned reads displayed.

AGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAAC  >  NC_000913/657434‑657632
                                                                                                   |                                                                                                    
aGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTaa                                                                                                      <  2:2119716/100‑1 (MQ=255)
aGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTaa                                                                                                      <  2:1150954/100‑1 (MQ=255)
  cgcgTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTCCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATc                                                                                                    <  2:1905338/100‑1 (MQ=255)
    cgTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAg                                                                                                  >  2:2070397/1‑100 (MQ=255)
     gcagatTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAg                                                                                                  <  1:219622/97‑1 (MQ=255)
     gTAGCGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGc                                                                                                 >  1:2608538/1‑100 (MQ=255)
     gTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGc                                                                                                 >  1:928656/1‑100 (MQ=255)
      tAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCg                                                                                                <  2:1366078/100‑1 (MQ=255)
      tAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCg                                                                                                <  1:2742248/100‑1 (MQ=255)
       agagTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCCTAATAAGCgg                                                                                               >  1:3788577/1‑100 (MQ=255)
        gagTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCggg                                                                                              >  2:700366/1‑100 (MQ=255)
          gTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCCACAATTCAAAACCCGCTCAATCCGCggg                                                                                              >  1:3596420/1‑98 (MQ=255)
          gTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAATAATTCAAAACCCGCTTAATCAGCGGGtt                                                                                           >  1:1089475/1‑100 (MQ=255)
          gTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGtt                                                                                           <  1:2015016/100‑1 (MQ=255)
          gTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGtt                                                                                           >  2:3723440/1‑100 (MQ=255)
          gTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGGtt                                                                                           <  2:3776152/100‑1 (MQ=255)
           ttCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGttt                                                                                          <  2:3770006/100‑1 (MQ=255)
           ttCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGGttt                                                                                          <  1:829220/100‑1 (MQ=255)
            tCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGGtttt                                                                                         <  1:2710149/100‑1 (MQ=255)
             cGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGttttt                                                                                        >  1:1122987/1‑100 (MQ=255)
             cGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGttat                                                                                         >  1:181964/1‑97 (MQ=255)
             cGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGGttttn                                                                                        <  1:1634114/100‑2 (MQ=255)
             cGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGGttttg                                                                                        <  2:2324498/100‑2 (MQ=39)
             cGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGGttttg                                                                                        <  1:129740/100‑2 (MQ=39)
             cGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGGttttg                                                                                        <  1:1633894/100‑2 (MQ=39)
              gAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGtttttt                                                                                       >  2:3494900/1‑100 (MQ=255)
              gAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGttttt                                                                                        >  2:3278707/1‑99 (MQ=255)
              gAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGGttttga                                                                                       <  1:43465/100‑3 (MQ=39)
              gAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGGttttga                                                                                       <  1:874705/100‑3 (MQ=39)
                aaTCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGtttttttt                                                                                     <  1:2240827/100‑1 (MQ=255)
                aaTCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGtttttttg                                                                                     >  2:3012578/1‑99 (MQ=255)
                aaTCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGGttttgaat                                                                                     <  2:3334517/100‑5 (MQ=38)
                 aTCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTg                                                                                    <  1:3687839/100‑1 (MQ=255)
                 aTCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTg                                                                                    <  2:2274516/100‑1 (MQ=255)
                 aTCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGGttttgaatt                                                                                    <  2:1118989/100‑6 (MQ=38)
                  tCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCCCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTgg                                                                                   <  2:3550017/100‑1 (MQ=255)
                     cTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCt                                                                                <  2:406085/100‑1 (MQ=255)
                       aaCGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCttt                                                                              <  2:1202674/100‑1 (MQ=255)
                         cGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAg                                                                            <  2:357960/100‑1 (MQ=255)
                         cGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAg                                                                            <  2:2358045/100‑1 (MQ=255)
                          ggTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTTGCTCTTTAg                                                                            <  2:1981507/100‑1 (MQ=255)
                              ccAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGc                                                                       >  1:298639/1‑100 (MQ=255)
                               cAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCg                                                                      >  2:1534463/1‑100 (MQ=255)
                                    gCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGa                                                                 >  1:552527/1‑100 (MQ=255)
                                     cccTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAg                                                                <  2:2485566/100‑1 (MQ=255)
                                      ccTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAgg                                                               <  1:1459097/100‑1 (MQ=255)
                                        ttCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGcc                                                             <  1:3427414/100‑1 (MQ=255)
                                           tgctgcAAAAGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAACCCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTATGCCGTTGAGGcc                                                             >  1:2055092/1‑97 (MQ=255)
                                            gctgcAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGaaa                                                         <  2:1938392/100‑1 (MQ=255)
                                            gctgcAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGaaa                                                         >  2:1712981/1‑100 (MQ=255)
                                            gctgcAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGaaa                                                         >  2:1327490/1‑100 (MQ=255)
                                              tgcAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGaaaa                                                        >  2:2908256/1‑99 (MQ=255)
                                              tgcAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAAc                                                       >  2:2551687/1‑100 (MQ=255)
                                                cAAACGTAATCGCTCTGTAAGTGACGTCAGCAAGAATTCAAAAACCGCTGATTCAGCGGGTTTTTTTTGATCTTTAGTGTGAGGTTGAGGCCGAAAACAg                                                     <  1:2996359/100‑1 (MQ=255)
                                                cAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGGTTTTTTTTTTTTTTATCTTGCGGTGGAGGCCGAATACAg                                                     >  1:3351166/1‑100 (MQ=255)
                                                cAAACGTAATCGCTCTGTAAGATAAGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAg                                                     <  1:2066813/100‑1 (MQ=255)
                                                   aCGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCa                                                  >  2:1389834/1‑100 (MQ=255)
                                                    cGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAg                                                 >  2:2414202/1‑100 (MQ=255)
                                                      tAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGaa                                               >  2:3592976/1‑100 (MQ=255)
                                                       aaTCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAAt                                              <  2:2969056/100‑1 (MQ=255)
                                                       aaTCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAAt                                              <  1:1148805/100‑1 (MQ=255)
                                                        aTCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGGATg                                             <  1:1008836/100‑1 (MQ=255)
                                                          cGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGcc                                           <  1:1915178/100‑1 (MQ=255)
                                                          cGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGcc                                           <  1:164800/100‑1 (MQ=255)
                                                            ctctGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAg                                         <  1:694598/100‑1 (MQ=255)
                                                                  aaGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGt                                   <  2:6689/100‑1 (MQ=255)
                                                                  aaGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGt                                   <  1:1438396/100‑1 (MQ=255)
                                                                  aaGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGt                                   <  1:1193295/100‑1 (MQ=255)
                                                                  aaGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGt                                   <  2:705332/100‑1 (MQ=255)
                                                                  aaGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGt                                   >  2:3055960/1‑100 (MQ=255)
                                                                   aGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTc                                  >  2:1014271/1‑100 (MQ=255)
                                                                    gATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGtga                                 >  2:3689821/1‑98 (MQ=255)
                                                                    gATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCa                                 <  2:2249887/100‑1 (MQ=255)
                                                                      tACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATg                               >  1:2723700/1‑100 (MQ=255)
                                                                       aCGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATgg                              <  1:488511/100‑1 (MQ=255)
                                                                        cGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGc                             <  2:906884/100‑1 (MQ=255)
                                                                        cGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGc                             <  2:2608553/100‑1 (MQ=255)
                                                                         gTCAGCAAGAATTCAAAACCCGCTTCATCAGAGGTTTTTTTTTGGTCTTTAGTGTGCCGTTTAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATgg                              <  1:3538900/98‑1 (MQ=255)
                                                                         gTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCa                            <  1:542879/100‑1 (MQ=255)
                                                                            aGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTTTTTTTTTTTGTGCGGCTGAGGTCGAAATCAGTCATTATGTCAGTTGGGTTATGGCacaa                         >  1:2421785/1‑100 (MQ=255)
                                                                             gCAAGAATTCAAAACCCGCTTAATAAGCGGGTTTTTTTTGTTCTTTAGTTTGAGGTTGATGCAGAAAACAGACAGAATGCCGTTTCTGTCATGGCaaa                          >  1:2417837/1‑98 (MQ=255)
                                                                               gaccaaaaaaaaCCCGCTGATTAAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGAc                      >  2:49131/9‑100 (MQ=37)
                                                                                  caaaaaaaaCCCGCTGATTAAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGAcccc                   >  2:346366/6‑100 (MQ=38)
                                                                                   aaaaaaaaCCCGCTGATTAAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGTCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCa                  >  1:3355688/5‑100 (MQ=38)
                                                                                   aTTCAAGACCCGCTTAAGCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGAcccc                   <  1:2392959/99‑1 (MQ=255)
                                                                                   aTTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCa                  >  2:3216400/1‑100 (MQ=255)
                                                                                    ttCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAg                 >  2:3031541/1‑100 (MQ=255)
                                                                                     tCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGa                <  1:115007/100‑1 (MQ=255)
                                                                                     tCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGa                <  1:532590/100‑1 (MQ=255)
                                                                                     tCAAAACACGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGa                >  1:1152555/1‑100 (MQ=255)
                                                                                      cAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGTTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTGATGGCAAAAGACCCCAGa                >  1:2563504/1‑99 (MQ=255)
                                                                                       aaaaCCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAg              >  2:300984/1‑100 (MQ=255)
                                                                                       aaaaCCCGCTGATTAAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAg              >  2:182832/1‑100 (MQ=255)
                                                                                         aaCCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGt            <  2:937158/100‑1 (MQ=255)
                                                                                          aCCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGtt           >  2:2136264/1‑100 (MQ=255)
                                                                                           cccGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTg          <  1:3488849/100‑1 (MQ=255)
                                                                                           cccGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTg          <  1:700362/100‑1 (MQ=255)
                                                                                            ccGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGa         <  2:716240/100‑1 (MQ=255)
                                                                                                   aTCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAAc  >  2:3085840/1‑100 (MQ=255)
                                                                                                   aTCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAAc  <  2:1954077/100‑1 (MQ=255)
                                                                                                   |                                                                                                    
AGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAAC  >  NC_000913/657434‑657632

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: