Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 922,316:1 +C 18.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.C18.5% 545.8 / 141.8 178intergenic (+21/+52)macB/cspDmacrolide ABC transporter peremase/ATPase/inhibitor of DNA replication, cold shock protein homolog
Reads supporting (aligned to +/- strand):  ref base . (73/72);  new base C (16/17);  total (89/89)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 185 total aligned reads displayed.

CACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑C‑CGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGT  >  NC_000913/922215‑922412
                                                                                                    |                                                                                                   
cACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAg                                                                                                      >  1:3892672/1‑100 (MQ=255)
cACCGGGATTTTATGTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAg                                                                                                      <  1:106911/100‑1 (MQ=255)
 aCCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑c                                                                                                    <  2:3801801/100‑1 (MQ=255)
  ccGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑c‑c                                                                                                  >  2:2288841/1‑100 (MQ=255)
    gggATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑C‑CGa                                                                                                >  2:744164/1‑100 (MQ=255)
      gATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑C‑CGa                                                                                                >  1:3448416/1‑98 (MQ=255)
       aTTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑C‑CGATc                                                                                              <  1:696963/99‑1 (MQ=255)
       aTTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑C‑CGATCg                                                                                             <  2:2209095/100‑1 (MQ=255)
          ttATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑C‑CGATCGGGc                                                                                          <  2:3495003/100‑1 (MQ=255)
             tttGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑C‑CGATCGGGCTgg                                                                                       >  1:108707/1‑100 (MQ=255)
             tttGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑C‑CGATCGGGCTgg                                                                                       <  1:1655931/100‑1 (MQ=255)
             tttGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑C‑CGATCGGGCTgg                                                                                       >  2:1568083/1‑100 (MQ=255)
             tttGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑C‑CGATCGGGCTgg                                                                                       >  1:3838279/1‑100 (MQ=255)
              tttgCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAATATGCCAG‑C‑CGATCGGGCTgg                                                                                       <  1:2695965/96‑1 (MQ=255)
              ttGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCC‑CGATC‑GGCTGGc                                                                                      <  2:3884416/100‑1 (MQ=255)
                ggCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCC‑CGATC‑GGCTGGCAt                                                                                    <  2:3324476/100‑1 (MQ=39)
                 gCTTGTTACCCTCACGTAATTCGGCACGACTGGATCTAGTAGATTCGCTTGCAAGAGATTAATTTTTTAGATAAAAATGTCAG‑c‑ctctcg                                                                                             <  1:3487893/90‑6 (MQ=255)
                 gCTGGTGACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCC‑CGATC‑GGATGGCAtt                                                                                   <  1:103956/100‑1 (MQ=255)
                   tGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑C‑CGATCGGGCTGGCAtttt                                                                                 <  1:1202185/100‑1 (MQ=255)
                    ggTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCC‑CGATC‑GGCTGGCAttttt                                                                                <  2:478505/100‑2 (MQ=38)
                    ggTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCC‑CGATC‑GGCTGGCAttttt                                                                                <  2:2130991/100‑2 (MQ=38)
                    ggTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑C‑CGATCGGGCTGGCATTTTg                                                                                <  1:4319757/100‑1 (MQ=255)
                     gTTACCCGCACGAAATGCGGCACGACTGTATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCC‑CGATC‑GGCTGGCAttttta                                                                               <  2:1468867/100‑3 (MQ=38)
                      ttACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑C‑CGATCGGGCTGGCATTTTGcc                                                                              <  2:3833329/100‑1 (MQ=255)
                      ttACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑C‑CGATCGGGCTGGCATTTTGcc                                                                              <  2:402963/100‑1 (MQ=255)
                      ttACCCGCACGAAATGCGGCACGACTGGATCCAGGAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑C‑CGATCGGGCTGGCATTTTGcc                                                                              <  2:2055883/100‑1 (MQ=255)
                      tgaCCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑C‑CGATCGGGCTGGCATTTTGcc                                                                              <  1:941963/98‑1 (MQ=255)
                           cgtaCGAAATGCCGTACGACTGCATCCAGTAGATGCTTTGCCACGAGAGTCATTTTTGAGATAAAAATGCCTG‑CACGATC‑GGCTGGCATTTTgatat                                                                            <  1:3544333/94‑5 (MQ=255)
                           cGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGGGAGTAATTTTTGAGATAAAAATGCCAG‑C‑CGATCGGGCTGGCATTTTGCCTTTAg                                                                         >  2:2545883/1‑100 (MQ=255)
                           cGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑C‑CGATCGGGCTGGCATTTTGCCTTTAg                                                                         >  2:1346828/1‑100 (MQ=255)
                           cGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCCGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑C‑CGATCGGGCTGGCATTTTGCCTTTAg                                                                         >  1:1522337/1‑100 (MQ=255)
                            gCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑C‑CGATCGGGCTGGCATTTTGCCTTTAgg                                                                        <  1:355641/100‑1 (MQ=255)
                             cACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCC‑CGATC‑GGCTGGCAtttttatctcaaaa                                                                       <  2:2952737/100‑11 (MQ=25)
                             cACGAAATGCGGCACGACTGGATCCAGGAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑C‑CGATCGGGCTGGCATTTTGCCTTTAGGa                                                                       <  2:3157181/100‑1 (MQ=255)
                              aCGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑C‑CGATCGGGCTGGCATTTTGCCTTTAGGAt                                                                      >  1:3930418/1‑100 (MQ=255)
                               cGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑C‑CGATCGGGCTGGCATTTTGCCTTTAgg                                                                        <  1:1821061/97‑1 (MQ=255)
                               cGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑C‑CGATCGGGCTGGCATTTTGCCTTTAGGa                                                                       >  1:1958890/1‑98 (MQ=255)
                                    tGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑C‑CGATCGGGCTGGCATTTTGCCTTTAGGATGTacac                                                                <  2:1251767/100‑1 (MQ=255)
                                     gCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑C‑CGATCGGGCTGGCATTTTGCCTTTAGGATGTAcaca                                                               <  1:413353/100‑1 (MQ=255)
                                      cGGCACGACTGGATCCAGTAGATGCTCTGGCACGATAGTAATTTTTGAGATAAAAATGCCAG‑C‑CGATCGGGCTGGCATTTTGCCTTTAGGATGTAcaca                                                               <  1:1986271/99‑1 (MQ=255)
                                        gCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑C‑CGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAAt                                                             >  1:2752846/1‑99 (MQ=255)
                                         cACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTTTTTTAATTATGCAAG‑C‑CGATCGGGCTGGCATTTTGCCTTTAGGATGTAGACAATga                                                           >  1:2115546/1‑100 (MQ=255)
                                         cACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑C‑CGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATga                                                           <  1:3397792/100‑1 (MQ=255)
                                         cACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑C‑CGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATga                                                           <  1:2933392/100‑1 (MQ=255)
                                                ggATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTTTTTTTTTATTGCCAG‑C‑CGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACag                                                      >  1:356705/1‑98 (MQ=255)
                                                ggATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑C‑CGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACagaa                                                    >  1:44609/1‑100 (MQ=255)
                                                 gATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑C‑CGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACagaag                                                   <  2:2092084/100‑1 (MQ=255)
                                                  aTCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑C‑CGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACagaaga                                                  >  1:2425485/1‑100 (MQ=255)
                                                        tAGATGCTCTGGCACGGGAGTAATTTTTGAGATAAAAATGCCAG‑C‑CGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATg                                            <  1:154809/100‑1 (MQ=255)
                                                        tAGATGCGCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑C‑CGATCGGGCTGGGATTTTGCCTTTACGATGTACACAATGAGACAGAAGAGCTATg                                            <  1:3228252/100‑1 (MQ=255)
                                                         aGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑C‑CGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGc                                           <  2:487743/100‑1 (MQ=255)
                                                            tGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑C‑CGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGAc                                        >  2:4087395/1‑100 (MQ=255)
                                                             gCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑C‑CGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACt                                       <  2:3632520/100‑1 (MQ=255)
                                                               tctGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑C‑CGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACt                                       >  1:930019/1‑98 (MQ=255)
                                                               tctGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑C‑CGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGc                                     >  1:4032089/1‑100 (MQ=255)
                                                                ctGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑C‑CGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGcc                                    >  2:2126923/1‑100 (MQ=255)
                                                                 tGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑C‑CGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCg                                   <  2:2344031/100‑1 (MQ=255)
                                                                   gCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑C‑CGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCt                                 <  1:2795100/100‑1 (MQ=255)
                                                                   gCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑C‑CGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCt                                 <  1:231982/100‑1 (MQ=255)
                                                                    cACGAGAGTAATTTTTGAGATAAAAATGCCAG‑C‑CGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCtt                                >  1:474036/1‑100 (MQ=255)
                                                                    cACGAGAGTAATTTTTGAGATAAAAATGCCAG‑C‑CGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCtt                                >  1:409511/1‑100 (MQ=255)
                                                                     aCGAGAGTAATTTTTGAGATAAAAATGCAAG‑C‑CGATCGGGGTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTc                               <  2:1746066/100‑1 (MQ=255)
                                                                      cGAGAGTAATTTTTGAGATAAAAATGCCAG‑C‑CGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCt                              <  1:122820/100‑1 (MQ=255)
                                                                       gagaGTAATTTTTGAGATAAAAATGCCAG‑C‑CGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTa                             >  2:604301/1‑100 (MQ=255)
                                                                        agagTAATTTTTGAGATAAAAATGCCAG‑C‑CGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTAc                            <  2:181717/100‑1 (MQ=255)
                                                                          agTAATTTTTGAGATAAAAATGCCAG‑C‑CGATCGTGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACtt                          <  2:473038/100‑1 (MQ=255)
                                                                          agTAATTTTTGAGATAAAAATGCCAG‑C‑CGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACtt                          >  2:2392882/1‑100 (MQ=255)
                                                                          agTAATTTTTGAGATAAAAATGCCAG‑C‑CGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACtt                          <  2:2157630/100‑1 (MQ=255)
                                                                           gTAATTTTTGAGATAAAAATGCCAG‑C‑CGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTCCTTc                         >  2:3772676/1‑100 (MQ=255)
                                                                           gTAATTTTTGAGATAAAAATGCCAG‑C‑CGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTc                         >  2:3989235/1‑100 (MQ=255)
                                                                               tttttGAGATAAAAATGCCAG‑C‑CGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACg                     >  2:3563221/1‑100 (MQ=255)
                                                                               tttttGAGATAAAAATGCCAG‑C‑CGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACg                     >  2:882073/1‑100 (MQ=255)
                                                                                ttttGATATAAAAATGCCAG‑C‑CGATCGGGCTCGCATTTTTCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACg                     <  1:599012/99‑1 (MQ=255)
                                                                                ttttGAGATAAAAATGCCAG‑C‑CGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACg                     <  1:2545870/99‑1 (MQ=255)
                                                                                  ttGAGATAAAAATGCCAG‑C‑CGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGc                  <  2:3395902/100‑1 (MQ=255)
                                                                                   tGAGATAAAAATGCCAG‑C‑CGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGca                 <  1:4224981/100‑1 (MQ=255)
                                                                                    gagaTAAAAATGCCAG‑C‑CGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGCCAGAAGAGCTATGCGACTGCCGCTTCTAATTCGACgg                    >  1:206429/1‑96 (MQ=255)
                                                                                    gagaTAAAAATGCCAG‑C‑CGATCGGGCTGGCATTATGCCTTTAGGCTCTACACACTGCGACAACACAGCTATGCGCCTGCCGCCTCTACTccgccgaccac                >  1:4352146/1‑95 (MQ=255)
                                                                                    gagaTAAAAATGCCAG‑C‑CGATAGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCAGCTTATACTTCGACGGGca                 >  1:3475383/1‑99 (MQ=255)
                                                                                     aaggcaaaaTGCCAGCC‑CGATC‑GGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGcaca               >  2:2362872/6‑100 (MQ=37)
                                                                                      gaTAAAAATGCCAG‑C‑CGATCGGGCTGGCATTTTGTCTATAAAATGTACACAATGAGACAGAAGACATATGCGACTGAAGCTTCTACttaaaacggc                  >  1:1284771/1‑88 (MQ=255)
                                                                                      gaTAAAAATGCCAG‑C‑CGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACaa              >  1:2860238/1‑100 (MQ=255)
                                                                                        tAAAAATGCCAG‑C‑CGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATa            >  1:777757/1‑100 (MQ=255)
                                                                                        tAAAAATGCCAG‑C‑CGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATa            >  1:1752796/1‑100 (MQ=255)
                                                                                        tAAAAATGCCAG‑C‑CGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATa            >  1:1793929/1‑100 (MQ=255)
                                                                                         caaaaTGCCAGCC‑CGATC‑GGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGGAGAGCTATGCGACTGCCGATTCTACTTCGACGGGCACAATaa           >  1:2043201/2‑100 (MQ=38)
                                                                                         caaaaTGCCAGCC‑CGATC‑GGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATaa           >  2:192932/2‑100 (MQ=38)
                                                                                         caaaaTGCCAGCC‑CGATC‑GGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATaa           >  2:484729/2‑100 (MQ=38)
                                                                                         aaaaaTGCCAG‑C‑CGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATaa           <  2:3662234/100‑1 (MQ=255)
                                                                                           aaaTGCCAGCC‑CGATC‑GGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAaca         >  2:22004/1‑100 (MQ=39)
                                                                                           aaaTGCCAG‑C‑CGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATa            <  1:3046866/97‑1 (MQ=255)
                                                                                            aaTGCCAGCC‑CGATC‑GGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAacac        >  2:1165064/1‑100 (MQ=39)
                                                                                            aaTGCCAGCC‑CGATC‑GGCTGGCATTTTGCCTTTAGGATATACACAATGAGACAGAAGAGCTATGCGACTACCGCTTCTACTTCGACGGGCACAAt             >  1:591045/1‑95 (MQ=38)
                                                                                             aTGCCAG‑C‑CGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATa            <  1:3058062/95‑1 (MQ=255)
                                                                                              tGCCAGCC‑CGATC‑GGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAacac        >  1:2085198/1‑98 (MQ=39)
                                                                                               gCCAGCC‑CGATC‑GGCTGGCATTTTGCCTTTAGGATGTACACCATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTgg     >  1:868155/1‑100 (MQ=255)
                                                                                               gCCAGCC‑CGATC‑GGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTgg     >  2:2700580/1‑100 (MQ=255)
                                                                                                 cAG‑C‑CGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCg   <  2:4399878/100‑1 (MQ=255)
                                                                                                 cAG‑C‑CGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCg   <  2:1786152/100‑1 (MQ=255)
                                                                                                  aG‑C‑CGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGt  >  2:1043054/1‑100 (MQ=255)
                                                                                                    |                                                                                                   
CACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑C‑CGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGT  >  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: ‑