Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 657,533 A→T 17.8% 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,5330AT17.8% 320.5 / 64.7 141intergenic (+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 (47/68);  new base T (9/16);  total (57/84)
Fisher's exact test for biased strand distribution p-value = 8.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.63e-02

Only 100 of 144 total aligned reads displayed.

GCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAAC  >  NC_000913/657435‑657632
                                                                                                  |                                                                                                    
gcgcGTAGAGTTCGAAATCACTAACGGTGCCAGAGGCCCTTCTGCTGCCAACGTAATCGCTCTGTAAGATACGTCAGCAAGTATTCAAAACCCGCTTAAt                                                                                                     >  1:1714532/1‑100 (MQ=255)
gcgcGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAAt                                                                                                     >  1:1022346/1‑100 (MQ=255)
 cgcgTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATc                                                                                                    <  2:644760/100‑1 (MQ=255)
 cgcgTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATc                                                                                                    >  2:3080561/1‑100 (MQ=255)
    gTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGc                                                                                                 <  1:1682542/100‑1 (MQ=255)
     tAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCg                                                                                                <  2:1466737/100‑1 (MQ=255)
     tAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCg                                                                                                <  1:127844/100‑1 (MQ=255)
     tAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCg                                                                                                <  2:1159038/100‑1 (MQ=255)
     tAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGATGCAAACGTAAACGCGCTGTAAGATACGTCAGTAAGAATTCAAAACTCGCTAAATCa                                                                                                   >  2:274283/1‑97 (MQ=255)
         gTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑tt                                                                                           <  1:10279/100‑1 (MQ=255)
         gTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑tt                                                                                           >  1:683554/1‑100 (MQ=255)
         gTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑tt                                                                                           >  1:3456325/1‑100 (MQ=255)
           tCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑tttt                                                                                         <  1:3799591/100‑1 (MQ=255)
           tCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGG‑tttt                                                                                         <  2:900531/100‑1 (MQ=255)
           tCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGG‑tttt                                                                                         <  1:797300/100‑1 (MQ=255)
           tCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGG‑tttt                                                                                         <  1:308396/100‑1 (MQ=255)
            cGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGNG‑ttttt                                                                                        >  1:3419226/1‑100 (MQ=255)
            cGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGG‑ttttg                                                                                        <  1:3910050/100‑2 (MQ=39)
            cGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCAATCTGTATGCTCCGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑ttttt                                                                                        >  2:1492627/1‑100 (MQ=255)
            cGAAATCACTAACGGTACCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGG‑ttttg                                                                                        <  2:1293572/100‑2 (MQ=38)
             gAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑tttttt                                                                                       <  1:1826265/100‑1 (MQ=255)
             gAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑ttattt                                                                                       >  2:2576607/1‑100 (MQ=255)
             gAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGG‑ttttga                                                                                       <  1:138547/100‑3 (MQ=38)
              aaaTCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑ttttttt                                                                                      <  1:3080581/100‑1 (MQ=255)
               aaTCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAGAACCCGCTTAATCAGCGGG‑tttttttt                                                                                     <  1:231812/100‑1 (MQ=255)
               aaTCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGG‑ttttgaat                                                                                     <  1:144876/100‑5 (MQ=38)
               aaTCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGG‑ttttgaat                                                                                     <  2:2007093/100‑5 (MQ=38)
                 tCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTgg                                                                                   <  1:2736469/100‑1 (MQ=255)
                 tCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTgg                                                                                   <  2:3942091/100‑1 (MQ=255)
                 tCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTgg                                                                                   <  2:2363236/100‑1 (MQ=255)
                 tCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTgg                                                                                   <  1:454296/100‑1 (MQ=255)
                   aCTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTc                                                                                 >  2:3048953/1‑100 (MQ=255)
                    atAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTACTCGCTCTGTAAGATACATCGGCCAGAACTCAGAACCCGCTTACTCAGCGGG‑TATTTTTCGGt                                                                                  >  2:708266/2‑98 (MQ=255)
                      aaCGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCttt                                                                              <  2:1415138/100‑1 (MQ=255)
                      aaCGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCttt                                                                              <  2:2252758/100‑1 (MQ=255)
                        cGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAg                                                                            <  1:3081155/100‑1 (MQ=255)
                          gTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAg                                                                            <  2:231179/98‑1 (MQ=255)
                            gCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGtgtg                                                                        <  2:321981/100‑1 (MQ=255)
                            gCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGtgtg                                                                        >  1:3832106/1‑100 (MQ=255)
                               aaaGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCgg                                                                     <  2:3700195/100‑1 (MQ=255)
                                aaGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGt                                                                    >  2:1718229/1‑100 (MQ=255)
                                 aGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGtt                                                                   <  2:348410/100‑1 (MQ=255)
                                  ggCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTg                                                                  <  1:1092857/100‑1 (MQ=255)
                                   gCCCTTCTGCTTCAAACGTAGTCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGt                                                                    <  2:2685327/97‑1 (MQ=255)
                                     ccTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAgg                                                               <  2:1430036/100‑1 (MQ=255)
                                      cTTCTGCTGCAAACGTAATTGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGc                                                              <  1:405086/100‑1 (MQ=255)
                                      cTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGc                                                              >  1:2344430/1‑100 (MQ=255)
                                      cTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGc                                                              >  1:3313245/1‑100 (MQ=255)
                                       ttCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGcc                                                             <  1:3604342/100‑1 (MQ=255)
                                            ctgcAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGaaaa                                                        <  2:381899/100‑1 (MQ=255)
                                            ctgcAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGaaa                                                         >  2:2711216/1‑99 (MQ=255)
                                               cAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGTG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACa                                                      <  2:3759926/99‑1 (MQ=255)
                                               cAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACa                                                      <  2:2472110/99‑1 (MQ=255)
                                                  aCGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCa                                                  <  2:104125/100‑1 (MQ=255)
                                                   cGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAg                                                 >  1:1001416/1‑100 (MQ=255)
                                                    gTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTATGGTCTTTAGTGTGCGGTTGAGGCTGAAAACAGCCAGa                                                <  2:2434401/100‑1 (MQ=255)
                                                     tAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGaa                                               >  1:3863623/1‑100 (MQ=255)
                                                           ctctGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAg                                         <  1:1672482/100‑1 (MQ=255)
                                                             ctGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTg                                       <  1:257284/100‑1 (MQ=255)
                                                             ctGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTg                                       >  1:1249082/1‑100 (MQ=255)
                                                               gTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTTTTTAGTGTGCGATTGAGGCCGAAAACAGCCAGAATGCCAGTGCg                                     >  2:887472/1‑100 (MQ=255)
                                                               gTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCg                                     >  1:2507299/1‑100 (MQ=255)
                                                                 aaGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGt                                   >  1:1409331/1‑100 (MQ=255)
                                                                 aaGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGt                                   <  2:2226118/100‑1 (MQ=255)
                                                                     tACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATg                               <  2:1016021/100‑1 (MQ=255)
                                                                     tACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATg                               <  1:710980/100‑1 (MQ=255)
                                                                       cGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATgg                              <  2:1079059/99‑1 (MQ=255)
                                                                       cGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGc                             <  1:3310506/100‑1 (MQ=255)
                                                                          cAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCaaa                          >  1:370833/1‑100 (MQ=255)
                                                                           aGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCaaaa                         >  1:2703600/1‑100 (MQ=255)
                                                                            gCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTTAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAg                        >  2:2921740/1‑100 (MQ=255)
                                                                            gCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAg                        >  1:3231258/1‑100 (MQ=255)
                                                                            gCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAg                        <  1:2251186/100‑1 (MQ=255)
                                                                            gCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAg                        <  2:1526719/100‑1 (MQ=255)
                                                                              aaGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGa                       >  1:121333/1‑100 (MQ=255)
                                                                              aaGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGAc                      <  2:3843642/100‑1 (MQ=255)
                                                                               aGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGAcc                     >  1:643805/1‑100 (MQ=255)
                                                                                 aaTTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGAcccc                   <  1:1044127/100‑1 (MQ=255)
                                                                                 aaTTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGAcccc                   >  2:2735151/1‑100 (MQ=255)
                                                                                  aaaaaaaaCCCGCTGATTAAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCATAATGCCAGGGCGGTCATGGCAAAAGACCCCa                  >  2:2657973/5‑100 (MQ=37)
                                                                                  aaaaaaaaCCCGCTGATTAAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCa                  >  1:1322015/5‑100 (MQ=38)
                                                                                  aaaaaaaaCCCGCTGATTAAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCa                  >  2:2978760/5‑100 (MQ=38)
                                                                                  aaaaaaaaCCCGCTGATTAAGCGGG‑TTTTGTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCa                  >  2:405052/5‑100 (MQ=38)
                                                                                  aTTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCa                  >  1:3364551/1‑100 (MQ=255)
                                                                                   aaaaaaaCCCGCTGATTAAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAg                 >  1:324261/4‑100 (MQ=38)
                                                                                     cagaccccGGTTGATCTGCGGG‑TGTTTTTTGGTCTTTACTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGc                             <  2:3537188/81‑1 (MQ=255)
                                                                                      aaaaCCCGCTGATTAAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAg              >  1:711891/1‑100 (MQ=255)
                                                                                       aaaCCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGAGGTTGAGGGAGAAAACAGTCAGAATTACAGTGCGGTCGTGGCAAAGTACCCCAg                 >  2:3389706/1‑96 (MQ=255)
                                                                                        aaCCCGCTTAATCAGCGGG‑TTTGTTTTGGTCTTTAGGGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGt            <  2:1274152/100‑1 (MQ=255)
                                                                                        aaCCCGCTGATTAAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCGAGTGCGGTCATGGCAAAAGACCCCAGa                >  2:3330221/1‑96 (MQ=255)
                                                                                         aCCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTCAGGCCGAAAGCAGCCAGAATGCCAGTGCGGTCATGGCGAAAGACCCCa                  >  2:2303551/1‑93 (MQ=255)
                                                                                          cccGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTg          <  2:1695921/100‑1 (MQ=255)
                                                                                           ccGCTTAATCAGCGGG‑GTTTTTTTGGGCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGt            <  2:2358724/97‑1 (MQ=255)
                                                                                             gCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACg       <  1:1718238/100‑1 (MQ=255)
                                                                                             gCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACg       >  1:3430832/1‑100 (MQ=255)
                                                                                             gCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACg       <  1:2191023/100‑1 (MQ=255)
                                                                                                tAATCAGCGGG‑TTTTTTTTTTTTTTTTTTTTGCGGTTTAGGCCGAAAACATCCAGAATGCTTGTGCGGGCATGTCAAAAGACCCTAGAAGGTTGACGaaa    >  2:2445211/1‑100 (MQ=255)
                                                                                                 aaTCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGaaaa   <  1:3353880/100‑1 (MQ=255)
                                                                                                  aTCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAAc  >  1:3847284/1‑100 (MQ=255)
                                                                                                  aTCAGCGGG‑TTTTTTTTAGTCTTTAGTGTGAGGTTGAGGACGCAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGATTGACGaaaa   >  2:369100/1‑99 (MQ=255)
                                                                                                  |                                                                                                    
GCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAAC  >  NC_000913/657435‑657632

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑