Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 239,062 T→A 5.3% L48Q (CTG→CAG)  hemD → uroporphyrinogen‑III synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947239,0620TA 5.3% 303.3 / 4.7 95L48Q (CTG→CAG) hemDuroporphyrinogen‑III synthase
Reads supporting (aligned to +/- strand):  ref base T (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 = 9.28e-01

ACCCGTGAGCCAATGGCGCCTGTTGCTGACCCGGCCTGCCGAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGT  >  NC_002947/238916‑239206
                                                                                                                                                  |                                                                                                                                                   
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                                                                                    >  5:186157/1‑149 (MQ=255)
                                                             cAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCt                                                                                    <  7:32157/149‑1 (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                                                                                <  1:278536/129‑1 (MQ=255)
                                                                                     ggCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTg                                                                                >  2:278536/1‑129 (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                                               >  1:395418/1‑149 (MQ=255)
                                                                                                 ctgccCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGACCAGCATCCAGcc                                               <  3:145776/149‑1 (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                                                                                                                      <  1:109340/68‑1 (MQ=255)
                                                                                                            ggAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAgg                                                                                                                      >  2:109340/1‑68 (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:185372/149‑1 (MQ=255)
                                                                                                                     aCCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCa                           >  7:339987/1‑149 (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                                                           >  5:93235/5‑110 (MQ=21)
                                                                                                                            gtcatGCAAGCCTTCCAGCGACAGC‑CGCTGTGGGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGAc                                                           <  6:93235/106‑1 (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                 <  5:320611/149‑1 (MQ=255)
                                                                                                                               gTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGa                 <  2:395418/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                                                                             <  4:360745/81‑1 (MQ=255)
                                                                                                                                        ccACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCtt                                                                             >  3:360745/1‑81 (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)
                                                                                                                                                  |                                                                                                                                                   
ACCCGTGAGCCAATGGCGCCTGTTGCTGACCCGGCCTGCCGAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGT  >  NC_002947/238916‑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: ‑