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.

CACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGT  >  NC_000913/922215‑922412
                                                                                                    |                                                                                                   
cACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAg                                                                                                      >  1:3892672/1‑100 (MQ=255)
cACCGGGATTTTATGTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAg                                                                                                      <  1:106911/100‑1 (MQ=255)
 aCCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGc                                                                                                    <  2:3801801/100‑1 (MQ=255)
  ccGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGcc                                                                                                  >  2:2288841/1‑100 (MQ=255)
    gggATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGa                                                                                                >  2:744164/1‑100 (MQ=255)
      gATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGa                                                                                                >  1:3448416/1‑98 (MQ=255)
       aTTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATc                                                                                              <  1:696963/99‑1 (MQ=255)
       aTTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCg                                                                                             <  2:2209095/100‑1 (MQ=255)
          ttATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGc                                                                                          <  2:3495003/100‑1 (MQ=255)
             tttGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTgg                                                                                       >  1:108707/1‑100 (MQ=255)
             tttGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTgg                                                                                       <  1:1655931/100‑1 (MQ=255)
             tttGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTgg                                                                                       >  2:1568083/1‑100 (MQ=255)
             tttGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTgg                                                                                       >  1:3838279/1‑100 (MQ=255)
              tttgCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAATATGCCAGCCGATCGGGCTgg                                                                                       <  1:2695965/96‑1 (MQ=255)
              ttGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCCGATCGGCTGGc                                                                                      <  2:3884416/100‑1 (MQ=255)
                ggCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCCGATCGGCTGGCAt                                                                                    <  2:3324476/100‑1 (MQ=39)
                 gCTTGTTACCCTCACGTAATTCGGCACGACTGGATCTAGTAGATTCGCTTGCAAGAGATTAATTTTTTAGATAAAAATGTCAGcctctcg                                                                                             <  1:3487893/90‑6 (MQ=255)
                 gCTGGTGACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCCGATCGGATGGCAtt                                                                                   <  1:103956/100‑1 (MQ=255)
                   tGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCAtttt                                                                                 <  1:1202185/100‑1 (MQ=255)
                    ggTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCCGATCGGCTGGCAttttt                                                                                <  2:478505/100‑2 (MQ=38)
                    ggTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCCGATCGGCTGGCAttttt                                                                                <  2:2130991/100‑2 (MQ=38)
                    ggTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTg                                                                                <  1:4319757/100‑1 (MQ=255)
                     gTTACCCGCACGAAATGCGGCACGACTGTATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCCGATCGGCTGGCAttttta                                                                               <  2:1468867/100‑3 (MQ=38)
                      ttACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGcc                                                                              <  2:3833329/100‑1 (MQ=255)
                      ttACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGcc                                                                              <  2:402963/100‑1 (MQ=255)
                      ttACCCGCACGAAATGCGGCACGACTGGATCCAGGAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGcc                                                                              <  2:2055883/100‑1 (MQ=255)
                      tgaCCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGcc                                                                              <  1:941963/98‑1 (MQ=255)
                           cgtaCGAAATGCCGTACGACTGCATCCAGTAGATGCTTTGCCACGAGAGTCATTTTTGAGATAAAAATGCCTGCACGATCGGCTGGCATTTTgatat                                                                            <  1:3544333/94‑5 (MQ=255)
                           cGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGGGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAg                                                                         >  2:2545883/1‑100 (MQ=255)
                           cGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAg                                                                         >  2:1346828/1‑100 (MQ=255)
                           cGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCCGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAg                                                                         >  1:1522337/1‑100 (MQ=255)
                            gCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAgg                                                                        <  1:355641/100‑1 (MQ=255)
                             cACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCCGATCGGCTGGCAtttttatctcaaaa                                                                       <  2:2952737/100‑11 (MQ=25)
                             cACGAAATGCGGCACGACTGGATCCAGGAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGa                                                                       <  2:3157181/100‑1 (MQ=255)
                              aCGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGAt                                                                      >  1:3930418/1‑100 (MQ=255)
                               cGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAgg                                                                        <  1:1821061/97‑1 (MQ=255)
                               cGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGa                                                                       >  1:1958890/1‑98 (MQ=255)
                                    tGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTacac                                                                <  2:1251767/100‑1 (MQ=255)
                                     gCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTAcaca                                                               <  1:413353/100‑1 (MQ=255)
                                      cGGCACGACTGGATCCAGTAGATGCTCTGGCACGATAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTAcaca                                                               <  1:1986271/99‑1 (MQ=255)
                                        gCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAAt                                                             >  1:2752846/1‑99 (MQ=255)
                                         cACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTTTTTTAATTATGCAAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTAGACAATga                                                           >  1:2115546/1‑100 (MQ=255)
                                         cACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATga                                                           <  1:3397792/100‑1 (MQ=255)
                                         cACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATga                                                           <  1:2933392/100‑1 (MQ=255)
                                                ggATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTTTTTTTTTATTGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACag                                                      >  1:356705/1‑98 (MQ=255)
                                                ggATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACagaa                                                    >  1:44609/1‑100 (MQ=255)
                                                 gATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACagaag                                                   <  2:2092084/100‑1 (MQ=255)
                                                  aTCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACagaaga                                                  >  1:2425485/1‑100 (MQ=255)
                                                        tAGATGCTCTGGCACGGGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATg                                            <  1:154809/100‑1 (MQ=255)
                                                        tAGATGCGCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGGATTTTGCCTTTACGATGTACACAATGAGACAGAAGAGCTATg                                            <  1:3228252/100‑1 (MQ=255)
                                                         aGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGc                                           <  2:487743/100‑1 (MQ=255)
                                                            tGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGAc                                        >  2:4087395/1‑100 (MQ=255)
                                                             gCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACt                                       <  2:3632520/100‑1 (MQ=255)
                                                               tctGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACt                                       >  1:930019/1‑98 (MQ=255)
                                                               tctGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGc                                     >  1:4032089/1‑100 (MQ=255)
                                                                ctGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGcc                                    >  2:2126923/1‑100 (MQ=255)
                                                                 tGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCg                                   <  2:2344031/100‑1 (MQ=255)
                                                                   gCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCt                                 <  1:2795100/100‑1 (MQ=255)
                                                                   gCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCt                                 <  1:231982/100‑1 (MQ=255)
                                                                    cACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCtt                                >  1:474036/1‑100 (MQ=255)
                                                                    cACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCtt                                >  1:409511/1‑100 (MQ=255)
                                                                     aCGAGAGTAATTTTTGAGATAAAAATGCAAGCCGATCGGGGTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTc                               <  2:1746066/100‑1 (MQ=255)
                                                                      cGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCt                              <  1:122820/100‑1 (MQ=255)
                                                                       gagaGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTa                             >  2:604301/1‑100 (MQ=255)
                                                                        agagTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTAc                            <  2:181717/100‑1 (MQ=255)
                                                                          agTAATTTTTGAGATAAAAATGCCAGCCGATCGTGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACtt                          <  2:473038/100‑1 (MQ=255)
                                                                          agTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACtt                          >  2:2392882/1‑100 (MQ=255)
                                                                          agTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACtt                          <  2:2157630/100‑1 (MQ=255)
                                                                           gTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTCCTTc                         >  2:3772676/1‑100 (MQ=255)
                                                                           gTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTc                         >  2:3989235/1‑100 (MQ=255)
                                                                               tttttGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACg                     >  2:3563221/1‑100 (MQ=255)
                                                                               tttttGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACg                     >  2:882073/1‑100 (MQ=255)
                                                                                ttttGATATAAAAATGCCAGCCGATCGGGCTCGCATTTTTCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACg                     <  1:599012/99‑1 (MQ=255)
                                                                                ttttGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACg                     <  1:2545870/99‑1 (MQ=255)
                                                                                  ttGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGc                  <  2:3395902/100‑1 (MQ=255)
                                                                                   tGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGca                 <  1:4224981/100‑1 (MQ=255)
                                                                                    gagaTAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGCCAGAAGAGCTATGCGACTGCCGCTTCTAATTCGACgg                    >  1:206429/1‑96 (MQ=255)
                                                                                    gagaTAAAAATGCCAGCCGATCGGGCTGGCATTATGCCTTTAGGCTCTACACACTGCGACAACACAGCTATGCGCCTGCCGCCTCTACTccgccgaccac                >  1:4352146/1‑95 (MQ=255)
                                                                                    gagaTAAAAATGCCAGCCGATAGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCAGCTTATACTTCGACGGGca                 >  1:3475383/1‑99 (MQ=255)
                                                                                     aaggcaaaaTGCCAGCCCGATCGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGcaca               >  2:2362872/6‑100 (MQ=37)
                                                                                      gaTAAAAATGCCAGCCGATCGGGCTGGCATTTTGTCTATAAAATGTACACAATGAGACAGAAGACATATGCGACTGAAGCTTCTACttaaaacggc                  >  1:1284771/1‑88 (MQ=255)
                                                                                      gaTAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACaa              >  1:2860238/1‑100 (MQ=255)
                                                                                        tAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATa            >  1:777757/1‑100 (MQ=255)
                                                                                        tAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATa            >  1:1752796/1‑100 (MQ=255)
                                                                                        tAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATa            >  1:1793929/1‑100 (MQ=255)
                                                                                         caaaaTGCCAGCCCGATCGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGGAGAGCTATGCGACTGCCGATTCTACTTCGACGGGCACAATaa           >  1:2043201/2‑100 (MQ=38)
                                                                                         caaaaTGCCAGCCCGATCGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATaa           >  2:192932/2‑100 (MQ=38)
                                                                                         caaaaTGCCAGCCCGATCGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATaa           >  2:484729/2‑100 (MQ=38)
                                                                                         aaaaaTGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATaa           <  2:3662234/100‑1 (MQ=255)
                                                                                           aaaTGCCAGCCCGATCGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAaca         >  2:22004/1‑100 (MQ=39)
                                                                                           aaaTGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATa            <  1:3046866/97‑1 (MQ=255)
                                                                                            aaTGCCAGCCCGATCGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAacac        >  2:1165064/1‑100 (MQ=39)
                                                                                            aaTGCCAGCCCGATCGGCTGGCATTTTGCCTTTAGGATATACACAATGAGACAGAAGAGCTATGCGACTACCGCTTCTACTTCGACGGGCACAAt             >  1:591045/1‑95 (MQ=38)
                                                                                             aTGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATa            <  1:3058062/95‑1 (MQ=255)
                                                                                              tGCCAGCCCGATCGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAacac        >  1:2085198/1‑98 (MQ=39)
                                                                                               gCCAGCCCGATCGGCTGGCATTTTGCCTTTAGGATGTACACCATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTgg     >  1:868155/1‑100 (MQ=255)
                                                                                               gCCAGCCCGATCGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTgg     >  2:2700580/1‑100 (MQ=255)
                                                                                                 cAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCg   <  2:4399878/100‑1 (MQ=255)
                                                                                                 cAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCg   <  2:1786152/100‑1 (MQ=255)
                                                                                                  aGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGt  >  2:1043054/1‑100 (MQ=255)
                                                                                                    |                                                                                                   
CACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGT  >  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: