Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 239,060 G→A 5.2% R47R (CGG→CGA)  hemD → uroporphyrinogen‑III synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947239,0600GA 5.2% 299.1 / 3.8 95R47R (CGG→CGA) hemDuroporphyrinogen‑III synthase
Reads supporting (aligned to +/- strand):  ref base G (42/48);  new base A (3/2);  total (45/50)
Fisher's exact test for biased strand distribution p-value = 6.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.69e-01

GGTCACCCGTGAGCCAATGGCGCCTGTTGCTGACCCGGCCTGCCGAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGT  >  NC_002947/238912‑239206
                                                                                                                                                    |                                                                                                                                                     
ggTCACCCGTGAGCCAATGGCGCCTGTTGCTGACCCGGCCTGCCGAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCgg                                                                                                                                                       >  7:210588/1‑149 (MQ=255)
ggTCACCCGTGAGCCAATGGCGCCTGTTGCTGACCCGGCCTGCCGAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCgg                                                                                                                                                       >  8:131335/1‑149 (MQ=255)
ggTCACCCGTGAGCCAATGGCGCCTGTTGCTGACCCGGCCTGCCGAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCgg                                                                                                                                                       <  8:210588/149‑1 (MQ=255)
    aCCCGTGAGCCAATGGCGCCTGTTGCTGACCCGGCCTGCCGAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGt                                                                                                                                                   <  7:107354/149‑1 (MQ=255)
     cccGTGAGCCAATGGCGCCTGTTGCTGACCCGGCCTGCCGAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑c                                                                                                                                                 <  2:391011/149‑1 (MQ=255)
       cGTGAGCCAATGGCGCCTGTTGCTGACCCGGCCTGCCGAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGc                                                                                                                                               >  8:280784/1‑149 (MQ=255)
             ccAATGGCGCCTGTTGCTGACCCGGCCTGCCGAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAg                                                                                                                                         >  7:112226/1‑149 (MQ=255)
             ccAATGGCGCCTGTTGCTGACCCGGCCTGCCGAGGACTGCGCGGCGCTGCCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCTTCCAGCGACAGC‑CGCTgtggg                                                                                                                                         >  3:242999/1‑145 (MQ=255)
                   gcgcCTGTTGCTGACCCGGCCTGCCGAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGc                                                                                                                                   <  6:48903/149‑1 (MQ=255)
                     gcCTGTTGCTGACCCGGCCTGCCGAGGACTGCGCGGCGCTGCCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCTTCCAGCGACAGC‑CGCTGTGGGGCTTGC‑ac                                                                                                                                <  4:242999/149‑2 (MQ=38)
                        tGTTGCTGACCCGGCCTGCCGAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGAc                                                                                                                             >  8:32157/1‑149 (MQ=255)
                          ttGCTGACCCGGCCTGCCGAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACtt                                                                                                                           <  3:210434/149‑1 (MQ=255)
                              tGACCCGGCCTGCCGAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGc                                                                                                                                   >  6:103467/1‑138 (MQ=255)
                              tGACCCGGCCTGCCGAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGc                                                                                                                                   <  5:103467/138‑1 (MQ=255)
                               gACCCGGCCTGCCGAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAgg                                                                                                                      >  8:289628/1‑149 (MQ=255)
                                  ccGGCCTGCCGAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGa                                                                                                                   <  7:131335/149‑1 (MQ=255)
                                  ccGGCCTGCCGAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGa                                                                                                                   <  2:263756/149‑1 (MQ=255)
                                    ggCCTGCCGAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGAt                                                                                                                  >  1:110023/1‑148 (MQ=255)
                                          ccGAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATtgtg                                                                                                           <  7:289628/149‑1 (MQ=255)
                                           cGAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTgg                                                                                                          >  3:114211/1‑149 (MQ=255)
                                            gAGGACTGCGCGGGGCTGGCGCATAGCCTGGCGGCAGCGGGAGTGGGCAGCTGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGt                                                                                                         <  2:22245/149‑1 (MQ=255)
                                            gAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGt                                                                                                         >  1:155196/1‑149 (MQ=255)
                                            gAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGt                                                                                                         <  3:194770/149‑1 (MQ=255)
                                            gAGGACTGCGCGGCGCTGCCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCAAGGCGATCATTGTGGt                                                                                                         >  8:375163/1‑149 (MQ=255)
                                                aCTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGc                                                                                                     <  7:379555/149‑1 (MQ=255)
                                                 cTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGc                                                                                                     >  3:76151/1‑148 (MQ=255)
                                                  tGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACt                                                                                                                            >  2:164605/1‑124 (MQ=255)
                                                  tGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACt                                                                                                                            <  1:164605/124‑1 (MQ=255)
                                                              gcgcAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCg                                                                                       >  4:269118/1‑149 (MQ=255)
                                                              gcgcAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCg                                                                                       >  2:1441/1‑149 (MQ=255)
                                                                gcAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCGACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGc                                                                                     >  1:128323/1‑149 (MQ=255)
                                                                gcAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGc                                                                                     >  6:405342/1‑149 (MQ=255)
                                                                gcAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGc                                                                                     >  4:221752/1‑149 (MQ=255)
                                                                 cAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCt                                                                                    <  7:32157/149‑1 (MQ=255)
                                                                 cAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCt                                                                                    >  5:186157/1‑149 (MQ=255)
                                                                   aaGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTTCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGt                                                                                  <  3:16022/149‑1 (MQ=255)
                                                                    aGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGACAGCAAGCCGGCCGCCCGGCTGtt                                                                                 <  6:251559/149‑1 (MQ=255)
                                                                     gCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTACGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGtt                                                                                 <  1:121380/149‑1 (MQ=255)
                                                                     gCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTg                                                                                <  2:155196/149‑1 (MQ=255)
                                                                         ggcggcAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGa                                                                                                                                           <  5:183445/87‑1 (MQ=255)
                                                                         ggcggcAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGa                                                                                                                                           >  6:183445/1‑87 (MQ=255)
                                                                          gcggcAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑Ga                                                                          <  4:114211/149‑1 (MQ=255)
                                                                          gcggcAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑Ga                                                                          <  5:271807/149‑1 (MQ=255)
                                                                             gcAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCa                                                                       >  5:185278/1‑149 (MQ=255)
                                                                                  gggAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTgg                                                                  <  4:76151/149‑1 (MQ=255)
                                                                                    gAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGcc                                                                >  4:80375/1‑149 (MQ=255)
                                                                                    gAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCAATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGcc                                                                >  6:198805/1‑149 (MQ=255)
                                                                                         ggCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTg                                                                                >  2:278536/1‑129 (MQ=255)
                                                                                         ggCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTg                                                                                <  1:278536/129‑1 (MQ=255)
                                                                                         ggCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCgaga                                                            <  3:221752/148‑1 (MQ=255)
                                                                                         ggCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGAc                                                           >  6:61325/1‑149 (MQ=255)
                                                                                         ggCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGAc                                                           <  6:185278/149‑1 (MQ=255)
                                                                                             gcagTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGACCAGc                                                       <  4:208996/149‑1 (MQ=255)
                                                                                                 ttGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGACCAGCATc                                                    >  3:204500/1‑148 (MQ=255)
                                                                                                    cctgccCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGACCAGCATCCAGc                                                >  5:59643/1‑149 (MQ=255)
                                                                                                     ctgccCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGAc                                                           <  6:28654/137‑1 (MQ=255)
                                                                                                     ctgccCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGACCAGCATCCAGcc                                               <  3:145776/149‑1 (MQ=255)
                                                                                                     ctgccCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGACCAGCATCCAGcc                                               >  1:395418/1‑149 (MQ=255)
                                                                                                       gccCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCAAAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAACCGGCCGCCCGGCTGTTGCTT‑GAGCGTCTGGCCGAGACCAGCATCCAGCCTc                                             >  7:109103/1‑149 (MQ=255)
                                                                                                        ccccTGCTGGAGAATGAACCCATTACCGAGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGACCAGCATCCAGCCTc                                             <  8:173995/148‑1 (MQ=255)
                                                                                                          ccTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTAGAGCATCTGGCCGAGACCAGCATCCAGCCTCCg                                           <  1:259152/149‑1 (MQ=255)
                                                                                                               tgGAGATTGAAGCCATTACCGTGCTAGCGGCTCAGCGGGTGT‑CGCTGGAAGGCTTGC‑ATGACTTCGAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGa                                     <  4:204500/149‑1 (MQ=255)
                                                                                                                ggAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAgg                                                                                                                      >  2:109340/1‑68 (MQ=255)
                                                                                                                ggAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAgg                                                                                                                      <  1:109340/68‑1 (MQ=255)
                                                                                                                ggAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAg                                    >  8:185372/1‑149 (MQ=255)
                                                                                                                 gagaTTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGa                                     <  5:380588/147‑1 (MQ=255)
                                                                                                                     ttGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTgg                               <  1:1441/149‑1 (MQ=255)
                                                                                                                      tGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCTATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTgg                               <  5:61325/148‑1 (MQ=255)
                                                                                                                      tGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGt                              >  6:320611/1‑149 (MQ=255)
                                                                                                                         aCCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCa                           >  7:339987/1‑149 (MQ=255)
                                                                                                                         aCCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCa                           <  7:185372/149‑1 (MQ=255)
                                                                                                                          cccATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCAc                          <  1:187715/149‑1 (MQ=255)
                                                                                                                           ccATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCAcc                         >  6:76931/1‑149 (MQ=255)
                                                                                                                             aTTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGc                                       <  7:276089/133‑1 (MQ=255)
                                                                                                                             aTTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGc                                       >  8:276089/1‑133 (MQ=255)
                                                                                                                             aTTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGt                       >  6:32285/1‑149 (MQ=255)
                                                                                                                             aTTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGACTGGTTc                            <  4:49335/144‑1 (MQ=255)
                                                                                                                             aTTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGACTGGTTc                            >  3:49335/1‑144 (MQ=255)
                                                                                                                              ttACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTg                      <  3:237394/149‑1 (MQ=255)
                                                                                                                                gtcatGCAAGCCTTCCAGCGACAGC‑CGCTGTGGGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGAc                                                           <  6:93235/106‑1 (MQ=21)
                                                                                                                                gtcatGCAAGCCTTCCAGCGACAGC‑CGCTGTGGGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGAc                                                           >  5:93235/5‑110 (MQ=21)
                                                                                                                                aCCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTgg                     <  6:59643/148‑1 (MQ=255)
                                                                                                                                 tcatGCAAGCCTTCCAGCGACAGC‑CGCTGTGGGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGt                   >  8:353268/4‑149 (MQ=37)
                                                                                                                                 ccGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGt                   >  2:58020/1‑149 (MQ=255)
                                                                                                                                   gTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGa                 <  2:395418/149‑1 (MQ=255)
                                                                                                                                   gTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGa                 <  5:320611/149‑1 (MQ=255)
                                                                                                                                   gTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGa                 <  2:128323/149‑1 (MQ=255)
                                                                                                                                    tGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGCGATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGt                   >  1:271607/1‑147 (MQ=255)
                                                                                                                                       aaGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGc              >  2:99763/1‑148 (MQ=255)
                                                                                                                                         gCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAAc           <  6:186157/149‑1 (MQ=255)
                                                                                                                                           cccACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCg         >  7:2175/1‑149 (MQ=255)
                                                                                                                                            ccACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCtt                                                                             >  3:360745/1‑81 (MQ=255)
                                                                                                                                            ccACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCtt                                                                             <  4:360745/81‑1 (MQ=255)
                                                                                                                                            ccACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCCGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGc        >  8:84352/1‑149 (MQ=255)
                                                                                                                                            ccACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGc        >  4:209527/1‑149 (MQ=255)
                                                                                                                                              acaGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTg      <  3:209527/149‑1 (MQ=255)
                                                                                                                                               caGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTgg     <  2:110023/149‑1 (MQ=255)
                                                                                                                                                  cGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGgtgt  <  3:269118/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GGTCACCCGTGAGCCAATGGCGCCTGTTGCTGACCCGGCCTGCCGAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGT  >  NC_002947/238912‑239206

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