Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,251,287 T→C 46.0% pseudogene (327/1323 nt) yjbI → pseudogene, SopA‑related, pentapeptide repeats‑containing

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,251,2870TC46.0% 43.9 / 254.3 191pseudogene (327/1323 nt)yjbIpseudogene, SopA‑related, pentapeptide repeats‑containing
Reads supporting (aligned to +/- strand):  ref base T (56/47);  new base C (42/46);  total (98/93)
Fisher's exact test for biased strand distribution p-value = 3.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 193 total aligned reads displayed.

CAATGTGATTTTTATCAAAGCGAGTGGTGAATGTGTCGATTTCAGCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAATTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATCTTTACGCTTCGCACTTCACCAGAGCAGAGTTTCTGTCTGCCAAAGAAATATCA  >  NC_000913/4251188‑4251386
                                                                                                   |                                                                                                   
caATGTGATTTTTATCAAAGCGAGTGGTGAATGTGTCGATTTCAGCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAAt                                                                                                     >  2:1469567/1‑100 (MQ=255)
    gtgATTTTTATCAAAGCGAGTGGTGAATGTGTCGATTTCAGCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAACTTTc                                                                                                 >  1:2692944/1‑100 (MQ=255)
     tgATTTTTATCAAAGCGAGTTGTGAATGTGTCGATTTCAGCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAATTTTcg                                                                                                <  2:214488/100‑1 (MQ=255)
     tgATTTTTATCAAAGCGAGTGGTGAATGTGTCGATTTCAGCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAACTTTcg                                                                                                <  2:1568512/100‑1 (MQ=255)
      gATTTTTATCAAAGCGAGTGGTGAATGTGTCGATTTCAGCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAACTTTcgc                                                                                               >  1:522596/1‑100 (MQ=255)
       aTTTTTATCAAAGCGAGTGGTGAATGTGTCGATTTCAGCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAATTTTcgcg                                                                                              <  1:899736/100‑1 (MQ=255)
        tttttATCAAAGCGAGTGGTGAATGTGTCGATTTCAGCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAATTTTCGCGa                                                                                             >  2:3426866/1‑100 (MQ=255)
        tttttATCAAAGCGAGTGGTGAATGTGTCGATTTCAGCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAACTTTCGCGa                                                                                             >  1:595312/1‑100 (MQ=255)
         ttttATCAAAGCGAGTGGTGAATGTGTCGATTTCAGCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAACTTTCGCGaa                                                                                            <  2:3671570/100‑1 (MQ=255)
           ttATCAAAGCGAGTTGTGAATGTGTCGATTACAGCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAACTTTCGCGAATg                                                                                          <  2:3553537/100‑1 (MQ=255)
           ttATCAAAGCGAGTGGTGAATGTGTCGATTTCAGCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAATTTTCGCGAATg                                                                                          <  2:3963525/100‑1 (MQ=255)
           ttATCAAAGCGAGTGGTGAATGTGTCGATTTCAGCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAACTTTTGCGAATg                                                                                          <  1:2168478/100‑1 (MQ=255)
             aTCAAAGCGAGTGGTGAATGTGTCGATTTCAGCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAATTTTCGCGAATGTc                                                                                        >  1:4012984/1‑100 (MQ=255)
               cAAAGCGAGTGGTGAATGTGTCGATTTCAGCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAATTTTCGCGAATGTCAg                                                                                      >  1:1023101/1‑100 (MQ=255)
               cAAAGCGAGTGGTGAATGTGTCGATTTCAGCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAATTTTCGCGAATGTCAg                                                                                      >  2:2371965/1‑100 (MQ=255)
               cAAAGCGAGTGGTGAATGTGTCGATTTCAGCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAATTTTCGCGAATGTCAg                                                                                      >  1:1483375/1‑100 (MQ=255)
               cAAAGCGAGTGGTGAATGTGTCGATTTCAGCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAACTTTCGCGAATGTCAg                                                                                      <  2:1420014/100‑1 (MQ=255)
                aaaGCGAGTGGTGAATGTGTCGATTTCAGCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAATTTTCGCGAATGTCAGa                                                                                     >  2:464777/1‑100 (MQ=255)
                aaaGCGAGTGGTGAATGTGTCGATTTCAGCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAATTTTCGCGAATGTCAGa                                                                                     >  2:2527557/1‑100 (MQ=255)
                aaaGCGAGTGGTGAATGTGTCGATTTCAGCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAACTTTCGCGAATGTCAGa                                                                                     >  2:90244/1‑100 (MQ=255)
                aaaGCGAGTGGTGAATGTGTCGATTTCAGCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAACTTTCGCGAATGTCAGa                                                                                     <  1:465251/100‑1 (MQ=255)
                    cGAGTGGTGAATGTGTCGATTTCAGCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAATTTTCGCGAATGTCAGATTAg                                                                                 >  2:1634757/1‑100 (MQ=255)
                         ggtgAATGTGTCGATTTCAGCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAACTTTCGCGAATGTCAGATTAGAAAtt                                                                            <  1:3037342/100‑1 (MQ=255)
                            gAATGTGTCGATTTCAGCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAACTTTCGCGAATGTCAGATTAGAAATTCaa                                                                         >  2:1964582/1‑100 (MQ=255)
                             aaTGTGTCGATTTCAGCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAACTTTCGCGAATGTCAGATTAGAAATTCaaa                                                                        >  1:2993954/1‑100 (MQ=255)
                              aTGTGTCGATTTCAGCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAATTTTCGCGAATGTCAGATTAGAAATTCAAAc                                                                       >  1:922208/1‑100 (MQ=255)
                              aTGTGTCGATTTCAGCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAATTTTCGCGAATGTCAGATTAGAAATTCAAAc                                                                       <  1:4128577/100‑1 (MQ=255)
                              aTGTGTCGATTTCAGCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAACTTTCGCGAATGTCAGATTAGAAATTCAAAc                                                                       >  2:3810659/1‑100 (MQ=255)
                              aTGTGTCGATTTCAGAAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAATTTTCGCGAATGTCAGATTAGAAATTCAAAc                                                                       <  2:3417986/100‑1 (MQ=255)
                               tgtgTCGATTTCAGCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAATTTTCGCGAATGTCAGATTAGAAATTCAAACt                                                                      <  1:710688/100‑1 (MQ=255)
                               tgtgTCGATTTCAGCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAACTTTCGCGAATGTCAGATTAGAAATTCAAACt                                                                      >  2:3193956/1‑100 (MQ=255)
                                  gtCGATTTCAGCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAATTTTCGCGAATGTCAGATTAGAAATTCAAACTTCg                                                                   <  1:3940912/100‑1 (MQ=255)
                                  gtCGATTTCAGCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAATTTTCGCGAATGTCAGATTAGAAATTCAAACTTCg                                                                   <  2:1820117/100‑1 (MQ=255)
                                  gtCGATTTCAGCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAACTTTCGCGAATGTCAGATTAGAAATTCAAACTTCg                                                                   >  2:4048837/1‑100 (MQ=255)
                                   tCGATTTCAGCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAATTTTCGCGAATGTCACATTAGAAATTCAAACTTCGa                                                                  >  2:1226428/1‑100 (MQ=255)
                                   tCGATTTCAGCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAACTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGa                                                                  <  1:1275267/100‑1 (MQ=255)
                                    cGATTTCAGCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAATTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGAt                                                                 >  1:3262986/1‑100 (MQ=255)
                                    cGATTTCAGCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAACTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGAt                                                                 >  1:1908970/1‑100 (MQ=255)
                                    cGATTTCAGCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAACTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGAt                                                                 >  1:3487455/1‑100 (MQ=255)
                                      aTTTCAGCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAATTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATaa                                                               >  1:644926/1‑100 (MQ=255)
                                      aTTTCAGCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAATTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATa                                                                >  2:3076910/1‑99 (MQ=255)
                                         tCAGCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAACTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTg                                                            <  1:3596251/100‑1 (MQ=255)
                                           aGCAAAAATATTCTTGATACAGTTTACTTCTCGCAGAGTCAACTTGGTCATAGTAATTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTg                                                            <  2:2767548/98‑1 (MQ=255)
                                           aGCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAACTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGtt                                                          >  1:3831565/1‑100 (MQ=255)
                                            gCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAATTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTa                                                         <  1:3228025/100‑1 (MQ=255)
                                            gCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAACTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTa                                                         >  1:3365637/1‑100 (MQ=255)
                                            gCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAACTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTa                                                         >  1:959197/1‑100 (MQ=255)
                                             cAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAATTTTCGCTAATGTCGGATTAGAAATTCAAACTTCGATAATTGTTAt                                                        >  1:3714703/1‑100 (MQ=255)
                                             cAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAACTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTAt                                                        >  2:3751840/1‑100 (MQ=255)
                                              aaaaaTATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAATTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATc                                                       <  1:2073212/100‑1 (MQ=255)
                                              aaaaaTATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAATTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATc                                                       <  1:3205684/100‑1 (MQ=255)
                                               aaaaTATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAATTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATc                                                       <  2:1766656/99‑1 (MQ=255)
                                               aaaaTATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAACTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATCt                                                      >  2:1702317/1‑100 (MQ=255)
                                                  atatTCTTGATACAGTTGACTTTTCGCAGAGTCATATTGGTCATAGTAACTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATCttt                                                    <  2:819843/99‑1 (MQ=255)
                                                  atatTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAACTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATCTTTa                                                   <  1:725367/100‑1 (MQ=255)
                                                     ttCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAACTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATCTTTACGc                                                <  2:727324/100‑1 (MQ=255)
                                                     ttCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAACTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATCTTTACGc                                                >  1:2890755/1‑100 (MQ=255)
                                                       cTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAATTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATCTTTACGCtt                                              >  1:1587037/1‑100 (MQ=255)
                                                          gATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAATTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATCTTTACGCTTCGc                                           >  2:2749731/1‑100 (MQ=255)
                                                          gATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAATTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATCTTTACGCTTCGc                                           <  2:2939085/100‑1 (MQ=255)
                                                          gATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAACTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATCTTTACGCTTCGc                                           <  2:1152355/100‑1 (MQ=255)
                                                          gATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAACTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATCTTTACGCTTCGc                                           <  1:1657877/100‑1 (MQ=255)
                                                          gATACAGTTGACTCCTCGCAGAGTCAAATTGGTCATAGTAACTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATCTTTACGCTTCGc                                           <  2:1104552/100‑1 (MQ=255)
                                                            tACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAATTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATCTTTACGCTTCGCAc                                         >  1:2965719/1‑100 (MQ=255)
                                                            tACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAACTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATCTTTACGCTTCGCAc                                         <  2:1392851/100‑1 (MQ=255)
                                                            tACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAACTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATCTTTACGCTTCGCAc                                         <  1:3578402/100‑1 (MQ=255)
                                                              cAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAATTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATCTTTACGCTTCGCACtt                                       >  2:352131/1‑100 (MQ=255)
                                                              cAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAATTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATCTTTACGCTTCGCACtt                                       >  1:3347700/1‑100 (MQ=255)
                                                              cAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAACTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATCTTTACGCTTCGCACtt                                       >  1:1436737/1‑100 (MQ=255)
                                                                gTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAACTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATCTTTACGCTTCGCACTTCa                                     >  2:3335923/1‑100 (MQ=255)
                                                                  tGACTTCTCGCAGAGTCAACTTGGTCATAGTAATTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATCTTTACGCTTCGCACTTCAcc                                   <  2:2171825/100‑1 (MQ=255)
                                                                  tGACTTCTCGCAGAGTCAACTTGGTCATAGTAATTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATCTTTACGCTTCGCACTTCAcc                                   <  2:3224998/100‑1 (MQ=255)
                                                                  tGACTTCTCGCAGAGTCAACTTGGTCATAGTAACTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATCTTTACGCTTCGCACTTCAcc                                   <  1:1968337/100‑1 (MQ=255)
                                                                    aCTTCTCGCAGAGTCAACTTGGTCATAGTAACTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATCTTTACGCTTCGCACTTCACcag                                 <  2:3284008/100‑1 (MQ=255)
                                                                     cTTCTCGCAGAGTCAACTTGGTCATAGTAATTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATCTTTACGCTTCGCACTTCACcaga                                >  2:1304152/1‑100 (MQ=255)
                                                                     cTTCTCGCAGAGTCAACTTGGTCATAGTAATTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATCTTTACGCTTCGCACTTCACcaga                                >  1:2800636/1‑100 (MQ=255)
                                                                      ttCTCGCAGAGTCAACTTGGTCATAGTAATTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATCTTTACGCTTCGCACTTCACcagag                               <  1:3172914/100‑1 (MQ=255)
                                                                       tctcGCAGAGTCAACTTGGTCATAGTAATTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATCTTTACGCTTCGCACTTCACcagagc                              <  2:1918509/100‑1 (MQ=255)
                                                                        ctcGCAGAGTCAACTTGGTCATATGAACTTTCTCGAATGTCAGATTAGAAATTCAAACTGCGATAATTGTTATCTTTACGCTTCGCACTTCACcagagca                             >  2:3449087/1‑100 (MQ=255)
                                                                        ctcGCAGAGTCAACTTGGTCATAGTAATTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATCTTTACGCTTCGCACTTCACcagagca                             >  1:1837574/1‑100 (MQ=255)
                                                                         tcGCAGAGTCAACTTGGTCATAGTAATTTTCGCGAATGTCAGATTAGAAATTCGAACTTCGATAATTGTTATCTTTACGCTTCGCACTTCACcagagcag                            <  2:2284616/100‑1 (MQ=255)
                                                                          cGCAGAGTCAACTTGGTCATAGTAATTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATCTTTACGCTTCGCACTTCACcagagcaga                           >  2:2299750/1‑100 (MQ=255)
                                                                          cGCAGAGTCAACTTGGTCATAGTAATTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATCTTTACGCTTCGCACTTCACcagagcag                            >  1:823238/1‑99 (MQ=255)
                                                                           gCAGAGTCAACTTGGTCATAGTAATTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATCTTTACGCTTCGCACTTCACcagagcagag                          <  1:13625/100‑1 (MQ=255)
                                                                              gagTCAACTTGGTCATAGTAATTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATCTTTACGCTTCGCACTTCACCAGAGCAGAGttt                       <  1:3543136/100‑1 (MQ=255)
                                                                               agTCAACTTGGTCATAGTAATTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATCTTTACGCTTCGCACTTCACCAGAGCAGAGTTtc                      <  1:1156563/100‑1 (MQ=255)
                                                                                  cAACTTGGTCATAGTAATTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATCTTTACGCTTCGCACTTCACCAGAGCAGAGTTtctgt                   >  2:1575616/1‑100 (MQ=255)
                                                                                  cAACTTGGTCATAGTAATTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATCTTTACGCTTCGCACTTCACCAGACCAGAGTTtctgt                   >  1:761759/1‑100 (MQ=255)
                                                                                  cAACTTGGTCATAGTAACTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATCTTTACGCTTCGCACTTCACCAGAGCAGAGTTtctgt                   >  1:2457873/1‑100 (MQ=255)
                                                                                   aaCTTGGTCATAGTAATTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATCTTTACGCTTCGCACTTCACCAGAGCAGAGTTtctgtc                  <  2:1675967/100‑1 (MQ=255)
                                                                                     cTTGGTCATAGTAATTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATCTTTACGCTTCGCACTTCACCAGAGCAGAGTTtctgtctg                >  2:2697502/1‑100 (MQ=255)
                                                                                         gTCATAGTAACTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATCTTTACGCTTCGCACTTCACCAGAGCAGAGTTTCTGTCTGCCaa            >  1:2317621/1‑100 (MQ=255)
                                                                                           cATAGTAATTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATCTTTACGCTTCGCACTTCACCAGAGCAGAGTTTCTGTCTGCCAAAg          >  2:659887/1‑100 (MQ=255)
                                                                                           cATAGTAACTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATCTTTACGCTTCGCACTTCACCAGAGCAGAGTTTCTGTCTGCCAAAg          >  1:4019205/1‑100 (MQ=255)
                                                                                             tAGTAATTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATCTTTACGGTTCGCAATTCACCAGAGCAGAGTTTCTGTCTGCCAAAGaa        <  2:1039578/100‑1 (MQ=255)
                                                                                             tAGTAATTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATCTTTACGCTTCGCACTTCACCAGAGCAGAGTTTCTGTCTGCCAAAGaa        <  2:969419/100‑1 (MQ=255)
                                                                                              aGTAATTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATCTTTACGCTTCGCACTTCACCAGAGCAGAGTTTCTGTCTGCCAAAGaaa       <  1:2559694/100‑1 (MQ=255)
                                                                                              aGTAACTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATCTTTACGCTTCGCACTTCACCAGAGCAGAGTTTCTGTCTGCCAAAGaaa       >  2:2053655/1‑100 (MQ=255)
                                                                                              aGTAACTTGCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATCTTTACGCTTCGCACTTCACCAGAGCAGAGTTTCTGTCTGCCAAAGa         >  2:223704/1‑98 (MQ=255)
                                                                                                tAATTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATCTTTACGCTTCGCACTTCACCAGAGCAGAGTTTCTGTCTGCCAAAGAAata     <  2:2800620/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
CAATGTGATTTTTATCAAAGCGAGTGGTGAATGTGTCGATTTCAGCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAATTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATCTTTACGCTTCGCACTTCACCAGAGCAGAGTTTCTGTCTGCCAAAGAAATATCA  >  NC_000913/4251188‑4251386

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