Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 657,531 T→G 15.4% intergenic (+30/+24) 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,5310TG15.4% 347.7 / 57.4 140intergenic (+30/+24)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 T (55/62);  new base G (5/17);  total (61/79)
Fisher's exact test for biased strand distribution p-value = 3.77e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.38e-01

Only 100 of 140 total aligned reads displayed.

CAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGA  >  NC_000913/657433‑657625
                                                                                                  |                                                                                               
cAGCGCGTAGAGTTCGAAATCACTAACGGTTCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTa                                                                                                <  2:1199867/100‑1 (MQ=255)
cAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTa                                                                                                >  1:3769976/1‑100 (MQ=255)
 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)
         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)
             tCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑tttt                                                                                  <  2:2815865/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                                                                                 <  1:129740/100‑2 (MQ=39)
              cGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGG‑ttttg                                                                                 <  2:2324498/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                                                                                <  2:1902935/100‑3 (MQ=38)
               gAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGG‑ttttga                                                                                <  1:874705/100‑3 (MQ=39)
               gAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGG‑ttttga                                                                                <  1:43465/100‑3 (MQ=39)
                aaaTCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑ttttttt                                                                               <  2:1528063/100‑1 (MQ=255)
                 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)
                               ccAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGc                                                                <  2:305363/100‑1 (MQ=255)
                                cAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCg                                                               >  2:1534463/1‑100 (MQ=255)
                                     gCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGa                                                          >  1:552527/1‑100 (MQ=255)
                                      cccTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAg                                                         <  2:1953134/100‑1 (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:1712981/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)
                                                     cGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAg                                          >  2:2414202/1‑100 (MQ=255)
                                                       tAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGaa                                        >  2:3592976/1‑100 (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                            <  1:1438396/100‑1 (MQ=255)
                                                                   aaGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGt                            <  2:6689/100‑1 (MQ=255)
                                                                   aaGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGt                            <  1:1193295/100‑1 (MQ=255)
                                                                   aaGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGt                            >  2:3055960/1‑100 (MQ=255)
                                                                   aaGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGt                            <  2:705332/100‑1 (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                        <  2:2979814/100‑1 (MQ=255)
                                                                       tACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATg                        >  1:2723700/1‑100 (MQ=255)
                                                                        aCGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATgg                       <  1:488511/100‑1 (MQ=255)
                                                                         cGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGc                      <  2:2608553/100‑1 (MQ=255)
                                                                         cGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGc                      <  2:906884/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)
                                                                                   aaTTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGAcccc            >  2:1348965/1‑100 (MQ=255)
                                                                                    aaaaaaaaCCCGCTGATTAAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGTCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCa           >  1:3355688/5‑100 (MQ=38)
                                                                                    aTTCAAGACCCGCTTAAGCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGAcccc            <  1:2392959/99‑1 (MQ=255)
                                                                                    aTTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCa           >  2:3216400/1‑100 (MQ=255)
                                                                                      tCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGa         <  1:532590/100‑1 (MQ=255)
                                                                                      tCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGa         <  1:115007/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)
                                                                                          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)
                                                                                                  |                                                                                               
CAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGA  >  NC_000913/657433‑657625

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: ‑