Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 454,856 G→A 15.6% 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.6% 419.5 / 53.6 175intergenic (+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 (64/82);  new base A (11/16);  total (76/99)
Fisher's exact test for biased strand distribution p-value = 8.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.03e-01

Only 100 of 187 total aligned reads displayed.

CCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATAACCG  >  NC_000913/454757‑454955
                                                                                                   |                                                                                                   
 ctccCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCgg                                                                                                    >  2:2202888/1‑100 (MQ=255)
  tccCTGTCGTGGAGCAGGAAGCATCGCGTAAAATCGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCggg                                                                                                   <  1:3185004/100‑1 (MQ=255)
  tccCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCggg                                                                                                   <  1:3734310/100‑1 (MQ=255)
  tccCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCggg                                                                                                   >  1:2714345/1‑100 (MQ=255)
     cTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATGTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCg                                                                                                >  1:2759199/1‑100 (MQ=255)
      tGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCgt                                                                                               <  2:240555/100‑1 (MQ=255)
       gTCGTGTAGCAGGAATCATATAGTAAAAACGCATTTGTAACTGTCGGCGCTTTTGCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtt                                                                                              <  1:2238156/100‑1 (MQ=255)
        tCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTATATATTTTATGAAAAACGGCCTGCGGGCCgtt                                                                                              <  2:1802457/99‑1 (MQ=255)
         cGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtttt                                                                                            <  1:3331827/100‑1 (MQ=255)
           tGGAGCAGGAAGCATCGCGTAAAAACGTATTTGCAACTGTCTGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtttttc                                                                                          <  1:2106029/100‑3 (MQ=255)
           tGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAATTGTTGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtttttc                                                                                          <  1:562061/100‑3 (MQ=255)
           tGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtttttc                                                                                          <  2:3411897/100‑3 (MQ=255)
            ggaggaggaAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTTTTTTGTTAAAAAAAAAAAAAAAAACGGCCTGCGGGcc                                                                                                 <  2:15258/92‑1 (MQ=255)
            ggAGCAGGAAGCATCGCGTAAAAACGTATTTGTAATTGTCGGCGTTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtttttca                                                                                         <  1:1098845/100‑4 (MQ=255)
            ggAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGTTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCgttttgtt                                                                                         <  1:3285708/100‑1 (MQ=255)
            ggAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCgttttgtt                                                                                         >  1:234141/1‑100 (MQ=255)
            ggAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtttttca                                                                                         <  2:3102680/100‑4 (MQ=255)
             gtgCAAGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTACGTTGCCAAAGATCTTATGAAAAGCGGCCTGCTGGCCggtttg                                                                                           >  2:1329982/3‑97 (MQ=255)
             gAGCAGGAAGCATCGCGTAAACTCGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCgttttgttt                                                                                        <  2:2540453/100‑1 (MQ=255)
             gAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtttttcat                                                                                        <  1:2027304/100‑5 (MQ=39)
             gAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtttttcat                                                                                        <  1:866600/100‑5 (MQ=39)
             gAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtttttca                                                                                         <  1:3562110/99‑4 (MQ=255)
             gAGCAAGAAGCATCGCTTAAAATTGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtttttcat                                                                                        <  1:3843865/100‑5 (MQ=37)
              aGCAGTAAGCATCGCGTAAAAACTTTTTTGTAACTGTCGGTGTTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCgttttgttt                                                                                        <  1:1735093/99‑1 (MQ=255)
              aGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCgttttgtttt                                                                                       <  2:1079670/100‑1 (MQ=255)
               gCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCgttttgttt                                                                                        <  2:457220/98‑1 (MQ=255)
               gCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGttttgttttg                                                                                      <  1:1784568/100‑1 (MQ=255)
               gCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGttttgttttg                                                                                      >  2:3511534/1‑100 (MQ=255)
               gCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGttttgttttg                                                                                      <  2:245405/100‑1 (MQ=255)
                cAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTtttgttttgt                                                                                     >  2:2375849/1‑100 (MQ=255)
                cAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTtttgttttgt                                                                                     <  2:2911869/100‑1 (MQ=255)
                cAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtttttcata                                                                                       <  2:3886601/98‑6 (MQ=38)
                  ggAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCt                                                                                   >  1:1649774/1‑100 (MQ=255)
                  ggAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACAGCCTGCGGGCAGTTTTGTTTTGTCt                                                                                   >  1:874234/1‑100 (MQ=255)
                   gAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTg                                                                                  <  1:1536023/100‑1 (MQ=255)
                     aGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCATTATGTTGCTAAAGATTTTATGAAAGACGGCCTGCGGGCCGTTTTGTTTTGTCTGGa                                                                                <  2:892955/100‑1 (MQ=255)
                        aTCGCGTAAAAACTCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGAt                                                                               <  1:3870331/98‑1 (MQ=255)
                         tCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGAtttt                                                                            >  2:435751/1‑100 (MQ=255)
                           gcgTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTgc                                                                          <  2:3285932/100‑1 (MQ=255)
                             gTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTgcgc                                                                        <  2:936841/100‑1 (MQ=255)
                                aaaaCGCATTTGCAACTGTCGGCGCTTTTCCAATATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGATTTGTTTTGTCTGGATTTTGCGCttt                                                                     >  1:3468939/1‑100 (MQ=255)
                                   aCGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTg                                                                  >  2:2906086/1‑100 (MQ=255)
                                    cGCATTTGCAACTGTCGGTGCTTTTCCAGTATGTTGTTTTTTTTTTTATGAAAAACGGCGTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTg                                                                  <  1:3209465/99‑1 (MQ=255)
                                           gCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCAt                                                          <  1:2690758/100‑1 (MQ=255)
                                            cAACTGTCGTCGCTTTTCCAGGAGTTTGCTAAAAATTTTATGAAAAACGTCCTGCGTGCCGTTTTGTTTTGTCTGGATTGTGCGCTTTTTGCCCAg                                                             <  2:226503/96‑1 (MQ=255)
                                            cAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCAtt                                                         >  2:2388408/1‑100 (MQ=255)
                                            cAACTGTCGGCGCTTTTCCAGCATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCAtt                                                         >  2:1725285/1‑100 (MQ=255)
                                             aaCTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTc                                                        <  1:2159301/100‑1 (MQ=255)
                                              aCTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATTAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCa                                                       <  1:1705483/100‑1 (MQ=255)
                                              aCTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCa                                                       <  2:1242913/100‑1 (MQ=255)
                                               cTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAg                                                      >  2:194652/1‑100 (MQ=255)
                                               cTGTCGGCGCTTTTCCAGTATGTTGCGAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTCCCGCTTTTTGCCCAGCATTCa                                                       <  2:3729304/99‑1 (MQ=255)
                                                   cGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGa                                                  >  1:3037957/1‑100 (MQ=255)
                                                    ggCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGaa                                                 >  2:3604733/1‑100 (MQ=255)
                                                    ggCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGaa                                                 >  2:1436459/1‑100 (MQ=255)
                                                     gcgcTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGaaa                                                >  2:3472946/1‑100 (MQ=255)
                                                          tttCCAGTATGTTGCTAAAGATTTTTTGAAAAACTGCCTGTGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGc                                           <  1:1348272/100‑1 (MQ=255)
                                                          tttCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGc                                           <  2:2892442/100‑1 (MQ=255)
                                                          tttCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTATATGGATTTTGCGATTTTTGCACAGAATTCAGACGAAAATTGc                                           >  1:606011/1‑100 (MQ=255)
                                                              cAGTATGTTGCTAAATATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGGCGAAAATTGCCCgg                                       <  1:2586131/100‑1 (MQ=255)
                                                               aGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCggg                                      <  1:158383/100‑1 (MQ=255)
                                                               aGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCggg                                      <  2:3374627/100‑1 (MQ=255)
                                                                gTATGTTGCTAAAGATTTTATAAAAAACGGCCTGCGGTCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACTAAAATTGccc                                         >  2:773939/1‑96 (MQ=255)
                                                                  aTGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCGGTTTTGTTTTTTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAAt                                   >  2:2078903/1‑100 (MQ=255)
                                                                   tGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGAGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAAtt                                  >  1:1981079/1‑100 (MQ=255)
                                                                    gTTTCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATtg                                 <  1:689067/100‑1 (MQ=255)
                                                                    gTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGTCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATtg                                 <  1:291275/100‑1 (MQ=255)
                                                                    gTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCTGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATtg                                 <  1:1879748/100‑1 (MQ=255)
                                                                    gTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATtg                                 <  2:1376665/100‑1 (MQ=255)
                                                                       gCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGa                              <  1:3930434/100‑1 (MQ=255)
                                                                       gCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGa                              >  1:3231094/1‑100 (MQ=255)
                                                                          aaaGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTAAGACGAAAATTGCCCGGGAATGGTGa                              <  2:2592155/97‑1 (MQ=255)
                                                                          aaaGATTTTATGAAAAACGGACTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGaaa                            >  2:1180985/1‑99 (MQ=255)
                                                                            aGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGaaaaaa                         <  2:3921690/100‑1 (MQ=255)
                                                                             cagacaaaacaaaaCGGCCCGCAGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAAt                        >  2:562381/11‑100 (MQ=37)
                                                                               ttttATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAACATTGCCCGGGAATTGTGAAAAAATAc                      <  1:2388430/100‑1 (MQ=255)
                                                                               ttttATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTTCCCGGGAATTGTGAAAAAATAc                      <  2:2647642/100‑1 (MQ=255)
                                                                               ttttATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCAAGCATTCAGACGAAAATTACCCGGGAATTGTGAAAAAATAc                      >  1:471656/1‑100 (MQ=255)
                                                                                tttATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATAcg                     <  2:2308058/100‑1 (MQ=255)
                                                                                 ttATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATAcgc                    <  2:1968716/100‑1 (MQ=255)
                                                                                 ttATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATAcgc                    >  1:980062/1‑100 (MQ=255)
                                                                                   aaacaaaaCGGCCCGCAGGCCGATTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGAACAAAATCCCCCGGGAATTGTGAAAAAATAc                      >  2:2264386/5‑96 (MQ=255)
                                                                                   aaacaaaaCGGCCCGCAGGACGATTTGTTTTGTATGGATTTTGCGATTTTTGCCCAGCATTCAGACGAAAATTGCACGGGAACTGTGAAAAAATAc                      >  2:3145874/5‑96 (MQ=255)
                                                                                    tGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGAc                 >  2:694754/1‑100 (MQ=255)
                                                                                    tGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGAc                 >  2:1220327/1‑100 (MQ=255)
                                                                                     gAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCTCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACa                >  1:1417347/1‑100 (MQ=255)
                                                                                     acaaaaCGGCCCGCAGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTTAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACa                >  2:2824782/3‑100 (MQ=255)
                                                                                      caaaaCGGCCCGCAGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACa                >  1:791473/2‑99 (MQ=255)
                                                                                      aaaaaCGGCCTGCGGGCCGTTTTGATTTGTTTGGATTATGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAg               >  2:1516354/1‑100 (MQ=255)
                                                                                       aaaaCGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAgc              >  1:2496349/1‑100 (MQ=255)
                                                                                        aaaCGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAgcg             <  1:3140302/100‑1 (MQ=255)
                                                                                          aCGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAgcgcg           >  1:400992/1‑100 (MQ=255)
                                                                                           cGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAgcgcgc          <  1:3053987/100‑1 (MQ=255)
                                                                                            ggCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCa         <  2:46047/100‑1 (MQ=255)
                                                                                            ggCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCAGGGAATTGTGAAAAAATACGCGACAGCGCGCa         >  1:552137/1‑100 (MQ=255)
                                                                                             gCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCaa        >  1:2371647/1‑100 (MQ=255)
                                                                                                tGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATaa     >  2:3542669/1‑100 (MQ=255)
                                                                                                 gCGGGCCGTTTTGTTTTGTCTGGATTTTGCGTTTTTTTTCCAGTATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATAAc    <  1:1096679/100‑1 (MQ=255)
                                                                                                 gCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTTTTTTTTTTTCAGACGAAAATTGCCCTAGAATTGTGAAAAAATACGCGACAGCACACAATAAc    >  2:510401/1‑100 (MQ=255)
                                                                                                   gggCCGTTGTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATAAcc   <  2:1876007/99‑1 (MQ=255)
                                                                                                   |                                                                                                   
CCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATAACCG  >  NC_000913/454757‑454955

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: