Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 239,163 Δ1 bp 6.2% 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.2% 259.8 / 9.9 80coding (244/768 nt)hemDuroporphyrinogen‑III synthase
Reads supporting (aligned to +/- strand):  ref base C (38/37);  new base . (2/3);  total (40/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.53e-01

GCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGTTTGCTTTCCGCCGCGTGGCGATGACAGCGAAGCCTTGCTGGCACTGCCGGCCCTGCGCCAGGCACTTACCGAGCCGGCAC  >  NC_002947/239015‑239310
                                                                                                                                                    |                                                                                                                                                     
gccCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTc                                                                                                                                                       >  6:172643/1‑149 (MQ=255)
 ccccTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAACCTcc                                                                                                                                                      <  5:240873/149‑1 (MQ=255)
    ctgctgGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCa                                                                                                                                                   <  7:14806/149‑1 (MQ=255)
       ctgGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGc                                                                                                                                                 >  2:328297/1‑148 (MQ=255)
           agaTTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGc                                                                                                                                           >  2:330140/1‑149 (MQ=255)
                 aaCCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGtt                                                                                                                                      >  8:195074/1‑148 (MQ=255)
                 aaCCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTc                                                                                                                                     <  5:211080/149‑1 (MQ=255)
                 aaCCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCGCTGCGGAGGCTggat                                                                                                                                      <  4:270477/148‑3 (MQ=255)
                    ccATTACCGTGCAAGCCCCACACAGGCTGTCGCTTGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGCCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCAcc                                                                                                                                  <  5:321237/149‑1 (MQ=255)
                          ccGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTgg                                                                                                                                        >  5:29474/1‑137 (MQ=255)
                          ccGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTgg                                                                                                                                        <  6:29474/137‑1 (MQ=255)
                          ccGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTggg                                                                                                                             <  7:195074/148‑1 (MQ=255)
                          ccGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGt                                                                                                                            <  7:140474/149‑1 (MQ=255)
                           cGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTg                                                                                                                           <  6:408786/149‑1 (MQ=255)
                            gTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGa                                                                                                                          <  6:245219/149‑1 (MQ=255)
                               cAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAgg                                                                                                                        <  6:107447/148‑1 (MQ=255)
                                aaGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGc                                                                                                                       >  4:350131/1‑148 (MQ=255)
                                       acaGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTg                                                                                                               >  2:94356/1‑149 (MQ=255)
                                        caGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTgg                                                                                                              <  5:172643/149‑1 (MQ=255)
                                         aGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGgt                                                                                                             >  8:159184/1‑149 (MQ=255)
                                                gTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGtgctg                                                                                                       >  2:408037/1‑148 (MQ=255)
                                                   gCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGtgct                                                                                                        <  3:15176/144‑1 (MQ=255)
                                                   gCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGtgct                                                                                                        >  4:15176/1‑144 (MQ=255)
                                                    cTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGc                                                                                                                 <  4:296057/134‑1 (MQ=255)
                                                    cTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGc                                                                                                                 >  3:296057/1‑134 (MQ=255)
                                                     tGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGt                                                                                                 >  2:56852/1‑149 (MQ=255)
                                                         aGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGTCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTg                                                                                             <  3:149342/149‑1 (MQ=255)
                                                          ggCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTgg                                                                                            <  3:350131/149‑1 (MQ=255)
                                                            cTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGc                                                                                          <  7:245697/149‑1 (MQ=255)
                                                             ttGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCt                                                                                         >  3:334219/1‑149 (MQ=255)
                                                                cATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAAc                                                                                                                    >  7:233259/1‑119 (MQ=255)
                                                                cATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAAc                                                                                                                    <  8:233259/119‑1 (MQ=255)
                                                                 aTGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAgg                                                                                                   <  6:194303/135‑1 (MQ=255)
                                                                 aTGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAgg                                                                                                   >  5:194303/1‑135 (MQ=255)
                                                                 aTGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAg                                                                                    >  8:265979/1‑149 (MQ=255)
                                                                   gACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGt                                                                                  >  4:158487/1‑149 (MQ=255)
                                                                    aCTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGtt                                                                                 <  1:330140/149‑1 (MQ=255)
                                                                     cTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGttt                                                                                <  1:408037/149‑1 (MQ=255)
                                                                      ttCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCGCTGCGGAGGCTgga                                                                                                                                       <  3:328791/94‑2 (MQ=255)
                                                                      ttCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCGCTGCGGAGGCTgga                                                                                                                                       >  4:328791/1‑93 (MQ=255)
                                                                        ccAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGTTTg                                                                               <  4:267898/147‑1 (MQ=255)
                                                                        ccAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGTTTg                                                                               >  3:267898/1‑147 (MQ=255)
                                                                            gCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGTTTGCTTTcc                                                                         >  1:384317/1‑149 (MQ=255)
                                                                                   ttGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGTTTGCTTTCCGCCGCGt                                                                  >  6:306740/1‑149 (MQ=255)
                                                                                   ttGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGGGTGCTGCAGGGTGCTGGGCTGAAGGTTTGCTTTCCGCCGCGt                                                                  >  7:91162/1‑149 (MQ=255)
                                                                                      tgGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGc                                                                                                                                           <  6:55967/74‑1 (MQ=255)
                                                                                      tgGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGc                                                                                                                                           >  5:55967/1‑74 (MQ=255)
                                                                                       ggTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGTTTGCTTTc                                                                          <  7:302898/137‑1 (MQ=255)
                                                                                       ggTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGTTTGCTTTc                                                                          >  8:302898/1‑137 (MQ=255)
                                                                                         tCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGTTTGCTTTCCGCCGCGTGGCGAt                                                            <  2:384317/149‑1 (MQ=255)
                                                                                           aGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGTTTGCTTTCCGCCGCGTGGCGATGa                                                          >  8:52400/1‑149 (MQ=255)
                                                                                           aGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGTTTGCTTTCCGCCGCGTGGCGATGa                                                          <  3:158487/149‑1 (MQ=255)
                                                                                            gCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAgg                                                                                                   <  8:161491/108‑1 (MQ=255)
                                                                                            gCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAgg                                                                                                   >  7:161491/1‑108 (MQ=255)
                                                                                             cAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGTTTGCTTTCCGCCGCGTGGCGATGa                                                          >  2:21306/1‑147 (MQ=255)
                                                                                             cAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGTTTGCTTTCCGCCGCGTGGCGATGACa                                                        <  1:94356/149‑1 (MQ=255)
                                                                                              aaGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGTTTGCTTTccgcc                                                                      <  1:123037/134‑1 (MQ=255)
                                                                                              aaGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGTTTGCTTTccgcc                                                                      >  2:123037/1‑134 (MQ=255)
                                                                                                   ggccgCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGtgtg                                                                                                          >  1:345835/1‑94 (MQ=255)
                                                                                                   ggccgCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGtgtg                                                                                                          <  2:345835/94‑1 (MQ=255)
                                                                                                        cccGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGTTTGCTTTCCGCCGCGTGGCGATGACAGCGAAGCCTTg                                             >  7:222602/1‑149 (MQ=255)
                                                                                                          cGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGTTTGCTTTCCGCCGCGTGGCGATGACAGCGAAGCCTTGCt                                           <  7:265979/149‑1 (MQ=255)
                                                                                                           ggCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGTTTGCTTTCCGCCGCGTGGCGATGACAGCGAAGCCTTGCTg                                          <  7:52400/149‑1 (MQ=255)
                                                                                                           ggCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGTTTGCTTTCCGCCGCGTGGCGATGACAGCGAAGCCTTGCTg                                          <  8:222602/149‑1 (MQ=255)
                                                                                                                ttGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGTTTGCTTTCCGCCGCGTGGCGATGACAGCGAAGCCTTGCTGGCACt                                     >  3:352659/1‑149 (MQ=255)
                                                                                                                     tGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGTTTGCTTTCCGCCGCGTGGCGATGACAGCGAAGCCTTGCTGGCACTGCCgg                                >  3:394893/1‑149 (MQ=255)
                                                                                                                              ggCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGTTTGCTTTCCGCCGCGTGGCGATGACAGCGAAGCCTTGCTGGCACTGCCGGCCCTgcgc                        >  2:362389/1‑148 (MQ=255)
                                                                                                                                 cGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGTTTGCTTTCCGCCGCGTGGCGATGACAGCGAAGCCTTGCTGGCACTGCCGGCCCTGCGCCAgg                    >  5:336272/1‑149 (MQ=255)
                                                                                                                                 cGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGTTTGCTTTCCGCCGCGTGGCGATGACAGCGAAGCCTTGCTGGCACTGCCGGCCCTGCGCCAgg                    >  5:266824/1‑149 (MQ=255)
                                                                                                                                  gagaCCCGCATCCCGCCTCCCCAGCCAGGCGGGTTCACCGGGGGTTAGGCAACCGCTGGGGTGCTGCAGGGGGCTGGGCTTAATCGTTTGCTTTCCGCCGCGTGGCGATGACAGCGAAGCCTTGCTGGCACTGCCGGCCCTGCGCCAgg                    <  1:328297/149‑1 (MQ=255)
                                                                                                                                     aCCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGTTTGCTTTCCGCCGCGTGGCGATGACAGCGAAGCCTTGCTGGCACTGCCGGCCCTGCGCCAGGCACt                >  5:357099/1‑149 (MQ=255)
                                                                                                                                       cAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGTTTGCTTTCCGCCGCGTGGCGATGACAGCGAAGCCTTGCTGGCACTGCCGGCCCTGCGCCAGGCa                  >  2:334987/1‑145 (MQ=255)
                                                                                                                                       cAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGTTTGCTTTCCGCCGCGTGGCGATGACAGCGAAGCCTTGCTGGCACTGCCGGCCCTGCGCCAGGCa                  <  1:334987/145‑1 (MQ=255)
                                                                                                                                        aGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGTTTGCTTTCCGCCGCGTGGCGATGACAGCGAAGCCTTGCTGGCACTGCCGGCCCTGCGCCAGGCACTTAc             <  7:159184/149‑1 (MQ=255)
                                                                                                                                        aGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGTTTGCTTTCCGCCGCGTGGCGATGACAGCGAAGCCTTGCTGGCACTGCCGGCCCTGCGCCAGGCACTTAc             <  4:352659/149‑1 (MQ=255)
                                                                                                                                            tCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGTTTGCTTTCCGCCGCGTGGCGATGACAGCGAAGCCTTGCTGGCACTGCCGGCCCTGCGCCAGGCACTTAcc            >  1:296469/1‑146 (MQ=255)
                                                                                                                                            tCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGTTTGCTTTCCGCCGCGTGGCGATGACAGCGAAGCCTTGCTGGCACTGCCGGCCCTGCGCCAGGCACTTACCGa          <  2:296469/148‑1 (MQ=255)
                                                                                                                                            tCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGGGCTGGGCTGAAGGTTTGCTTTCCGCCGCGTGGCGATGACAGCGAAGCCTTGCTGGCACTGCCGGCCCTGCGCCAGGCACTTACCGAg         <  8:91162/149‑1 (MQ=255)
                                                                                                                                             ccAGCCTCGCTGCGGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAgg                                                                                                   <  6:31551/59‑1 (MQ=25)
                                                                                                                                             ccAGCCTCGCTGCGGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAgg                                                                                                   >  5:31551/1‑59 (MQ=25)
                                                                                                                                                gCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGTTTGCTTTCCGCCGCGTGGCGATGACAGCGAAGCCTTGCTGGCACTGCCGGCCCTGCGCCAGGCACTTACCGAGCCgg     >  5:36284/1‑149 (MQ=255)
                                                                                                                                                   tCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGTTTGCTTTCCGCCGCGTGGCGATGACAGCGAAGCCTTGCTGGCACTGCCGGCCCTGCGCCAGGCACTTACCGAGCCGGCAc  >  7:187966/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTCGCTGGAAGGCTTGCATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTGAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGCTGGTGTGCTGCAGGGTGCTGGGCTGAAGGTTTGCTTTCCGCCGCGTGGCGATGACAGCGAAGCCTTGCTGGCACTGCCGGCCCTGCGCCAGGCACTTACCGAGCCGGCAC  >  NC_002947/239015‑239310

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: