Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 239,163 Δ1 bp 6.0% coding (244/768 nt) hemD → uroporphyrinogen‑III synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947239,1630C. 6.0% 609.1 / 33.1 183coding (244/768 nt)hemDuroporphyrinogen‑III synthase
Reads supporting (aligned to +/- strand):  ref base C (89/83);  new base . (5/6);  total (94/89)
Fisher's exact test for biased strand distribution p-value = 7.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

Only 100 of 185 total aligned reads displayed.

TGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGTTTGCTTTCCGCCGCGTGGCGATGACAGCGAAGCCTTGCTGGCACTGCCGGCCCTGCGCCAGGCACTTACCGAGCCGGCAC  >  NC_002947/239014‑239310
                                                                                                                                                     |                                                                                                                                                    
 gccCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTcc                                                                                                                                                     <  2:2483724/150‑1 (MQ=255)
    ccTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCa                                                                                                                                                  <  2:2351754/150‑1 (MQ=255)
     ctgctgGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAg                                                                                                                                                 <  1:878319/150‑1 (MQ=255)
      tgctgGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGc                                                                                                                                                >  2:4703932/1‑150 (MQ=255)
        ctgGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGa                                                                                                                                             <  2:230977/150‑1 (MQ=255)
           gagaTTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGc                                                                                                                                          >  1:4881979/1‑150 (MQ=255)
              aTTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTgg                                                                                                                                       <  2:2580334/150‑1 (MQ=255)
              aTTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTgg                                                                                                                                       <  1:1318627/150‑1 (MQ=255)
               ttGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCGCTGCGGAGGCTgga                                                                                                                                      <  1:3096106/150‑2 (MQ=255)
               ttGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCGCTGCGGAGGCTgga                                                                                                                                      <  1:218142/150‑2 (MQ=255)
               ttGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCAGCCAGCCTCCGCCGCGCGGCTGGt                                                                                                                                      >  2:3619984/1‑150 (MQ=255)
                   aCCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTCGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCAc                                                                                                                                  >  1:1748653/1‑150 (MQ=255)
                     ccATTACCGTGCAAGCCCCATAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCg                                                                                                                                >  1:4929358/1‑150 (MQ=255)
                       aTTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTg                                                                                                                              >  2:4990758/1‑150 (MQ=255)
                       aTTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTg                                                                                                                              >  2:2220303/1‑150 (MQ=255)
                       aTTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTg                                                                                                                              >  2:4052813/1‑150 (MQ=255)
                        ttACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTgg                                                                                                                             >  2:527897/1‑150 (MQ=255)
                        ttACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTgg                                                                                                                             <  1:3858802/150‑1 (MQ=255)
                        tgaCCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTgg                                                                                                                             <  2:4444732/148‑1 (MQ=255)
                          aCCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGt                                                                                                                           <  1:925950/150‑1 (MQ=255)
                          aCCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGt                                                                                                                           <  1:933911/150‑1 (MQ=255)
                            cgagCAAGCCCCACAGCTGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTTTGGACAGCAACCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCAGAGACCACCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGa                                                                                                                         <  2:1593299/147‑1 (MQ=255)
                            cGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGa                                                                                                                         >  1:617982/1‑150 (MQ=255)
                            cGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGa                                                                                                                         >  2:4618725/1‑150 (MQ=255)
                              tGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAgg                                                                                                                       <  2:985608/150‑1 (MQ=255)
                               gCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCCGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGc                                                                                                                      >  1:1459287/1‑150 (MQ=255)
                                    ccccACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCg                                                                                                                 >  2:2727989/1‑150 (MQ=255)
                                       cacaGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTg                                                                                                              <  1:3524656/150‑1 (MQ=255)
                                       cacaGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGAACATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTg                                                                                                              <  1:2742569/150‑1 (MQ=255)
                                        acaGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTgg                                                                                                             >  1:4860171/1‑150 (MQ=255)
                                         caGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCACGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGgt                                                                                                            <  2:3758160/150‑1 (MQ=255)
                                           gCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGgtgt                                                                                                          >  1:2200882/1‑150 (MQ=255)
                                             ggCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGtgc                                                                                                        <  2:3595853/150‑1 (MQ=255)
                                                tGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGtgctgc                                                                                                     >  2:2066351/1‑150 (MQ=255)
                                                  tCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAg                                                                                                   <  1:1971658/150‑1 (MQ=255)
                                                  tCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAg                                                                                                   <  2:2628635/150‑1 (MQ=255)
                                                  tCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCGCTGCGGAGGCTgga                                                                                                                                      <  1:2485371/115‑2 (MQ=255)
                                                   cGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAgg                                                                                                  <  2:219760/150‑1 (MQ=255)
                                                   cGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAgg                                                                                                  <  1:1701926/150‑1 (MQ=255)
                                                    gCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAggg                                                                                                 >  2:538926/1‑150 (MQ=255)
                                                     cTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGt                                                                                                <  1:298475/150‑1 (MQ=255)
                                                           ggCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTggg                                                                                          <  1:727462/150‑1 (MQ=255)
                                                           ggCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTggg                                                                                          <  2:3914746/150‑1 (MQ=255)
                                                              ttGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTgg                                                                                                                             <  2:3666225/112‑1 (MQ=255)
                                                              ttGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTg                                                                                       >  1:3451578/1‑150 (MQ=255)
                                                                gCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGaa                                                                                     <  2:2573219/150‑1 (MQ=255)
                                                                gCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGaa                                                                                     >  2:1393895/1‑150 (MQ=255)
                                                                gCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGaa                                                                                     <  2:1652841/150‑1 (MQ=255)
                                                                 cATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAg                                                                                    <  1:3379701/150‑1 (MQ=255)
                                                                 cATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAg                                                                                    >  2:64974/1‑150 (MQ=255)
                                                                    gACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGtt                                                                                 <  2:4594957/150‑1 (MQ=255)
                                                                    gACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGtt                                                                                 <  1:2331134/150‑1 (MQ=255)
                                                                     aCTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGttt                                                                                <  1:186338/150‑1 (MQ=255)
                                                                     aCTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGttt                                                                                <  2:1916939/150‑1 (MQ=255)
                                                                     aCTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGttt                                                                                <  1:4376526/150‑1 (MQ=255)
                                                                     aCTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCGCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGttt                                                                                <  1:1090863/150‑1 (MQ=255)
                                                                      attCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGTTTg                                                                               <  1:2220303/149‑1 (MQ=255)
                                                                       ttCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGTTTGc                                                                              >  1:960324/1‑150 (MQ=255)
                                                                       ttCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGTTTGc                                                                              <  2:1286658/150‑1 (MQ=255)
                                                                          cAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGTTTGCttt                                                                           <  2:1448840/150‑1 (MQ=255)
                                                                          cAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGTTTGCttt                                                                           <  2:2019531/150‑1 (MQ=255)
                                                                           aGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGTTTGCTTTc                                                                          >  2:4194836/1‑150 (MQ=255)
                                                                            ggCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGTTTGCTTTcc                                                                         <  1:2063620/150‑1 (MQ=255)
                                                                            ggCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGTTTGCTTTcc                                                                         <  1:4431465/150‑1 (MQ=255)
                                                                             gCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGTTTGCTTTccg                                                                        >  1:4123531/1‑150 (MQ=255)
                                                                                   aTTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGTTTGCTTTCCGCCGCGt                                                                  <  2:1149035/150‑1 (MQ=255)
                                                                                   aTTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGTTTGCTTTCCGCCGCGt                                                                  >  1:3145821/1‑150 (MQ=255)
                                                                                    ttGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGTTTGCTTTCCGCCGCGTg                                                                 <  2:1805893/150‑1 (MQ=255)
                                                                                      gtgGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCTGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGTTTGCTTTCCGCCGCGTGGc                                                               <  2:4857472/150‑1 (MQ=255)
                                                                                         gTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGTTTGCTTTCCGCCGCGTGGCGAt                                                            <  2:879731/150‑1 (MQ=255)
                                                                                           cAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGTTTGCTTTCCGCCGCGTGGCGATGa                                                          >  2:2019433/1‑150 (MQ=255)
                                                                                                 gccggccgCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGTTTGCTTTCCGCCGCGTGGCGATGACAGCGa                                                    >  2:2919071/1‑150 (MQ=255)
                                                                                                  ccggccgCCCGGCTGTTGCTTGAGCCTCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGTTTGCTTTCCGCCGCGTGGCGATGACAGCGaa                                                   <  1:998107/150‑1 (MQ=255)
                                                                                                   cggccgCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGTTTGCTTTCCGCCGCGTGGCGATGACAGCGAAg                                                  >  2:3396052/1‑150 (MQ=255)
                                                                                                       cgccCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGTTTGCTTTCCGCCGCGTGGCGATGACAGCGAAGCCtt                                              <  1:2210311/150‑1 (MQ=255)
                                                                                                         cccGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGTTTGCTTTCCGCCGCGTGGCGATGACAGCGAAGCCTTGc                                            >  1:3663246/1‑150 (MQ=255)
                                                                                                         cccGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGTTTGCTTTCCGCCGCGTGGCGATGACAGCGAAGCCTTGc                                            <  2:3236472/150‑1 (MQ=255)
                                                                                                          ccGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGTTTGCTTTCCGCCGCGTGGCGATGACAGCGAAGCCTTGCt                                           <  2:612757/150‑1 (MQ=255)
                                                                                                              cTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGTTTGCTTTCCGCCGCGTGGCGATGACAGCGAAGCCTTGCTGGCa                                       <  1:3294148/150‑1 (MQ=255)
                                                                                                                    cTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGTTTGCTTTCCGCCGCGTGGCGATGACAGCGAAGCCTTGCTGGCACTGCCg                                 >  1:1210283/1‑150 (MQ=255)
                                                                                                                     ttGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTg                                                                                                              <  2:5035957/72‑1 (MQ=255)
                                                                                                                     ttGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGTTTGCTTTCCGCCGCGTGGCGATGACAGCGAAGCCTTGCTGGCACTGCCgg                                >  1:2962399/1‑150 (MQ=255)
                                                                                                                        aGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGTTTGCTTTCCGCCGCGTGGCGATGACAGCGAAGCCTTGCTGGCACTGCCGGccc                             <  2:642752/150‑1 (MQ=255)
                                                                                                                        aGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGTTTGCTTTCCGCCGCGTGGCGATGACAGCGAAGCCTTGCTGGCACTGCCGGccc                             <  1:4703932/150‑1 (MQ=255)
                                                                                                                         gCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGTTTGCTTTCCGCCGCGTGGCGATGACAGCGAAGCCTTGCCGGCAGTGCCAGcaag                            <  1:527897/150‑4 (MQ=255)
                                                                                                                              tGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGtgct                                                                                                       <  1:3730607/70‑1 (MQ=255)
                                                                                                                              tGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGTTTGCTTTCCGCCGCGTGGCGATGACAGCGAAGCCTTGCTGGCACTGCCGGCCCTGCGcc                       >  1:3641818/1‑150 (MQ=255)
                                                                                                                               ggCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGTTTGCTTTCCGCCGCGTGGCGATGACAGCGAAGCCTTGCTGGCACTGCCGGCCCTGCGCCa                      <  1:2233333/150‑1 (MQ=255)
                                                                                                                                   gagaCCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGTTTGCTTTCCGCCGCGTGGCGATGACAGCGAAGCCTTGCTGGCACTGCCGGCCCTGCGCCAGGCa                  >  1:458350/1‑150 (MQ=255)
                                                                                                                                      aCCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGTTTGCTTTCCGCCGCGTGGCGATGACAGCGAAGCCTTGCTGGCACTGCCGGCCCTGCGCCAGGCACtt               >  2:1213228/1‑150 (MQ=255)
                                                                                                                                        cAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGTTTGCTTTCCGCCGCGTGGCGATGACAGCGAAGCCTTGCTGGCACTGCCGGCCCTGCGCCAGGCACTTAc             >  1:1646008/1‑150 (MQ=255)
                                                                                                                                        cAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGTTTGCTTTCCGCCGCGTGGCGATGACAGCGAAGCCTTGCTGGCACTGCCGGCCCTGCGCCAGGCACTTAc             <  2:4728989/150‑1 (MQ=255)
                                                                                                                                        cAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGTTTGCTTTCCGCCGCGTGGCGATGACAGCGAAGCCTTGCTGGCACTGCCGGCCCTGCGCCAGGCACTTAc             >  1:2365629/1‑150 (MQ=255)
                                                                                                                                         aGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGTTTGCTTTCCGCCGCGTGGCGATGACAGCGAAGCCTTGCTGGCACTGCCGGCCCTGCGCCAGGCACTTAcc            >  2:5106982/1‑150 (MQ=255)
                                                                                                                                           cATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGTTTGCTTTCCGCCGCGTGGCGATGACAGCGAAGCCTTGCTGGCACTGCCGGCCCTGCGCCAGGCACTTACCGa          <  2:1748653/150‑1 (MQ=255)
                                                                                                                                            aTCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGTTTGCTTTCCGCCGCGTGGCGATGACAGCGAAGCCTTGCTGGCACTGCCGGCCCTGCGCCAGGCACTTACCGAg         >  1:3995913/1‑150 (MQ=255)
                                                                                                                                              ccAGCCTCGCTGCGGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGTTTGCTTTCCGCCGCGTGGCGATGACAGCGAAGCCTTGCTGGCACTGCCGGCCCTGCGCCAGGCACTTACCGAGcc       >  1:3473276/1‑150 (MQ=255)
                                                                                                                                              ccAGCCTCGCTGCGGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGTTTGCTTTCCGCCGCGTGGCGATGACAGCGAAGCCTTGCTGGCACTGCCGGCCCTGCGCCAGGCACTTACCGAGcc       >  2:2856332/1‑150 (MQ=255)
                                                                                                                                               cAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGTTTGCTTTCCGCCGCGTGGCGATGACAGCGAAGCCTTGCTGGCACTGCCGGCCCTGCGCCAGGCACTTACCGAGCCg      >  2:2974960/1‑150 (MQ=255)
                                                                                                                                                   cTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGTTTGCTTTCCGCCGCGTGGCGATGACAGCGAAGCCTTGCTGGCACTGCCGGCCCTGCGCCAGGCACTTACCGAGCCGGCAc  >  1:1006196/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGTTTGCTTTCCGCCGCGTGGCGATGACAGCGAAGCCTTGCTGGCACTGCCGGCCCTGCGCCAGGCACTTACCGAGCCGGCAC  >  NC_002947/239014‑239310

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: