Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 657,531 T→G 18.2% 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,5310TG18.2% 324.8 / 74.3 139intergenic (+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 (45/68);  new base G (9/16);  total (55/84)
Fisher's exact test for biased strand distribution p-value = 8.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.59e-01

Only 100 of 144 total aligned reads displayed.

TCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAA  >  NC_000913/657432‑657630
                                                                                                   |                                                                                                    
tCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCtt                                                                                                      >  2:1226163/1‑100 (MQ=255)
   gcgcGTAGAGTTCGAAATCACTAACGGTGCCAGAGGCCCTTCTGCTGCCAACGTAATCGCTCTGTAAGATACGTCAGCAAGTATTCAAAACCCGCTTAAt                                                                                                   >  1:1714532/1‑100 (MQ=255)
   gcgcGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAAt                                                                                                   >  1:1022346/1‑100 (MQ=255)
    cgcgTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATc                                                                                                  <  2:644760/100‑1 (MQ=255)
    cgcgTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATc                                                                                                  >  2:3080561/1‑100 (MQ=255)
    cgcgTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATc                                                                                                  <  2:2733550/100‑1 (MQ=255)
       gTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGc                                                                                               <  1:1682542/100‑1 (MQ=255)
        tAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCg                                                                                              <  2:1159038/100‑1 (MQ=255)
        tAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCg                                                                                              <  1:127844/100‑1 (MQ=255)
          gagTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCggg                                                                                            <  1:3789634/100‑1 (MQ=255)
            gTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCggg                                                                                            >  2:1865071/1‑98 (MQ=255)
            gTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑tt                                                                                         >  1:683554/1‑100 (MQ=255)
            gTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑tt                                                                                         <  1:10279/100‑1 (MQ=255)
            gTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑tt                                                                                         >  1:3456325/1‑100 (MQ=255)
              tCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑tttt                                                                                       <  1:3799591/100‑1 (MQ=255)
              tCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGG‑tttt                                                                                       <  2:900531/100‑1 (MQ=255)
              tCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGG‑tttt                                                                                       <  1:797300/100‑1 (MQ=255)
              tCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGG‑tttt                                                                                       <  1:308396/100‑1 (MQ=255)
               cGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGNG‑ttttt                                                                                      >  1:3419226/1‑100 (MQ=255)
               cGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGG‑ttttg                                                                                      <  1:3910050/100‑2 (MQ=39)
               cGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGG‑ttttg                                                                                      <  2:2238136/100‑2 (MQ=39)
               cGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCAATCTGTATGCTCCGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑ttttt                                                                                      >  2:1492627/1‑100 (MQ=255)
                gAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑tttttt                                                                                     <  1:1826265/100‑1 (MQ=255)
                gAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGG‑ttttga                                                                                     <  1:138547/100‑3 (MQ=38)
                gAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGG‑ttttga                                                                                     <  1:322017/100‑3 (MQ=38)
                 aaaTCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑ttttttt                                                                                    <  1:3080581/100‑1 (MQ=255)
                 aaaTCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGG‑ttttgaa                                                                                    <  2:1472467/100‑4 (MQ=38)
                  aaTCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAGAACCCGCTTAATCAGCGGG‑tttttttt                                                                                   <  1:231812/100‑1 (MQ=255)
                  aaTCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGG‑ttttgaat                                                                                   <  1:144876/100‑5 (MQ=38)
                  aaTCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGG‑ttttgaat                                                                                   <  2:2007093/100‑5 (MQ=38)
                   aTCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTg                                                                                  <  2:2823376/100‑1 (MQ=255)
                    tCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTgg                                                                                 <  1:454296/100‑1 (MQ=255)
                    tCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTgg                                                                                 <  1:2736469/100‑1 (MQ=255)
                    tCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTgg                                                                                 <  2:3942091/100‑1 (MQ=255)
                      aCTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTc                                                                               >  2:3048953/1‑100 (MQ=255)
                       atAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTACTCGCTCTGTAAGATACATCGGCCAGAACTCAGAACCCGCTTACTCAGCGGG‑TATTTTTCGGt                                                                                >  2:708266/2‑98 (MQ=255)
                         aaCGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCttt                                                                            <  2:2252758/100‑1 (MQ=255)
                         aaCGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCttt                                                                            <  2:1415138/100‑1 (MQ=255)
                           cGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAg                                                                          <  1:3081155/100‑1 (MQ=255)
                                  aaaGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCgg                                                                   <  2:3700195/100‑1 (MQ=255)
                                   aaGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGt                                                                  <  1:3322684/100‑1 (MQ=255)
                                    aGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGtt                                                                 <  2:348410/100‑1 (MQ=255)
                                     ggCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTg                                                                <  1:1092857/100‑1 (MQ=255)
                                      gCCCTTCTGCTTCAAACGTAGTCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGt                                                                  <  2:2685327/97‑1 (MQ=255)
                                        ccTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAgg                                                             <  2:1430036/100‑1 (MQ=255)
                                         cTTCTGCTGCAAACGTAATTGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGc                                                            <  1:405086/100‑1 (MQ=255)
                                         cTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGc                                                            >  1:2344430/1‑100 (MQ=255)
                                          ttCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTTTTTAGTGTGCGGTTGAGGcc                                                           <  2:2412312/100‑1 (MQ=255)
                                          ttCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGcc                                                           <  1:3604342/100‑1 (MQ=255)
                                             tgctgcAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGaa                                                        <  2:2448450/100‑1 (MQ=255)
                                              gctgcAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGaaa                                                       <  2:1596893/100‑1 (MQ=255)
                                               ctgcAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGaaaa                                                      <  2:381899/100‑1 (MQ=255)
                                                 gcAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGAGGG‑TTGTTTTTGTTCTTTAGTGTGCGGTTGAGGCAGAAAAc                                                     >  2:3658763/1‑99 (MQ=255)
                                                  cAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGTG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACa                                                    <  2:3759926/99‑1 (MQ=255)
                                                  cAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACa                                                    <  2:2472110/99‑1 (MQ=255)
                                                     aCGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCa                                                <  2:104125/100‑1 (MQ=255)
                                                      cGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAg                                               >  1:1001416/1‑100 (MQ=255)
                                                       gTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTATGGTCTTTAGTGTGCGGTTGAGGCTGAAAACAGCCAGa                                              <  2:2434401/100‑1 (MQ=255)
                                                        tAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGaa                                             >  1:3863623/1‑100 (MQ=255)
                                                            cGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGcc                                         >  2:2339452/1‑100 (MQ=255)
                                                              ctctGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAg                                       <  1:1672482/100‑1 (MQ=255)
                                                              ctctGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAg                                       <  1:3846883/100‑1 (MQ=255)
                                                                ctGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTg                                     >  1:1249082/1‑100 (MQ=255)
                                                                ctGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTg                                     <  1:257284/100‑1 (MQ=255)
                                                                  gTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTTTTTAGTGTGCGATTGAGGCCGAAAACAGCCAGAATGCCAGTGCg                                   >  2:887472/1‑100 (MQ=255)
                                                                  gTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCg                                   >  1:2507299/1‑100 (MQ=255)
                                                                    aaGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGt                                 <  2:3456302/100‑1 (MQ=255)
                                                                    aaGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGt                                 >  1:1409331/1‑100 (MQ=255)
                                                                        tACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATg                             <  2:1016021/100‑1 (MQ=255)
                                                                        tACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATg                             <  1:710980/100‑1 (MQ=255)
                                                                          cGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATgg                            <  2:1079059/99‑1 (MQ=255)
                                                                          cGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGc                           <  1:3310506/100‑1 (MQ=255)
                                                                           gTCAGCAAGAATTCAAAACCCGCTTAATCAGCGAG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCa                          <  1:3290491/100‑1 (MQ=255)
                                                                             cAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCaaa                        >  1:370833/1‑100 (MQ=255)
                                                                              aGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCaaaa                       >  1:2703600/1‑100 (MQ=255)
                                                                               gCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTTAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAg                      >  2:2921740/1‑100 (MQ=255)
                                                                               gCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAg                      >  1:3231258/1‑100 (MQ=255)
                                                                               gCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAg                      <  1:2251186/100‑1 (MQ=255)
                                                                                 aaGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGa                     >  1:121333/1‑100 (MQ=255)
                                                                                 aaGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGAc                    <  2:3843642/100‑1 (MQ=255)
                                                                                  aGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGAcc                   >  1:643805/1‑100 (MQ=255)
                                                                                   ccaaaaaaaaCCCGCTGATTAAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGAccc                  >  2:137942/7‑100 (MQ=37)
                                                                                    aaTTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGAcccc                 <  1:1044127/100‑1 (MQ=255)
                                                                                    aaTTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGAcccc                 >  2:2735151/1‑100 (MQ=255)
                                                                                     aaaaaaaaCCCGCTGATTAAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCATAATGCCAGGGCGGTCATGGCAAAAGACCCCa                >  2:2657973/5‑100 (MQ=37)
                                                                                     aaaaaaaaCCCGCTGATTAAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCa                >  1:1322015/5‑100 (MQ=38)
                                                                                     aaaaaaaaCCCGCTGATTAAGCGGG‑TTTTGTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCa                >  2:405052/5‑100 (MQ=38)
                                                                                     aTTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTTTTTTTTTTTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGAccc                  >  2:3008354/1‑98 (MQ=255)
                                                                                     aTTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCa                >  1:3364551/1‑100 (MQ=255)
                                                                                         aaaaCCCGCTGATTAAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAg            >  1:711891/1‑100 (MQ=255)
                                                                                           aaCCCGCTTAATCAGCGGG‑TTTGTTTTGGTCTTTAGGGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGt          <  2:1274152/100‑1 (MQ=255)
                                                                                           aaCCCGCTGATTAAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCGAGTGCGGTCATGGCAAAAGACCCCAGa              >  2:3330221/1‑96 (MQ=255)
                                                                                            aCCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGtt         >  2:2639573/1‑100 (MQ=255)
                                                                                            aCCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTCAGGCCGAAAGCAGCCAGAATGCCAGTGCGGTCATGGCGAAAGACCCCa                >  2:2303551/1‑93 (MQ=255)
                                                                                             cccGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTg        <  2:1695921/100‑1 (MQ=255)
                                                                                              ccGCTTAATCAGCGGG‑GTTTTTTTGGGCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGt          <  2:2358724/97‑1 (MQ=255)
                                                                                                gCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACg     >  1:3430832/1‑100 (MQ=255)
                                                                                                gCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACg     <  1:2191023/100‑1 (MQ=255)
                                                                                                gCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACg     <  1:1718238/100‑1 (MQ=255)
                                                                                                   tAATCAGCGGG‑TTTTTTTTTGTCTTTAGTTTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTCACGaaa  >  2:3264292/1‑100 (MQ=255)
                                                                                                   |                                                                                                    
TCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGG‑TTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAA  >  NC_000913/657432‑657630

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: ‑