Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 657,533 A→T 6.4% 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,5330AT 6.4% 300.1 / 5.8 105intergenic (+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 (43/54);  new base T (3/4);  total (46/59)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.29e-01

Only 100 of 106 total aligned reads displayed.

CACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAACACCACT  >  NC_000913/657385‑657681
                                                                                                                                                    |                                                                                                                                                      
cacTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAAt                                                                                                                                                       >  2:79788/1‑150 (MQ=255)
 acTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTAAAAACCCCCTTAAta                                                                                                                                                      >  2:244805/1‑149 (MQ=255)
  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                                                                                                                                         <  2:591899/150‑2 (MQ=255)
            aTCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGG‑‑ttttg                                                                                                                                         <  1:1679641/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:26488/150‑1 (MQ=255)
                aGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑‑TTTTTTTTg                                                                                                                                     <  1:709187/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)
                                        cTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCg                                                                                                             >  1:263614/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)
                                                    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                                                                         <  1:501977/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)
                                                                                                            gCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAAt                                         >  1:458168/1‑150 (MQ=255)
                                                                                                            gCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAACGTTGACGAAAACGTTCAGTAAt                                         >  1:1635501/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)
                                                                                                                                 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:694628/150‑1 (MQ=255)
                                                                                                                                        aaaaCCCGCTTAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACa             <  2:926145/150‑1 (MQ=255)
                                                                                                                                          aaCCCGCTTAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACaaa           >  1:1640327/1‑150 (MQ=255)
                                                                                                                                           aCCCGCTTAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGTCCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCGAGCCAAAGCGGCCCTCTTGTAACaaaa          <  1:443426/150‑1 (MQ=255)
                                                                                                                                            cccGCTGATTAAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACaaaaa         >  2:1003757/1‑150 (MQ=255)
                                                                                                                                              cGCTTAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAaca       >  1:1326090/1‑150 (MQ=255)
                                                                                                                                               gCTTAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAacac      >  1:730389/1‑150 (MQ=255)
                                                                                                                                               gCTTAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAacac      <  1:1700792/150‑1 (MQ=255)
                                                                                                                                                  tAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAAcaccac   >  1:925751/1‑150 (MQ=255)
                                                                                                                                                  tAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAAcaccac   >  2:4143/1‑150 (MQ=255)
                                                                                                                                                   aaTCAGCGGG‑‑TTTTTTTTGGTCTATAGTGTGCGGGTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAACACCACt  <  1:715591/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
CACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAACACCACT  >  NC_000913/657385‑657681

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: