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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947239,0530CT 5.2% 300.5 / 3.6 96P45L (CCA→CTA) ‡hemDuroporphyrinogen‑III synthase
Reads supporting (aligned to +/- strand):  ref base C (48/43);  new base T (2/3);  total (50/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 = 4.34e-01

Only 100 of 102 total aligned reads displayed.

GAGGCCGGTCACCCGTGAGCCAATGGCGCCTGTTGCTGACCCGGCCTGCCGAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGC  >  NC_002947/238906‑239200
                                                                                                                                                      |                                                                                                                                                    
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)
                      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)
                             cTGTTGCTGACCCGGCCTGCCGAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGa                                                                                                                        <  3:195647/149‑1 (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)
                                             cTGCCGAGGACTGCGCAGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGc                                                                                                               >  6:276530/1‑142 (MQ=255)
                                                ccGAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATtgtg                                                                                                     <  1:22008/149‑1 (MQ=255)
                                                ccGAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATtgt                                                                                                      >  6:94750/1‑148 (MQ=255)
                                                ccGAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATtgt                                                                                                      <  5:94750/148‑1 (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                                                                                                                       >  3:222396/1‑109 (MQ=255)
                                                                       gcAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGAc                                                                                                                       <  4:222396/109‑1 (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                                                                     <  2:168843/149‑1 (MQ=255)
                                                                                  gcggcAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGa                                                                     <  6:198025/149‑1 (MQ=255)
                                                                                       aGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATc                                                                >  2:4211/1‑149 (MQ=255)
                                                                                        gCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCtt                                                                       >  6:64113/1‑141 (MQ=255)
                                                                                        gCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCtt                                                                       <  5:64113/141‑1 (MQ=255)
                                                                                            gAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGcc                                                           >  2:136408/1‑149 (MQ=255)
                                                                                            gAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGcc                                                           <  7:89204/149‑1 (MQ=255)
                                                                                                gggCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCgaga                                                       >  3:103565/1‑149 (MQ=255)
                                                                                                  gcagcaGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCt                                                                        >  7:250423/1‑130 (MQ=255)
                                                                                                  gcagcaGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCt                                                                        <  8:250423/130‑1 (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                                                                           <  6:11722/117‑4 (MQ=255)
                                                                                                            cctgccCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGGCGGCCGGCttgc                                                                           >  5:11722/1‑114 (MQ=255)
                                                                                                               gacccTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTc                                       >  1:91203/3‑149 (MQ=255)
                                                                                                                ccccTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTcc                                      >  7:8233/1‑149 (MQ=255)
                                                                                                                   ctgctgGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCa                                   >  8:294561/1‑149 (MQ=255)
                                                                                                                         gagaTTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGa                               <  1:64697/147‑1 (MQ=255)
                                                                                                                         gagaTTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGa                               <  5:202357/147‑1 (MQ=255)
                                                                                                                         gagaTTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAgg                             >  2:184194/1‑149 (MQ=255)
                                                                                                                         gagaTTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAgg                             <  1:247100/149‑1 (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                                                                                        <  2:374131/80‑1 (MQ=255)
                                                                                                                                    cATTACCATGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAgccgg                                                                                        >  1:374131/1‑80 (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  <  8:8233/149‑1 (MQ=255)
                                                                                                                                                     ccACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGc  >  3:276782/1‑149 (MQ=255)
                                                                                                                                                     ccACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGc  >  5:317736/1‑149 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
GAGGCCGGTCACCCGTGAGCCAATGGCGCCTGTTGCTGACCCGGCCTGCCGAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGC  >  NC_002947/238906‑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: