Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 657,535 C→A 18.7% 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,5350CA18.7% 347.1 / 70.3 139intergenic (+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 (49/64);  new base A (5/21);  total (54/85)
Fisher's exact test for biased strand distribution p-value = 2.61e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.44e-01

Only 100 of 140 total aligned reads displayed.

CGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAAC  >  NC_000913/657436‑657632
                                                                                                   |                                                                                                 
ctcgTAGAGTTCTAAATCATTAACAGTGCCAAAGTCCATTCTGCTGCAAACGTAATCGCTCTGTAAAACACGTCAGCAAGAATTCAAAACCCGCTTAATc                                                                                                   <  1:647798/98‑1 (MQ=255)
  cgTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAg                                                                                                 <  1:1790227/100‑1 (MQ=255)
   gTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGc                                                                                                >  2:3530068/1‑100 (MQ=255)
      gagTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCggg                                                                                             <  1:2648339/100‑1 (MQ=255)
      gagTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCggg                                                                                             <  1:2026934/100‑1 (MQ=255)
      gagTTCAAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCggg                                                                                             >  2:3505727/1‑100 (MQ=255)
       agTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGGt                                                                                            <  2:3575674/100‑1 (MQ=255)
        gTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGtt                                                                                           >  2:549176/1‑100 (MQ=255)
        gTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGtt                                                                                           >  1:3247116/1‑100 (MQ=255)
        gTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGtt                                                                                           >  2:2132664/1‑100 (MQ=255)
        gTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGt                                                                                            >  1:1275216/1‑99 (MQ=255)
         ttCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGGttt                                                                                          <  2:715704/100‑1 (MQ=255)
         ttCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGGttt                                                                                          <  1:2909926/100‑1 (MQ=255)
          tCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGGtttt                                                                                         <  1:3811660/100‑1 (MQ=255)
          tCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGGtttt                                                                                         <  2:1059452/100‑1 (MQ=255)
           cGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGGttttg                                                                                        <  1:1621246/100‑2 (MQ=39)
           cGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGGttttg                                                                                        <  2:3627918/100‑2 (MQ=39)
            gAAATCACTAACGGTGCCAAAGGCCCTTCTGCTTCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGGttttga                                                                                       <  1:87426/100‑3 (MQ=38)
            gAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGGttttga                                                                                       <  2:2714487/100‑3 (MQ=38)
             aaaaCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACATCAGCAAGAATTCAAAACCCGCTGATTAAGCGGGttttgaa                                                                                      <  1:1199033/100‑4 (MQ=38)
             aaaTCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGGttttgaa                                                                                      <  2:2250295/100‑4 (MQ=38)
             aaaTCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGGttttgaa                                                                                      <  2:665183/100‑4 (MQ=39)
              aaTCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGtttttttt                                                                                     <  1:3003211/100‑1 (MQ=255)
              aaTCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGtttttttt                                                                                     <  1:1805708/100‑1 (MQ=255)
              aaTCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGtttttttt                                                                                     <  2:1629091/100‑1 (MQ=255)
              aaTCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGGttttgaat                                                                                     <  1:1537282/100‑5 (MQ=38)
              aaTCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGGttttgaat                                                                                     <  1:3423952/100‑5 (MQ=38)
              aaTCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGGttttgaat                                                                                     <  2:2195036/100‑5 (MQ=38)
               aTCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTg                                                                                    <  1:2867589/100‑1 (MQ=255)
                tCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTgg                                                                                   <  2:1275223/100‑1 (MQ=255)
                tCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGAATCAGCGGGTTTTTTTTgg                                                                                   <  1:1050419/100‑1 (MQ=255)
                 cACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGt                                                                                  >  1:271663/1‑100 (MQ=255)
                  aCTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTc                                                                                 >  2:1395589/1‑100 (MQ=255)
                  aCTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTc                                                                                 <  2:1472983/100‑1 (MQ=255)
                  aCTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGCTTTTTTTGGTc                                                                                 >  2:2853132/1‑100 (MQ=255)
                   cTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCt                                                                                >  1:1008449/1‑100 (MQ=255)
                   cTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCt                                                                                <  1:786858/100‑1 (MQ=255)
                     aaCGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCttt                                                                              <  1:476306/100‑1 (MQ=255)
                     aaCGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCttt                                                                              <  2:1398764/100‑1 (MQ=255)
                     aaCGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCttt                                                                              <  1:3267193/100‑1 (MQ=255)
                      aCGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTa                                                                             <  1:405539/100‑1 (MQ=255)
                           gCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGtgtg                                                                        <  2:2093273/100‑1 (MQ=255)
                           gCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGtgtg                                                                        >  2:866190/1‑100 (MQ=255)
                            ccAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGc                                                                       <  2:719578/100‑1 (MQ=255)
                             cAAAGGCCCTTCTGATGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGtttttttttttttt                                                                               >  1:2247349/1‑85 (MQ=255)
                               aaGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGt                                                                    <  1:2885512/100‑1 (MQ=255)
                                 ggCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTg                                                                  <  2:2470001/100‑1 (MQ=255)
                                 ggCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTg                                                                  <  2:398391/100‑1 (MQ=255)
                                   cccTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAg                                                                >  1:923040/1‑100 (MQ=255)
                                   cccTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAg                                                                <  1:836563/100‑1 (MQ=255)
                                   cccTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAg                                                                <  2:3090889/100‑1 (MQ=255)
                                    ccTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAgg                                                               <  1:1296404/100‑1 (MQ=255)
                                    ccTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAgg                                                               <  1:1296363/100‑1 (MQ=255)
                                     cTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGa                                                                 <  1:1317512/97‑1 (MQ=255)
                                          gctgcAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGaaa                                                         >  2:561706/1‑100 (MQ=255)
                                          gctgcAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGaaa                                                         >  1:3588324/1‑100 (MQ=255)
                                          gctgcAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGaaa                                                         <  1:561699/100‑1 (MQ=255)
                                            tgcAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAAc                                                       >  2:938896/1‑100 (MQ=255)
                                            tgcAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAAc                                                       >  2:3080887/1‑100 (MQ=255)
                                             gcAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACa                                                      <  2:3365818/100‑1 (MQ=255)
                                              cAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAg                                                     <  2:3809176/100‑1 (MQ=255)
                                               aaaCGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACACGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGc                                                    >  1:2440114/1‑100 (MQ=255)
                                                aaCGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGcc                                                   <  1:468011/100‑1 (MQ=255)
                                                   gTAATCGCTCTGTAAGATATGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGa                                                <  1:89355/100‑1 (MQ=255)
                                                        cGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGcc                                           <  1:866183/100‑1 (MQ=255)
                                                          ctctGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTGTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAg                                         <  1:1859508/100‑1 (MQ=255)
                                                          ctctGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAg                                         <  2:633417/100‑1 (MQ=255)
                                                           tctGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTGTTTGTTGGTCGTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGt                                        >  2:2075429/1‑100 (MQ=255)
                                                             tGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGc                                      <  2:2491991/100‑1 (MQ=255)
                                                             tGTAAGATACGTCAGCAAGAATTAAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTAAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGc                                      >  2:2685182/1‑100 (MQ=255)
                                                              gTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCg                                     >  2:2424131/1‑100 (MQ=255)
                                                               tAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCgg                                    <  1:74562/100‑1 (MQ=255)
                                                                aaGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGt                                   >  1:3594549/1‑100 (MQ=255)
                                                                      cGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGc                             >  2:1525833/1‑100 (MQ=255)
                                                                       gatctcAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCa                            <  1:608439/95‑1 (MQ=255)
                                                                       gTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCa                            >  1:3055137/1‑100 (MQ=255)
                                                                        tCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCaa                           >  2:1357887/1‑100 (MQ=255)
                                                                        tCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCaa                           <  1:450811/100‑1 (MQ=255)
                                                                          aGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGTTCTTTAGTGTGCGGTTGAAGCCGAAAACAGCCAGCATGGCAGTGCGGTCATGGCacaa                         >  1:1406659/1‑100 (MQ=255)
                                                                           gCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAg                        <  2:2618439/100‑1 (MQ=255)
                                                                           gCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAg                        <  2:364616/100‑1 (MQ=255)
                                                                           gCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAg                        >  2:1641018/1‑100 (MQ=255)
                                                                            cAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCaagatc                       >  2:608447/1‑96 (MQ=255)
                                                                              aGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGAcc                     >  2:934993/1‑100 (MQ=255)
                                                                               gAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGAccc                    <  2:1498635/100‑1 (MQ=255)
                                                                                aaTTCAAAACCCGCTTAATCAGCGGGTTTTTTTTTTTCTTTTTTTTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGAcccc                   >  1:2140340/1‑100 (MQ=255)
                                                                                aaTTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGAcccc                   <  1:1027332/100‑1 (MQ=255)
                                                                                 aTTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCa                  <  2:3588346/100‑1 (MQ=255)
                                                                                 aTTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCa                  <  2:3543233/100‑1 (MQ=255)
                                                                                  aaaaaaaCCCGCTGATTAAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAg                 >  2:589705/4‑100 (MQ=38)
                                                                                   aaaaaaCCCGCTGATTAAGCGGGTTTTTTTTGGTCTTTAGTGTTCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACACCAGa                >  1:3042732/3‑100 (MQ=38)
                                                                                    cAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGaa               >  2:3313705/1‑100 (MQ=255)
                                                                                     aaaaCCCGCTGATTAAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAg              >  2:3275517/1‑100 (MQ=255)
                                                                                        aCCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGt            <  2:500262/99‑1 (MQ=255)
                                                                                        aCCCGCTGATTAAGCGGGTTTTTTTTGGTTTTTAGTGCTCGGTTGAGGCCGAAAACAGCAAGAATGTCAGTGCGGTCATGGCCCAAGACCCCAGaaag             >  1:408990/1‑96 (MQ=255)
                                                                                         cccGCTTAATCAGCGGGTTTTTTTTTTTTTTTTTTTTTCTGTTGAGGCCGAAATCAGCCAGAATGCCATTGCGGTCATGGCAAAAGACCCCAGAAGGtt           >  1:2720117/1‑99 (MQ=255)
                                                                                          ccGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGa         <  2:387621/100‑1 (MQ=255)
                                                                                           cGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGAc        <  1:3817368/100‑1 (MQ=255)
                                                                                           cGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGAc        <  1:1463186/100‑1 (MQ=255)
                                                                                                 aTCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTAATGGCAAAAGACCCCAGAAGGTTGACGAAAAc  >  2:2382558/1‑100 (MQ=255)
                                                                                                   |                                                                                                 
CGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAAC  >  NC_000913/657436‑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: