Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 657,535 C→A 15.1% intergenic (+34/+20) cspE → / ← flc constitutive cold shock family transcription antitermination protein; negative regulator of cspA transcription; RNA melting protein; ssDNA‑binding protein/fluoride efflux channel, dual topology membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913657,5350CA15.1% 369.4 / 46.0 133intergenic (+34/+20)cspE/flcconstitutive cold shock family transcription antitermination protein; negative regulator of cspA transcription; RNA melting protein; ssDNA‑binding protein/fluoride efflux channel, dual topology membrane protein
Reads supporting (aligned to +/- strand):  ref base C (48/64);  new base A (7/13);  total (56/77)
Fisher's exact test for biased strand distribution p-value = 6.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.79e-01

Only 100 of 134 total aligned reads displayed.

GCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGT  >  NC_000913/657437‑657634
                                                                                                  |                                                                                                    
gcgTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCa                                                                                                     <  1:10460/100‑1 (MQ=255)
 cgTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAg                                                                                                    >  1:3192795/1‑100 (MQ=255)
 cgTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAg                                                                                                    >  2:1787562/1‑100 (MQ=255)
   tAGAGTTCGAAATCACTAACGGTGCCAAAGGCCATTGTGCAGCAAACGTAATGGCTCTGTAAGATACGTCAGCAAGAAAGCAAAACACGCTTAATCAGCg                                                                                                  >  2:669762/1‑100 (MQ=255)
   tAGAGTTCGAAATCACTAACGGTACCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGACTTCAAAACCCGCTTAATCAGCg                                                                                                  <  2:2432729/100‑1 (MQ=255)
      agTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGG‑t                                                                                              <  1:639386/100‑1 (MQ=255)
       gTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCggg                                                                                                >  2:4179144/1‑98 (MQ=255)
       gTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑tt                                                                                             >  2:92980/1‑100 (MQ=255)
        ttCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGG‑ttt                                                                                            <  1:3654542/100‑1 (MQ=255)
         tCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGG‑tttt                                                                                           <  2:1644466/100‑1 (MQ=255)
         tCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGG‑tttt                                                                                           <  2:2373247/100‑1 (MQ=255)
         tCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACCTCAGCAAGAATTCAAAACCCGCTGGTTAAGCGGG‑tttt                                                                                           <  2:2589014/100‑1 (MQ=255)
          cGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGG‑ttttg                                                                                          <  2:725382/100‑2 (MQ=39)
          cGAAATCACTAACGGTGCCAAAGGCCCTTCTCCTTCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGG‑ttttc                                                                                          <  2:700374/100‑2 (MQ=255)
           gAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑tttttt                                                                                         >  1:3724696/1‑100 (MQ=255)
            aaaTCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑ttttttt                                                                                        <  1:2174104/100‑1 (MQ=255)
            aaaTCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑ttttttt                                                                                        <  1:3170998/100‑1 (MQ=255)
             aaTCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑tttttttt                                                                                       <  2:4100992/100‑1 (MQ=255)
             aaTCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑tttttttt                                                                                       >  1:4123768/1‑100 (MQ=255)
             aaTCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAACGGG‑tttttttt                                                                                       <  1:2750948/100‑1 (MQ=255)
             aaTCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGG‑ttttgaat                                                                                       <  1:1953049/100‑5 (MQ=38)
              aTCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTg                                                                                      >  2:2490035/1‑100 (MQ=255)
              aTCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTg                                                                                      <  1:984868/100‑1 (MQ=255)
              aTCACTAACGGTGCCAAAGGCCCTGCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGG‑ttttgaat                                                                                       <  2:566285/99‑5 (MQ=38)
               tcttccgatctGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTgg                                                                                     <  2:3871195/90‑1 (MQ=255)
               tCACTAATGGTGCGAAAGGCGCTTGTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGAAAGAATTCAAAACCCGCTGATTAAGCGGG‑ttttgaat                                                                                       <  2:1539450/98‑5 (MQ=25)
               tCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTgg                                                                                     <  2:254482/100‑1 (MQ=255)
                 aCTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTc                                                                                   >  1:2427224/1‑100 (MQ=255)
                   tAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCtt                                                                                 <  1:1716786/100‑1 (MQ=255)
                   tAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCtt                                                                                 <  2:828047/100‑1 (MQ=255)
                    aaCGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCttt                                                                                >  1:1043054/1‑100 (MQ=255)
                        gTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAgtg                                                                            <  1:1451738/100‑1 (MQ=255)
                        gTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAgtg                                                                            <  1:2616295/100‑1 (MQ=255)
                          gCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGtgtg                                                                          >  1:2648320/1‑100 (MQ=255)
                                ggCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTg                                                                    >  1:27501/1‑100 (MQ=255)
                                ggCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCATAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTg                                                                    <  2:2576492/100‑1 (MQ=255)
                                 gCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGa                                                                   >  2:362060/1‑100 (MQ=255)
                                  cctttCTGCGGCAAACGTAACCGTTTTGTAAGTTACGCCAGAAAGAGTTCAAAACCCGCTTAATAAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAg                                                                  <  1:3258547/97‑1 (MQ=255)
                                  cccTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAg                                                                  <  1:337995/100‑1 (MQ=255)
                                   ccTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAgg                                                                 <  2:2991675/100‑1 (MQ=255)
                                    cttcttCTGCAACCGTATTGGCTCTGTAAGCTACGTGATCAAGAATTTAAGACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGttttacg                                                                 <  2:399631/99‑5 (MQ=255)
                                     ttCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGc                                                                <  2:1112630/100‑1 (MQ=255)
                                     ttCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGcc                                                               <  1:3271588/100‑1 (MQ=255)
                                     ttCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGcc                                                               <  1:1137651/100‑1 (MQ=255)
                                         gctgaAAACGGAATCCCTCTGTAAGATACGTCAGCAAGGATTCAAAACCGCCTTAGTCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGaaa                                                           <  2:285317/100‑1 (MQ=255)
                                           tgcAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAAc                                                         <  1:2884770/100‑1 (MQ=255)
                                            gcAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCATAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACa                                                        >  2:2295716/1‑100 (MQ=255)
                                            gcAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACa                                                        <  1:883345/100‑1 (MQ=255)
                                            gcAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACa                                                        >  1:4088922/1‑100 (MQ=255)
                                            gcAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACa                                                        <  1:2290320/100‑1 (MQ=255)
                                            gcAAACGTAATCGCACTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCAGG‑TTTTTTTATTTTTTTATTGTGCGGTTGAGGCAGaaaa                                                          >  2:1444479/1‑98 (MQ=255)
                                             cAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAg                                                       <  1:3132812/100‑1 (MQ=255)
                                               aaCGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACa                                                        >  2:2145472/1‑97 (MQ=255)
                                               aaCGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGcc                                                     >  1:2955824/1‑100 (MQ=255)
                                                aCGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCa                                                    >  1:3454281/1‑100 (MQ=255)
                                                aCGCAATCCCTCTGTAAGATACGTCACCAAGAGTTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTCGTCTATAGTGTGCGGTTGAGGCCCAAAACAGcc                                                     <  2:3540517/99‑1 (MQ=255)
                                                 cGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAg                                                   >  2:145218/1‑100 (MQ=255)
                                                    aaTCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAAt                                                >  2:708248/1‑100 (MQ=255)
                                                     aTCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATg                                               >  2:1624110/1‑100 (MQ=255)
                                                     aTCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATg                                               <  2:3373639/100‑1 (MQ=255)
                                                       cGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGcc                                             <  2:1258158/100‑1 (MQ=255)
                                                       cGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGcc                                             <  2:2465933/100‑1 (MQ=255)
                                                       cGCACTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGcc                                             <  2:220604/100‑1 (MQ=255)
                                                            tGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGc                                        <  2:4064920/100‑1 (MQ=255)
                                                              tAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCgg                                      <  1:341129/100‑1 (MQ=255)
                                                               aaGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCGAGAATGGCAGTGCgg                                      >  2:504690/1‑99 (MQ=255)
                                                               aaGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATACCAGTGCGGt                                     >  2:1735715/1‑100 (MQ=255)
                                                                 gATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCa                                   <  1:2139083/100‑1 (MQ=255)
                                                                 gATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCa                                   <  1:3833896/100‑1 (MQ=255)
                                                                    aCGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATgg                                <  1:3407430/100‑1 (MQ=255)
                                                                     cGTCAGCAAGAATTCAAAGCCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGc                               <  2:817918/100‑1 (MQ=255)
                                                                     cGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGc                               <  2:3214765/100‑1 (MQ=255)
                                                                      gTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCa                              >  1:3432684/1‑100 (MQ=255)
                                                                      gTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCa                              >  1:2107111/1‑100 (MQ=255)
                                                                         aGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCaaaa                           >  2:2790500/1‑100 (MQ=255)
                                                                          gCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTTTTTTTTACTGTGCGGTTTAGGCCGAATACATCCTGCATGCCAGTGCGCTCTTGGCaaaa                           >  2:2785177/1‑99 (MQ=255)
                                                                          gCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAg                          <  2:1382761/100‑1 (MQ=255)
                                                                          gCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAg                          <  1:3133320/100‑1 (MQ=255)
                                                                          gCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAg                          <  1:1014435/100‑1 (MQ=255)
                                                                           cAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGa                         >  1:2971267/1‑100 (MQ=255)
                                                                           cAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGa                         <  2:1795264/100‑1 (MQ=255)
                                                                              ccaaaaaaaaCCCGCTGATTAAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGAccc                      >  1:2759386/7‑100 (MQ=37)
                                                                               caaaaaaaaCCCGCTGATTAAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGAcccc                     >  1:3547562/6‑100 (MQ=38)
                                                                                aaaaaaaaCCCGCTGATTAAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCa                    >  2:3472011/5‑100 (MQ=38)
                                                                                aaaaaaaaCCCGCTGATTAAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCa                    >  1:367503/5‑100 (MQ=38)
                                                                                aTTCAAAACCCGCTTGATGAGCGGG‑TGTTATTAGGTCTTTAGCGTGCGGTTGAGGGGGAAAACAGCCAAAATGCAAGAGAGTTAATGGGAAAAGACCCCa                    >  2:2105620/1‑100 (MQ=255)
                                                                                aTTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCa                    >  1:2367643/1‑100 (MQ=255)
                                                                                  tCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGa                  >  1:156314/1‑100 (MQ=255)
                                                                                   cAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGaa                 >  1:572420/1‑100 (MQ=255)
                                                                                   cAAAACCCGCTTAAGCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGaa                 <  2:796342/100‑1 (MQ=255)
                                                                                    aaaaCCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAg                <  1:3670263/100‑1 (MQ=255)
                                                                                     aaaCCCGCTGATTAAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAgg               >  1:899223/1‑100 (MQ=255)
                                                                                      aaCCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTCTGCGGTTGAGGCCGAAAACCACCACAATTCCAGTGCGGTCATGGCCAAACACCCCAGAAgg               >  2:643676/1‑99 (MQ=255)
                                                                                        cccGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGtt             <  2:372244/99‑1 (MQ=255)
                                                                                        cccGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGtt             <  2:911969/99‑1 (MQ=255)
                                                                                        cccGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTg            >  1:1937878/1‑100 (MQ=255)
                                                                                         ccGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGa           <  1:2619812/100‑1 (MQ=255)
                                                                                          cGCTTAATCAGCGGG‑TTGTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGAc          <  2:3160821/100‑1 (MQ=255)
                                                                                           gCTTAGTCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAAGGCCAGTGCGGTCTTGGCAAAAGACCCCAGAAGGTTGACg         <  2:1861959/100‑1 (MQ=255)
                                                                                                  cAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGt  <  2:4130134/100‑1 (MQ=255)
                                                                                                  |                                                                                                    
GCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGT  >  NC_000913/657437‑657634

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