Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 31,028 T→C 7.7% W200R (TGG→CGG)  PP_0026 → cobalt/cadmium/zinc exporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294731,0280TC 7.7% 178.8 / 4.6 65W200R (TGG→CGG) PP_0026cobalt/cadmium/zinc exporter
Reads supporting (aligned to +/- strand):  ref base T (31/29);  new base C (2/3);  total (33/32)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.48e-01

AAAGGAGCCTACCTCGAAGTTTGGAGCGACATGCTCGGCTCGCTTGGCGTAATCATCGCGGCACTTATCATCCGCTTCACCGGCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑AGGGGGCTCGACGTTACGGCAATCGAGGCCACCATCCTCGGAGTAGATGGCGTGACGGACGTTCATGACTTGCACGTGTGGGCTGTCAGCAGTGG  >  NC_002947/30884‑31170
                                                                                                                                                 |                                                                                                                                               
aaaGGAGCCTACCTCGAAGTTTGGAGCGACATGCTCGGCTCGCTTGGCGTAATCATCGCGGCACTTATCATCCGCTTCACCGGCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCa                                                                                                                                             <  5:39641/149‑1 (MQ=255)
   ggAGCCTACCTCGAAGTTTGGAGCGACATGCTCGGCTCGCTTGGCGTAATCATCGCGGCACTTATCATCCGCTTCACCGGCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCACGTCCGCGGCAgaa                                                                                                                                          <  3:56424/149‑3 (MQ=255)
    gAGCCTACCTCGAAGTTTGGAGCGACATGCTCGGCTCGCTTGGCGTAATCATCGCGGCACTTATCATCCGCTTCACCGGCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCACGTCCGCGGCAgaac                                                                                                                                         <  7:47106/149‑4 (MQ=255)
     aGCCTACCTCGAAGTTTGGAGCGACATGCTCGGCTCGCTTGGCGTAATCATCGCGGCACTTATCATCCGCTTCACCGGCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGc                                                                                                                                        <  8:150851/149‑1 (MQ=255)
       ccTACCTCGAAGTTTGGAGCGACATGCTCGGCTCGCTTGGCGTAATCATCGCGGCACTTATCATCCGCTTCACCGGCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCACGTCCGCGGCAgaacc                                                                                                                                        <  1:157349/147‑5 (MQ=255)
       ccTACCTCGAAGTTTGGAGCGACATGCTCGGCTCGCTTGGCGTAATCATCGCGGCACTTATCATCCGCTTCACCGGCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCACGTCCGCGGCAgaacc                                                                                                                                        >  2:157349/1‑143 (MQ=255)
        cTACCTCGAAGTTTGGAGCGACATGCTCGGCTCGCTTGGCGTAATCATCGCGGCACTTATCATCCGCTTCACCGGCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTc                                                                                                                                     <  4:213827/149‑1 (MQ=255)
         tACCTCGAAGTTTGGAGCGACATGCTCGGCTCGCTTGGCGTAATCATCGCGGCACTTATCATCCGCTTCACCGGCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTcg                                                                                                                                    >  1:160979/1‑149 (MQ=255)
            cTCGAAGTTTGGAGCGACATGCTCGGCTCGCTTGGCGTAATCATCGCGGCACTTATCATCCGCTTCACCGGCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTcgcg                                                                                                                                  >  3:155782/1‑148 (MQ=255)
            cTCGAAGTTTGGAGCGACATGCTCGGCTCGCTTGGCGTAATCATCGCGGCACTTATCATCCGCTTCACCGGCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTcgcg                                                                                                                                  >  4:308703/1‑148 (MQ=255)
             tCGAAGTTTGGAGCGACATGCTCGGCTCGCTTGGCGTAATCATCGCGGCACTTATCATCCGCTTCACCGGCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGt                                                                                                                                           <  2:174425/138‑1 (MQ=255)
             tCGAAGTTTGGAGCGACATGCTCGGCTCGCTTGGCGTAATCATCGCGGCACTTATCATCCGCTTCACCGGCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGt                                                                                                                                           >  1:174425/1‑138 (MQ=255)
                  gTTTGGAGCGACATGCTCGGCTCGCTTGGCGTAATCATCGCGGCACTTATCATCCGCTTCACCGGCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCt                                                                                                                           >  1:213219/1‑149 (MQ=255)
                   tttGGAGCGACATGCTCGGCTCGCTTGGCGTAATCATCGCGGCACTTATCATCCGCTTCACCGGCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTg                                                                                                                          >  1:167660/1‑149 (MQ=255)
                     tGGAGCGACATGCTCGGCTCGCTTGGCGTAATCATCGCGGCACTTATCATCCGCTTCACCGGCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTgg                                                                                                                         <  6:96064/148‑1 (MQ=255)
                               tGCTCGGCTCGCTTGGCGTAATCATCGCGGCACTTATCATCCGCTTCACCGGCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATg                                                                                                              >  3:142931/1‑149 (MQ=255)
                                   cggctcgCTTGGCGTAATCATCGCGGCACTTATCATCCGCTTCACCGGCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAgg                                                                                                          >  5:267550/1‑149 (MQ=255)
                                    ggctcgCTTGGCGTAATCATCGCGGCACTTATCATCCGCTTCACCGGCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAggg                                                                                                         >  4:136004/1‑149 (MQ=255)
                                    ggctcgCTTGGCGTAATCATCGCGGCACTTATCATCCGCTTCACCGGCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAg                                                                                                           <  8:83528/147‑1 (MQ=255)
                                          cTTGGCGTAATCATCGCGGCACTTATCATCCGCTTCACCGGCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTccc                                                                                                   <  3:308703/149‑1 (MQ=255)
                                                gTAATCATCGCGGCACTTATCATCCGCTTCACCGGCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATcc                                                                                                                   >  1:32592/1‑127 (MQ=255)
                                                gTAATCATCGCGGCACTTATCATCCGCTTCACCGGCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATcc                                                                                                                   <  2:32592/127‑1 (MQ=255)
                                                gTAATCATCGCGGCACTTATCATCCGCTTCACCGGCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑Agggg                                                                                            >  2:123347/1‑149 (MQ=255)
                                                    tcatcGCGGCACTTATCATCCGCTTCACCGGCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑AGGGGGCTc                                                                                        >  8:58731/1‑149 (MQ=255)
                                                    tcatcGCGGCACTTATCAACCGCTTCACCGGCTGGACCTGGGACG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTTGCATTTGCTTCGCGAAAGACTGGGGATCCTCATTGAGGGTGTCCCG‑AGGGGGCTc                                                                                        >  5:290661/1‑149 (MQ=255)
                                                            gCACTTATCATCCGCTTCACCGGCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCTGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑AGGGGGCTCGACGTTAc                                                                                <  5:63920/149‑1 (MQ=255)
                                                             cACTTATCATCCGCTTCACCGGCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑AGGGGGCTCGACGTTACg                                                                               <  2:167660/149‑1 (MQ=255)
                                                              aCTTATCATCCGCTTCACCGGCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑AGGGGGCTCGACGTTACgg                                                                              <  6:122972/149‑1 (MQ=255)
                                                                ttATCATCCGCTTCACCGGCTGGACCTGGGTCGTACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑AGGGGGCTCGACGTTACgg                                                                              >  1:44373/1‑148 (MQ=255)
                                                                  atcatcCGCTTCACCGGCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑AGGGGGCTCGACGTTACGGCAAt                                                                          >  3:126393/1‑149 (MQ=255)
                                                                  atcatcCGCTTCACCGGCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑AGGGGGCTCGACGTTACGGCAAt                                                                          >  4:240577/1‑149 (MQ=255)
                                                                   tcatcCGCTTCACCGGCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑AGGGGGCTCGACGTTACGGCAAt                                                                          <  1:123347/148‑1 (MQ=255)
                                                                     atcCGCTTCACCGGCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑AGGGGGCTCGACGTTACGGCAATc                                                                         >  3:218995/1‑147 (MQ=255)
                                                                     atcCGCTTCACCGGCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑AGGGGGCTCGACGTTACGGCAATCGa                                                                       <  2:213219/149‑1 (MQ=255)
                                                                        cGCTTCACCGGCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑AGGGGGCTCGACGTTACGGCAATCGAGGc                                                                    <  2:160979/149‑1 (MQ=255)
                                                                               ccGGCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑Aggg                                                                                             >  6:202580/1‑117 (MQ=255)
                                                                               ccGGCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑Aggg                                                                                             <  5:202580/117‑1 (MQ=255)
                                                                               ccGGCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑AGGGGGCTCGACGTTACGGCAATCGAGGCCACCATc                                                             >  2:320664/1‑149 (MQ=255)
                                                                                  gCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑AGGGGGCTCGACGTTACGGCAATCGAGGCCACCATCCTc                                                          >  1:277166/1‑149 (MQ=255)
                                                                                   cTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑AGGGGGCTCGACGTTACGGCAATCGAGGCCACCATCCTCg                                                         <  4:218995/149‑1 (MQ=255)
                                                                                    tGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCAGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑AGGGGGCTCGACGTTACGGCAATCGAGGCCACCATCCTCg                                                         >  6:73502/1‑148 (MQ=255)
                                                                                      gACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑AGGGGGCTCGACGTTACGGCAATCGAGGCCACCATCCTCGGAg                                                      <  4:126393/149‑1 (MQ=255)
                                                                                       aCCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑AGGGGGCTCGACGTTACGGCAATCGAGGCCACCATCCTCGGAGt                                                     <  7:58731/149‑1 (MQ=255)
                                                                                        ccTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑AGGGGGCTCGACGTTACGGCAATCGAGGCCACCATCCTCGGAg                                                      >  2:232093/1‑147 (MQ=255)
                                                                                           gggTCG‑ACACGATTGTCGCTGTGGAGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑AGGGGGCTCGACGTTACGGCAATCGAGGCCACCATCCTCGGAGTAGAt                                                 >  1:324335/1‑149 (MQ=255)
                                                                                            ggTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCGTAGGGGGCTCGACGTTACGGCAATCGAGGCCACCATCCTCGGAGTAGa                                                  <  1:72304/148‑1 (MQ=255)
                                                                                                  acacGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑AGGGGGCTCGACGTTACGGCAATCGAGGCCACCATCCTCGGAGTAGATGGCGt                                            >  1:146615/1‑148 (MQ=255)
                                                                                                   cacGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑AGGGGGCTCGACGTTACGGCAATCGAGGCCACCATCCTCGGAGTAGATGGc                                              >  8:65592/1‑145 (MQ=255)
                                                                                                   cacGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑AGGGGGCTCGACGTTACGGCAATCGAGGCCACCATCCTCGGAGTAGATGGc                                              <  7:65592/145‑1 (MQ=255)
                                                                                                    acGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑AGGGGGCTCGACGTTACGGCAATCGAGGc                                                                    >  5:54627/1‑122 (MQ=255)
                                                                                                    acGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑AGGGGGCTCGACGTTACGGCAATCGAGGc                                                                    <  6:54627/122‑1 (MQ=255)
                                                                                                    acGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑AGGGGGCTCGACGTTACGGCAATCGAGGCCACCATCCTCGGAGTAGATGGCGTgac                                         <  3:136004/149‑1 (MQ=255)
                                                                                                      gATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑AGGGGGCTCGACGTTACGGCAATCGAGGCCACCATCCTCGGAGTAGATGGCGTgacg                                        >  1:129770/1‑148 (MQ=255)
                                                                                                              cTGTGGCGATCGGCCTTTGGGTTCTGCGGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑AGGGGGCTCGACGTTACGGCAATCGAGGCCACCATCCTCGGAGTAGATGGCGTGACGGACGTTCAt                               <  5:73502/149‑1 (MQ=255)
                                                                                                              cTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑AGGGGGCTCGACGTTACGGCAATCGAGGCCACCATCCTCGGAGTAGATGGCGTGACGGACGTTCAt                               >  7:298902/1‑149 (MQ=255)
                                                                                                                 tgGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑AGGGGGCTCGACGTTACGGCAATCGAGGCCACCATCCTCGGAGTAGATGGCGTGACGGACGTTCATGAc                            <  6:267550/149‑1 (MQ=255)
                                                                                                                   gCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑AGGGGGCTCGACGTTACGGCAATCGAGGCCACCATCCTCGGAGTAGATGGCGTGACGGACGTTCATGACtt                          <  3:240577/149‑1 (MQ=255)
                                                                                                                   gCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑AGGGGGCTCGACGTTACGGCAATCGAGGCCACCATCCTCGGAGTAGATGGCGTGACGGACGTTCATGACtt                          >  2:242130/1‑149 (MQ=255)
                                                                                                                             ttttGGTTCTGCCGCGTAGGTGGCAGTTGCTTCGCGATAGCCTGGGGATCCTATTGGAGGGTGTCCCG‑TGGTGGCTCTACGTTAAGGCAATCGAGGCCACCATCCTCGGAGTAGATGGCGTGACGGACGTTCATGACTTGCACGTGTgg                <  6:290661/149‑1 (MQ=255)
                                                                                                                             tttGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑AGGGGGCTCGACGTTACGGCAATCGAGGCCACCATCCTCGGAGTAGATGGCGTGACGGACGTTCATGACTTGCACGTGTgg                >  2:113937/1‑149 (MQ=255)
                                                                                                                              gaagcaacTGCCACGTCCGCGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑AGGGGGCTCGACGTTACGGCAATCGAGGCCACCATCCTCGGAGTAGATGGCGTGACGGACGTTCATGACTTGCACGTGTggg               >  5:45373/8‑149 (MQ=255)
                                                                                                                               tGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑AGGGGGCTCGACGTTACGGCAATCGAGGCCACCATCCTCGGAGTAGATGGCGTGACGGACGTTCATGACTTGCACGTGTggg               <  2:277166/148‑1 (MQ=255)
                                                                                                                                         cgcgGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑AGGGGGCTCGACGTTACGGCAATCGAGGCCACCATCCTCGGAGTAGATGGCGTGACGGACGTTCATGACTTGCACGTGTGGGCTGt           >  6:87130/1‑142 (MQ=255)
                                                                                                                                         cgcgGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑AGGGGGCTCGACGTTACGGCAATCGAGGCCACCATCCTCGGAGTAGATGGCGTGACGGACGTTCATGACTTGCACGTGTGGGCTGt           <  5:87130/142‑1 (MQ=255)
                                                                                                                                           cgGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑AGGGGGCTCGACGTTACGGCAATCGAGGCCACCATCCTCGGAGTAGATGCCGTGACGGACGTTCATGACTTGCACGTGTGGGCTGTCAGCAGTgg  <  1:232093/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
AAAGGAGCCTACCTCGAAGTTTGGAGCGACATGCTCGGCTCGCTTGGCGTAATCATCGCGGCACTTATCATCCGCTTCACCGGCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑AGGGGGCTCGACGTTACGGCAATCGAGGCCACCATCCTCGGAGTAGATGGCGTGACGGACGTTCATGACTTGCACGTGTGGGCTGTCAGCAGTGG  >  NC_002947/30884‑31170

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: ‑