Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,808,189 G→C 14.6% intergenic (‑298/‑257) 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,1890GC14.6% 786.1 / 115.7 308intergenic (‑298/‑257)PP_4234/dsbD‑IIbifunctional aldehyde oxidase/xanthine dehydrogenase/thiol/disulfide interchange protein 2
Reads supporting (aligned to +/- strand):  ref base G (96/167);  new base C (29/16);  total (125/183)
Fisher's exact test for biased strand distribution p-value = 5.31e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.19e-03

Only 100 of 446 total aligned reads displayed.

GGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACA  >  NC_002947/4808041‑4808322
                                                                                                                                                    |                                                                                                                                     
gggccagttattcctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCcgcg                                                                                                                                       <  8:52684‑M1/16‑1 (MQ=255)
 ggccagttattcctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCcgcgc                                                                                                                                      <  3:141393‑M1/17‑1 (MQ=255)
   ccagttattcctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCgcgg                                                                                                                                    <  8:75227‑M1/19‑1 (MQ=255)
   ccagttattcctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCgcgg                                                                                                                                    >  4:276072‑M1/131‑149 (MQ=255)
    cagttattcctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttcCCCCGGAGCCGCCGCgcggg                                                                                                                                   <  3:205999‑M1/20‑1 (MQ=255)
    cagttattcctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCgcggg                                                                                                                                   <  1:141944‑M1/20‑1 (MQ=255)
    cagttattcctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCgcggg                                                                                                                                   <  5:360933‑M1/20‑1 (MQ=255)
    cagttattcctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcaggtttaCCCGCGAAGCGCCGCgcggg                                                                                                                                   <  7:50495‑M1/20‑1 (MQ=255)
     agttattcctgtgggacgcgcacgcgagccacctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCgcgggc                                                                                                                                  <  8:172411‑M1/21‑1 (MQ=255)
      gttattcctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacatacggaggagcgggtttgCCCGCGAAGCGCCGCgcgggcg                                                                                                                                 <  5:263750‑M1/22‑1 (MQ=255)
      gttattcctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCgcgggcg                                                                                                                                 <  2:267019‑M1/22‑1 (MQ=255)
      gttattcctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCgcgggcg                                                                                                                                 >  6:205835‑M1/128‑149 (MQ=255)
      gttattcctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCgcgggcg                                                                                                                                 <  7:190960‑M1/22‑1 (MQ=255)
       ttattcctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCccgcgcgg                                                                                                                                <  5:86062‑M1/23‑1 (MQ=255)
        tattcctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGggcggc                                                                                                                               <  2:134935‑M1/24‑1 (MQ=255)
        tattcctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgcgtttaCCCGCGAAGCGCCGCgcgggcgg                                                                                                                                >  1:202376‑M1/126‑148 (MQ=255)
         attcctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCgcggg                                                                                                                                   >  1:132520‑M1/125‑144 (MQ=255)
         attcctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGgcggcg                                                                                                                              <  8:338908‑M1/25‑1 (MQ=255)
         attcctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGgcggcg                                                                                                                              <  2:21121‑M1/25‑1 (MQ=255)
         attcctgtgggacgcgcacgcgagccccctggcccacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCcgcgc                                                                                                                                      >  5:314149‑M1/125‑141 (MQ=255)
          ttcctgtgggacgcgcacgcgagccccctcgcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGgcgc                                                                                                                             >  5:254210‑M1/124‑149 (MQ=255)
           tcctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCcgcgcg                                                                                                                                     >  4:258327‑M1/123‑140 (MQ=255)
           tcctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCt                                                                                                                            <  3:186638‑M1/27‑1 (MQ=255)
           tcctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCt                                                                                                                            <  1:180277‑M1/27‑1 (MQ=255)
            cctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggtggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGgcgc                                                                                                                             >  3:253802‑M1/122‑147 (MQ=255)
            cctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCCCGCGCGGCGCt                                                                                                                            >  6:86062‑M1/122‑148 (MQ=255)
            cctgtgggacgcgcacgcgagccccctggcccacgcctctatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCt                                                                                                                            >  3:226620‑M1/122‑148 (MQ=255)
             ctgttggactcgcacgcgagccacctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaaaaaacgtaggagccggtttaCCCGCGAAGCGCCGAGCGGGCGGCGCTc                                                                                                                           <  4:224591‑M1/28‑1 (MQ=255)
             ctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTc                                                                                                                           <  2:67873‑M1/28‑1 (MQ=255)
             ctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTc                                                                                                                           <  2:148082‑M1/28‑1 (MQ=255)
             ctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTc                                                                                                                           >  4:174892‑M1/121‑148 (MQ=255)
             ctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTc                                                                                                                           >  6:315110‑M1/121‑148 (MQ=255)
             ctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTc                                                                                                                           <  8:246674‑M1/28‑1 (MQ=255)
             ctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTc                                                                                                                           <  2:253454‑M1/28‑1 (MQ=255)
             ctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGCGCGGCGCTc                                                                                                                           <  6:360965‑M1/28‑1 (MQ=255)
             ctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggacgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTc                                                                                                                           <  2:98898‑M1/28‑1 (MQ=255)
             ctgtgggacgcgcacgagagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagccggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTc                                                                                                                           <  4:99013‑M1/28‑1 (MQ=255)
              tgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGggcggcgcccat                                                                                                                         <  3:38228‑M1/30‑5 (MQ=255)
              tgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGa                                                                                                                         <  6:156318‑M1/30‑1 (MQ=255)
              tgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGa                                                                                                                         <  7:173468‑M1/30‑1 (MQ=255)
              tgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGa                                                                                                                         <  4:205417‑M1/30‑1 (MQ=255)
              tgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGa                                                                                                                         <  3:28886‑M1/30‑1 (MQ=255)
               gtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCcgcgcg                                                                                                                                     <  5:149693‑M1/18‑1 (MQ=255)
                tgggacgcgcccgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCCCGCGCggcgcttcgcg                                                                                                                       >  5:274895‑M1/118‑144 (MQ=255)
                tgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGGAGCGCCGCCCGCGCggcgcttcgcg                                                                                                                       >  8:318631‑M1/118‑144 (MQ=255)
                tgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGa                                                                                                                         >  5:217791‑M1/118‑147 (MQ=255)
                tgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGAtc                                                                                                                       >  4:308270‑M1/118‑149 (MQ=255)
                tgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGAtc                                                                                                                       <  6:177797‑M1/32‑1 (MQ=255)
                tgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGCGGCGCTCGAt                                                                                                                        >  5:299678‑M1/118‑147 (MQ=255)
                ggggacgcgcacgcgagccccctggcgcacgcctcgatctcgcgtcgtacaaacaaggcgggcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTcgacc                                                                                                                       >  4:123904‑M1/118‑147 (MQ=255)
                 gggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGcggcgcttcgcgg                                                                                                                      >  1:21121‑M1/117‑143 (MQ=255)
                  ggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagtgattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGAtc                                                                                                                       <  6:302706‑M1/32‑1 (MQ=255)
                  ggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGggcggc                                                                                                                               <  7:278963‑M1/24‑1 (MQ=255)
                  ggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGAtctc                                                                                                                     >  4:28886‑M1/116‑149 (MQ=255)
                  ggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGAtct                                                                                                                      <  4:66077‑M1/33‑1 (MQ=255)
                  ggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGAtct                                                                                                                      <  4:146182‑M1/33‑1 (MQ=255)
                  ggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGAtct                                                                                                                      <  2:196356‑M1/33‑1 (MQ=255)
                  ggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGAtct                                                                                                                      <  1:289521‑M1/33‑1 (MQ=255)
                  ggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGAtct                                                                                                                      <  1:105782‑M1/33‑1 (MQ=255)
                  ggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCCCGCGcggcgcttcgcgg                                                                                                                      <  6:288340‑M1/33‑7 (MQ=255)
                  ggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtgtaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGAtctc                                                                                                                     <  7:78010‑M1/34‑1 (MQ=255)
                  ggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcggggttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGAtct                                                                                                                      <  7:310501‑M1/33‑1 (MQ=255)
                  ggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattgaccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGAtct                                                                                                                      <  3:312168‑M1/33‑1 (MQ=255)
                  gacgcgcacgcgagccccctggcgcacgcctcggatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGAtct                                                                                                                      <  1:225058‑M1/33‑1 (MQ=255)
                   gacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTca                                                                                                                    >  5:225361‑M1/115‑149 (MQ=255)
                   gacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAGCGCCGCCCGCGcggcgcttcgcgg                                                                                                                      <  7:115934‑M1/32‑7 (MQ=255)
                   gacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAGCGCCGCCCGCGcggcgcttcgcgg                                                                                                                      >  8:115934‑M1/115‑140 (MQ=255)
                    acgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCgcgggcgg                                                                                                                                >  6:238683‑M1/114‑136 (MQ=255)
                    acgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGgcgc                                                                                                                             >  6:351974‑M1/114‑139 (MQ=255)
                    acgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGgcgc                                                                                                                             <  5:351974‑M1/26‑1 (MQ=255)
                     cgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCgcgggcg                                                                                                                                 >  6:136156‑M1/113‑134 (MQ=255)
                      gcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCgcgg                                                                                                                                    <  5:32671‑M1/19‑1 (MQ=255)
                      gcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTc                                                                                                                           <  8:199017‑M1/28‑1 (MQ=255)
                       cgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGAt                                                                                                                        >  3:58729‑M1/111‑141 (MQ=255)
                          acgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCCCGCGCGgcgcttcgc                                                                                                                        <  6:325549‑M1/31‑5 (MQ=255)
                           gcgcgaccgccttggggcatgtctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgactggtgcaggagggactggccgaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCt                                                                                                            <  1:40355‑M1/43‑1 (MQ=255)
                           cgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGAt                                                                                                                        <  5:133941‑M1/31‑1 (MQ=255)
                                 ccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCCCGCGCGGcgcttcg                                                                                                                         <  2:37961‑M1/30‑4 (MQ=255)
                                    cctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGccgcgcgggtaaa                                                                                                                               >  3:344116‑M1/98‑117 (MQ=255)
                                    cctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGccgcgcgggtaaa                                                                                                                               <  4:344116‑M1/24‑5 (MQ=255)
                                       ggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgccttgcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCgcgg                                                                                                                                    >  6:56505‑M1/95‑113 (MQ=255)
                                         cgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGAtc                                                                                                                       >  6:226783‑M1/93‑124 (MQ=255)
                                               cctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCgcggg                                                                                                                                   >  1:342652‑M1/87‑106 (MQ=255)
                                                tctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgactggtgcaggagggactggccgaaacaaacgtaggagcgggattaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCcggcaaagttc                                                                                       <  2:321314‑M1/64‑5 (MQ=255)
                                                ctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCCCGCGCGgcgcttcgc                                                                                                                        <  7:331790‑M1/31‑5 (MQ=255)
                                                 tcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGGAGCGCCGCCCGCGCGGcgcttcg                                                                                                                         <  6:25675‑M1/30‑4 (MQ=255)
                                                   gatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGa                                                                                                                         >  1:324038‑M1/83‑112 (MQ=255)
                                                                 acaaacaaggcggtcgcgcgcgactggtgcaggagggactggccgaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCATCCGGCAAACACCTCAAAATTCCAACGCCg                                                                      <  8:291034‑M1/81‑1 (MQ=255)
                                                                  aaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGgcgc                                                                                                                             <  5:231720‑M1/26‑1 (MQ=255)
                                                                  aaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGAtct                                                                                                                      <  7:223691‑M1/33‑1 (MQ=255)
                                                                  aaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGAtct                                                                                                                      >  8:223691‑M1/68‑100 (MQ=255)
                                                                    acaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTc                                                                                                                           >  5:46871‑M1/66‑93 (MQ=255)
                                                                         gcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGgcgc                                                                                                                             <  6:119507‑M1/26‑1 (MQ=255)
                                                                                 gcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTc                                                                                                                           >  7:9859‑M1/53‑80 (MQ=255)
                                                                                     gcctggcacaggagggattggccagaaacaaacgtcggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTca                                                                                                                    >  7:144898‑M1/49‑83 (MQ=255)
                                                                                     gcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTca                                                                                                                    <  8:144898‑M1/35‑1 (MQ=255)
                                                                                      gactggtgcaggagggactggccgaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTcc                                                       >  5:20553‑M1/48‑143 (MQ=255)
                                                                                                          gccgaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAAc                              >  2:331644‑M1/28‑148 (MQ=255)
                                                                                                                  aaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCACCGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTcc                     <  3:265619‑M1/130‑1 (MQ=255)
                                                                                                                       taggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCg                >  1:905‑M1/15‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                     
GGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACA  >  NC_002947/4808041‑4808322

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: