Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 239,052 C→T 5.2% P45S (CCA→TCA) ‡ hemD → uroporphyrinogen‑III synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947239,0520CT 5.2% 300.5 / 3.6 97P45S (CCA→TCA) ‡hemDuroporphyrinogen‑III synthase
Reads supporting (aligned to +/- strand):  ref base C (48/43);  new base T (2/3);  total (51/46)
Fisher's exact test for biased strand distribution p-value = 6.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.17e-01

Only 100 of 103 total aligned reads displayed.

GCGAGGCCGGTCACCCGTGAGCCAATGGCGCCTGTTGCTGACCCGGCCTGCCGAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGC  >  NC_002947/238904‑239200
                                                                                                                                                       |                                                                                                                                                     
gCGAGGCCGGTCACCCGTGAGCCAATGGCGCCTGTTGCTGACCCGGCCTGCCGAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGccc                                                                                                                                                       <  3:185083/149‑1 (MQ=255)
  gAGGCCGGTCACCCGTGAGCCAATGGCGCCTGTTGCTGACCCGGCCTGCCGAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCca                                                                                                                                                     >  8:14216/1‑149 (MQ=255)
        ggTCACCCGTGAGCCAATGGCGCCTGTTGCTGACCCGGCCTGCCGAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCgg                                                                                                                                               >  1:206845/1‑149 (MQ=255)
        ggTCACCCGTGAGCCAATGGCGCCTGTTGCTGACCCGGCCTGCCGAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCgg                                                                                                                                               >  2:247100/1‑149 (MQ=255)
          tCACCCGTGAGCCAATGGCGCCTGTTGCTGACCCGGCCTGCCGAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCt                                                                                                                                             <  2:206845/149‑1 (MQ=255)
              ccGTGAGCCAATGGCGCCTGTTGCTGACCCGGCCTGCCGAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCg                                                                                                                                         >  4:195647/1‑149 (MQ=255)
                 tGAGCCAATGGCGCCTGTTGCTGACCCGGCCTGCCGAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTg                                                                                                                                      >  4:84296/1‑149 (MQ=255)
                   aGCCAATGGCGCCTGTTGCTGACCCGGCCTGCCGAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTgg                                                                                                                                     <  7:216804/148‑1 (MQ=255)
                     ccAATGGCGCCTGTTGCTGACCCGGCCTGCCGAGGACTGCGCGGCGCTGGCGCGAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTgg                                                                                                                                     >  1:283912/1‑147 (MQ=255)
                     ccAATGGCGCCTGTTGCTGACCCGGCCTGCCGAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAg                                                                                                                                  >  4:303484/1‑149 (MQ=255)
                        aTGGCGCCTGTTGCTGACCCGGCCTGCCGAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCTTCCAGCGACAGCCGCTgtggggct                                                                                                                               >  8:127964/1‑142 (MQ=255)
                           gcgcCTGTTGCTGACCCGGCCTGCCGAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTg                                                                                                                             >  8:303382/1‑148 (MQ=255)
                           gcgcCTGTTGCTGACCCGGCCTGCCGAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGc                                                                                                                            <  8:344259/149‑1 (MQ=255)
                           gcgcCTGTTGCTGACCCGGCCTGCCGAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGc                                                                                                                            <  4:156233/149‑1 (MQ=255)
                           gcgcCTGTTGCTGACCCGGCCTGCCGAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGc                                                                                                                            <  7:170778/149‑1 (MQ=255)
                            cgcCTGTTGCTGACCCGGCCTGCCGAGGACTGCGCGGCGGTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAATGCCCCACAGCGGCTGTCGCTGGAAGGCTTGc                                                                                                                            <  1:38458/149‑1 (MQ=255)
                            cgcCTGTTGCTGACCCGGCCTGCCGAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGATCCCATTACCGTGCAAGCCTTCCAGCGACAGCCGCTGTGGGGCTTGCa                                                                                                                           <  6:385403/149‑1 (MQ=38)
                            cgcCTGTTGCTGACCCGGCCTGCCGAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCTTCCAGCGACAGCCGCTGTGGGGCTTGCa                                                                                                                           <  7:127964/149‑1 (MQ=38)
                             gcCTGTTGCTGACCCGGCCTGCCGAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCAt                                                                                                                          <  4:356452/149‑1 (MQ=255)
                              ccTGTTGCTGACCCGGCCTGCCGAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCAt                                                                                                                          >  1:168843/1‑148 (MQ=255)
                                  ttGCTGACCCGGCCTGCCGAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCTTCCAGCGACAGCCGCTGTGGGGCTTGCacggt                                                                                                                       <  8:121578/147‑5 (MQ=37)
                                  ttGCTGACCCGGCCTGCCGAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGAc                                                                                                                       <  3:56994/147‑1 (MQ=255)
                                     cTGACCCGGCCTGCCGAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACGGCGGCTGTCGCTGGAAGTCATGCAAGCCTTCCa                                                                                                                  <  4:26010/149‑1 (MQ=255)
                                     cTGACCCGGCCTGCCGAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCa                                                                                                                  <  4:174303/149‑1 (MQ=255)
                                        aCCCGGCCTGCCGAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAgg                                                                                                                <  7:14216/148‑1 (MQ=255)
                                        aCCCGGCCTGCCGAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGc                                                                                                               <  1:56972/149‑1 (MQ=255)
                                          ccGGCCTGCCGAGGACTGCGCGGCGCATGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGc                                                                                                                            >  1:193736/1‑135 (MQ=255)
                                          ccGGCCTGCCGAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGc                                                                                                                            <  2:193736/134‑1 (MQ=255)
                                           cGGCCTGCCGAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGAt                                                                                                            >  5:198025/1‑149 (MQ=255)
                                           cGGCCTGCCGAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGAt                                                                                                            >  1:35914/1‑149 (MQ=255)
                                            ggCCTGCCGAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATc                                                                                                           <  3:303484/149‑1 (MQ=255)
                                               cTGCCGAGGACTGCGCAGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGc                                                                                                               <  5:276530/142‑1 (MQ=255)
                                                  ccGAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATtgtg                                                                                                     <  1:22008/149‑1 (MQ=255)
                                                  ccGAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATtgt                                                                                                      <  5:94750/148‑1 (MQ=255)
                                                  ccGAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATtgt                                                                                                      >  6:94750/1‑148 (MQ=255)
                                                    gAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGt                                                                                                   <  3:84296/149‑1 (MQ=255)
                                                     aGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTc                                                                                                  <  3:248597/149‑1 (MQ=255)
                                                          tGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACt                                                                                                                      >  2:66186/1‑124 (MQ=255)
                                                          tGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACt                                                                                                                      <  1:66186/124‑1 (MQ=255)
                                                          tGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTc                                                                                                  <  7:38146/144‑1 (MQ=255)
                                                          tGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTc                                                                                                  >  8:38146/1‑144 (MQ=255)
                                                             gcggcgCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCGCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCCCAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAgcc                                                                                          >  1:268185/1‑149 (MQ=255)
                                                             gcggcgCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAgcc                                                                                          >  7:182100/1‑149 (MQ=255)
                                                                  gcTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAgccgggcg                                                                                     >  6:193338/1‑149 (MQ=255)
                                                                  gcTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAgccggccg                                                                                     >  4:113867/1‑149 (MQ=255)
                                                                         gcAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGAc                                                                                                                       <  4:222396/109‑1 (MQ=255)
                                                                         gcAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGAc                                                                                                                       >  3:222396/1‑109 (MQ=255)
                                                                            aaaGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCGCTGCTGGAGATTGAACCCATTACCGTGCACGCCCCCCAGCGGCTGTCGCTGGAAGGGTTGCATGACTTCCCGGCGATCATGGTGGTCAGCAAGCCGGCCGCCCGGcgg                                                                             >  8:194336/1‑147 (MQ=255)
                                                                              aGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGtt                                                                           >  8:74646/1‑149 (MQ=255)
                                                                                 cTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCCTTCCCGTGCAAGCCCCCCCGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGGCAGCAAGCCGGGCGGCCGGCttgctga                                                                        <  2:347316/149‑7 (MQ=255)
                                                                                 cTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGGCGGCCGGCttgctga                                                                        <  3:200806/149‑7 (MQ=255)
                                                                                 cTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCt                                                                        >  4:362212/1‑149 (MQ=255)
                                                                                    gcggcAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGa                                                                     <  6:198025/149‑1 (MQ=255)
                                                                                    gcggcAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGa                                                                     <  2:168843/149‑1 (MQ=255)
                                                                                         aGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATc                                                                >  2:4211/1‑149 (MQ=255)
                                                                                          gCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCtt                                                                       <  5:64113/141‑1 (MQ=255)
                                                                                          gCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCtt                                                                       >  6:64113/1‑141 (MQ=255)
                                                                                              gAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGcc                                                           >  2:136408/1‑149 (MQ=255)
                                                                                              gAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGcc                                                           <  7:89204/149‑1 (MQ=255)
                                                                                                  gggCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCgaga                                                       >  3:103565/1‑149 (MQ=255)
                                                                                                    gcagcaGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCt                                                                        <  8:250423/130‑1 (MQ=255)
                                                                                                    gcagcaGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCt                                                                        >  7:250423/1‑130 (MQ=255)
                                                                                                    gcagcaGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGAc                                                      >  2:64697/1‑148 (MQ=255)
                                                                                                     cagcagTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCGGTTGCTTGAGCATCTGGCCGAGAc                                                      >  6:48555/1‑147 (MQ=255)
                                                                                                        cagTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCa                                                >  6:202357/1‑149 (MQ=255)
                                                                                                             gcctgccCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCa                                            <  6:989/148‑1 (MQ=255)
                                                                                                              cctgccCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGGCGGCCGGCttgc                                                                           >  5:11722/1‑114 (MQ=255)
                                                                                                              cctgccCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGGCGGCCGGCttgc                                                                           <  6:11722/117‑4 (MQ=255)
                                                                                                                 gacccTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTc                                       >  1:91203/3‑149 (MQ=255)
                                                                                                                  ccccTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTcc                                      >  7:8233/1‑149 (MQ=255)
                                                                                                                     ctgctgGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCa                                   >  8:294561/1‑149 (MQ=255)
                                                                                                                           gagaTTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGa                               <  5:202357/147‑1 (MQ=255)
                                                                                                                           gagaTTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGa                               <  1:64697/147‑1 (MQ=255)
                                                                                                                           gagaTTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAgg                             <  1:247100/149‑1 (MQ=255)
                                                                                                                           gagaTTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAgg                             >  2:184194/1‑149 (MQ=255)
                                                                                                                            agaTTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGc                            <  8:182100/149‑1 (MQ=255)
                                                                                                                             gaTTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCt                           >  4:206121/1‑149 (MQ=255)
                                                                                                                               ttGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTgg                         <  7:303382/149‑1 (MQ=255)
                                                                                                                                tGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGt                        >  2:29479/1‑149 (MQ=255)
                                                                                                                                 gAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGtt                       >  6:61674/1‑149 (MQ=255)
                                                                                                                                  aaCCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCCATCCAGCCTCCGCAGCGAGGCTGGtt                       <  1:4211/149‑1 (MQ=255)
                                                                                                                                  aaCCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGtt                       >  6:277388/1‑148 (MQ=255)
                                                                                                                                    cccATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCAc                    <  2:283912/149‑1 (MQ=255)
                                                                                                                                     ccATTACCGTGCAAGCCCCACAGCGGCTGTCGCTAGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCAcc                   >  1:49372/1‑149 (MQ=255)
                                                                                                                                      cATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCg                  <  3:113867/149‑1 (MQ=255)
                                                                                                                                      cATTACCATGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAgccgg                                                                                        >  1:374131/1‑80 (MQ=255)
                                                                                                                                      cATTACCATGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAgccgg                                                                                        <  2:374131/80‑1 (MQ=255)
                                                                                                                                       aTTACCGTGCAAGCCACACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGACGGCCGACCGGCTGTTGCTTGAGAATCTCGCAGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGt                 >  8:353903/1‑149 (MQ=255)
                                                                                                                                          aCCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTgg               <  7:74646/148‑1 (MQ=255)
                                                                                                                                          aCCGTGCAAGCCCCACAGCGGCTGTAGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGTTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGAAGCGAGGCTGGTTCACCGTggg              >  7:282643/1‑149 (MQ=255)
                                                                                                                                             gTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGa           <  1:136408/149‑1 (MQ=255)
                                                                                                                                              tGCAAGCCTTCCAGCGACAGCCGCTGTGGGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAg          >  1:220823/1‑149 (MQ=38)
                                                                                                                                               gCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAgg         >  3:248350/1‑149 (MQ=255)
                                                                                                                                                cAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAgg         >  7:265571/1‑148 (MQ=255)
                                                                                                                                                cAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAgg         >  4:161915/1‑148 (MQ=255)
                                                                                                                                                     ccccACAGCGGCTGTCGCTGGCAGGCTTGCAAGACTGCCAGGCGATCATTCAGGTCAGCAAGACGGCCGCCCGGCTGTTGCTAGAGCATCTGGCCGAGCCCAGCCTCCGGCCGCCCCAGCGAGGCGCGTACCCCGGGGGAGAGGCAccc    >  6:213271/1‑149 (MQ=255)
                                                                                                                                                      cccACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCg   >  5:297915/1‑149 (MQ=255)
                                                                                                                                                       ccACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGc  >  3:276782/1‑149 (MQ=255)
                                                                                                                                                       ccACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGc  >  5:317736/1‑149 (MQ=255)
                                                                                                                                                       ccACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGc  <  8:8233/149‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                     
GCGAGGCCGGTCACCCGTGAGCCAATGGCGCCTGTTGCTGACCCGGCCTGCCGAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGC  >  NC_002947/238904‑239200

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: