Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,808,192 G→C 15.2% intergenic (‑301/‑254) PP_4234 ← / → dsbD‑II bifunctional aldehyde oxidase/xanthine dehydrogenase/thiol/disulfide interchange protein 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,808,1920GC15.2% 753.6 / 114.6 295intergenic (‑301/‑254)PP_4234/dsbD‑IIbifunctional aldehyde oxidase/xanthine dehydrogenase/thiol/disulfide interchange protein 2
Reads supporting (aligned to +/- strand):  ref base G (91/159);  new base C (28/17);  total (119/176)
Fisher's exact test for biased strand distribution p-value = 1.58e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.43e-03

Only 100 of 398 total aligned reads displayed.

GGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACA  >  NC_002947/4808044‑4808322
                                                                                                                                                    |                                                                                                                                  
ccagttattcctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCgcgg                                                                                                                                    <  4:96321‑M1/19‑1 (MQ=255)
 cagttattcctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCgcggg                                                                                                                                   >  5:337372‑M1/130‑149 (MQ=255)
 cagttattcctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCgcggg                                                                                                                                   <  5:258056‑M1/20‑1 (MQ=255)
 cagttattcctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcaggtttaCCCGCGAAGCGCCGCgcggg                                                                                                                                   <  7:50495‑M1/20‑1 (MQ=255)
  agttattcctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCgcgggc                                                                                                                                  <  1:139055‑M1/21‑1 (MQ=255)
  agttattcctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCgcgggc                                                                                                                                  >  6:274497‑M1/129‑149 (MQ=255)
   gttattcctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCgcgggcg                                                                                                                                 <  5:269259‑M1/22‑1 (MQ=255)
   gttattcctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCgcgggcg                                                                                                                                 <  6:18552‑M1/22‑1 (MQ=255)
   gttattcctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcacgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCgcgggcg                                                                                                                                 <  8:161423‑M1/22‑1 (MQ=255)
    ttattcctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCccgcgcgg                                                                                                                                <  5:86062‑M1/23‑1 (MQ=255)
      attcctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCgcggg                                                                                                                                   <  2:132520‑M1/20‑1 (MQ=255)
      attcctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCgcgg                                                                                                                                    <  4:272662‑M1/19‑1 (MQ=255)
      attcctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGgcggcg                                                                                                                              <  8:338908‑M1/25‑1 (MQ=255)
      attcctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGgcggcg                                                                                                                              <  7:189502‑M1/25‑1 (MQ=255)
      attcctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGgcggcg                                                                                                                              <  6:138595‑M1/25‑1 (MQ=255)
       ttcctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGgcgc                                                                                                                             <  3:202351‑M1/26‑1 (MQ=255)
       ttcctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGgcgc                                                                                                                             >  3:216323‑M1/124‑149 (MQ=255)
         cctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCt                                                                                                                            >  3:48965‑M1/122‑148 (MQ=255)
         cctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTc                                                                                                                           <  8:21272‑M1/28‑1 (MQ=255)
         cctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCCCGCGCGGCGCt                                                                                                                            >  6:86062‑M1/122‑148 (MQ=255)
         cctgggggacgcgctcgcgagccccctggcgcacgcctcgatatcgcggcgtacaatcaaggcgggcgcgcgcgcctggcacaggaggtattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGgcggcg                                                                                                                              >  6:232941‑M1/122‑146 (MQ=255)
         cctctgggacgcgcacgcgcgccccctggcgcacgcctcgatatcccgtcgtacaaccaaggcggtcgcgcgcgcctggcacacgagggatagcccagaaacaaacgtaggagccggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTc                                                                                                                           <  6:119554‑M1/28‑1 (MQ=255)
          ctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGggcggc                                                                                                                               <  4:334450‑M1/24‑1 (MQ=255)
          ctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGTGGCGCTc                                                                                                                           <  5:195685‑M1/28‑1 (MQ=255)
          ctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTc                                                                                                                           <  3:174892‑M1/28‑1 (MQ=255)
          ctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTc                                                                                                                           <  1:208057‑M1/28‑1 (MQ=255)
          ctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTc                                                                                                                           >  6:315110‑M1/121‑148 (MQ=255)
          ctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTc                                                                                                                           <  8:246674‑M1/28‑1 (MQ=255)
          ctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTc                                                                                                                           <  2:148082‑M1/28‑1 (MQ=255)
          ctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTc                                                                                                                           <  5:143653‑M1/28‑1 (MQ=255)
          ctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTc                                                                                                                           <  2:253454‑M1/28‑1 (MQ=255)
          ctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTc                                                                                                                           <  2:339795‑M1/28‑1 (MQ=255)
          ctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTc                                                                                                                           <  2:67873‑M1/28‑1 (MQ=255)
          ctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCg                                                                                                                          >  4:186638‑M1/121‑149 (MQ=255)
          ctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCg                                                                                                                          >  5:353388‑M1/121‑149 (MQ=255)
          ctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGCGCGGCGCTc                                                                                                                           <  6:360965‑M1/28‑1 (MQ=255)
          ctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGCGAGGCACTc                                                                                                                           <  6:251554‑M1/28‑1 (MQ=255)
          ctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCCCGCGCGGCgcttc                                                                                                                          >  3:288454‑M1/121‑147 (MQ=255)
          ctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggacgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTc                                                                                                                           <  2:98898‑M1/28‑1 (MQ=255)
           tgtgggacgcgcacgcgagccccctggcgcacgcctcgatctcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGa                                                                                                                         <  8:81644‑M1/30‑1 (MQ=255)
           tgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGa                                                                                                                         <  7:173468‑M1/30‑1 (MQ=255)
            gtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGAt                                                                                                                        <  7:117353‑M1/31‑1 (MQ=255)
             tgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGGAGCGCCGCCCGCGCggcgcttcgcg                                                                                                                       >  8:318631‑M1/118‑144 (MQ=255)
             tgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGGAGCGCCGCCCGCGCggcgcttcgcg                                                                                                                       >  5:214298‑M1/118‑144 (MQ=255)
             tgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGAtc                                                                                                                       <  4:307992‑M1/32‑1 (MQ=255)
             tgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGAtc                                                                                                                       >  4:133547‑M1/118‑149 (MQ=255)
             tgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGAtc                                                                                                                       <  7:268457‑M1/32‑1 (MQ=255)
               ggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagtgattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGAtc                                                                                                                       <  6:302706‑M1/32‑1 (MQ=255)
               ggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGGAGCGCCGCCCGCGcggcgcttcgcgg                                                                                                                      >  2:231004‑M1/116‑142 (MQ=255)
               ggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCgcgggcgg                                                                                                                                <  1:165359‑M1/23‑1 (MQ=255)
               ggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGAtctc                                                                                                                     >  4:28886‑M1/116‑149 (MQ=255)
               ggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGAtct                                                                                                                      <  2:196356‑M1/33‑1 (MQ=255)
               ggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGAtct                                                                                                                      <  3:34046‑M1/33‑1 (MQ=255)
               ggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGAtct                                                                                                                      <  1:289521‑M1/33‑1 (MQ=255)
               ggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGAtct                                                                                                                      <  5:280669‑M1/33‑1 (MQ=255)
               ggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGAtct                                                                                                                      <  1:219492‑M1/33‑1 (MQ=255)
               ggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGAtct                                                                                                                      <  5:129330‑M1/33‑1 (MQ=255)
               ggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGAtct                                                                                                                      <  1:105782‑M1/33‑1 (MQ=255)
               ggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCCCAAtct                                                                                                                      >  1:189210‑M1/116‑148 (MQ=255)
               ggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCACGCGGGCGGCGCTCGAtct                                                                                                                      <  1:117843‑M1/33‑1 (MQ=255)
               ggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcggggttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGAtct                                                                                                                      <  7:310501‑M1/33‑1 (MQ=255)
               ggacgcgcacgcgagccccctggagcacgcctcgatatctcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaACCGCGGAGCGCCGCCCGCGcggcgcttcgcgg                                                                                                                      >  8:237080‑M1/116‑142 (MQ=255)
                gacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTca                                                                                                                    <  7:98965‑M1/35‑1 (MQ=255)
                 acgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCgcgggcgg                                                                                                                                <  5:238683‑M1/23‑1 (MQ=255)
                 acgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGgcgc                                                                                                                             >  6:351974‑M1/114‑139 (MQ=255)
                 acgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGgcgc                                                                                                                             >  8:123897‑M1/114‑139 (MQ=255)
                   gcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCgcgg                                                                                                                                    >  6:32671‑M1/112‑130 (MQ=255)
                     gcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacgaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTc                                                                                                                           <  4:213608‑M1/28‑1 (MQ=255)
                     gcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCgcgggc                                                                                                                                  <  1:236888‑M1/21‑1 (MQ=255)
                       acgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGa                                                                                                                         <  5:224298‑M1/30‑1 (MQ=255)
                       acgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCCCGCGCGgcgcttcgc                                                                                                                        >  5:325549‑M1/108‑134 (MQ=255)
                        gcgcgaccgccttggggcatgtctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgactggtgcaggagggactggccgaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCt                                                                                                            <  1:40355‑M1/43‑1 (MQ=255)
                              ccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCCCGcgcggc                                                                                                                               <  5:303797‑M1/24‑1 (MQ=255)
                              ccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCCCGCGCGGcgcttcg                                                                                                                         <  2:37961‑M1/30‑4 (MQ=255)
                                ccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcgatcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGgcggcg                                                                                                                              <  3:295742‑M1/25‑1 (MQ=255)
                                    ggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgccttgcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCgcgg                                                                                                                                    >  6:56505‑M1/95‑113 (MQ=255)
                                    ggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctcgcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCgcgg                                                                                                                                    <  5:56505‑M1/19‑1 (MQ=255)
                                     gcgcacgcctcgatatcgcgtcgtacacacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGAtc                                                                                                                       >  3:35384‑M1/94‑125 (MQ=255)
                                            cctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCgcggg                                                                                                                                   >  1:342652‑M1/87‑106 (MQ=255)
                                             ctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTca                                                                                                                    >  2:322410‑M1/86‑120 (MQ=255)
                                             ctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCCCGCGCGgcgcttcgc                                                                                                                        <  7:331790‑M1/31‑5 (MQ=255)
                                               tcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgactggtgcaggagggactggccgaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCt                                                                                      <  4:243232‑M1/65‑1 (MQ=255)
                                                gatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGa                                                                                                                         >  1:324038‑M1/83‑112 (MQ=255)
                                                       cgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGggcggc                                                                                                                               >  4:114469‑M1/76‑99 (MQ=255)
                                                       cgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCCCGCGCGgcgcttcgc                                                                                                                        <  6:161028‑M1/31‑5 (MQ=255)
                                                           gtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCgcggg                                                                                                                                   >  7:82650‑M1/72‑91 (MQ=255)
                                                            tacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCgcgggc                                                                                                                                  >  6:22662‑M1/71‑91 (MQ=255)
                                                            tacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCgcgggc                                                                                                                                  <  5:22662‑M1/21‑1 (MQ=255)
                                                            tacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCgcgggc                                                                                                                                  >  8:49600‑M1/71‑91 (MQ=255)
                                                              acaaacaaggcggtcgcgcgcgactggtgcaggagggactggccgaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCATCCGGCAAACACCTCAAAATTCCAACGCCg                                                                      <  8:291034‑M1/81‑1 (MQ=255)
                                                               aaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGAtct                                                                                                                      <  7:223691‑M1/33‑1 (MQ=255)
                                                                 acaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCt                                                                                                                            >  2:272411‑M1/66‑92 (MQ=255)
                                                                   caaggcggtcgcgcgcgactggtgcaggagggactggccgaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCAc                                                                 >  1:65552‑M1/64‑149 (MQ=255)
                                                                              gcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGa                                                                                                                         <  8:179434‑M1/30‑1 (MQ=255)
                                                                                  gcctggcacaggagggattggccagaaacaaacgtcggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTca                                                                                                                    >  7:144898‑M1/49‑83 (MQ=255)
                                                                                                  gactggccgaaacaaacgtaggagcgggtttaCCCGCGTAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAATATTCCAACGCCGTGCACTCAAAATTCCAACGCCGCCCCCTCAAAATTc                                  >  4:166776‑M1/33‑149 (MQ=255)
                                                                                                           aaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAAcgccgc                         >  3:243232‑M1/24‑149 (MQ=255)
                                                                                                               aaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCACCGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTcc                     <  3:265619‑M1/130‑1 (MQ=255)
                                                                                                                                  CCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCGCCGCTCAATTCCAACGCCGCCTGCGCTCGAAAATCcgacaac   <  5:366874‑M1/148‑1 (MQ=255)
                                                                                                                                   CCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCcaacaaca  <  4:51259‑M1/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                  
GGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACA  >  NC_002947/4808044‑4808322

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