Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 922,316:1 +C 15.0% 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.C15.0% 469.4 / 89.4 133intergenic (+21/+52)macB/cspDmacrolide ABC transporter peremase/ATPase/inhibitor of DNA replication, cold shock protein homolog
Reads supporting (aligned to +/- strand):  ref base . (52/61);  new base C (11/9);  total (63/70)
Fisher's exact test for biased strand distribution p-value = 4.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.06e-01

Only 100 of 138 total aligned reads displayed.

CACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTG  >  NC_000913/922215‑922413
                                                                                                    |                                                                                                   
cACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAg                                                                                                      <  1:697145/100‑1 (MQ=255)
cACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATACGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAg                                                                                                      >  1:1891264/1‑100 (MQ=255)
      gATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATc                                                                                               <  2:2549982/100‑1 (MQ=255)
            aTTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTg                                                                                         >  2:3345302/1‑100 (MQ=255)
             tttGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTgg                                                                                        <  2:2973849/100‑1 (MQ=255)
             tttGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTgg                                                                                        <  2:1116922/100‑1 (MQ=255)
              ttGGCTGGTTACCCGCACGAAATGCGTCACGACTGGTTCCAGTAGATTCTCTGGCACGAGAGGAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGc                                                                                       <  1:872471/100‑1 (MQ=255)
              ttGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGc                                                                                       <  1:575253/100‑1 (MQ=255)
              ttGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGc                                                                                       <  2:2497704/100‑1 (MQ=255)
                ggCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCCGATC‑GGCTGGCAt                                                                                     <  2:1248234/100‑1 (MQ=39)
                 gCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCCGATC‑GGCTGGCAtt                                                                                    <  1:1963354/100‑1 (MQ=39)
                 gCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCCGATC‑GGCTGGCAtt                                                                                    <  2:2452265/100‑1 (MQ=39)
                  cTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCAttt                                                                                   >  2:327108/1‑100 (MQ=255)
                  cTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCAttt                                                                                   >  1:3364190/1‑100 (MQ=255)
                  cTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAAGGCCAG‑CCGATCGGGCTGTCAtgt                                                                                   <  1:3323004/100‑3 (MQ=255)
                   tGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCAtttt                                                                                  <  1:2547825/100‑1 (MQ=255)
                    ggTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCCGATCGGG‑TGGCAtttt                                                                                  >  2:888755/1‑99 (MQ=39)
                    ggTTACCCGCACGAAATGCCGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGCAAATTCTGAGATAACAATGCCAG‑CCGAGCGGGCCGGCATTTTg                                                                                 >  1:991461/1‑100 (MQ=255)
                     gTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCCGATC‑GGCTGGCAttttta                                                                                <  1:1880006/100‑3 (MQ=38)
                     gTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGc                                                                                <  2:3359470/100‑1 (MQ=255)
                      ttACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGcc                                                                               <  1:1514972/100‑1 (MQ=255)
                      ttACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGcc                                                                               <  2:651286/100‑1 (MQ=255)
                       tACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCCGATC‑GGCTGGCAtttttatc                                                                              <  2:2345341/100‑5 (MQ=38)
                       tACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAAAGTAATTTTTGAGATAAAAATGCCAGCCCGATC‑TGCTGGCAtttttatc                                                                              <  2:2024382/100‑5 (MQ=37)
                        aCCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCCGATC‑GGCTGGCAtttttatct                                                                             <  2:3533851/100‑6 (MQ=37)
                          ccGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTa                                                                           >  1:2871388/1‑100 (MQ=255)
                          ccGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTa                                                                           >  1:1535649/1‑100 (MQ=255)
                           cGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAg                                                                          >  2:715682/1‑100 (MQ=255)
                            gctcGAAGTGGGGCACGACTGGATCCAGTAGATGCTCTAGCACGAGAGTAATTTCTGAGATAAAAATGCCAGCCCGATC‑GGCTGGCAtttttatctc                                                                            <  1:3235418/94‑7 (MQ=21)
                            cttccgatctcGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAgg                                                                         <  2:1107621/90‑1 (MQ=255)
                                gAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGt                                                                     >  1:2339524/1‑100 (MQ=255)
                                   aTGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGCAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTaca                                                                  >  1:1078494/1‑100 (MQ=255)
                                    tGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGGGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTacac                                                                 >  1:2779424/1‑100 (MQ=255)
                                    tGCGGCACGACTGGATACAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCGAC‑ACGATCGGGCTGGCATTTTGCCTTTAGGATGTacac                                                                 >  2:714665/1‑100 (MQ=255)
                                     gCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTAcaca                                                                >  2:1278225/1‑100 (MQ=255)
                                      cGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGagatcggaag                                                               >  1:1107614/1‑91 (MQ=255)
                                      cGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACaa                                                               <  2:2122369/100‑1 (MQ=255)
                                      cGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACaa                                                               >  2:80502/1‑100 (MQ=255)
                                       ggCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAAt                                                              <  1:180476/100‑1 (MQ=255)
                                        gCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATg                                                             >  1:61368/1‑100 (MQ=255)
                                        gCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATg                                                             >  2:1026146/1‑100 (MQ=255)
                                        gCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATg                                                             <  1:2111208/100‑1 (MQ=255)
                                         cACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATga                                                            >  2:2365073/1‑100 (MQ=255)
                                            gACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGAc                                                         <  2:1579434/100‑1 (MQ=255)
                                             aCTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGGCTTTAGGATGTACACAATGAGACa                                                        >  2:3014157/1‑100 (MQ=255)
                                             aCTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACa                                                        >  1:1478359/1‑100 (MQ=255)
                                              cTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACag                                                       <  1:1876214/100‑1 (MQ=255)
                                                ggATCCAGTAGATGCTCTGGCACGATAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTCGGATGTACACAATGAGAc                                                         <  1:890231/96‑1 (MQ=255)
                                                ggATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACagaa                                                     >  2:490899/1‑100 (MQ=255)
                                                  aTCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCGG‑CCGCTCGGGCTGGCATTTTGCTTTTAGGATGTACACAATGAGCCagaag                                                    >  1:2727549/1‑99 (MQ=255)
                                                  aTCCAGTAGATGCTCTGGCACGAGAGTAATGTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTATGATGTACACAATGAGACagaag                                                    >  1:2508812/1‑99 (MQ=255)
                                                   tCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAagag                                                  >  2:1800913/1‑100 (MQ=255)
                                                     cAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGcn                                                <  1:1640058/100‑2 (MQ=255)
                                                     cAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCt                                                >  2:39911/1‑100 (MQ=255)
                                                      aGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTa                                               <  2:37123/100‑1 (MQ=255)
                                                      aGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTa                                               >  1:1429315/1‑100 (MQ=255)
                                                         aGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATg                                             <  1:2013552/99‑1 (MQ=255)
                                                          gATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑TCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGCGCTATGCg                                           >  2:2778801/1‑100 (MQ=255)
                                                           aTGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCgn                                          <  1:2703365/100‑2 (MQ=255)
                                                            tGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGcaga                                         >  2:155954/1‑97 (MQ=255)
                                                              ctctGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACt                                        >  1:1229805/1‑99 (MQ=255)
                                                              ctctGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTg                                       <  2:325539/100‑1 (MQ=255)
                                                              ctctGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTg                                       >  1:504479/1‑100 (MQ=255)
                                                              ctctGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTg                                       <  1:52210/100‑1 (MQ=255)
                                                              ctctGGCACGAGAGTAATCTTTGCGATAAAAATGCCAG‑CCGATCGGGCTGGCCTTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGa                                          >  1:2627845/1‑97 (MQ=255)
                                                                  ggCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGc                                   >  2:475310/1‑100 (MQ=255)
                                                                  ggCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGc                                   <  2:2182681/100‑1 (MQ=255)
                                                                   gCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCt                                  <  2:1774281/100‑1 (MQ=255)
                                                                    cACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCtt                                 >  1:1754214/1‑100 (MQ=255)
                                                                    cACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCtt                                 >  2:325976/1‑100 (MQ=255)
                                                                      cGAGAGTAATTTTTTATTTATAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACAAAATGAGACAGAAGAGCTATGCGACTGCCGGTTCt                               >  1:1375626/1‑100 (MQ=255)
                                                                      cGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCt                               <  1:2307764/100‑1 (MQ=255)
                                                                      cGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCt                               <  1:2990045/100‑1 (MQ=255)
                                                                       gagaGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTa                              <  2:1543335/100‑1 (MQ=255)
                                                                        agagTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGTCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCtt                                 <  1:2655301/96‑1 (MQ=255)
                                                                         gagTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACt                            >  2:105892/1‑100 (MQ=255)
                                                                           gTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTc                          <  2:1136644/100‑1 (MQ=255)
                                                                            tAATTTTTGAGATAAAAATGCCAG‑TCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCg                         <  2:1491223/100‑1 (MQ=255)
                                                                            tAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCg                         <  1:3389916/100‑1 (MQ=255)
                                                                                  ttGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGc                   <  2:2809800/100‑1 (MQ=255)
                                                                                    gagaTAGAAATACCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGcac                 <  1:708217/100‑1 (MQ=255)
                                                                                    gagaTAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGcac                 <  1:399425/100‑1 (MQ=255)
                                                                                     agaTAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGcaca                >  1:2229082/1‑100 (MQ=255)
                                                                                       aTAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAAt              <  2:1429322/100‑1 (MQ=255)
                                                                                        gcaaaaTGCCAGCCCGATC‑GGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATa             >  1:531485/3‑100 (MQ=38)
                                                                                          aaaaTGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAac           >  2:2918690/1‑100 (MQ=255)
                                                                                           aaaTGCCAGCCCGATC‑GGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAaca          >  2:666519/1‑100 (MQ=39)
                                                                                            aaTGCCAGCCCGATC‑GGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAacac         >  2:717093/1‑100 (MQ=39)
                                                                                            aaTGCCAGCCCGATC‑GGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAacac         >  1:1509974/1‑100 (MQ=39)
                                                                                            aaTGCCAGCCCGATC‑GGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAacac         >  2:2135567/1‑100 (MQ=39)
                                                                                            aaTGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAacac         <  2:1304120/100‑1 (MQ=255)
                                                                                               gCCAGCCCGATC‑GGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTgg      >  2:48957/1‑100 (MQ=255)
                                                                                               gCCAGCCCGATC‑GGCTGGCATTGTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTg       >  1:2065314/1‑99 (MQ=255)
                                                                                                ccAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGc     >  2:2023534/1‑100 (MQ=255)
                                                                                                ccAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGc     >  2:1956200/1‑100 (MQ=255)
                                                                                                 cAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCg    <  1:1348937/100‑1 (MQ=255)
                                                                                                 cAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCg    <  1:1996473/100‑1 (MQ=255)
                                                                                                  aG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCg    <  1:670464/99‑1 (MQ=255)
                                                                                                  aG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGt   >  1:1339523/1‑100 (MQ=255)
                                                                                                   g‑CCGATCGGGCTGGCATTTTGGCTTTGGGGTGTACTCAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTg  <  1:243310/100‑1 (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: ‑