Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 657,531 T→G 19.0% intergenic (+30/+24) cspE → / ← flc constitutive cold shock family transcription antitermination protein; negative regulator of cspA transcription; RNA melting protein; ssDNA‑binding protein/fluoride efflux channel, dual topology membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913657,5310TG19.0% 309.9 / 77.0 142intergenic (+30/+24)cspE/flcconstitutive cold shock family transcription antitermination protein; negative regulator of cspA transcription; RNA melting protein; ssDNA‑binding protein/fluoride efflux channel, dual topology membrane protein
Reads supporting (aligned to +/- strand):  ref base T (49/66);  new base G (5/22);  total (54/88)
Fisher's exact test for biased strand distribution p-value = 2.67e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.99e-01

Only 100 of 143 total aligned reads displayed.

TCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGAC  >  NC_000913/657432‑657626
                                                                                                   |                                                                                               
   gcgcGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAAt                                                                                              >  1:2261709/1‑100 (MQ=255)
   gcgcGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAAt                                                                                              <  2:215933/100‑1 (MQ=255)
    ctcgTAGAGTTCTAAATCATTAACAGTGCCAAAGTCCATTCTGCTGCAAACGTAATCGCTCTGTAAAACACGTCAGCAAGAATTCAAAACCCGCTTAATc                                                                                             <  1:647798/98‑1 (MQ=255)
      cgTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAg                                                                                           <  1:1790227/100‑1 (MQ=255)
          gagTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCggg                                                                                       <  1:2026934/100‑1 (MQ=255)
          gagTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCggg                                                                                       >  1:1793256/1‑100 (MQ=255)
          gagTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCggg                                                                                       <  1:2648339/100‑1 (MQ=255)
          gagTTCAAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCggg                                                                                       >  2:3505727/1‑100 (MQ=255)
            gTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGtt                                                                                     >  1:3247116/1‑100 (MQ=255)
            gTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGtt                                                                                     >  2:2933175/1‑100 (MQ=255)
            gTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGtt                                                                                     >  2:549176/1‑100 (MQ=255)
            gTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGt                                                                                      >  1:1275216/1‑99 (MQ=255)
            gTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCCCTGTAAGCTACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCggg                                                                                       >  1:460098/1‑98 (MQ=255)
             ttCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGGttt                                                                                    <  1:2909926/100‑1 (MQ=255)
              tCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGGtttt                                                                                   <  1:3811660/100‑1 (MQ=255)
              tCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGGtttt                                                                                   <  2:3265742/100‑1 (MQ=255)
               cGAAATCACTAACGGTGCCAAAGGTCCTTCTGCTCCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGGttttg                                                                                  <  2:2632871/100‑2 (MQ=38)
               cGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGttttt                                                                                  >  2:3715936/1‑100 (MQ=255)
               cGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGGttttg                                                                                  <  1:1621246/100‑2 (MQ=39)
                gAAATCACTAACGGTGCCAAAGGCCCTTCTGCTTCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGGttttga                                                                                 <  1:87426/100‑3 (MQ=38)
                gAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGtttttt                                                                                 >  1:3699965/1‑100 (MQ=255)
                gAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGGttttga                                                                                 <  2:2714487/100‑3 (MQ=38)
                 aaaaCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACATCAGCAAGAATTCAAAACCCGCTGATTAAGCGGGttttgaa                                                                                <  1:1199033/100‑4 (MQ=38)
                 aaaTCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGGttttgaa                                                                                <  2:665183/100‑4 (MQ=39)
                  aaTCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGtttttttt                                                                               <  2:1629091/100‑1 (MQ=255)
                  aaTCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGtttttttt                                                                               <  1:3003211/100‑1 (MQ=255)
                  aaTCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGGttttgaat                                                                               <  1:3423952/100‑5 (MQ=38)
                  aaTCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGGttttgaat                                                                               <  2:2195036/100‑5 (MQ=38)
                  aaTCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGGttttgaat                                                                               <  1:1537282/100‑5 (MQ=38)
                  aaTCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGGttttgaat                                                                               <  2:2190544/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)
                       cTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCt                                                                          >  1:1008449/1‑100 (MQ=255)
                         aaCGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCttt                                                                        <  1:3267193/100‑1 (MQ=255)
                         aaCGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCttt                                                                        <  2:1398764/100‑1 (MQ=255)
                          aCGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTa                                                                       <  1:405539/100‑1 (MQ=255)
                             gTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAgtg                                                                    <  1:486860/100‑1 (MQ=255)
                               gCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGtgtg                                                                  <  2:2093273/100‑1 (MQ=255)
                               gCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGtgtg                                                                  >  2:866190/1‑100 (MQ=255)
                               gCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGtgtg                                                                  >  2:3599707/1‑100 (MQ=255)
                                ccAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGc                                                                 <  2:719578/100‑1 (MQ=255)
                                 cAAAGGCCCTTCTGATGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGtttttttttttttt                                                                         >  1:2247349/1‑85 (MQ=255)
                                  aaaGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCgg                                                               <  1:810635/100‑1 (MQ=255)
                                   aaGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGt                                                              <  1:2885512/100‑1 (MQ=255)
                                     ggCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTg                                                            <  2:398391/100‑1 (MQ=255)
                                      gCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGa                                                           <  2:2476063/100‑1 (MQ=255)
                                       cccTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAg                                                          >  1:923040/1‑100 (MQ=255)
                                       cccTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAg                                                          <  2:623006/100‑1 (MQ=255)
                                       cccTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAg                                                          <  1:836563/100‑1 (MQ=255)
                                        ccTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAgg                                                         <  1:1296404/100‑1 (MQ=255)
                                        ccTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAgg                                                         <  1:1296363/100‑1 (MQ=255)
                                          ttCTGCTGCAAACGTAACCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGcc                                                       >  1:731126/1‑100 (MQ=255)
                                              gctgcAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGaaa                                                   >  1:3588324/1‑100 (MQ=255)
                                              gctgcAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGaaa                                                   <  1:561699/100‑1 (MQ=255)
                                              gctgcAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCAGGTTTTTTTTGGTCTTTAGTGTGCGGATGAGGTCGaaa                                                   >  1:3227507/1‑100 (MQ=255)
                                                tgcAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAAc                                                 >  2:938896/1‑100 (MQ=255)
                                                tgcAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAAc                                                 >  2:3080887/1‑100 (MQ=255)
                                                 gcAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACa                                                <  1:1897541/100‑1 (MQ=255)
                                                 gcAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACa                                                <  2:3365818/100‑1 (MQ=255)
                                                  cAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAg                                               <  2:3809176/100‑1 (MQ=255)
                                                   aaaCGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACACGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGc                                              >  1:2440114/1‑100 (MQ=255)
                                                    aaCGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGcc                                             <  2:3458742/100‑1 (MQ=255)
                                                    aaCGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGcc                                             <  1:468011/100‑1 (MQ=255)
                                                       gTAATCGCTGTGTAAGATACGTCAGCAAGAAGGGAAAGCCCGCTTGATCAGCGGGTTTTGTTTGGGCTTTAGTGTGCGGTTGAGGCCGGAAACAGCCAg                                           <  1:888232/99‑1 (MQ=255)
                                                       gTAATCGCTCTGTAAGATATGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGa                                          <  1:89355/100‑1 (MQ=255)
                                                            cGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGcc                                     <  1:3818462/100‑1 (MQ=255)
                                                              ctctGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTGTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAg                                   <  1:1859508/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)
                                                                     aGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTc                            <  2:549342/100‑1 (MQ=255)
                                                                          cGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGc                       >  2:1525833/1‑100 (MQ=255)
                                                                            tCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCaa                     >  2:1357887/1‑100 (MQ=255)
                                                                             cAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCaaa                    <  1:1357879/100‑1 (MQ=255)
                                                                              aGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGTTCTTTAGTGTGCGGTTGAAGCCGAAAACAGCCAGCATGGCAGTGCGGTCATGGCacaa                   >  1:1406659/1‑100 (MQ=255)
                                                                              aGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGTCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCaaaa                   >  2:452077/1‑100 (MQ=255)
                                                                               gCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAg                  <  2:364616/100‑1 (MQ=255)
                                                                               gCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAg                  >  2:1641018/1‑100 (MQ=255)
                                                                               gCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAg                  <  2:1051219/100‑1 (MQ=255)
                                                                                cAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCaagatc                 >  2:608447/1‑96 (MQ=255)
                                                                                  aGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGAcc               >  2:934993/1‑100 (MQ=255)
                                                                                   gAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGAccc              <  2:1498635/100‑1 (MQ=255)
                                                                                    aaTTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGAcccc             <  1:1027332/100‑1 (MQ=255)
                                                                                     aaaaaaaaCCCGCTGATTAAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCa            >  2:3543053/5‑100 (MQ=38)
                                                                                     aTTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCa            <  2:3588346/100‑1 (MQ=255)
                                                                                     aTTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCa            <  2:3543233/100‑1 (MQ=255)
                                                                                      ttCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAg           <  2:2343159/100‑1 (MQ=255)
                                                                                      aaaaaaaCCCGCTGATTAAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAg           >  2:589705/4‑100 (MQ=38)
                                                                                       aaaaaaCCCGCTGATTAAGCGGGTTTTTTTTGGTCTTTAGTGTTCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACACCAGa          >  1:3042732/3‑100 (MQ=38)
                                                                                         aaaaCCCGCTGATTAAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAg        >  2:3275517/1‑100 (MQ=255)
                                                                                            aCCCGCTGATTAAGCGGGTTTTTTTTGGTTTTTAGTGCTCGGTTGAGGCCGAAAACAGCAAGAATGTCAGTGCGGTCATGGCCCAAGACCCCAGaaag       >  1:408990/1‑96 (MQ=255)
                                                                                             cccGCTTAATCAGCGGGTTTTTTTTTTTTTTTTTTTTTCTGTTGAGGCCGAAATCAGCCAGAATGCCATTGCGGTCATGGCAAAAGACCCCAGAAGGtt     >  1:2720117/1‑99 (MQ=255)
                                                                                              ccGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGa   <  2:68153/100‑1 (MQ=255)
                                                                                               cGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGAc  <  1:1463186/100‑1 (MQ=255)
                                                                                                   |                                                                                               
TCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGAC  >  NC_000913/657432‑657626

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: