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.

ATTGGTTTTTCACCGTTGGCGCTGCTGCTGGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAA‑TTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCT  >  NC_000913/922165‑922462
                                                                                                                                                       |                                                                                                                                                    
aTTGGTCTTTCACCGTTGGCGCTGCTGCTGGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAA‑TTTTTGAGATAAAAATGCCAg                                                                                                                                                       <  2:996590/150‑1 (MQ=255)
  tGGGTTTTCACCGTTGGCGCTGCTGCTGGCGGTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAAGGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAA‑TTTTTGAGATAAAAATGCCAG‑cc                                                                                                                                                    <  1:1101646/150‑1 (MQ=255)
    gTTTTTCACCGTTGGCGCTGCTGCTGGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAA‑TTTTTGAGATAAAAATGCCAG‑CCGa                                                                                                                                                  >  1:703567/1‑150 (MQ=255)
        ttCACCGTTGGCGCTGCTGCTGGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAA‑TTTTTGAGATAAAAATGCCAG‑CCGATCgg                                                                                                                                              >  1:566886/1‑150 (MQ=255)
          cACCGTTGGCGCTGCTGCTGGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAA‑TTTTTGAGATAAAAATGCCAGCCCGATCggc                                                                                                                                             <  2:383388/150‑2 (MQ=255)
          cACCGTTGGCGCTGCTGCGGGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAA‑TTTTTGAGATAAAAATGCCAG‑CCGATCGGGc                                                                                                                                            <  1:767810/150‑1 (MQ=255)
           aCCGTTGGCGCTGCTGCTGGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAA‑TTTTTGAGATAAAAATGCCAG‑CCGATCGGGCt                                                                                                                                           >  2:391540/1‑150 (MQ=255)
           aCCGTTGGCGCTGCTGCTGGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAA‑TTTTTGAGATAAAAATGCCAG‑CCGATCGGGCt                                                                                                                                           >  2:1182989/1‑150 (MQ=255)
            ccGTTGGCGCTGCTGCTGGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGGTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAA‑TTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTg                                                                                                                                          <  2:239893/150‑1 (MQ=255)
             cGTTGGCGCTGCTGCTGGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAA‑TTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTgg                                                                                                                                         <  1:169410/150‑1 (MQ=255)
              gTTGGCGCTGCTGCTGGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAA‑TTTTTGAGATAAAAATGCCAGCCCGATC‑GGCTGGc                                                                                                                                        <  1:741614/150‑1 (MQ=255)
                tGGCGCTGCTGCTGGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAA‑TTTTTGAGATAAAAATGC‑‑‑‑CCGATC‑GGCTGGCAttttt                                                                                                                                  <  2:536878/150‑2 (MQ=255)
                     gtgctgctgGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGTATCCAGTAGATGCTCTGGCACGAGAGTAA‑TTTTTGAGATAAAAATGCCAGCCCGATC‑GGCTGGCAttttta                                                                                                                                 <  2:1288063/149‑3 (MQ=255)
                     ctgctgctgGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAA‑TTTTTGAGATAAAAATGCCAGCCCGATC‑GGCTGGCAttttta                                                                                                                                 <  1:962925/150‑3 (MQ=255)
                     ctgctgctgGCGTTTCTCTGCTCGACGGTCACCAGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAA‑TTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTg                                                                                                                                  >  2:662719/1‑149 (MQ=255)
                      tgctgctgGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAA‑TTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGcc                                                                                                                                <  1:786311/150‑1 (MQ=255)
                        ctgctgGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAA‑TTTTTGAGATAAAAATGCCAGCCCGATC‑GGCTGGCAtttttatct                                                                                                                              <  1:694504/150‑6 (MQ=255)
                          gctgGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAA‑TTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTa                                                                                                                            >  2:904068/1‑150 (MQ=255)
                           ctgGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAA‑TTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAg                                                                                                                           >  2:242198/1‑150 (MQ=255)
                           ctgGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAA‑TTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAg                                                                                                                           >  1:1487486/1‑150 (MQ=255)
                           ctgGCGTTTCTCTGCTCGACGGTCACCGGGATTTTAATTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAA‑TTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAg                                                                                                                           >  2:695072/1‑150 (MQ=255)
                            tgGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAA‑TTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAgg                                                                                                                          <  1:1231547/150‑1 (MQ=255)
                            tgGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAA‑TTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAgg                                                                                                                          <  1:774643/150‑1 (MQ=255)
                            gggCGTTTCCCTGCTCGCCGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAA‑TTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAgg                                                                                                                          <  1:813401/149‑1 (MQ=255)
                               cGTTTGTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAA‑TTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATg                                                                                                                       <  2:752696/150‑1 (MQ=255)
                               cGTTTCTCTGCTCGACGGTCACCGGTATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAA‑TTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATg                                                                                                                       <  1:262320/150‑1 (MQ=255)
                                gTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAA‑TTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCATTTAGGATGt                                                                                                                      >  2:1030951/1‑150 (MQ=255)
                                gTTTCTCTGCTCGACGGGCACCGGGATTTTATTTGGCTGCTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAA‑TTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGt                                                                                                                      <  1:291791/150‑1 (MQ=255)
                                   tctcTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAA‑TTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTaca                                                                                                                   >  2:184327/1‑150 (MQ=255)
                                    ctctGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAA‑TTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTacac                                                                                                                  >  2:757608/1‑150 (MQ=255)
                                        gCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAA‑TTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATg                                                                                                              >  2:1589683/1‑150 (MQ=255)
                                        gCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAA‑TTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATg                                                                                                              >  1:699993/1‑150 (MQ=255)
                                          tCGACGGTCACCGGGTTTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAA‑TTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATgag                                                                                                            <  1:305942/150‑1 (MQ=255)
                                          tCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGGCTGGATCCAGTAGATGCTCTGGCACGAGAGTAA‑TTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATgag                                                                                                            >  1:628500/1‑150 (MQ=255)
                                           cGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAA‑TTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATgaga                                                                                                           >  2:1364046/1‑150 (MQ=255)
                                             aCGGTCACCGTGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAA‑TTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACa                                                                                                         >  2:973104/1‑150 (MQ=255)
                                             aCGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAA‑TTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACa                                                                                                         >  2:953033/1‑150 (MQ=255)
                                              cGGTCTCGGGGATTTTATTTGGCTGGTTACCCGCACGAATTGCGGAACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAA‑TTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACag                                                                                                        <  1:1504019/150‑1 (MQ=255)
                                               ggTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAA‑TTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACaga                                                                                                       <  1:189218/150‑1 (MQ=255)
                                                   aCCGGGATTTTATTTGGCTGGTTACCCGCACGAACTGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAA‑TTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAagag                                                                                                   <  1:153539/150‑1 (MQ=255)
                                                    ccGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAA‑TTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGc                                                                                                  >  2:61783/1‑150 (MQ=255)
                                                         aTTTTATTTGGCTGGTTACCCGCACGAAATGCGACACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAA‑TTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGc                                                                                             >  2:457471/1‑150 (MQ=255)
                                                           tttATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAA‑TTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGa                                                                                           >  1:544965/1‑150 (MQ=255)
                                                            ttATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAA‑TTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCAtttt                                                                                                                                   <  2:1043772/109‑1 (MQ=255)
                                                            ttATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAA‑TTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCAtttt                                                                                                                                   >  1:1043772/1‑109 (MQ=255)
                                                             tATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAA‑TTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACt                                                                                         >  1:1723050/1‑150 (MQ=255)
                                                                ttGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAA‑TTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGcc                                                                                      <  1:242198/150‑1 (MQ=255)
                                                                  ggCTGGTTCCCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAA‑TTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGc                                                                                    >  2:537824/1‑150 (MQ=255)
                                                                  ggCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAA‑TTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGc                                                                                    >  1:183401/1‑150 (MQ=255)
                                                                   gCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAA‑TTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCt                                                                                   >  1:12023/1‑150 (MQ=255)
                                                                    cTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAA‑TTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCtt                                                                                  <  1:494157/150‑1 (MQ=255)
                                                                    cTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAA‑TTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCtt                                                                                  >  2:265030/1‑150 (MQ=255)
                                                                     gggTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAA‑TTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTc                                                                                 <  2:444155/149‑1 (MQ=255)
                                                                      ggTTACCCGCACGAAATGTGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAA‑TTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATCAGACAGAAGAGCTATGCGACTGCCGCTTCt                                                                                >  2:354871/1‑150 (MQ=255)
                                                                      ggTTACCCGCACGAAATGCGGCACGACTGGAGCCAGTAGATGCTCTGGCACGAGAGTAA‑TTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCt                                                                                <  2:1053118/150‑1 (MQ=255)
                                                                       ttttcccGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTTA‑TTTTTGAGATAATAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGTTGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTa                                                                               <  1:125131/146‑1 (MQ=255)
                                                                       gTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAA‑TTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTa                                                                               >  1:544610/1‑150 (MQ=255)
                                                                        ttACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAA‑TTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTAc                                                                              <  1:1589683/150‑1 (MQ=255)
                                                                          aCCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAA‑TTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATgaga                                                                                                           <  1:1646050/119‑1 (MQ=255)
                                                                          aCCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAA‑TTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATgaga                                                                                                           >  2:1646050/1‑119 (MQ=255)
                                                                           cccGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAA‑TTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTc                                                                           >  2:1151638/1‑150 (MQ=255)
                                                                            ccGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGAGGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTc                                                                           <  1:32967/150‑1 (MQ=255)
                                                                             cGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAA‑TTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGa                                                                         >  2:1711976/1‑150 (MQ=255)
                                                                                    aaTGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAA‑TTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGcac                                                                  <  1:1173862/150‑1 (MQ=255)
                                                                                    aaTGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAA‑TTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGcac                                                                  <  2:1463641/150‑1 (MQ=255)
                                                                                     aTGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAA‑TTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGcaca                                                                 <  1:630789/150‑1 (MQ=255)
                                                                                      ttcGGCACGACGGGCTCCAGTAGTTGCTCGGGTACTTGATTAA‑TTTTTGATATAAAAATGCCAG‑CCGATCTGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTAAGCGACTGCCGCATCTACTTCGACGGGCACaa                                                                <  2:158256/148‑1 (MQ=255)
                                                                                        cGGCACGCCTGGATCCAGTAGATGCTCTGGCACGAGAGTAA‑TTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCTCAATa                                                              >  1:90594/1‑150 (MQ=255)
                                                                                                cTGGATCCAGTAGATGCTCTGGCACGAGAGTAA‑TTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGc                                                      >  1:444388/1‑150 (MQ=255)
                                                                                                cTGGATCCAGTAGATGCTCTGGCAAGAGAGTAA‑TTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGCAGAGCTATGCGACTGCCGCTTCTCCTTCGACGGGCACAATAACACTGGc                                                      <  2:1556100/150‑1 (MQ=255)
                                                                                                   gATCCAGTAGATGCTCTGGCACGAGAGTAA‑TTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTg                                                   >  2:1287571/1‑150 (MQ=255)
                                                                                                   gATCCAGTAGATGCTCTGGCACGAGAGTAA‑TTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTg                                                   <  2:198937/150‑1 (MQ=255)
                                                                                                   gATCCAGTAGATGCTCTGGCACGAGAGTAA‑TTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTg                                                   >  2:604358/1‑150 (MQ=255)
                                                                                                     tCCAGTAGATGCTCTGGCACGAGAGTAA‑TTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGAt                                                 <  2:628500/150‑1 (MQ=255)
                                                                                                       cAGTAGATGCTCTGGCACGAGAGTAA‑TTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTg                                               >  1:570819/1‑150 (MQ=255)
                                                                                                              tGCTCTGTGACGAGAGTAA‑TTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTg                                        <  2:168723/150‑1 (MQ=255)
                                                                                                              tGCTCTGGCACGAGAGTAA‑TTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTg                                        <  1:1341843/150‑1 (MQ=255)
                                                                                                                 tctGGCACGAGAGTAA‑TTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGcc                                     <  1:614140/150‑1 (MQ=255)
                                                                                                                  ctGGCACGAGAGTAA‑TTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGccc                                    >  1:583836/1‑150 (MQ=255)
                                                                                                                         gagaGTAA‑TTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCtggtgg                             <  2:1423927/150‑1 (MQ=255)
                                                                                                                            agTAA‑TTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTAAACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTATTGGCCCCTGGTGGACa                          >  2:1350843/1‑150 (MQ=255)
                                                                                                                             gTAA‑TTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACAt                         >  2:546399/1‑150 (MQ=255)
                                                                                                                                a‑TTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACAATGGCGTGATTGCCTTTTGGCCACTGGTGGACATCac                      >  1:1292588/1‑149 (MQ=255)
                                                                                                                                   ttttGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCGTGGTGGACATCAAAc                    >  1:1436298/1‑150 (MQ=255)
                                                                                                                                    tttGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACt                   <  1:646893/150‑1 (MQ=255)
                                                                                                                                     ttGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTg                  <  1:1255613/150‑1 (MQ=255)
                                                                                                                                      tGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGa                 <  1:457471/150‑1 (MQ=255)
                                                                                                                                        agaTAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAAc               >  1:1614779/1‑150 (MQ=255)
                                                                                                                                        aaggcaaaaTGCCAGCCCGAAC‑GGCTGGAATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAAc               >  2:1262030/6‑150 (MQ=255)
                                                                                                                                         gaTAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACg              <  1:1711976/150‑1 (MQ=255)
                                                                                                                                         aggcaaaaTGCCAGCCCGATC‑GGCTGGCATTTTGCCTTTAGGTTTTACACAATGAGACAGAAGAGCTATGCGACTGCCGGTTCTACTTCGACGGTCACAATAACACTGGCGTGCTTGCCGTTTGGCCCCTGTGGGACATCAAACTGACCg              >  1:85710/5‑150 (MQ=255)
                                                                                                                                           gcaaaaTGCCAGCCCGATC‑GGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGa            >  2:1594059/3‑150 (MQ=255)
                                                                                                                                            caaaaTGCCAGCCCGATC‑GGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGAt           >  1:1660328/2‑150 (MQ=255)
                                                                                                                                            aaaaaTGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGAt           <  1:1143156/150‑1 (MQ=255)
                                                                                                                                               aaTGCCAGCCCGATC‑GGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGATATCAAACTGAACGGATTGt        >  2:1366332/1‑150 (MQ=255)
                                                                                                                                                aTGCCAGCCCGATC‑GGCTGGCATTTTGCCTTTAGGAGGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTc       >  1:937325/1‑150 (MQ=255)
                                                                                                                                                 tgcacg‑CCGATTGGGCTTGCATTTTGCCTATAGGATTTACTAAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTcc      <  2:320490/145‑1 (MQ=255)
                                                                                                                                                  gCCAGCCCGATC‑GGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCa     >  2:1672045/1‑150 (MQ=255)
                                                                                                                                                    cAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGCCATCAAACTGAACGGATTGTCCAGc   >  2:635333/1‑150 (MQ=255)
                                                                                                                                                     aG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCt  <  1:953033/150‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                    
ATTGGTTTTTCACCGTTGGCGCTGCTGCTGGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAA‑TTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCT  >  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: ‑