Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 454,856 G→A 15.8% intergenic (+67/‑277) bolA → / → tig stationary‑phase morphogene, transcriptional repressor for mreB; also regulator for dacA, dacC, and ampC/peptidyl‑prolyl cis/trans isomerase (trigger factor)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913454,8560GA15.8% 442.8 / 57.0 181intergenic (+67/‑277)bolA/tigstationary‑phase morphogene, transcriptional repressor for mreB; also regulator for dacA, dacC, and ampC/peptidyl‑prolyl cis/trans isomerase (trigger factor)
Reads supporting (aligned to +/- strand):  ref base G (75/76);  new base A (5/24);  total (81/100)
Fisher's exact test for biased strand distribution p-value = 1.81e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

Only 100 of 183 total aligned reads displayed.

CTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATAACCG  >  NC_000913/454758‑454955
                                                                                                  |                                                                                                   
ctccCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCgg                                                                                                    <  1:2778863/100‑1 (MQ=255)
    cTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCg                                                                                                >  2:2841997/1‑100 (MQ=255)
     tGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgt                                                                                               <  1:3596963/100‑1 (MQ=255)
     tGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTAGGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCgt                                                                                               >  2:905659/1‑100 (MQ=255)
      gTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtt                                                                                              <  1:1161429/100‑1 (MQ=255)
       tCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCgttt                                                                                             <  1:2146747/100‑1 (MQ=255)
        ggTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtttt                                                                                            <  2:2740656/99‑1 (MQ=255)
         gTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgttttt                                                                                           <  1:396343/100‑2 (MQ=255)
         gTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgttttt                                                                                           <  1:468384/100‑2 (MQ=255)
          tGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtttttc                                                                                          <  1:1801104/100‑3 (MQ=255)
          tGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtttttc                                                                                          <  1:2978893/100‑3 (MQ=255)
           ggAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCgttttgtt                                                                                         >  1:3802263/1‑100 (MQ=255)
           ggAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCgttttgtt                                                                                         <  1:3346792/100‑1 (MQ=255)
           ggAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCgttttgtt                                                                                         <  2:244325/100‑1 (MQ=255)
           ggAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCgttttgt                                                                                          >  1:2080906/1‑99 (MQ=255)
           ggAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtttttca                                                                                         <  2:358788/100‑4 (MQ=255)
           ggAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtttttca                                                                                         <  1:302235/100‑4 (MQ=255)
           ggAGCAGGAAGCATCGCGTAAAAACGCATTTGAAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtttttca                                                                                         <  2:881400/100‑4 (MQ=255)
            gAGCAGGAAGCATCGCGTAAAAACGTATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTACAGATTTTATGAAAAACGGCCCGCAGGCCgtttttcat                                                                                        <  2:2229271/100‑5 (MQ=38)
            gAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCgttttgttt                                                                                        <  2:549221/100‑1 (MQ=255)
            gAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtttttcat                                                                                        <  2:3609655/100‑5 (MQ=39)
            gAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtttttcat                                                                                        <  1:1968696/100‑5 (MQ=39)
             aGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCgttttgtttt                                                                                       <  2:1383607/100‑1 (MQ=255)
              gCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGttttgttttg                                                                                      <  1:4191743/100‑1 (MQ=255)
               cAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTttcgttttgt                                                                                     <  1:2203562/100‑1 (MQ=255)
                aGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTc                                                                                    <  1:3124081/100‑1 (MQ=255)
                aGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTc                                                                                    >  2:1085613/1‑100 (MQ=255)
                aGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGAATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTtttgttttgt                                                                                     >  2:2578800/1‑99 (MQ=255)
                 ggAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTGTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCt                                                                                   <  2:3062161/100‑1 (MQ=255)
                  gAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTg                                                                                  <  2:427432/100‑1 (MQ=255)
                  gAAGCATCGCGTAAAAACGCAGTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTg                                                                                  <  1:249719/100‑1 (MQ=255)
                   aaGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTgg                                                                                 <  1:2723500/100‑1 (MQ=255)
                     gCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTgg                                                                                 >  1:433419/1‑98 (MQ=255)
                       aTCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGAttt                                                                             >  2:3137633/1‑100 (MQ=255)
                       aTCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGAttt                                                                             <  2:3390334/100‑1 (MQ=255)
                         cgcgTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTg                                                                           <  1:2134904/100‑1 (MQ=255)
                         cgcgTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTg                                                                           <  1:3645973/100‑1 (MQ=255)
                          gcgTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGACAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTgc                                                                          >  2:851095/1‑100 (MQ=255)
                          gcgTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTgc                                                                          >  1:3493716/1‑100 (MQ=255)
                            gTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTgcgc                                                                        <  1:349577/100‑1 (MQ=255)
                             tAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCt                                                                       >  2:3068600/1‑100 (MQ=255)
                               aaaaCGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCttt                                                                     <  2:96940/100‑1 (MQ=255)
                                aaaCGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCtttt                                                                    >  1:399872/1‑100 (MQ=255)
                                 aaCGCATTTGCAACTGTCGGCGCGTTTCCAGTATGTTGCTATAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGTATTTTGCGCtttt                                                                    <  2:491000/99‑1 (MQ=255)
                                   cGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGc                                                                 <  1:22186/100‑1 (MQ=255)
                                   cGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGc                                                                 >  2:4069100/1‑100 (MQ=255)
                                   cGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGc                                                                 >  1:1679699/1‑100 (MQ=255)
                                   cGCATGTACAACCGTCGGCGACTCTCCAGTATGCTGATAAAGATTGTATAAAAATCGGACTACAGCCAGTTTTATTTTGTCTGGATTTTGCGCTTTTAGc                                                                 <  2:3477006/100‑1 (MQ=255)
                                     cATTTGCAACTGTCGGCGCTTTTCCAGTATGTTTCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGccc                                                               <  2:2469447/100‑1 (MQ=255)
                                     cATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGccc                                                               <  2:1718365/100‑1 (MQ=255)
                                       tttGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAg                                                             <  1:1466456/100‑1 (MQ=255)
                                        ttGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGc                                                            <  2:293058/100‑1 (MQ=255)
                                          gCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCAt                                                          <  2:643328/100‑1 (MQ=255)
                                             aCTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCa                                                       >  1:338485/1‑100 (MQ=255)
                                             aCTGTCGGCGATTTTCCAGTATTTTGTTAATGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTTCGCTTTTTGCCCAGCATTCa                                                       <  2:3647437/100‑1 (MQ=255)
                                              cTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAg                                                      >  2:1589875/1‑100 (MQ=255)
                                                 tCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACg                                                   <  1:3750604/100‑1 (MQ=255)
                                                  cGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGa                                                  >  1:2639435/1‑100 (MQ=255)
                                                   ggCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTTTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGaa                                                 >  2:1661537/1‑100 (MQ=255)
                                                   ggCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGGTTTTTCGCTTTTTTCCCAGCATTCAAACGaa                                                 >  2:804942/1‑100 (MQ=255)
                                                   ggCGCTTTTCCAGAATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGaa                                                 >  2:2596816/1‑100 (MQ=255)
                                                        ttttCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTg                                            <  2:543557/100‑1 (MQ=255)
                                                        ttttCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTg                                            >  1:2899366/1‑100 (MQ=255)
                                                         tttCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGc                                           >  2:4178566/1‑100 (MQ=255)
                                                          ttCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGcc                                          >  1:1835941/1‑100 (MQ=255)
                                                           tCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGccc                                         <  1:201540/100‑1 (MQ=255)
                                                           tCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGccc                                         <  1:269568/100‑1 (MQ=255)
                                                           tCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGccc                                         <  1:195776/100‑1 (MQ=255)
                                                            caaGTATGTTGCTAAAGATTTTATGAAAAACCGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCg                                        >  1:2826666/3‑100 (MQ=255)
                                                             cAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCgg                                       >  1:1629963/1‑100 (MQ=255)
                                                                  tGTTGCTAACGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAAtt                                  >  2:1514787/1‑100 (MQ=255)
                                                                   gTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATtg                                 <  2:2113160/100‑1 (MQ=255)
                                                                    ttGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATtgt                                >  1:3759654/1‑100 (MQ=255)
                                                                     tGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATtgtg                               <  1:1139120/100‑1 (MQ=255)
                                                                      gCTAAAGATTTTATGAAAAACTGCCTGCGTGCCGTATTGTCTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCCGACGAAAATTGCCCGGGAATTGTGa                              >  2:38240/1‑100 (MQ=255)
                                                                      gCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGa                              <  1:696950/100‑1 (MQ=255)
                                                                       cTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGaa                             >  1:1471/1‑100 (MQ=255)
                                                                         aaaGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGaaaa                           >  2:2970353/1‑100 (MQ=255)
                                                                          aaGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGaaaaa                          <  1:1264555/100‑1 (MQ=255)
                                                                              ttttATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATAc                      <  2:119649/100‑1 (MQ=255)
                                                                               tttATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAACAAATAc                      >  1:3886521/1‑99 (MQ=255)
                                                                                  aaacaaaaCGGCCCGCAGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGa                  >  2:3970776/5‑100 (MQ=39)
                                                                                  aaacaaaaCGGCCCGCAGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGa                  >  1:968721/5‑100 (MQ=39)
                                                                                   tGAAAGACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGAc                 >  1:532687/1‑100 (MQ=255)
                                                                                    gAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGTTTTTTCGCTTTTGGCCCAGCACTCCAACGAAAAATGCAAGGGACTTGTGaaaccatcc                      >  2:2057147/1‑88 (MQ=255)
                                                                                    gAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACa                >  1:1859811/1‑100 (MQ=255)
                                                                                      aaaaCGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACa                >  1:284425/1‑98 (MQ=255)
                                                                                      aaaaCGGCCCGCAGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCCGCATTCAGACGAAAATTGCACGGGGATTGTGAAAAAATACGCGACAgc              >  2:795601/1‑100 (MQ=255)
                                                                                      aaaaCGGCCCGCAGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAgc              >  2:1852230/1‑100 (MQ=255)
                                                                                        aaCGGCCCGCAGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAgcgc            >  1:2043040/1‑100 (MQ=255)
                                                                                          cGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAgcgcgc          <  2:1171917/100‑1 (MQ=255)
                                                                                           ggCCTGCTGGCCGTTTTGATTTGTCTGGCTTGTGCGCTTTTTGCCCAGCGTACAGACGAAGATGGGCCGGGTGTTGTCATAAAATACGCGACAgc              >  2:2497311/1‑95 (MQ=255)
                                                                                           ggCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGGATTGTGAAAAAATACGCGACAGCGCGCa         >  2:1034788/1‑100 (MQ=255)
                                                                                           ggCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCa         <  1:7562/100‑1 (MQ=255)
                                                                                             ccTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAAt       <  2:3442604/100‑1 (MQ=255)
                                                                                             ccTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAAt       <  1:2568276/100‑1 (MQ=255)
                                                                                                gCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATAAc    <  1:3285389/100‑1 (MQ=255)
                                                                                                 cGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATAAcc   >  2:1007447/1‑100 (MQ=255)
                                                                                                 cGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATAAcc   <  2:1482251/100‑1 (MQ=255)
                                                                                                  gggCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATAACCg  <  1:230713/100‑1 (MQ=255)
                                                                                                  |                                                                                                   
CTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATAACCG  >  NC_000913/454758‑454955

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: