Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 922,314:1 +C 7.2% 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.C 7.2% 402.5 / 27.9 97intergenic (+21/+52)macB/cspDmacrolide ABC transporter peremase/ATPase/inhibitor of DNA replication, cold shock protein homolog
Reads supporting (aligned to +/- strand):  ref base . (42/48);  new base C (2/5);  total (44/53)
Fisher's exact test for biased strand distribution p-value = 4.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.96e-01

TTTTTCACCGTTGGCGCTGCTGCTGGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAG  >  NC_000913/922170‑922460
                                                                                                                                                 |                                                                                                                                                  
tttttCACCGTTGGCGCTGCTGCTGGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGcccg                                                                                                                                                >  4:8624/1‑147 (MQ=255)
   ttCACCGTTGGCGCTGCTGCTGGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCg                                                                                                                                             <  3:779888/149‑1 (MQ=255)
    tCACCGTTGGCGCTGCTGCTGGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCgg                                                                                                                                            <  1:281727/149‑1 (MQ=255)
      aCCGTTGGCGCTGCTGCTGGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGc                                                                                                                                          >  3:1285/1‑149 (MQ=255)
       ccGTTGGCGCTGCTGCTGGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCt                                                                                                                                         <  1:882832/149‑1 (MQ=255)
       ccGTTGGCGCTGCTGCTGGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCt                                                                                                                                         <  1:775113/149‑1 (MQ=255)
         gTTGGCGCTGCTGCTGGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTgg                                                                                                                                       <  4:54158/149‑1 (MQ=255)
            ggCGCTGCTGCTGGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCAt                                                                                                                                    <  3:453075/149‑1 (MQ=255)
             gcgcTGCTGCTGGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCAtt                                                                                                                                   >  4:725254/1‑149 (MQ=255)
             gcgcTGCTGCTGGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCCCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCCGATCGGCTGGCAtt                                                                                                                                   <  3:352096/149‑1 (MQ=255)
              cgcTGCTGCTGGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCCGATCGGCTGGCAttt                                                                                                                                  <  3:671252/149‑1 (MQ=255)
               gctgctgctGGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCAtttt                                                                                                                                 <  2:481961/149‑1 (MQ=255)
               gctgctgctGGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCAtttt                                                                                                                                 >  4:496205/1‑149 (MQ=255)
                ctgctgctgGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCCGATCGGCTGGCAttttt                                                                                                                                <  3:8624/149‑2 (MQ=255)
                 tgctgctgGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGc                                                                                                                               <  2:210095/149‑1 (MQ=255)
                    tgctgGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCCGATCGGCTGGCAtttttatct                                                                                                                            <  4:556016/149‑6 (MQ=255)
                     gctgGCGTTTCTCTGCTCGACGGTCACCGTGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCttt                                                                                                                           <  4:43020/149‑1 (MQ=255)
                     gctgGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCttt                                                                                                                           >  1:767342/1‑149 (MQ=255)
                      ctgGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTa                                                                                                                          <  3:252187/149‑1 (MQ=255)
                       tgGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAg                                                                                                                         <  4:559292/149‑1 (MQ=255)
                         gCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGa                                                                                                                       <  4:337285/149‑1 (MQ=255)
                         gCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGa                                                                                                                       >  2:687115/1‑149 (MQ=255)
                           gTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATg                                                                                                                     >  2:834255/1‑149 (MQ=255)
                                tctGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCCGATCGGCTGGCAtt                                                                                                                                   <  4:953758/130‑1 (MQ=255)
                                tctGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCCGATCGGCTGGCAtt                                                                                                                                   >  3:953758/1‑130 (MQ=255)
                                tctGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTacac                                                                                                                >  4:836746/1‑149 (MQ=255)
                                tctGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTacac                                                                                                                >  4:238007/1‑149 (MQ=255)
                                tctGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTacac                                                                                                                >  4:477294/1‑149 (MQ=255)
                                  tGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACaa                                                                                                              <  1:371107/149‑1 (MQ=255)
                                       gACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATgaga                                                                                                         <  1:697436/149‑1 (MQ=255)
                                       gACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATgaga                                                                                                         <  1:687115/149‑1 (MQ=255)
                                           gTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCAtttt                                                                                                                                 <  1:671789/121‑1 (MQ=255)
                                           gTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCAtttt                                                                                                                                 >  2:671789/1‑121 (MQ=255)
                                              aCCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACagaaga                                                                                                  >  4:716205/1‑149 (MQ=255)
                                                cGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGc                                                                                                <  1:1038456/149‑1 (MQ=255)
                                                cGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGc                                                                                                <  3:1015932/149‑1 (MQ=255)
                                                 gggATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCt                                                                                               <  4:276861/149‑1 (MQ=255)
                                                 gggATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCt                                                                                               >  4:910852/1‑149 (MQ=255)
                                                   gATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTAt                                                                                             >  3:859271/1‑149 (MQ=255)
                                                             ggCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCg                                                                                   >  4:306096/1‑149 (MQ=255)
                                                             ggCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCg                                                                                   >  4:169105/1‑149 (MQ=255)
                                                              gCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACACTGAGACAGAAGAGCTATGCGACTGCCGc                                                                                  <  3:73907/149‑1 (MQ=255)
                                                              gCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGc                                                                                  <  1:138974/149‑1 (MQ=255)
                                                               cTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCt                                                                                 >  3:715456/1‑149 (MQ=255)
                                                                tGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCAttt                                                                                                                                  <  3:118085/99‑1 (MQ=255)
                                                                tGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCAttt                                                                                                                                  >  4:118085/1‑99 (MQ=255)
                                                                 ggTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTc                                                                               <  2:122554/149‑1 (MQ=255)
                                                                 ggTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTc                                                                               <  4:632538/149‑1 (MQ=255)
                                                                       ccGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAg                                                                                                                         >  4:682224/1‑101 (MQ=255)
                                                                       ccGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAg                                                                                                                         <  3:682224/101‑1 (MQ=255)
                                                                       ccGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTc                                                                         <  2:39015/149‑1 (MQ=255)
                                                                       ccGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTc                                                                         >  3:387037/1‑149 (MQ=255)
                                                                          cACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGCTGTACCCAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGAc                                                                      <  3:348077/149‑1 (MQ=255)
                                                                            cGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACgg                                                                    <  1:998985/149‑1 (MQ=255)
                                                                             gAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACggg                                                                   >  3:66734/1‑149 (MQ=255)
                                                                                  gCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACaa                                                              <  2:767342/149‑1 (MQ=255)
                                                                                    ggCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATa                                                            >  2:815992/1‑149 (MQ=255)
                                                                                     gCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATaa                                                           >  3:165779/1‑149 (MQ=255)
                                                                                     gCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATaa                                                           >  4:417933/1‑149 (MQ=255)
                                                                                        cGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAATAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAacac                                                        <  3:890977/149‑1 (MQ=255)
                                                                                        cGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAacac                                                        <  1:1093297/149‑1 (MQ=255)
                                                                                         gACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACt                                                       <  4:859271/149‑1 (MQ=255)
                                                                                           cTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTgg                                                     >  4:376844/1‑149 (MQ=255)
                                                                                                tCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGCGa                                                >  4:147741/1‑149 (MQ=255)
                                                                                                   aGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTg                                             >  4:149320/1‑149 (MQ=255)
                                                                                                    gTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCtt                                                                                <  4:795898/113‑1 (MQ=255)
                                                                                                    gTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCtt                                                                                >  3:795898/1‑113 (MQ=255)
                                                                                                      aGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAAGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCt                                          >  1:592682/1‑149 (MQ=255)
                                                                                                        aTGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCCTTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCttt                                        <  1:815992/149‑1 (MQ=255)
                                                                                                         tGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCtttt                                       >  1:330916/1‑149 (MQ=255)
                                                                                                           ctctGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTgg                                     >  4:377854/1‑149 (MQ=255)
                                                                                                            tctGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGc                                    <  2:1076586/149‑1 (MQ=255)
                                                                                                            tctGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGc                                    >  4:72823/1‑149 (MQ=255)
                                                                                                            tctGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGc                                    >  1:762768/1‑149 (MQ=255)
                                                                                                                  aCGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCtgg                              >  1:178608/1‑149 (MQ=255)
                                                                                                                   cGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCtggt                             <  3:155675/149‑1 (MQ=255)
                                                                                                                      gagTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGa                          >  4:211470/1‑149 (MQ=255)
                                                                                                                       agTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGAc                         <  2:96398/149‑1 (MQ=255)
                                                                                                                        gTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGCGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACa                        >  4:654338/1‑149 (MQ=255)
                                                                                                                        gTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGTCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACa                        >  3:758080/1‑149 (MQ=255)
                                                                                                                          aaTTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATc                      <  3:417933/149‑1 (MQ=255)
                                                                                                                             ttttGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCaaa                   <  1:959018/149‑1 (MQ=255)
                                                                                                                               ttGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTCCCTTTTGGCCCCTGGTGGACATCAAACt                 <  3:72823/149‑1 (MQ=255)
                                                                                                                                    aTAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACg            >  2:25497/1‑149 (MQ=255)
                                                                                                                                      aaaaaTGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGa          <  3:725254/149‑1 (MQ=255)
                                                                                                                                       aaaaTGCCAGCCCGGTCGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGAt         >  3:83625/1‑149 (MQ=255)
                                                                                                                                       aaaaTGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGAt         <  3:238007/149‑1 (MQ=255)
                                                                                                                                        aaaTGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCt                                <  2:582174/125‑1 (MQ=255)
                                                                                                                                        aaaTGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCt                                >  1:582174/1‑125 (MQ=255)
                                                                                                                                        aaaTGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGAtt        >  2:631533/1‑149 (MQ=255)
                                                                                                                                             ccAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCa   <  3:169105/149‑1 (MQ=255)
                                                                                                                                             ccAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCa   <  3:306096/149‑1 (MQ=255)
                                                                                                                                              cAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGc                                                    >  1:242117/1‑99 (MQ=255)
                                                                                                                                              cAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGc                                                    <  2:242117/99‑1 (MQ=255)
                                                                                                                                              cAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAg  >  2:149729/1‑149 (MQ=255)
                                                                                                                                              cAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAg  <  4:1031904/149‑1 (MQ=255)
                                                                                                                                               aGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACtt                                                                          <  4:58886/76‑1 (MQ=255)
                                                                                                                                               aGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACtt                                                                          >  3:58886/1‑76 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
TTTTTCACCGTTGGCGCTGCTGCTGGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAG  >  NC_000913/922170‑922460

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: