Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 922,316:1 +C 17.5% intergenic (+23/+50) macB → / ← cspD macrolide ABC transporter peremase/ATPase/inhibitor of DNA replication, cold shock protein homolog

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913922,3141.C17.5% 368.0 / 84.2 114intergenic (+21/+52)macB/cspDmacrolide ABC transporter peremase/ATPase/inhibitor of DNA replication, cold shock protein homolog
Reads supporting (aligned to +/- strand):  ref base . (46/48);  new base C (9/11);  total (55/59)
Fisher's exact test for biased strand distribution p-value = 8.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.89e-01

Only 100 of 121 total aligned reads displayed.

CACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTG  >  NC_000913/922215‑922413
                                                                                                    |                                                                                                   
cACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAg                                                                                                      <  1:1996416/100‑1 (MQ=255)
cACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAg                                                                                                      >  1:2657139/1‑100 (MQ=255)
  ccGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑cc                                                                                                   <  2:33933/100‑1 (MQ=255)
   cGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGTATCCAGTAGATGCGCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑c                                                                                                    <  2:685016/98‑1 (MQ=255)
   cGGGATTTTATTTGGCGGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑cc                                                                                                   <  2:2064885/99‑1 (MQ=255)
     ggATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGAt                                                                                                >  1:233029/1‑100 (MQ=255)
      gATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATc                                                                                               <  2:1603478/100‑1 (MQ=255)
       aTTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCg                                                                                              <  2:42789/100‑1 (MQ=255)
        ttttttttGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑GCGATCg                                                                                              <  2:2983098/99‑1 (MQ=255)
        ttttATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCgg                                                                                             <  2:2160205/100‑1 (MQ=255)
        ttttATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCgg                                                                                             >  2:2426863/1‑100 (MQ=255)
         tttATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCgg                                                                                             >  1:3401585/1‑99 (MQ=255)
          ttATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGc                                                                                           >  2:2876919/1‑100 (MQ=255)
            aTTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTg                                                                                         <  1:3596098/100‑1 (MQ=255)
             tttGGCTGGTTACCCGCACGCAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTgg                                                                                        <  2:1372913/100‑1 (MQ=255)
              ttGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGc                                                                                       <  1:2980358/100‑1 (MQ=255)
               tGGCTGGTTCCCCGCACGAAATGCGGCACTATTGGCTCCAGAAGGTGCGCTGGCACGAGAGGAATTTTTGAGATAAAAATGCCAG‑CCGATCggagtgg                                                                                        >  2:1155396/1‑93 (MQ=255)
               tGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCa                                                                                      >  1:2776701/1‑100 (MQ=255)
               tGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCa                                                                                      >  1:3387998/1‑100 (MQ=255)
               tGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCa                                                                                      >  1:407916/1‑100 (MQ=255)
               tGGCTAGTTACCCGCACGAAATGCGGCACTACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCCGATC‑GGCTGGCa                                                                                      <  1:3233729/100‑1 (MQ=255)
                ggCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCCGATC‑GGCTGGCAt                                                                                     <  2:3163614/100‑1 (MQ=39)
                 gCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCCGATC‑GGCTGGCAtt                                                                                    <  1:2121665/100‑1 (MQ=39)
                 gCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCAtt                                                                                    <  1:633804/100‑1 (MQ=255)
                   tGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCAtttt                                                                                  <  1:3021816/100‑1 (MQ=255)
                   tGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCAtttt                                                                                  <  1:1772256/100‑1 (MQ=255)
                    ggTTACCCGCACGAAATGCGGCACTACTGTATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCCGATC‑GGCTGGCAttttt                                                                                 <  2:3701782/100‑2 (MQ=38)
                     gTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGc                                                                                >  2:1563316/1‑100 (MQ=255)
                       tACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCCGATC‑GGCTGGCAtttttatc                                                                              <  2:1917382/100‑5 (MQ=38)
                       tACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCCGATC‑GGCTGGCAtttttatc                                                                              <  1:3590213/100‑5 (MQ=38)
                         cccGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTTTGTTTTTATTGCCAG‑CCGATCGGGCTGGCATTTTGCCttt                                                                            >  2:455203/1‑100 (MQ=255)
                         cccGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCCGATC‑GGCTGGCAtttttatct                                                                             <  1:1616531/99‑6 (MQ=37)
                         cccGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCttt                                                                            >  1:400462/1‑100 (MQ=255)
                         cccGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCttt                                                                            <  1:3715564/100‑1 (MQ=255)
                           cGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCGtttttatctcaa                                                                          >  1:3298851/1‑92 (MQ=255)
                           cGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAg                                                                          >  2:459995/1‑100 (MQ=255)
                           cGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAg                                                                          <  2:2070130/100‑1 (MQ=255)
                           cGCACGAAATGCCGGACGACTGTATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGTGATAAAAATGCCAGCCCGATC‑GGCTGGCAtttttatc                                                                              <  2:745490/96‑5 (MQ=25)
                            gCAGGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCttt                                                                            <  2:3035540/97‑1 (MQ=255)
                            gCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGCTCGGGCTGGCATTTTGCCTTTAgg                                                                         >  2:253250/1‑100 (MQ=255)
                            gCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAgg                                                                         <  2:3518734/100‑1 (MQ=255)
                            gCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAgg                                                                         >  1:3205033/1‑100 (MQ=255)
                            gCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAgg                                                                         <  2:751063/100‑1 (MQ=255)
                             cACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGa                                                                        <  1:172521/100‑1 (MQ=255)
                                gAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATg                                                                      <  2:3662316/99‑1 (MQ=255)
                                     gCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTAcaca                                                                >  1:137607/1‑100 (MQ=255)
                                       ggCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAAt                                                              >  1:2813456/1‑100 (MQ=255)
                                        gCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATg                                                             >  2:1945664/1‑100 (MQ=255)
                                        gCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATg                                                             <  2:468557/100‑1 (MQ=255)
                                          aCGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATgag                                                           <  2:2949390/100‑1 (MQ=255)
                                           cGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATgaga                                                          >  1:3435585/1‑100 (MQ=255)
                                           cGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATgaga                                                          >  1:2137547/1‑100 (MQ=255)
                                            gACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGAc                                                         <  1:1774662/100‑1 (MQ=255)
                                            gACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGAc                                                         <  1:425088/100‑1 (MQ=255)
                                            gACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGAc                                                         <  1:2762315/100‑1 (MQ=255)
                                             aCTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTGTAGGATGTACACAATGAGACa                                                        >  1:1484362/1‑100 (MQ=255)
                                              cTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACag                                                       >  2:1294692/1‑100 (MQ=255)
                                                ggATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACagaa                                                     >  2:1954802/1‑100 (MQ=255)
                                                 gATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACagaag                                                    >  1:1013716/1‑100 (MQ=255)
                                                  aTCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACagaaga                                                   <  1:98816/100‑1 (MQ=255)
                                                  aTCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACagaaga                                                   <  1:2121060/100‑1 (MQ=255)
                                                     cAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCt                                                >  1:152044/1‑100 (MQ=255)
                                                     cAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCt                                                >  2:1295482/1‑100 (MQ=255)
                                                          gATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCg                                           <  2:1861635/100‑1 (MQ=255)
                                                            tGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGAc                                         <  2:2711499/100‑1 (MQ=255)
                                                               tctGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGc                                      <  2:847506/100‑1 (MQ=255)
                                                                  ggCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGcc                                     >  2:888990/1‑98 (MQ=255)
                                                                   gCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCt                                  <  1:6680/100‑1 (MQ=255)
                                                                    cACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCtt                                 <  1:2756814/100‑1 (MQ=255)
                                                                      cGAGAGTAATTTTTTAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCAGCtttt                               >  2:867234/1‑98 (MQ=255)
                                                                      cGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCt                               <  1:2490765/100‑1 (MQ=255)
                                                                      cGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCt                               <  2:2794972/100‑1 (MQ=255)
                                                                      cGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCt                               <  1:2923099/100‑1 (MQ=255)
                                                                      cGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCt                               >  2:174927/1‑100 (MQ=255)
                                                                          agTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGTGACTGCCGCTTCTACtt                           >  1:1810065/1‑100 (MQ=255)
                                                                           gTAATTTTTGAGATAAAAATGCCAG‑CCGATCTGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTc                          <  1:2540984/100‑1 (MQ=255)
                                                                           gTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTc                          <  1:2118229/100‑1 (MQ=255)
                                                                              aTTTTTGTGATAAATATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGAc                       >  2:620628/1‑100 (MQ=255)
                                                                              aTTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGAc                       <  1:1019264/100‑1 (MQ=255)
                                                                                  ctaaaggcaaaaTGCCAGCCCGATC‑GGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGc                   >  1:312855/9‑100 (MQ=25)
                                                                                      gaTAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCGACTTCGACGGGCACaa               >  2:2226259/1‑100 (MQ=255)
                                                                                      gaTAAAAATGCCAG‑CCGATCGGGCTGGCATTTTACCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACaa               >  1:1574818/1‑100 (MQ=255)
                                                                                        tAAAAATGCCAG‑TCGAGCGGGCTGGCATTTTGCTTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGATTCTACTTCGACGGGCACAATa             <  2:164884/100‑1 (MQ=255)
                                                                                        tAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATa             >  2:3167574/1‑100 (MQ=255)
                                                                                        tAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATa             >  1:1882744/1‑100 (MQ=255)
                                                                                         caaaaTGCCAGCCCGATC‑GGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATaa            >  1:3547905/2‑100 (MQ=38)
                                                                                         caaaaTGCCAGCCCGATC‑GGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATaa            >  2:3025354/2‑100 (MQ=38)
                                                                                         caaaaTGCCAGCCCGAAC‑GGCTGGCATTTTACCTTTAGGATGTACAAAATGAGACACAAGAGCTATGCGACTGCCGCGTCTACTTCGACGGGGACAATaa            >  2:3303227/2‑100 (MQ=25)
                                                                                         aataaTGCCAG‑CCGATCGGGCTGGCATTGTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAAt              <  2:1040954/95‑1 (MQ=255)
                                                                                         aaaaaTGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATaa            >  2:822661/1‑100 (MQ=255)
                                                                                          aaaaTGCCAGCCCGATC‑GGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAac           >  1:2443010/1‑100 (MQ=39)
                                                                                          aaaaTGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAac           <  2:2244905/100‑1 (MQ=255)
                                                                                          aaaaTGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAac           <  2:928678/100‑1 (MQ=255)
                                                                                          aaaaTGCCAG‑CAGATCGGGCTGGCATTTTGCCTTTAGGATATACACAATGAGACAGAAGAGCTATGCGACAGCCGCTTCTACTTCGACAGACACAATa             >  2:518615/1‑98 (MQ=255)
                                                                                            aaTGCCAGCCCGATC‑GGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGCCGGGCACAATAaca          >  2:3059799/1‑99 (MQ=38)
                                                                                            aaTGCCAGCCCGATC‑GGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAacac         >  2:2639149/1‑100 (MQ=39)
                                                                                            aaTGCCAGCCCGATC‑GGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAacac         >  1:167002/1‑100 (MQ=39)
                                                                                             aTGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACAAAATGAGACAGAAGAGCGATGCGACTGCCGCTTGTACTTCGACGGGCACAATAACACt        >  1:1874832/1‑100 (MQ=255)
                                                                                                   g‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTg  <  1:1312492/100‑1 (MQ=255)
                                                                                                   g‑CCGATCGGGCTGCCATTTTGCCTTTAGGATGTACACAGTGAGACAGAAGAGCTATGCGACTGCCCCTTCTACTTCGACGGGCACAATAac           >  2:2356175/1‑91 (MQ=255)
                                                                                                    |                                                                                                   
CACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTG  >  NC_000913/922215‑922413

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