Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 657,531 T→G 5.9% 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,5310TG 5.9% 277.4 / 7.6 103intergenic (+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 (41/56);  new base G (3/3);  total (44/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 = 3.19e-01

Only 100 of 107 total aligned reads displayed.

GTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAACACCAC  >  NC_000913/657382‑657680
                                                                                                                                                     |                                                                                                                                                       
gTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCtt                                                                                                                                                         <  1:415152/150‑1 (MQ=255)
   cacTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAAt                                                                                                                                                      >  2:79788/1‑150 (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: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)
                                     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)
                                                              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                                                                        <  1:501977/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)
                                                                                                               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)
                                                                                                                        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                                 >  2:339453/1‑141 (MQ=255)
                                                                                                                               cAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGc                                 <  1:339453/141‑1 (MQ=255)
                                                                                                                                aGCAATAATTAAAAACCCGCTTAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGGCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGcc                       <  1:854712/150‑1 (MQ=255)
                                                                                                                                 gCAAGAATTCAAAACCCGCTTAATCAGCGGG‑‑TGTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAACGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCCGCCAAAGCGGccc                      >  1:1096635/1‑150 (MQ=255)
                                                                                                                                  cAAGAATTCAAAACCCGCTTAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCct                     >  1:1536121/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)
                                                                                                                                                 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  >  2:4143/1‑150 (MQ=255)
                                                                                                                                                     tAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAAcaccac  >  1:925751/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                       
GTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAACACCAC  >  NC_000913/657382‑657680

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: