Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 657,533 A→T 17.0% 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.0% 398.8 / 64.8 159intergenic (+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 (63/68);  new base T (12/15);  total (75/84)
Fisher's exact test for biased strand distribution p-value = 8.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

Only 100 of 162 total aligned reads displayed.

GCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAAC  >  NC_000913/657435‑657632
                                                                                                  |                                                                                                    
gcgcGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAAt                                                                                                     <  1:1611329/100‑1 (MQ=255)
gcgcGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTATGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAAt                                                                                                     >  1:2393540/1‑100 (MQ=255)
 cgcgTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATc                                                                                                    <  2:1611522/100‑1 (MQ=255)
 cgcgTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATc                                                                                                    <  1:1401149/100‑1 (MQ=255)
   cgTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCan                                                                                                  <  1:3840830/100‑2 (MQ=255)
    gTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCATTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGc                                                                                                 <  1:2171711/100‑1 (MQ=255)
     tAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCg                                                                                                <  1:1538879/100‑1 (MQ=255)
       gagTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCggg                                                                                              <  2:806815/100‑1 (MQ=255)
         gTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑tt                                                                                           >  2:807012/1‑100 (MQ=255)
         gTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGG‑tt                                                                                           <  1:1505475/100‑1 (MQ=255)
          ttCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGG‑ttt                                                                                          <  2:2004477/100‑1 (MQ=255)
          ttCGAAATCACTAACGGTGCCAAAGGCCCTACTGCTACAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGA‑ttt                                                                                          >  1:1466646/1‑100 (MQ=255)
           tCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGG‑tttt                                                                                         <  1:2736224/100‑1 (MQ=255)
           tCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGG‑tttt                                                                                         <  2:2389911/100‑1 (MQ=255)
           tCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGG‑tttt                                                                                         <  2:3839691/100‑1 (MQ=255)
            cGAAATCACTAACGGTGCCAAAGGCCCTTCTGTTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGG‑ttttg                                                                                        <  1:1684981/100‑2 (MQ=38)
            cGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑ttttt                                                                                        >  1:92806/1‑100 (MQ=255)
            cGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGG‑ttttg                                                                                        <  1:4121478/100‑2 (MQ=39)
            cGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGG‑ttttg                                                                                        <  1:4118960/100‑2 (MQ=39)
                aTCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTg                                                                                    <  2:1252809/100‑1 (MQ=255)
                aTCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTg                                                                                    <  2:205065/100‑1 (MQ=255)
                 tCACTAACGGTGCCAAAGGCCCTTCTGCTGCATACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTgg                                                                                   <  1:1736638/100‑1 (MQ=255)
                 tCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTgg                                                                                   <  1:4278067/100‑1 (MQ=255)
                 tCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTgg                                                                                   <  2:1934923/100‑1 (MQ=255)
                 tCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTgg                                                                                   <  2:636577/100‑1 (MQ=255)
                  cACTAACTGTGCCAAAGGCCCGTCTGCTTCAAACGTAATCGCGCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGt                                                                                  <  2:3017698/100‑1 (MQ=255)
                   aCTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTc                                                                                 >  1:2957179/1‑100 (MQ=255)
                    cTAACGGTGCCAAAGGCCCTTCTGCTGCAAGCGGAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTCTTTTTGGTCt                                                                                >  1:2897171/1‑100 (MQ=255)
                      aaCGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCttt                                                                              >  1:3421974/1‑100 (MQ=255)
                       aCGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTa                                                                             >  1:1745317/1‑100 (MQ=255)
                        cGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTTTTCttt                                                                              >  1:3304384/1‑98 (MQ=255)
                        cGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAg                                                                            <  2:2501271/100‑1 (MQ=255)
                        cGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAg                                                                            >  1:3559159/1‑100 (MQ=255)
                        cGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAg                                                                            <  1:3422776/100‑1 (MQ=255)
                          gTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAgtg                                                                          <  2:561306/100‑1 (MQ=255)
                           tGCCAAAGGCGCTTCTGCTGCAAACGTAATCGCTGTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAg                                                                            <  1:1404417/97‑1 (MQ=255)
                           tGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAgtgt                                                                         <  2:1876337/100‑1 (MQ=255)
                              cAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCg                                                                      >  1:3579704/1‑100 (MQ=255)
                              cAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCg                                                                      >  2:161996/1‑100 (MQ=255)
                              cAAAGGCCCTTCTCCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGtgtg                                                                        <  2:4352457/98‑1 (MQ=255)
                                   gCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGCTGa                                                                 >  2:3554570/1‑100 (MQ=255)
                                    cccTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAg                                                                >  1:3923959/1‑100 (MQ=255)
                                     ctttCTGTTTCAAACGTAATTTCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAgg                                                               <  2:3678693/98‑1 (MQ=255)
                                     ccTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAgg                                                               <  1:4041600/100‑1 (MQ=255)
                                     ccTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAgg                                                               <  2:1009677/100‑1 (MQ=255)
                                       ttCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAg                                                                >  1:1311619/1‑97 (MQ=255)
                                       ttCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGcc                                                             <  2:2957196/100‑1 (MQ=255)
                                         ctgctgCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGa                                                           >  1:2387044/1‑100 (MQ=255)
                                          tgctgcAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGa                                                           <  1:255264/100‑1 (MQ=255)
                                            ctgcAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGaaaa                                                        <  1:4231869/100‑1 (MQ=255)
                                            ctgcAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGaaaa                                                        >  2:3006768/1‑100 (MQ=255)
                                             tgcAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTTGTCTTTAGTGTGCGGTTGAGTCCGAAAAc                                                       >  2:653150/1‑100 (MQ=255)
                                             tgcAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAAc                                                       >  2:765231/1‑100 (MQ=255)
                                                aaaCGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAg                                                     <  2:192503/100‑1 (MQ=255)
                                                aaaCGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTAAGTGTGCGGTTGAGGCAGAAAACAGc                                                    >  1:3270562/1‑100 (MQ=255)
                                                     tAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGaa                                               >  2:1089673/1‑100 (MQ=255)
                                                     tAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGaa                                               >  1:129085/1‑100 (MQ=255)
                                                      aaTCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAAt                                              <  2:2521028/100‑1 (MQ=255)
                                                         cGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGTTCTTTAGTTTTCGGTTGAGGACGAAAACAGCCAGAAt                                              >  1:2587109/1‑97 (MQ=255)
                                                          gCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTTTGCGGTTGAGGCCGAAAACAGCCAGAATGCCa                                          >  1:3669917/1‑100 (MQ=255)
                                                          gCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCa                                          >  1:3914808/1‑100 (MQ=255)
                                                           ctctGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCa                                          >  1:120637/1‑99 (MQ=255)
                                                            tctGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATgagatc                                        >  1:1173676/1‑95 (MQ=255)
                                                             ctGTAATATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGC‑TTGTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAg                                         <  1:1130293/98‑1 (MQ=255)
                                                             ctGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAg                                         >  1:1519233/1‑98 (MQ=255)
                                                                 aaGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGt                                   >  2:2323843/1‑100 (MQ=255)
                                                                   gATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCa                                 >  2:4239550/1‑100 (MQ=255)
                                                                     tACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATg                               <  2:490480/100‑1 (MQ=255)
                                                                      aCGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATgg                              <  1:3476868/100‑1 (MQ=255)
                                                                       cGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGc                             >  2:2226200/1‑100 (MQ=255)
                                                                        ttcatcaAGAATTCAACACCCGCTTAATCAGCGAT‑TTTTTTTTGTGCTTTAGTTTGCGGTTGAGGACAAAAACACCCAGAATGCCAGTGCGGACATCGCa                            <  1:717417/99‑1 (MQ=255)
                                                                        gTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCa                            >  2:2518162/1‑100 (MQ=255)
                                                                          cAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCaaa                          <  2:120692/100‑1 (MQ=255)
                                                                            gCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAg                        >  2:2247595/1‑100 (MQ=255)
                                                                              aaGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGAc                      <  2:1671236/100‑1 (MQ=255)
                                                                              aaGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGAc                      <  2:1489711/100‑1 (MQ=255)
                                                                               aGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTTGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGa                       >  1:1034333/1‑98 (MQ=255)
                                                                                gAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGAccc                    <  2:3405690/100‑1 (MQ=255)
                                                                                ccaaaaaaaaCCCGCTGATTAAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGAccc                    >  2:1013481/7‑100 (MQ=38)
                                                                                ccaaaaaaaaCCCGCTGATTAAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGAccc                    >  1:586785/7‑100 (MQ=38)
                                                                                caaaaaaaaaCCCGCTGATTAAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTACGATCATGGCAAAAGAccc                    >  2:3671576/7‑100 (MQ=37)
                                                                                 caaaaaaaaCCCGCTGATTAAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGAcccc                   >  1:1805931/6‑100 (MQ=38)
                                                                                 aagtCAAATCCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGAcccc                   <  2:2284372/97‑1 (MQ=255)
                                                                                 aaTTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGAcccc                   <  2:3579729/100‑1 (MQ=255)
                                                                                  aaaaaaaaCCCGCTGATTAAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCa                  >  2:3882930/5‑100 (MQ=38)
                                                                                  aaaaaaaaCCCGCTGATTAAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCa                  >  2:2060659/5‑100 (MQ=38)
                                                                                  aTTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCa                  >  2:3267118/1‑100 (MQ=255)
                                                                                  aTTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCa                  <  1:192275/100‑1 (MQ=255)
                                                                                  aTTCAAAACCCGCTTAATCAGCGGG‑TATTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCa                  >  2:2677737/1‑100 (MQ=255)
                                                                                     cAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGaa               >  2:2754326/1‑100 (MQ=255)
                                                                                      aaaaCCCGCTGATTAAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAg              >  2:878131/1‑100 (MQ=255)
                                                                                       aaaCCCGCTGATTAAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAgg             >  2:4160419/1‑100 (MQ=255)
                                                                                        aaCCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGt            <  1:2226190/100‑1 (MQ=255)
                                                                                          cccGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTg          >  2:556759/1‑100 (MQ=255)
                                                                                           ccGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGa         <  2:3168914/100‑1 (MQ=255)
                                                                                            cGCTTAATCAGCGGG‑TTTTTTTTTTTTTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAACGTTGAc        >  1:3755628/1‑100 (MQ=255)
                                                                                              cTTAATCAGCGGG‑TTTTTTTTGTTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGa      >  2:3599437/1‑100 (MQ=255)
                                                                                              cTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGa      <  2:1690976/100‑1 (MQ=255)
                                                                                               ttAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGaa     <  1:3337937/100‑1 (MQ=255)
                                                                                                  aTCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACATCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAAc  >  1:1259573/1‑100 (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: ‑