Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 657,535 C→A 5.8% 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,5350CA 5.8% 309.7 / 6.7 104intergenic (+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 (41/57);  new base A (3/3);  total (44/60)
Fisher's exact test for biased strand distribution p-value = 6.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.75e-01

Only 100 of 105 total aligned reads displayed.

CTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAACACCACTTCT  >  NC_000913/657387‑657684
                                                                                                                                                    |                                                                                                                                                       
cTTCTCTGCAATCCAGACTAATGGTGTTAAAACTCTTGCTGAATGTCAGCGCGTAGAGTTCTAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCa                                                                                                                                                        <  1:865757/150‑1 (MQ=255)
  tctcTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGc                                                                                                                                                      <  1:1321421/150‑1 (MQ=255)
   ctctGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCg                                                                                                                                                     >  1:361724/1‑150 (MQ=255)
    tctGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAAAGGTGCCAAAGGACCTTCTGCTGCAAACGTAATCGCTATGTAAGAAACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCgg                                                                                                                                                    >  1:100639/1‑150 (MQ=255)
     ctGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCggg                                                                                                                                                   >  1:816155/1‑150 (MQ=255)
         aaTCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGG‑‑tttt                                                                                                                                             <  2:649651/150‑1 (MQ=255)
          aTCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGG‑‑ttttg                                                                                                                                            <  1:1679641/150‑2 (MQ=255)
          aTCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGG‑‑ttttg                                                                                                                                            <  2:591899/150‑2 (MQ=255)
            ccAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCTCTTAATCAGCGGG‑‑ttttttt                                                                                                                                          >  2:674804/1‑150 (MQ=255)
            ccAGACTAATGGTTTAAAAACTCATGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAATGCCCTTCTGCTGCAAACGTAATCGCTCAGTAAGATACGTCAGAAAGAATTCAAAACCCTCTTAATCAGCGGG‑‑ttttttt                                                                                                                                          >  2:854712/1‑150 (MQ=255)
             cAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑‑tttttttt                                                                                                                                         >  1:536181/1‑150 (MQ=255)
              aGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑‑TTTTTTTTg                                                                                                                                        <  1:709187/150‑1 (MQ=255)
              aGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑‑TTTTTTTTg                                                                                                                                        <  1:26488/150‑1 (MQ=255)
               taCTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAATCGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑‑TTTTTTTTgg                                                                                                                                       <  1:197788/149‑1 (MQ=255)
               gACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTCCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑‑TTTTTTTTgg                                                                                                                                       <  2:208016/150‑1 (MQ=255)
               gACTAATGGTTTTAAAACTCTTGCTGAAGGGCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑‑TTTTTTTTgg                                                                                                                                       <  1:1478823/150‑1 (MQ=255)
                      ggTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑‑TTTTTTTTGGTCTTTAg                                                                                                                                <  2:573757/150‑1 (MQ=255)
                            aaaaCTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTCACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCg                                                                                                                          <  1:864252/150‑1 (MQ=255)
                            aaaaCTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCg                                                                                                                          <  1:1188282/150‑1 (MQ=255)
                              aaCTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGt                                                                                                                        <  1:244805/150‑1 (MQ=255)
                               aCTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTATGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGtt                                                                                                                       <  2:1623260/150‑1 (MQ=255)
                                ctctTGCTGAAGGTCAGCGCGAAGAGTTCGAAATCACTACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGTATTCAAAACCCGCTTAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGa                                                                                                                     <  1:239034/150‑1 (MQ=255)
                                atctTGCTTAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGAGCCAAAGGCCCTTCTTCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTg                                                                                                                      <  2:1490116/149‑1 (MQ=255)
                                      cTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCg                                                                                                                >  1:263614/1‑150 (MQ=255)
                                       tGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGa                                                                                                               >  1:815066/1‑150 (MQ=255)
                                        gAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGaa                                                                                                              >  1:914338/1‑150 (MQ=255)
                                         aaGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGaaa                                                                                                             <  1:744181/150‑1 (MQ=255)
                                          aGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGaaaa                                                                                                            <  1:1153626/150‑1 (MQ=255)
                                             tCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAg                                                                                                         <  2:813523/150‑1 (MQ=255)
                                               aGCGCTTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTATGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTATTCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGcc                                                                                                       <  2:1072550/150‑1 (MQ=255)
                                               aGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCAGCTTAATCAGCTGG‑‑TTTTTTTTGGTCTTTAGTGTGAGGTTGAGGCAGAAAACAGcc                                                                                                       >  1:1098702/1‑150 (MQ=255)
                                               aGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCGCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑‑ATTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGcc                                                                                                       >  1:1104825/1‑150 (MQ=255)
                                                 ggcgTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAg                                                                                                     <  1:428692/149‑1 (MQ=255)
                                                 cgcgTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAg                                                                                                     <  1:410902/150‑1 (MQ=255)
                                                  gcgTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCTCTCTGTAAGATACGTCAGCAAGATTTCAAAACCCGCTTAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCAGAAAACAGCCAGa                                                                                                    >  1:1623008/1‑150 (MQ=255)
                                                    gTAGAGTTCGAAATCACTAACGTTGCCACAGGCCCTTCTGCTGCAACCGTACTCGCTCTGAAAGATACGTCTGCAAGAATTCAAAACCCGCTTAATCTGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAAt                                                                                                  <  1:281064/150‑1 (MQ=255)
                                                      agagTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGc                                                                                                <  1:62330/150‑1 (MQ=255)
                                                         gTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAg                                                                                             <  1:79788/150‑1 (MQ=255)
                                                            cGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGc                                                                                          <  1:1570402/150‑1 (MQ=255)
                                                            cGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGAAAGAATTAAAAACCCGCTTTATCTGGGGGGTTTTTTTTTGGTCTTTTGTTGCCGGTTGAGGCCCAAAACCGCCAGAATGccgga                                                                                            >  2:85231/1‑147 (MQ=255)
                                                               aaTCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCAGt                                                                                       >  2:759566/1‑150 (MQ=255)
                                                                   aCTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATg                                                                                   <  2:827768/150‑1 (MQ=255)
                                                                    cTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATgg                                                                                  >  1:800890/1‑150 (MQ=255)
                                                                     tAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGc                                                                                 <  2:1378309/150‑1 (MQ=255)
                                                                     tAACGGTGCCAAACGCCCTTCTGCTGCAAACGTAATCGCTCTGAAAGATACGTCAGCAAGACTTCACAACCCGCTTACTCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGc                                                                                 <  2:113852/150‑1 (MQ=255)
                                                                         ggTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCaaaa                                                                             >  1:1716387/1‑150 (MQ=255)
                                                                          gTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCATCAAGAATTCAAAACCCGCTTAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAg                                                                            <  2:1273327/150‑1 (MQ=255)
                                                                          gTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAg                                                                            <  2:361724/150‑1 (MQ=255)
                                                                          gTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAg                                                                            <  2:206936/150‑1 (MQ=255)
                                                                            gCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGAc                                                                          >  1:314904/1‑150 (MQ=255)
                                                                              cAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGTTACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCTGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGAccc                                                                        <  1:1008765/150‑1 (MQ=255)
                                                                              cAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGAccc                                                                        >  2:480842/1‑150 (MQ=255)
                                                                                 aGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAg                                                                     >  2:64810/1‑150 (MQ=255)
                                                                                    cccTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAAGCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAg                                                                  <  2:263614/150‑1 (MQ=255)
                                                                                      cTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGt                                                                >  2:662402/1‑150 (MQ=255)
                                                                                        tCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTg                                                              <  1:607646/150‑1 (MQ=255)
                                                                                         ctgctgCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGa                                                             <  2:1495984/150‑1 (MQ=255)
                                                                                         ctgctgCAAACGTAATCGATCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTTGGGCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTg                                                              <  2:824267/150‑1 (MQ=255)
                                                                                           gctgcAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACg                                                           >  2:534063/1‑150 (MQ=255)
                                                                                            ctgcAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGa                                                          >  1:926145/1‑150 (MQ=255)
                                                                                                 aaCGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACg                                                     <  2:421948/150‑1 (MQ=255)
                                                                                                    gTAATCGCTTTGTAAGATACGTAAGCAAGAATTCAAAACCCGCTTAATCATCGTG‑‑TTTTTTTTGGTCTTTAGTTTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTc                                                  <  1:1016227/150‑1 (MQ=255)
                                                                                                    gTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTc                                                  <  2:1168782/150‑1 (MQ=255)
                                                                                                    gTAATCGCTCTGTAAGATACGTCAGAAAGATTTCAAAACCCGCTTAATCAGCGGG‑‑TTTTTTTTTGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTc                                                  <  2:342261/150‑1 (MQ=255)
                                                                                                         cGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTaa                                             >  2:536416/1‑150 (MQ=255)
                                                                                                          gCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAAt                                            >  1:458168/1‑150 (MQ=255)
                                                                                                           ctctGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATg                                           <  2:380430/150‑1 (MQ=255)
                                                                                                            tctGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGc                                          <  2:911372/150‑1 (MQ=255)
                                                                                                             ctGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑‑GTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAACGACCACAGAAGGTTGACGAAAACGTTCAGTAATGcc                                         >  1:1160711/1‑150 (MQ=255)
                                                                                                              tGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGccc                                        <  2:536181/150‑1 (MQ=255)
                                                                                                                 aaGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGc                                     >  1:1354120/1‑150 (MQ=255)
                                                                                                                   gATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCa                                   >  1:518568/1‑150 (MQ=255)
                                                                                                                     tACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGAGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTACAGTAATGCCCAGCCaaa                                 >  1:427383/1‑150 (MQ=255)
                                                                                                                      aCGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGACAAAAGACCCCAGAAGGTTGACGAACACGTTCAGTAATGCCCAGCCAAAg                                <  2:374875/150‑1 (MQ=255)
                                                                                                                       cGTCAGCAAGAATTCAAAACCCGCATAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGc                               <  1:674804/150‑1 (MQ=255)
                                                                                                                          cAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGc                                     <  1:339453/141‑1 (MQ=255)
                                                                                                                          cAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGc                                     >  2:339453/1‑141 (MQ=255)
                                                                                                                           aGCAATAATTAAAAACCCGCTTAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGGCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGcc                           <  1:854712/150‑1 (MQ=255)
                                                                                                                            gCAAGAATTCAAAACCCGCTTAATCAGCGGG‑‑TGTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAACGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCCGCCAAAGCGGccc                          >  1:1096635/1‑150 (MQ=255)
                                                                                                                             cAAGAATTCAAAACCCGCTTAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCct                         >  1:1303057/1‑150 (MQ=255)
                                                                                                                             cAAGAATTCAAAACCCGCTTAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCct                         >  1:1536121/1‑150 (MQ=255)
                                                                                                                               aGAATTCAAAACCCGCTTAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAACGCGGCCctct                       >  1:1331718/1‑150 (MQ=255)
                                                                                                                                 caaaaaaaaCCCGCTGATTAAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGAGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTg                     >  1:1707165/6‑150 (MQ=255)
                                                                                                                                 aaTTCAAAACCCGCTTAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTTATGCCCAGCCAAAGCGGCCCTCTTg                     <  2:1716387/150‑1 (MQ=255)
                                                                                                                                 aaTTCAAAACCCGCTTAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTg                     >  2:1138964/1‑150 (MQ=255)
                                                                                                                                   ttCAAAACCCGCTTAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTa                   <  2:1473255/150‑1 (MQ=255)
                                                                                                                                   ttCAAAACCCGCTTAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTa                   <  2:1102181/150‑1 (MQ=255)
                                                                                                                                    tCAAAACCCGCTTAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTaa                  >  1:1074838/1‑150 (MQ=255)
                                                                                                                                     cAAAACCCGCTTAATCAGCGGG‑‑TTTTTTTTGGGCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAAc                 <  1:1015121/150‑1 (MQ=255)
                                                                                                                                     aaaaaCCCGCTGATTAAGCGGG‑‑TATTTTTTGGTCTTTAGTGTACGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCACAGAAGGTTGACGAAAACGTTCAGTAATGCCAAGCCAAAGCGGCCCTCTTGTAAc                 >  1:845736/2‑150 (MQ=255)
                                                                                                                                      aaaaCCCGCTTAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACa                <  2:926145/150‑1 (MQ=255)
                                                                                                                                      aaaaCCCGCTTAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACa                <  2:694628/150‑1 (MQ=255)
                                                                                                                                        aaCCCGCTTAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACaaa              >  1:1640327/1‑150 (MQ=255)
                                                                                                                                         aCCCGCTTAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGTCCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCGAGCCAAAGCGGCCCTCTTGTAACaaaa             <  1:443426/150‑1 (MQ=255)
                                                                                                                                            cGCTTAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAaca          >  1:1326090/1‑150 (MQ=255)
                                                                                                                                             gCTTAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAacac         <  1:1700792/150‑1 (MQ=255)
                                                                                                                                             gCTTAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAacac         >  1:730389/1‑150 (MQ=255)
                                                                                                                                                tAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAAcaccac      >  1:925751/1‑150 (MQ=255)
                                                                                                                                                 aaTCAGCGGG‑‑TTTTTTTTGGTCTATAGTGTGCGGGTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAACACCACt     <  1:715591/150‑1 (MQ=255)
                                                                                                                                                    cAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAACACCACTTCt  >  1:396223/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                       
CTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAACACCACTTCT  >  NC_000913/657387‑657684

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: