Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 922,314:1 +C 11.7% intergenic (+21/+52) 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.C11.7% 425.4 / 55.8 120intergenic (+21/+52)macB/cspDmacrolide ABC transporter peremase/ATPase/inhibitor of DNA replication, cold shock protein homolog
Reads supporting (aligned to +/- strand):  ref base . (53/53);  new base C (8/6);  total (61/59)
Fisher's exact test for biased strand distribution p-value = 7.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.04e-01

Only 100 of 123 total aligned reads displayed.

ATTGGTTTTTCACCGTTGGCGCTGCTGCTGGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCT  >  NC_000913/922165‑922462
                                                                                                                                                       |                                                                                                                                                    
aTTGGTCTTTCACCGTTGGCGCTGCTGCTGGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAg                                                                                                                                                       <  2:996590/150‑1 (MQ=255)
  tGGGTTTTCACCGTTGGCGCTGCTGCTGGCGGTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAAGGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGcc                                                                                                                                                    <  1:1101646/150‑1 (MQ=255)
    gTTTTTCACCGTTGGCGCTGCTGCTGGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGa                                                                                                                                                  >  1:703567/1‑150 (MQ=255)
        ttCACCGTTGGCGCTGCTGCTGGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCgg                                                                                                                                              >  1:566886/1‑150 (MQ=255)
          cACCGTTGGCGCTGCTGCTGGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCCGATCggc                                                                                                                                             <  2:383388/150‑2 (MQ=255)
          cACCGTTGGCGCTGCTGCGGGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGc                                                                                                                                            <  1:767810/150‑1 (MQ=255)
           aCCGTTGGCGCTGCTGCTGGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCt                                                                                                                                           >  2:391540/1‑150 (MQ=255)
           aCCGTTGGCGCTGCTGCTGGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCt                                                                                                                                           >  2:1182989/1‑150 (MQ=255)
            ccGTTGGCGCTGCTGCTGGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGGTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTg                                                                                                                                          <  2:239893/150‑1 (MQ=255)
             cGTTGGCGCTGCTGCTGGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTgg                                                                                                                                         <  1:169410/150‑1 (MQ=255)
              gTTGGCGCTGCTGCTGGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCCGATCGGCTGGc                                                                                                                                        <  1:741614/150‑1 (MQ=255)
                tGGCGCTGCTGCTGGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGC‑‑‑CCGATCGGCTGGCAttttt                                                                                                                                  <  2:536878/150‑2 (MQ=255)
                     gtgctgctgGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGTATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCCGATCGGCTGGCAttttta                                                                                                                                 <  2:1288063/149‑3 (MQ=255)
                     ctgctgctgGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCCGATCGGCTGGCAttttta                                                                                                                                 <  1:962925/150‑3 (MQ=255)
                     ctgctgctgGCGTTTCTCTGCTCGACGGTCACCAGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTg                                                                                                                                  >  2:662719/1‑149 (MQ=255)
                      tgctgctgGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGcc                                                                                                                                <  1:786311/150‑1 (MQ=255)
                        ctgctgGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCCGATCGGCTGGCAtttttatct                                                                                                                              <  1:694504/150‑6 (MQ=255)
                          gctgGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTa                                                                                                                            >  2:904068/1‑150 (MQ=255)
                           ctgGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAg                                                                                                                           >  2:242198/1‑150 (MQ=255)
                           ctgGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAg                                                                                                                           >  1:1487486/1‑150 (MQ=255)
                           ctgGCGTTTCTCTGCTCGACGGTCACCGGGATTTTAATTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAg                                                                                                                           >  2:695072/1‑150 (MQ=255)
                            tgGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAgg                                                                                                                          <  1:1231547/150‑1 (MQ=255)
                            tgGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAgg                                                                                                                          <  1:774643/150‑1 (MQ=255)
                            gggCGTTTCCCTGCTCGCCGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAgg                                                                                                                          <  1:813401/149‑1 (MQ=255)
                               cGTTTGTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATg                                                                                                                       <  2:752696/150‑1 (MQ=255)
                               cGTTTCTCTGCTCGACGGTCACCGGTATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATg                                                                                                                       <  1:262320/150‑1 (MQ=255)
                                gTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCATTTAGGATGt                                                                                                                      >  2:1030951/1‑150 (MQ=255)
                                gTTTCTCTGCTCGACGGGCACCGGGATTTTATTTGGCTGCTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGt                                                                                                                      <  1:291791/150‑1 (MQ=255)
                                   tctcTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTaca                                                                                                                   >  2:184327/1‑150 (MQ=255)
                                    ctctGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTacac                                                                                                                  >  2:757608/1‑150 (MQ=255)
                                        gCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATg                                                                                                              >  2:1589683/1‑150 (MQ=255)
                                        gCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATg                                                                                                              >  1:699993/1‑150 (MQ=255)
                                          tCGACGGTCACCGGGTTTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATgag                                                                                                            <  1:305942/150‑1 (MQ=255)
                                          tCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGGCTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATgag                                                                                                            >  1:628500/1‑150 (MQ=255)
                                           cGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATgaga                                                                                                           >  2:1364046/1‑150 (MQ=255)
                                             aCGGTCACCGTGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACa                                                                                                         >  2:973104/1‑150 (MQ=255)
                                             aCGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACa                                                                                                         >  2:953033/1‑150 (MQ=255)
                                              cGGTCTCGGGGATTTTATTTGGCTGGTTACCCGCACGAATTGCGGAACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACag                                                                                                        <  1:1504019/150‑1 (MQ=255)
                                               ggTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACaga                                                                                                       <  1:189218/150‑1 (MQ=255)
                                                   aCCGGGATTTTATTTGGCTGGTTACCCGCACGAACTGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAagag                                                                                                   <  1:153539/150‑1 (MQ=255)
                                                    ccGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGc                                                                                                  >  2:61783/1‑150 (MQ=255)
                                                         aTTTTATTTGGCTGGTTACCCGCACGAAATGCGACACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGc                                                                                             >  2:457471/1‑150 (MQ=255)
                                                           tttATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGa                                                                                           >  1:544965/1‑150 (MQ=255)
                                                            ttATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCAtttt                                                                                                                                   <  2:1043772/109‑1 (MQ=255)
                                                            ttATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCAtttt                                                                                                                                   >  1:1043772/1‑109 (MQ=255)
                                                             tATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACt                                                                                         >  1:1723050/1‑150 (MQ=255)
                                                                ttGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGcc                                                                                      <  1:242198/150‑1 (MQ=255)
                                                                  ggCTGGTTCCCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGc                                                                                    >  2:537824/1‑150 (MQ=255)
                                                                  ggCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGc                                                                                    >  1:183401/1‑150 (MQ=255)
                                                                   gCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCt                                                                                   >  1:12023/1‑150 (MQ=255)
                                                                    cTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCtt                                                                                  <  1:494157/150‑1 (MQ=255)
                                                                    cTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCtt                                                                                  >  2:265030/1‑150 (MQ=255)
                                                                     gggTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTc                                                                                 <  2:444155/149‑1 (MQ=255)
                                                                      ggTTACCCGCACGAAATGTGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATCAGACAGAAGAGCTATGCGACTGCCGCTTCt                                                                                >  2:354871/1‑150 (MQ=255)
                                                                      ggTTACCCGCACGAAATGCGGCACGACTGGAGCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCt                                                                                <  2:1053118/150‑1 (MQ=255)
                                                                       ttttcccGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTTATTTTTGAGATAATAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGTTGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTa                                                                               <  1:125131/146‑1 (MQ=255)
                                                                       gTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTa                                                                               >  1:544610/1‑150 (MQ=255)
                                                                        ttACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTAc                                                                              <  1:1589683/150‑1 (MQ=255)
                                                                          aCCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATgaga                                                                                                           <  1:1646050/119‑1 (MQ=255)
                                                                          aCCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATgaga                                                                                                           >  2:1646050/1‑119 (MQ=255)
                                                                           cccGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTc                                                                           >  2:1151638/1‑150 (MQ=255)
                                                                            ccGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGAGGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTc                                                                           <  1:32967/150‑1 (MQ=255)
                                                                             cGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGa                                                                         >  2:1711976/1‑150 (MQ=255)
                                                                                    aaTGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGcac                                                                  <  1:1173862/150‑1 (MQ=255)
                                                                                    aaTGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGcac                                                                  <  2:1463641/150‑1 (MQ=255)
                                                                                     aTGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGcaca                                                                 <  1:630789/150‑1 (MQ=255)
                                                                                      ttcGGCACGACGGGCTCCAGTAGTTGCTCGGGTACTTGATTAATTTTTGATATAAAAATGCCAGCCGATCTGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTAAGCGACTGCCGCATCTACTTCGACGGGCACaa                                                                <  2:158256/148‑1 (MQ=255)
                                                                                        cGGCACGCCTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCTCAATa                                                              >  1:90594/1‑150 (MQ=255)
                                                                                                cTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGc                                                      >  1:444388/1‑150 (MQ=255)
                                                                                                cTGGATCCAGTAGATGCTCTGGCAAGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGCAGAGCTATGCGACTGCCGCTTCTCCTTCGACGGGCACAATAACACTGGc                                                      <  2:1556100/150‑1 (MQ=255)
                                                                                                   gATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTg                                                   >  2:1287571/1‑150 (MQ=255)
                                                                                                   gATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTg                                                   <  2:198937/150‑1 (MQ=255)
                                                                                                   gATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTg                                                   >  2:604358/1‑150 (MQ=255)
                                                                                                     tCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGAt                                                 <  2:628500/150‑1 (MQ=255)
                                                                                                       cAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTg                                               >  1:570819/1‑150 (MQ=255)
                                                                                                              tGCTCTGTGACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTg                                        <  2:168723/150‑1 (MQ=255)
                                                                                                              tGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTg                                        <  1:1341843/150‑1 (MQ=255)
                                                                                                                 tctGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGcc                                     <  1:614140/150‑1 (MQ=255)
                                                                                                                  ctGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGccc                                    >  1:583836/1‑150 (MQ=255)
                                                                                                                         gagaGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCtggtgg                             <  2:1423927/150‑1 (MQ=255)
                                                                                                                            agTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTAAACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTATTGGCCCCTGGTGGACa                          >  2:1350843/1‑150 (MQ=255)
                                                                                                                             gTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACAt                         >  2:546399/1‑150 (MQ=255)
                                                                                                                                aTTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACAATGGCGTGATTGCCTTTTGGCCACTGGTGGACATCac                      >  1:1292588/1‑149 (MQ=255)
                                                                                                                                   ttttGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCGTGGTGGACATCAAAc                    >  1:1436298/1‑150 (MQ=255)
                                                                                                                                    tttGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACt                   <  1:646893/150‑1 (MQ=255)
                                                                                                                                     ttGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTg                  <  1:1255613/150‑1 (MQ=255)
                                                                                                                                      tGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGa                 <  1:457471/150‑1 (MQ=255)
                                                                                                                                        agaTAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAAc               >  1:1614779/1‑150 (MQ=255)
                                                                                                                                        aaggcaaaaTGCCAGCCCGAACGGCTGGAATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAAc               >  2:1262030/6‑150 (MQ=255)
                                                                                                                                         gaTAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACg              <  1:1711976/150‑1 (MQ=255)
                                                                                                                                         aggcaaaaTGCCAGCCCGATCGGCTGGCATTTTGCCTTTAGGTTTTACACAATGAGACAGAAGAGCTATGCGACTGCCGGTTCTACTTCGACGGTCACAATAACACTGGCGTGCTTGCCGTTTGGCCCCTGTGGGACATCAAACTGACCg              >  1:85710/5‑150 (MQ=255)
                                                                                                                                           gcaaaaTGCCAGCCCGATCGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGa            >  2:1594059/3‑150 (MQ=255)
                                                                                                                                            caaaaTGCCAGCCCGATCGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGAt           >  1:1660328/2‑150 (MQ=255)
                                                                                                                                            aaaaaTGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGAt           <  1:1143156/150‑1 (MQ=255)
                                                                                                                                               aaTGCCAGCCCGATCGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGATATCAAACTGAACGGATTGt        >  2:1366332/1‑150 (MQ=255)
                                                                                                                                                aTGCCAGCCCGATCGGCTGGCATTTTGCCTTTAGGAGGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTc       >  1:937325/1‑150 (MQ=255)
                                                                                                                                                 tgcacgCCGATTGGGCTTGCATTTTGCCTATAGGATTTACTAAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTcc      <  2:320490/145‑1 (MQ=255)
                                                                                                                                                  gCCAGCCCGATCGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCa     >  2:1672045/1‑150 (MQ=255)
                                                                                                                                                    cAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGCCATCAAACTGAACGGATTGTCCAGc   >  2:635333/1‑150 (MQ=255)
                                                                                                                                                     aGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCt  <  1:953033/150‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                    
ATTGGTTTTTCACCGTTGGCGCTGCTGCTGGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCT  >  NC_000913/922165‑922462

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: