Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,808,234 G→C 8.4% 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,2340GC 8.4% 221.5 / 9.7 71intergenic (‑298/‑257)PP_4234/dsbD‑IIbifunctional aldehyde oxidase/xanthine dehydrogenase/thiol/disulfide interchange protein 2
Reads supporting (aligned to +/- strand):  ref base G (26/39);  new base C (3/3);  total (29/42)
Fisher's exact test for biased strand distribution p-value = 6.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01

GGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGACCCT  >  NC_002947/4808086‑4808375
                                                                                                                                                    |                                                                                                                                             
gggccagttattcctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCcgcgc                                                                                                                                              <  2:1644849‑M1/17‑1 (MQ=255)
  gccagttattcctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCgcgg                                                                                                                                            <  1:1394342‑M1/19‑1 (MQ=255)
   ccagttattcctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCgcggg                                                                                                                                           <  1:559718‑M1/20‑1 (MQ=255)
     agttattcctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCgcgggcg                                                                                                                                         >  2:559718‑M1/129‑150 (MQ=255)
            cctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCt                                                                                                                                    <  1:336238‑M1/27‑1 (MQ=255)
            cctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCt                                                                                                                                    >  2:336238‑M1/122‑148 (MQ=255)
             ctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggggctggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGa                                                                                                                                 <  1:2240716‑M1/30‑1 (MQ=255)
             ctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGa                                                                                                                                 <  2:2111020‑M1/30‑1 (MQ=255)
               gtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGAtc                                                                                                                               <  1:464185‑M1/32‑1 (MQ=255)
                                    cgttggggcagggcccgatatcgcgtcgtacacacaaggcggtcgcgcgcgactggtgcaggagggactggacgaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAAtcctc                                                                                                          >  2:1073386‑M1/98‑150 (MQ=255)
                                    cctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTca                                                                                                                            >  1:1242989‑M1/98‑132 (MQ=255)
                                    cctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTca                                                                                                                            <  2:1242989‑M1/35‑1 (MQ=255)
                                                tctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgactggtgcaggagggactggccgaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCt                                                                                              <  1:1957022‑M1/65‑1 (MQ=255)
                                                tctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgactggtgcaggagggactggccgaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCt                                                                                              >  1:1570413‑M1/86‑150 (MQ=255)
                                                   cgatatcgcgtcgtacaaacaaggcggtcgcgcgcgactggtgcaggagggactggccgaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAAtcctc                                                                                                          <  2:791117‑M1/53‑1 (MQ=255)
                                                   cgatatcgcgtcgtacaaacaaggcggtcgcgcgcgactggtgcaggagggactggccgaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAAtcctc                                                                                                          >  1:791117‑M1/83‑135 (MQ=255)
                                                   cgatatcgcgtcgtacaaacaaggcggtcgcgcgcgactggtgcaggagggactggccgaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAcacaccccca                                                                                           <  2:373919‑M1/68‑5 (MQ=255)
                                                    gatatcgcgtcgtacaaacaaggcggtcgcgcgcgactggtgcaggagggactggccgaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCaaa                                                                                          <  2:715366‑M1/69‑1 (MQ=255)
                                                               tacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGGAGCGCCGCCCGCGCGGCgcttc                                                                                                                                  >  2:2004622‑M1/71‑97 (MQ=255)
                                                               tacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGGAGCGCCGCCCGCGCGGCgcttc                                                                                                                                  <  1:2004622‑M1/29‑3 (MQ=255)
                                                               gtacaaacaaggcggtcgcgcgcgactggtgcaggagggactggccgaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGgcgc                                                                                                                     <  2:1608394‑M1/42‑1 (MQ=255)
                                                               gtacaaacaaggcggtcgcgcgcgactggtgcaggagggactggccgaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGgcgc                                                                                                                     >  1:1608394‑M1/71‑112 (MQ=255)
                                                                 acaaacaaggcggtcgcgcgcgactggtgcaggagggactggccgaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGAt                                                                                                                                >  1:1245690‑M1/69‑99 (MQ=255)
                                                                 acaaacaaggcggtcgcgcgcgactggtgcaggagggactggccgaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGAt                                                                                                                                <  2:1245690‑M1/31‑1 (MQ=255)
                                                                 acaaacaaggcggtcgcgcgcgactggtgcaggagggactggccgaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGt                                                                             <  1:2060640‑M1/82‑1 (MQ=255)
                                                                   aaacaaggcggtcgcgcgcgactggtgcaggagggactggccgaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTcc                                                                           <  2:1645700‑M1/84‑1 (MQ=255)
                                                                     acaaggcggtcgcgcgcgactggtgcaggagggactggccgaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCAGTCCAc                                                                         <  2:285889‑M1/86‑1 (MQ=255)
                                                                        ggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTc                                                                                                                                   >  1:1932824‑M1/62‑89 (MQ=255)
                                                                        ggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTc                                                                                                                                   <  2:1932824‑M1/28‑1 (MQ=255)
                                                                         gcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGgcgc                                                                                                                                     >  2:1732302‑M1/61‑86 (MQ=255)
                                                                         gcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGgcgc                                                                                                                                     <  1:1732302‑M1/26‑1 (MQ=255)
                                                                              tcgcgcgcgactggtgcaggagggactggccgaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGaaaa                                                                                                               <  1:1876344‑M1/48‑1 (MQ=255)
                                                                              tcgcgcgcgactggtgcaggagggactggccgaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGaaaa                                                                                                               >  2:1876344‑M1/56‑103 (MQ=255)
                                                                                  gcgcgactggtgcaggagggactggccgaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAAc                                                            <  2:1223700‑M1/99‑1 (MQ=255)
                                                                                        ctggtgcaggagggactggccgaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACgccgcc                                                      >  2:1957022‑M1/46‑150 (MQ=255)
                                                                                           gtgcaggagggactggccgaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTc                                                                       <  2:1872214‑M1/88‑1 (MQ=255)
                                                                                           gtgcaggagggactggccgaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTc                                                                       >  1:1872214‑M1/43‑130 (MQ=255)
                                                                                           gtgcaggagggactggccgaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGccccc                                                   <  2:854288‑M1/108‑1 (MQ=255)
                                                                                                 gagggactggccgaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTcc                                                               <  1:2222121‑M1/96‑1 (MQ=255)
                                                                                                 gagggactggccgaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTcc                                                               >  2:2222121‑M1/37‑132 (MQ=255)
                                                                                                 gagggactggccgaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCaaaa                                             <  1:297018‑M1/114‑1 (MQ=255)
                                                                                                        tggccgaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAAc                                      <  1:2094582‑M1/121‑1 (MQ=255)
                                                                                                          gccgaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAAcgc                                    <  1:770470‑M1/123‑1 (MQ=255)
                                                                                                          gccgaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGAGCTAAATTTCACAGGCGCTGGAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAAcgc                                    <  1:1455165‑M1/123‑1 (MQ=255)
                                                                                                           ccgaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAAcgcc                                   <  2:1040297‑M1/124‑1 (MQ=255)
                                                                                                                  aaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTccc                            >  1:852160‑M1/20‑150 (MQ=255)
                                                                                                                  aaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGCAAATCCTCCGGCAAACACATCAAAATGCCAACGCCGTCCACGCAGAAGTCCATCGCCGCCCCCTCAAAATCCCAACGCCGCCAccc                            <  2:2015582‑M1/131‑1 (MQ=255)
                                                                                                                    acgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCt                          <  2:1634077‑M1/133‑1 (MQ=255)
                                                                                                                    acgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCt                          >  1:1185401‑M1/18‑150 (MQ=255)
                                                                                                                    acgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCt                          <  2:1130267‑M1/133‑1 (MQ=255)
                                                                                                                          gagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCCCTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGaaaa                    >  1:1645700‑M1/12‑150 (MQ=255)
                                                                                                                          gagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGaaaa                    >  1:1223700‑M1/12‑150 (MQ=255)
                                                                                                                               ggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGCTCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCcaa               >  1:373919‑M1/7‑150 (MQ=255)
                                                                                                                               ggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGccccc                                                   >  2:303016‑M1/7‑114 (MQ=255)
                                                                                                                               ggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGccccc                                                   <  1:303016‑M1/108‑1 (MQ=255)
                                                                                                                               ggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCACCGCCGCCTCCCCTCGAAAATCcaa               >  1:548696‑M1/7‑150 (MQ=255)
                                                                                                                               ggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCcaa               >  1:1040297‑M1/7‑150 (MQ=255)
                                                                                                                                gtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAGACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCcaac              <  2:1220498‑M1/145‑1 (MQ=255)
                                                                                                                                gtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCcaac              <  1:1641534‑M1/145‑1 (MQ=255)
                                                                                                                                     CCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCTCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCcaacaacaa         >  1:1338602‑M1/1‑150 (MQ=255)
                                                                                                                                     CCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCTCCCCCTCAAAATTCCAACGCCGCCTCCCCTCAAAAATCcaacaacaa         <  1:382426‑M1/150‑1 (MQ=255)
                                                                                                                                      CCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCt                          >  2:1273020‑M1/1‑132 (MQ=255)
                                                                                                                                      CCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCt                          <  1:1273020‑M1/132‑1 (MQ=255)
                                                                                                                                      CCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGaaaa                    >  2:82695‑M1/1‑138 (MQ=255)
                                                                                                                                      CCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGaaaa                    <  1:82695‑M1/138‑1 (MQ=255)
                                                                                                                                      CCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACaaa        <  1:174797‑M1/150‑1 (MQ=255)
                                                                                                                                      CCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACAAATCAAAAAACAAACGTGGTGCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACaaa        <  1:2035109‑M1/150‑1 (MQ=255)
                                                                                                                                       CGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCaa                                                             >  2:1865878‑M1/1‑96 (MQ=255)
                                                                                                                                       CGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCaa                                                             <  1:1865878‑M1/96‑1 (MQ=255)
                                                                                                                                       CGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAg       <  2:948365‑M1/150‑1 (MQ=255)
                                                                                                                                        atcgAGCGCCGCCCGCGCGGCGCTTGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCaaaa                                             >  2:67452‑M1/5‑111 (MQ=255)
                                                                                                                                        atcgAGCGCCGCCCGCGCGGCGCTTGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCaaaa                                             <  1:67452‑M1/107‑1 (MQ=255)
                                                                                                                                            AGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCGCAATTCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGACCCt  >  2:187611‑M1/1‑150 (MQ=255)
                                                                                                                                            AGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAAt                   >  1:1870721‑M1/1‑133 (MQ=255)
                                                                                                                                            AGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAAt                   <  2:1870721‑M1/133‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                             
GGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGACCCT  >  NC_002947/4808086‑4808375

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