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.

AGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAAC  >  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)
          gTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAATAATTCAAAACCCGCTTAATCAGCGGG‑tt                                                                                           >  1:1089475/1‑100 (MQ=255)
          gTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑tt                                                                                           <  1:2015016/100‑1 (MQ=255)
          gTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑tt                                                                                           >  2:3723440/1‑100 (MQ=255)
          gTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGG‑tt                                                                                           <  2:3776152/100‑1 (MQ=255)
           ttCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑ttt                                                                                          <  2:3770006/100‑1 (MQ=255)
           ttCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGG‑ttt                                                                                          <  1:829220/100‑1 (MQ=255)
            tCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGG‑tttt                                                                                         <  1:2710149/100‑1 (MQ=255)
             cGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑ttttt                                                                                        >  1:1122987/1‑100 (MQ=255)
             cGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑ttat                                                                                         >  1:181964/1‑97 (MQ=255)
             cGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGG‑ttttn                                                                                        <  1:1634114/100‑2 (MQ=255)
             cGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGG‑ttttg                                                                                        <  2:2324498/100‑2 (MQ=39)
             cGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGG‑ttttg                                                                                        <  1:129740/100‑2 (MQ=39)
             cGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGG‑ttttg                                                                                        <  1:1633894/100‑2 (MQ=39)
              gAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑tttttt                                                                                       >  2:3494900/1‑100 (MQ=255)
              gAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑ttttt                                                                                        >  2:3278707/1‑99 (MQ=255)
              gAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGG‑ttttga                                                                                       <  1:43465/100‑3 (MQ=39)
              gAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGG‑ttttga                                                                                       <  1:874705/100‑3 (MQ=39)
                aaTCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑tttttttt                                                                                     <  1:2240827/100‑1 (MQ=255)
                aaTCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑tttttttg                                                                                     >  2:3012578/1‑99 (MQ=255)
                aaTCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGG‑ttttgaat                                                                                     <  2:3334517/100‑5 (MQ=38)
                 aTCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTg                                                                                    <  1:3687839/100‑1 (MQ=255)
                 aTCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTg                                                                                    <  2:2274516/100‑1 (MQ=255)
                 aTCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGG‑ttttgaatt                                                                                    <  2:1118989/100‑6 (MQ=38)
                  tCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCCCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTgg                                                                                   <  2:3550017/100‑1 (MQ=255)
                     cTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCt                                                                                <  2:406085/100‑1 (MQ=255)
                       aaCGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCttt                                                                              <  2:1202674/100‑1 (MQ=255)
                         cGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAg                                                                            <  2:357960/100‑1 (MQ=255)
                         cGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAg                                                                            <  2:2358045/100‑1 (MQ=255)
                          ggTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTTGCTCTTTAg                                                                            <  2:1981507/100‑1 (MQ=255)
                              ccAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGc                                                                       >  1:298639/1‑100 (MQ=255)
                               cAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCg                                                                      >  2:1534463/1‑100 (MQ=255)
                                    gCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGa                                                                 >  1:552527/1‑100 (MQ=255)
                                     cccTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAg                                                                <  2:2485566/100‑1 (MQ=255)
                                      ccTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAgg                                                               <  1:1459097/100‑1 (MQ=255)
                                        ttCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGcc                                                             <  1:3427414/100‑1 (MQ=255)
                                           tgctgcAAAAGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAACCCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTATGCCGTTGAGGcc                                                             >  1:2055092/1‑97 (MQ=255)
                                            gctgcAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGaaa                                                         <  2:1938392/100‑1 (MQ=255)
                                            gctgcAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGaaa                                                         >  2:1712981/1‑100 (MQ=255)
                                            gctgcAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGaaa                                                         >  2:1327490/1‑100 (MQ=255)
                                              tgcAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGaaaa                                                        >  2:2908256/1‑99 (MQ=255)
                                              tgcAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAAc                                                       >  2:2551687/1‑100 (MQ=255)
                                                cAAACGTAATCGCTCTGTAAGTGACGTCAGCAAGAATTCAAAAACCGCTGATTCAGCGGG‑TTTTTTTTGATCTTTAGTGTGAGGTTGAGGCCGAAAACAg                                                     <  1:2996359/100‑1 (MQ=255)
                                                cAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑GTTTTTTTTTTTTTTATCTTGCGGTGGAGGCCGAATACAg                                                     >  1:3351166/1‑100 (MQ=255)
                                                cAAACGTAATCGCTCTGTAAGATAAGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAg                                                     <  1:2066813/100‑1 (MQ=255)
                                                   aCGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCa                                                  >  2:1389834/1‑100 (MQ=255)
                                                    cGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAg                                                 >  2:2414202/1‑100 (MQ=255)
                                                      tAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGaa                                               >  2:3592976/1‑100 (MQ=255)
                                                       aaTCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAAt                                              <  2:2969056/100‑1 (MQ=255)
                                                       aaTCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAAt                                              <  1:1148805/100‑1 (MQ=255)
                                                        aTCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGGATg                                             <  1:1008836/100‑1 (MQ=255)
                                                          cGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGcc                                           <  1:1915178/100‑1 (MQ=255)
                                                          cGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGcc                                           <  1:164800/100‑1 (MQ=255)
                                                            ctctGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAg                                         <  1:694598/100‑1 (MQ=255)
                                                                  aaGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGt                                   <  2:6689/100‑1 (MQ=255)
                                                                  aaGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGt                                   <  1:1438396/100‑1 (MQ=255)
                                                                  aaGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGt                                   <  1:1193295/100‑1 (MQ=255)
                                                                  aaGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGt                                   <  2:705332/100‑1 (MQ=255)
                                                                  aaGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGt                                   >  2:3055960/1‑100 (MQ=255)
                                                                   aGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTc                                  >  2:1014271/1‑100 (MQ=255)
                                                                    gATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGtga                                 >  2:3689821/1‑98 (MQ=255)
                                                                    gATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCa                                 <  2:2249887/100‑1 (MQ=255)
                                                                      tACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATg                               >  1:2723700/1‑100 (MQ=255)
                                                                       aCGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATgg                              <  1:488511/100‑1 (MQ=255)
                                                                        cGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGc                             <  2:906884/100‑1 (MQ=255)
                                                                        cGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGc                             <  2:2608553/100‑1 (MQ=255)
                                                                         gTCAGCAAGAATTCAAAACCCGCTTCATCAGAGGT‑TTTTTTTTGGTCTTTAGTGTGCCGTTTAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATgg                              <  1:3538900/98‑1 (MQ=255)
                                                                         gTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCa                            <  1:542879/100‑1 (MQ=255)
                                                                            aGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTTTTTTTTTTTGTGCGGCTGAGGTCGAAATCAGTCATTATGTCAGTTGGGTTATGGCacaa                         >  1:2421785/1‑100 (MQ=255)
                                                                             gCAAGAATTCAAAACCCGCTTAATAAGCGGG‑TTTTTTTTGTTCTTTAGTTTGAGGTTGATGCAGAAAACAGACAGAATGCCGTTTCTGTCATGGCaaa                          >  1:2417837/1‑98 (MQ=255)
                                                                               gaccaaaaaaaaCCCGCTGATTAAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGAc                      >  2:49131/9‑100 (MQ=37)
                                                                                  caaaaaaaaCCCGCTGATTAAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGAcccc                   >  2:346366/6‑100 (MQ=38)
                                                                                   aaaaaaaaCCCGCTGATTAAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGTCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCa                  >  1:3355688/5‑100 (MQ=38)
                                                                                   aTTCAAGACCCGCTTAAGCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGAcccc                   <  1:2392959/99‑1 (MQ=255)
                                                                                   aTTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCa                  >  2:3216400/1‑100 (MQ=255)
                                                                                    ttCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAg                 >  2:3031541/1‑100 (MQ=255)
                                                                                     tCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGa                <  1:115007/100‑1 (MQ=255)
                                                                                     tCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGa                <  1:532590/100‑1 (MQ=255)
                                                                                     tCAAAACACGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGa                >  1:1152555/1‑100 (MQ=255)
                                                                                      cAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGTTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTGATGGCAAAAGACCCCAGa                >  1:2563504/1‑99 (MQ=255)
                                                                                       aaaaCCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAg              >  2:300984/1‑100 (MQ=255)
                                                                                       aaaaCCCGCTGATTAAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAg              >  2:182832/1‑100 (MQ=255)
                                                                                         aaCCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGt            <  2:937158/100‑1 (MQ=255)
                                                                                          aCCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGtt           >  2:2136264/1‑100 (MQ=255)
                                                                                           cccGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTg          <  1:3488849/100‑1 (MQ=255)
                                                                                           cccGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTg          <  1:700362/100‑1 (MQ=255)
                                                                                            ccGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGa         <  2:716240/100‑1 (MQ=255)
                                                                                                   aTCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAAc  >  2:3085840/1‑100 (MQ=255)
                                                                                                   aTCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAAc  <  2:1954077/100‑1 (MQ=255)
                                                                                                   |                                                                                                    
AGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAAC  >  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: ‑