Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 657,535 C→A 15.7% intergenic (+34/+20) 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,5350CA15.7% 378.8 / 52.5 137intergenic (+34/+20)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 C (53/62);  new base A (7/15);  total (60/77)
Fisher's exact test for biased strand distribution p-value = 2.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.20e-01

Only 100 of 138 total aligned reads displayed.

CGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAAC  >  NC_000913/657436‑657632
                                                                                                   |                                                                                                  
cgcgTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTCCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATc                                                                                                    <  2:1905338/100‑1 (MQ=255)
  cgTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAg                                                                                                  >  2:2070397/1‑100 (MQ=255)
   gTAGCGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGc                                                                                                 >  1:2608538/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)
      gagTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGTTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCggg                                                                                              >  2:2336324/1‑100 (MQ=255)
      gagTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCggg                                                                                              >  2:700366/1‑100 (MQ=255)
      gagTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCggg                                                                                              <  1:689353/100‑1 (MQ=255)
        gTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCCACAATTCAAAACCCGCTCAATCCGCggg                                                                                              >  1:3596420/1‑98 (MQ=255)
        gTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAATAATTCAAAACCCGCTTAATCAGCGGGtt                                                                                           >  1:1089475/1‑100 (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)
          tCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGtttt                                                                                         <  2:2815865/100‑1 (MQ=255)
          tCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGGtttt                                                                                         <  1:2710149/100‑1 (MQ=255)
           cGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGttttt                                                                                        >  1:1122987/1‑100 (MQ=255)
           cGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGGttttn                                                                                        <  1:1634114/100‑2 (MQ=255)
           cGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGGttttg                                                                                        <  1:129740/100‑2 (MQ=39)
           cGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGGttttg                                                                                        <  2:2324498/100‑2 (MQ=39)
           cGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGGtttt                                                                                         <  2:364234/99‑1 (MQ=39)
            gAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGtttttt                                                                                       >  2:3494900/1‑100 (MQ=255)
            gAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGttttt                                                                                        >  2:3278707/1‑99 (MQ=255)
            gAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGGttttga                                                                                       <  1:874705/100‑3 (MQ=39)
            gAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGGttttga                                                                                       <  1:43465/100‑3 (MQ=39)
             aaaTCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGttttttt                                                                                      <  2:1528063/100‑1 (MQ=255)
              aaTCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGtttttttt                                                                                     <  1:2240827/100‑1 (MQ=255)
              aaTCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGtttttttg                                                                                     >  2:3012578/1‑99 (MQ=255)
               aTCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTg                                                                                    >  2:2072035/1‑100 (MQ=255)
                tCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTgg                                                                                   <  2:3454691/100‑1 (MQ=255)
                tCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCCCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTgg                                                                                   <  2:3550017/100‑1 (MQ=255)
                 cACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGt                                                                                  <  1:1740366/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)
                        ggTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTttttt                                                                             >  1:2437249/1‑97 (MQ=255)
                            ccAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGc                                                                       <  2:305363/100‑1 (MQ=255)
                                   cccTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAg                                                                <  2:1953134/100‑1 (MQ=255)
                                   cccTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAg                                                                <  2:2485566/100‑1 (MQ=255)
                                    ccTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAgg                                                               <  1:1459097/100‑1 (MQ=255)
                                      ttCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGcc                                                             <  1:3427414/100‑1 (MQ=255)
                                      ttCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGcc                                                             <  1:77287/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)
                                            tgcAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAAc                                                       <  1:1470466/100‑1 (MQ=255)
                                             gcAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACa                                                      >  2:2235083/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)
                                                   ggaaTCGCTGTGTAAGATACGTCAGCAAGAATTCAAGACCCGATTAATCAGCGGGTTTTGTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAg                                                 <  1:3226431/97‑1 (MQ=255)
                                                     aaTCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAAt                                              <  1:1148805/100‑1 (MQ=255)
                                                      aTCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGGATg                                             <  1:1008836/100‑1 (MQ=255)
                                                        cGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGcc                                           <  1:164800/100‑1 (MQ=255)
                                                        cGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGcc                                           <  1:1915178/100‑1 (MQ=255)
                                                         gCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGATTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCa                                          >  1:343865/1‑100 (MQ=255)
                                                          ctctGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAg                                         <  1:2274930/100‑1 (MQ=255)
                                                          ctctGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAg                                         <  1:694598/100‑1 (MQ=255)
                                                                aaGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGt                                   <  1:1438396/100‑1 (MQ=255)
                                                                aaGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGt                                   <  1:1193295/100‑1 (MQ=255)
                                                                aaGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGt                                   <  2:6689/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)
                                                                  gATACGCCAGCAAGAATTCAAACCCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTc                                  <  2:1122992/99‑1 (MQ=255)
                                                                    tACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATg                               <  2:2979814/100‑1 (MQ=255)
                                                                      cGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGc                             <  2:906884/100‑1 (MQ=255)
                                                                      cGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGGTTTTTTAGGTCTTTAATTTGCGGTTTAGGACGAAAACAGCCAGAATGCCAGTGCGATCATgac                             >  1:2631142/1‑98 (MQ=255)
                                                                       gTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTTGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCa                            <  1:194203/100‑1 (MQ=255)
                                                                       gTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCa                            <  1:542879/100‑1 (MQ=255)
                                                                          aGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTTTTTTTTTTTGTGCGGCTGAGGTCGAAATCAGTCATTATGTCAGTTGGGTTATGGCacaa                         >  1:2421785/1‑100 (MQ=255)
                                                                          aGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTTTTCTTTAGTGTTCGTTTGAGGCCGAAAACAGCCATAATGCCAGTGCGGTCAt                                >  1:502012/1‑93 (MQ=255)
                                                                          aGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCaaaa                         >  1:2841677/1‑100 (MQ=255)
                                                                            cAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCCTGGCaaaa                         >  1:3565870/1‑98 (MQ=255)
                                                                             gaccaaaaaaaaCCCGCTGATTAAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGAc                      >  2:49131/9‑100 (MQ=37)
                                                                              aGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGAcc                     >  2:3089563/1‑100 (MQ=255)
                                                                               gAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGAccc                    >  1:1151386/1‑100 (MQ=255)
                                                                               ccaaaaaaaaCCCGCTGATTAAGCGGGTTTTTTTTGGTCTTTAGTGTGCGTTTGAGGACGAAAACAGCCAGAATACCACTGCCGTCATTGCAAAAGAccc                    >  1:219133/7‑100 (MQ=21)
                                                                                caaaaaaaaCCCGCTGATTAAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGAcccc                   >  2:346366/6‑100 (MQ=38)
                                                                                aaTTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGAcccc                   <  2:2419256/100‑1 (MQ=255)
                                                                                 aTTCAAGACCCGCTTAAGCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGAcccc                   <  1:2392959/99‑1 (MQ=255)
                                                                                 aTTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCa                  >  2:3216400/1‑100 (MQ=255)
                                                                                   tCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGa                <  1:115007/100‑1 (MQ=255)
                                                                                   tCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGa                <  1:532590/100‑1 (MQ=255)
                                                                                     aaaaCCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAg              >  2:300984/1‑100 (MQ=255)
                                                                                     aaaaCCCGCTGATTAAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAg              >  2:182832/1‑100 (MQ=255)
                                                                                       aaCCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGt            <  2:937158/100‑1 (MQ=255)
                                                                                         cccGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTg          <  1:3488849/100‑1 (MQ=255)
                                                                                         cccGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTg          <  2:1007030/100‑1 (MQ=255)
                                                                                          ccGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGa         <  2:716240/100‑1 (MQ=255)
                                                                                                 aTCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAAc  <  1:392197/100‑1 (MQ=255)
                                                                                                   |                                                                                                  
CGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAAC  >  NC_000913/657436‑657632

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