Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 922,316:1 +C 18.1% 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.C18.1% 500.7 / 128.0 160intergenic (+21/+52)macB/cspDmacrolide ABC transporter peremase/ATPase/inhibitor of DNA replication, cold shock protein homolog
Reads supporting (aligned to +/- strand):  ref base . (62/69);  new base C (17/12);  total (79/81)
Fisher's exact test for biased strand distribution p-value = 3.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.30e-01

Only 100 of 167 total aligned reads displayed.

CACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGT  >  NC_000913/922215‑922412
                                                                                                    |                                                                                                  
 aCCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑c                                                                                                   >  1:2550857/1‑100 (MQ=255)
  ccGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑cc                                                                                                  >  2:3240273/1‑100 (MQ=255)
  ccGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑cc                                                                                                  <  1:2242933/100‑1 (MQ=255)
   cGGGATTTTATTTGGCTGGTTACCAGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑c                                                                                                   <  1:2138397/98‑1 (MQ=255)
    gggATTTTAGGTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑cc                                                                                                  <  2:3420587/98‑1 (MQ=255)
     ggATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGAt                                                                                               <  2:1102974/100‑1 (MQ=255)
      gATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATc                                                                                              >  2:1715610/1‑100 (MQ=255)
       aTTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCg                                                                                             >  1:1931528/1‑100 (MQ=255)
       aTTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCCGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCg                                                                                             <  2:3804874/100‑1 (MQ=255)
        ttttATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCgg                                                                                            >  1:2032508/1‑100 (MQ=255)
        ttttATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCgg                                                                                            >  2:1058926/1‑100 (MQ=255)
          ttATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGc                                                                                          <  2:3577760/100‑1 (MQ=255)
          ttATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGc                                                                                          >  1:2769175/1‑100 (MQ=255)
           tATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCt                                                                                         >  2:2691083/1‑100 (MQ=255)
           tATTTGGCTGGGTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCt                                                                                         <  1:3313348/100‑1 (MQ=255)
             tttGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTgg                                                                                       <  2:699934/100‑1 (MQ=255)
             tttGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTgg                                                                                       <  1:1652914/100‑1 (MQ=255)
             tttGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTgg                                                                                       <  2:3411725/100‑1 (MQ=255)
             tttGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTgg                                                                                       <  2:746214/100‑1 (MQ=255)
              ttGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGc                                                                                      <  1:2218420/100‑1 (MQ=255)
               tGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CTGATCGGGCTGGCa                                                                                     <  2:2441720/100‑1 (MQ=255)
               tGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCa                                                                                     <  2:2261281/100‑1 (MQ=255)
                ggCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCCGATC‑GGCTGGCAt                                                                                    <  1:1538062/100‑1 (MQ=39)
                 gCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCCGATC‑GGCTGGCAtt                                                                                   <  1:3007425/100‑1 (MQ=39)
                 gCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCCGATC‑GGCTGGCAtt                                                                                   <  1:472132/100‑1 (MQ=39)
                 gCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCCGATC‑GGCTGGCAtt                                                                                   <  1:356717/100‑1 (MQ=39)
                 gCTGGTTAACCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCCGATC‑GGCTGGCAtt                                                                                   <  2:887213/100‑1 (MQ=38)
                  cTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCAttt                                                                                  >  1:3731746/1‑100 (MQ=255)
                   tGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCCGATC‑GGCTGGCAtttt                                                                                 <  1:1467168/100‑1 (MQ=39)
                    ggTTACCCGCACGAACTGCAGCACGATTGGATCCAGTAGATGCTCTGGCACGATAATAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTg                                                                                <  2:415732/100‑1 (MQ=255)
                    ggTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCCGATC‑GGCTGGCAttttt                                                                                <  1:295981/100‑2 (MQ=38)
                       tACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCt                                                                             <  1:1123862/100‑1 (MQ=255)
                       tACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCt                                                                             <  1:2170222/100‑1 (MQ=255)
                       tACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCt                                                                             <  1:1715619/100‑1 (MQ=255)
                       tACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCt                                                                             <  1:2604099/100‑1 (MQ=255)
                         cccGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCttt                                                                           >  2:1297301/1‑100 (MQ=255)
                          ccGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTa                                                                          <  2:1407257/100‑1 (MQ=255)
                           cGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAg                                                                         >  1:1037662/1‑100 (MQ=255)
                           cGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAg                                                                         >  1:2722316/1‑100 (MQ=255)
                             cACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAgg                                                                        <  2:1764841/99‑1 (MQ=255)
                             cACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGa                                                                       <  2:2308996/100‑1 (MQ=255)
                                gAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGt                                                                    >  2:135422/1‑100 (MQ=255)
                                   aTGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTaca                                                                 >  1:1412990/1‑100 (MQ=255)
                                     gCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTAcaca                                                               <  1:2691100/100‑1 (MQ=255)
                                      cGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACaa                                                              >  1:1322876/1‑100 (MQ=255)
                                       ggCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAAt                                                             <  2:249153/100‑1 (MQ=255)
                                       ggCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAAt                                                             <  1:232605/100‑1 (MQ=255)
                                        catctaCTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATg                                                            >  1:1228408/6‑100 (MQ=255)
                                         cACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTCCACAATga                                                           >  1:2342981/1‑100 (MQ=255)
                                            gACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGAc                                                        <  2:3695100/100‑1 (MQ=255)
                                            gACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGAc                                                        >  1:1552986/1‑100 (MQ=255)
                                            gACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGAc                                                        <  2:2237939/100‑1 (MQ=255)
                                              cTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACag                                                      >  1:3127635/1‑100 (MQ=255)
                                                ggATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACagaa                                                    >  1:3491973/1‑100 (MQ=255)
                                                  aTCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACagaaga                                                  >  2:2002232/1‑100 (MQ=255)
                                                   tCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAagag                                                 >  2:1869315/1‑100 (MQ=255)
                                                    ccAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGTATGTACACAATGAGACAGAAGAGc                                                <  1:544369/100‑1 (MQ=255)
                                                    ccAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGc                                                <  1:1824230/100‑1 (MQ=255)
                                                         aGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGc                                           >  2:1235012/1‑100 (MQ=255)
                                                         aGATGCTCTGGAACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATg                                            <  2:1928354/99‑1 (MQ=255)
                                                            tGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGAc                                        <  2:3752387/100‑1 (MQ=255)
                                                            tGCGCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGAc                                        <  2:3705772/100‑1 (MQ=255)
                                                             gCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACt                                       >  1:972351/1‑100 (MQ=255)
                                                             gCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTATAGGAGGTACACAAGGAGACAGAAGAGCTATGCGACt                                       >  2:2659732/1‑100 (MQ=255)
                                                              ctctGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTg                                      <  2:2902502/100‑1 (MQ=255)
                                                                ctGGCANGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGcc                                    >  2:3257178/1‑100 (MQ=255)
                                                                ctGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTGTGCGACTGcc                                    <  2:2083218/100‑1 (MQ=255)
                                                                  ggCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGc                                  >  2:2051036/1‑100 (MQ=255)
                                                                   gCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCt                                 <  1:3050133/100‑1 (MQ=255)
                                                                      cGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCt                              >  2:2719643/1‑100 (MQ=255)
                                                                      cGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCt                              >  1:767575/1‑100 (MQ=255)
                                                                      cGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCt                              <  1:1297310/100‑1 (MQ=255)
                                                                       gagaGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCGTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTa                             <  2:1943445/100‑1 (MQ=255)
                                                                         gagTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACt                           >  1:1446384/1‑100 (MQ=255)
                                                                          agTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACtt                          <  1:2528474/100‑1 (MQ=255)
                                                                            tAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCg                        >  1:2787326/1‑100 (MQ=255)
                                                                               tttttGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGCAGAGCTATGCGACTGCCGCTTCTACTTCGACg                     >  1:982271/1‑100 (MQ=255)
                                                                                 tttGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACggg                   >  1:877148/1‑100 (MQ=255)
                                                                                 cctaaaggcaaaaTGCCAGCCCGATC‑GGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACggg                   >  2:2601272/10‑100 (MQ=25)
                                                                                  ttGGGATAAAAATGCCAG‑CCGATCGGACTAGCATTTGGCCTTTACGATGTGAACAATGAGACACCAGCGCTATGGGACTGCCGCCTCTACTTAGACGAGc                  >  2:3217447/1‑100 (MQ=255)
                                                                                  ctaaaggcaaaaTGCCAGCCCGATC‑GGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGc                  >  1:3636744/9‑100 (MQ=25)
                                                                                    gagaTAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGcac                >  1:3463855/1‑100 (MQ=255)
                                                                                     agaTAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGcaca               <  1:3325502/100‑1 (MQ=255)
                                                                                      gaTAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACaa              >  2:1044111/1‑100 (MQ=255)
                                                                                        tAAAAATTCCAG‑CCGATCGGGCTGGCATTTGGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGcac                <  2:951647/96‑1 (MQ=255)
                                                                                         caaaaTGCCAGCCCGATC‑GGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATaa           >  2:3097837/2‑100 (MQ=38)
                                                                                         caaaaTGCCAGCCCGATC‑GGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATaa           >  2:1576293/2‑100 (MQ=38)
                                                                                         caaaaTGCCAGCCCGATC‑GGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATaa           >  1:2974207/2‑100 (MQ=38)
                                                                                         caaaaTGCCAGCCCGATC‑GGCTGGCATTTTGCCTTTAGGATGTAAACAATGAGACAGAAGGGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATaa           >  1:3929164/2‑100 (MQ=38)
                                                                                          aaaaTGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAac          >  2:3797039/1‑100 (MQ=255)
                                                                                          aaaaTGCCAG‑CCGATCGGGCTGGCACTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAac          <  2:1718458/100‑1 (MQ=255)
                                                                                            aaTGCCAGCCCGATC‑GGCTGGCAGTTTGCCTTTAGGATGTACACAATGAGGCGGAAGAGCTATGCGACTGCCGCTTCTGCTTCTACGGGCACAATAacac        >  2:3331820/1‑100 (MQ=255)
                                                                                              tGCCAGCCCGATC‑GGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTg      >  2:1535488/1‑100 (MQ=255)
                                                                                              tGCCAGCCCGATC‑GGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTg      >  1:1732837/1‑100 (MQ=255)
                                                                                              tGCCAGCCCGATC‑GGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTg      >  1:1078032/1‑100 (MQ=255)
                                                                                              tGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTg      >  1:436300/1‑100 (MQ=255)
                                                                                               gCCAGCCCGATC‑GGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTgg     >  2:1528416/1‑100 (MQ=255)
                                                                                               gCCAGCCCGATC‑GGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTgg     >  1:3491770/1‑100 (MQ=255)
                                                                                                 cAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCg   <  1:3533260/100‑1 (MQ=255)
                                                                                                  aG‑CCGATAGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGt  >  1:2067238/1‑100 (MQ=255)
                                                                                                    |                                                                                                  
CACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGT  >  NC_000913/922215‑922412

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