Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 240,531 G→C 5.4% A276P (GCC→CCC)  PP_0188 → uroporphyrin‑III C‑methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947240,5310GC 5.4% 488.1 / 13.6 148A276P (GCC→CCC) PP_0188uroporphyrin‑III C‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (64/75);  new base C (4/4);  total (69/79)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

Only 100 of 155 total aligned reads displayed.

CCCGAGTTCGACAGTAACGCCGATGCGCTTGGCGCCCTGACCGCCGACGGTGATGGTGCCAGCAGGCTGTCGCAGTGGTGGGCAGAAATATCCAAGTACTTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGTGATGCCAAGGTTTACACGCAGGCGCTGGACGACGCCCGCAGCGTACTGCTGGCCAACTTCAATGCAGACAAC  >  NC_002947/240387‑240680
                                                                                                                                                |                                                                                                                                                     
cccGAGTTCGACAGTAACGCCGATGCGCTTGGCGCCCTGACCGCCGACGGTGATGGTGCCAGCAGGCTGTCGCAGTGGTGGGCAGAAATATCCAAGTACTTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCggg                                                                                                                                                  >  2:4790099/1‑150 (MQ=255)
cccGAGTTCGACAGTAACGCCGATGCGCTTGGCGCCCTGACCGCCGACGGTGATGGTGCCAGCAGGCTGTCGCAGTGGTGGGCAGAAATATCCAAGTACTTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCggg                                                                                                                                                  <  2:4540262/150‑1 (MQ=255)
  cGAGTTCGACAGTAACGCCGATGCGCTTGGCGCCCTGACCGCCGACGGTGATGGTGCCAGCAGGCTGTCGCAGTGGTGGGCAGAAATATCCAAGTACTTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGgca                                                                                                                                                >  1:2374730/1‑150 (MQ=255)
      ttCGACAGTAACGCCGATGCGCTTGGCGCCCTGACCGCCGACGGTGATGGTGCCAGCAGGCTGTCGCAGTGGTGGGCAGAAATATCCAAGTACTTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGCCCGGCCAGCaa                                                                                                                                            <  2:4177862/150‑1 (MQ=255)
       tCGACAGTAACGCCGATGCGCTTGGCGCCCTGACCGCCGACGGTGATGGTGCCAGCAGGCTGTCGCAGTGGTGGGCAGAAATATCCAAGTACTTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAAc                                                                                                                                           <  1:1248263/150‑1 (MQ=255)
       tCGACAGTAACGCCGATGCGCTTGGCGCCCTGACCGCCGACGGTGATGGTGCCAGCAGGCTGTCGCAGTGGTGGGCAGAAATATCCAAGTACTTCCAGATCGACTTCAATCCCGATGAAAACGTGCGTCCGTTGCTGCCCGGCCAGCAAc                                                                                                                                           <  2:4070955/150‑1 (MQ=255)
        cGACAGTAACGCCGATGCGCTTGGCGCCCTGACCGCCGACGGTGATGGTGCCAGCAGGCTGTCGCAGTGGTGGGCAGAAATATCCAAGTACTTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACt                                                                                                                                          <  1:4886057/150‑1 (MQ=255)
         gACAGTAACGCCGATGCGCTTGGCGCCCTGACCGCCGACGGTGATGGTGCCAGCAGGCTGTCGCAGTGGTGGGCAGAAATATCCAAGTACTTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTg                                                                                                                                         >  2:625659/1‑150 (MQ=255)
         gACAGTAACGCCGATGCGCTTGGCGCCCTGACCGCCGACGGTGATGGTGCCAGCAGGCTGTCGCAGTGGTGGGCAGAAATATCCAAGTACTTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTg                                                                                                                                         >  2:330523/1‑150 (MQ=255)
           cAGTAACGCCGATGCGCTTGGCGCCCTGACCGCCGACGGTGATGGTGCCAGCAGGCTGTCGCAGTGGTGGGCAGAAATATCCAAGTACTTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCTGTTGCTGCCCGGCCAGCAACGGac                                                                                                                                       <  2:1300140/150‑2 (MQ=255)
           cAGTAACGCCGATGCGCTTGGCGCCCTGACCGCCGACGGTGATGGTGCCAGCAGGCTGTCGCAGTGGTGGGCAGAAATATCCAAGTACTTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGaa                                                                                                                                       <  1:3541255/150‑1 (MQ=255)
           cAGTAACGCCGATGCGCTTGGCGCCCTGACCGCCGACGGTGATGGTGCCAGCAGGCTGTCGCAGTGGTGGGCAGAAATATCCAAGTACTTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGaa                                                                                                                                       >  2:1400351/1‑150 (MQ=255)
            aGTAACGCCGATGCGCTTGGCGCCCTGACCGCCGACGGTGATGGTGCCAGCAGGCTGTCGCAGTGGTGGGCAGAAATATCCAAGTACTTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAAc                                                                                                                                      <  1:4307175/150‑1 (MQ=255)
               aaCGCCGATGCGCTTGGCGCCCTGACCGCCGACGGTGATGGTGCCAGCAGGCTGTCGCAGTGGTGGGCATAAATATCCAAGTACTTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAg                                                                                                                                   <  1:452332/150‑1 (MQ=255)
                aCGCCGATGCGCTTGGCGCCCTGACCGCCGACGGTGATGGTGCCAGCAGGCTGTCGCAGTGGTGGGCAGAAATATCCAAGTACTTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGc                                                                                                                                  >  1:2678761/1‑150 (MQ=255)
                 cGCCGATGCGCTTGGCGCCCTGACCGCCGACGGTGATGGTGCCAGCAGGCTGTCGCAGTGGTGGGCAGAAATATCCAAGTACTTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCt                                                                                                                                 >  2:4148950/1‑150 (MQ=255)
                  gCCGATGCGCTTGGCGCCCTGACCGCCGACGGTGATGGTGCCAGCAGGCTGTCGCAGTGGTGGGCAGAAATATCCAAGTACTTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTg                                                                                                                                >  2:4999984/1‑150 (MQ=255)
                           cTTGGCGCCCTGACCGCCGACGGTGATGGTGCCAGCAGGCTGTCGCAGTGGTGGGCAGAAATATCCAAGTACTTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGcc                                                                                                                       <  2:4743030/150‑1 (MQ=255)
                            ttGGCGCCCTGACCGCCGACGGTGATGGTGCCAGCAGGCTGTCGCAGTGGTGGGCAGAAATATCCAAGTACTTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGccc                                                                                                                      >  1:897708/1‑150 (MQ=255)
                            ttGGCGCCCTGACCGCCGACGGTGATGGTGCCAGCAGGCTGTCGCAGTGGTGGGCAGAAATATCCAAGTACTTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGACCGGGCAGCAACTGAACCAGCTGCGCCTGGccc                                                                                                                      >  2:1650781/1‑150 (MQ=255)
                             tGGCGCCCTGACCGCCGACGGTGATGGTGCCAGCAGTCTGTCGCAGTGGTGGCCAGAAATATCCAAGTAATTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCt                                                                                                                     <  2:3106483/150‑1 (MQ=255)
                               gcgcCCTGACCGCCGACGGTGATGGTGCCAGCAGGCTGTCGCAGTGGTGGGCAGAAATATCCAAGTACTTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGa                                                                                                                   >  1:524048/1‑150 (MQ=255)
                                cgcCCTGACCGCCGACGGTGATGGTGCCAGCAGGCTGTCGCAGTGGTGGGCAGAAATATCCAAGTACTTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCt                                                                                                                     <  2:3530991/147‑1 (MQ=255)
                                cgcCCTGACCGCCGACGGTGATGGTGCCAGCAGGCTGTCGCAGTGGTGGGCAGAAATATCCAAGTACTTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCt                                                                                                                     >  1:3530991/1‑147 (MQ=255)
                                     tGACCGCCGACGGTGATGGTGCCAGCAGGCTGTCGCAGTGGTGGGCAGAAATATCCAAGTACTTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGa                                                                                                             >  2:1222552/1‑150 (MQ=255)
                                      gACCGCCGACGGTGATGGTGCCAGCAGGCTGTCGCAGTGGTGGGCAGAAATATCCAAGTACTTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGAc                                                                                                            >  1:4405157/1‑150 (MQ=255)
                                        ccgccgACGGTGATGGTGCCAGCAGGCTGTCGCAGTGGTGGGCAGAAATATCCAAGTACTTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTa                                                                                                          >  1:151257/1‑150 (MQ=255)
                                            cgACGGTGATGGTGCCAGCAGGCTGTCGCAGTGGTGGGCAGAAATATCCAAGTACTTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGa                                                                                                      <  1:1928013/150‑1 (MQ=255)
                                              aCGGTGATGGTGCCAGCAGGCTGTCGCAGTGGTGGGCAGAAATATCCAAGTACTTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGc                                                                                                    <  2:2161049/150‑1 (MQ=255)
                                               cGGTGTTGGTGCCAGCAGGCTGTCGCAGTGGTGGGCAGAAATATCCAAGTACTTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCa                                                                                                   <  1:736953/150‑1 (MQ=255)
                                               cGGTGATGGTGCCAGCAGGCTGTCGCAGTGGTGGGCAGAAATATCCAAGTACTTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCa                                                                                                   <  2:3373772/150‑1 (MQ=255)
                                                   gATGGTGCCAGCAGGCTGTCGCAGTGGTGGGCAGAAATATCCAAGTACTTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGcc                                                                                               <  1:1713969/150‑1 (MQ=255)
                                                           cagcagGCTGTCGCAGTGGTGGGCAGAAATATCCAAGTACTTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGggc                                                                                       >  2:4465756/1‑150 (MQ=255)
                                                            agcagGCTGTCGCAGTGGTGGGCAGAAATATCCAAGTACTTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGggcg                                                                                      <  1:2784233/150‑1 (MQ=255)
                                                            agcagGCTGTCGCAGTGGTGGGCAGAAATATCCAAGTACTTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGggcg                                                                                      <  2:701318/150‑1 (MQ=255)
                                                              cagGCTGTCGCAGTGGTGGGCAGAAATATCCAAGTACTTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGggcggc                                                                                    >  2:903746/1‑150 (MQ=255)
                                                                 gCTGTCGCAGTGGTGGGCAGAAATATCCAAGTACTTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCt                                                                                 <  1:884513/150‑1 (MQ=255)
                                                                   tGTCGCAGTGGTGGGCAGAAATATCCAAGTACTTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCa                                                                               >  1:2502521/1‑150 (MQ=255)
                                                                    gTCGCAGTGGTGGGCAGAAATATCCAAGTACTTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCATCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCaa                                                                              <  2:4405931/150‑1 (MQ=255)
                                                                      cGCAGTGGTGGGCAGAAATATCCAAGTACTTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATg                                                                            <  2:2261379/150‑1 (MQ=255)
                                                                        cAGTGGTGGGCAGACATATCCAAGTACTTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGt                                                                          <  2:4593495/150‑1 (MQ=255)
                                                                           tggtggGCAGAAATATCCAAGTACTTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGTGAt                                                                       <  1:4465756/150‑1 (MQ=255)
                                                                           tggtggGCAGAAATATCCAAGTACTTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGTGAt                                                                       <  2:601708/150‑1 (MQ=255)
                                                                             gtggGCAGAAATATCCAAGTACTTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGTGATGc                                                                     <  1:1222552/150‑1 (MQ=255)
                                                                                ggCAGAAATATCCAAGTACTTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGTGATGCCaa                                                                  <  1:4159780/150‑1 (MQ=255)
                                                                                ggCAGAAATATCCAAGTACTTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGTGATGCCaa                                                                  <  2:2122305/150‑1 (MQ=255)
                                                                                 gCAGAAATATCCAAGTACTTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGTGATGCCAAg                                                                 >  2:2005454/1‑150 (MQ=255)
                                                                                 gCAGAAATATCCAAGTACTTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGTGATGCCAAg                                                                 <  2:1646387/150‑1 (MQ=255)
                                                                                   aGAAATATCCAAGTACTTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGTGATGCCAAGGt                                                               >  1:4647260/1‑150 (MQ=255)
                                                                                   aGAAATATCCAAGTACTTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGTGATGCCAAGGt                                                               <  2:1893431/150‑1 (MQ=255)
                                                                                       atatCCAAGTACTTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGTGATGCCAAGGTTTac                                                           <  1:3404537/150‑1 (MQ=255)
                                                                                        tatCCAAGTACTTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGTGATGCCAAGGTTTaca                                                          <  1:3156347/150‑1 (MQ=255)
                                                                                        tatCCAAGTACTTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGTGATGCCAAGGTTTaca                                                          <  1:3156325/150‑1 (MQ=255)
                                                                                        tatCCAAGTACTTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGTGATGCCAAGGTTTaca                                                          <  1:481306/150‑1 (MQ=255)
                                                                                        tatCCAAGTACTTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGTGATGCCAAGGTTTaca                                                          <  1:1495941/150‑1 (MQ=255)
                                                                                         atCCAAGTACTTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGTGATGCCAAGGTTTacac                                                         >  2:4672962/1‑150 (MQ=255)
                                                                                            cAAGTACTTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGTGATGCCAAGGTTTACACGCa                                                      >  1:2745945/1‑150 (MQ=255)
                                                                                            cAAGTACTTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGTGATGCCAAGGTTTACACGCa                                                      >  1:2359436/1‑150 (MQ=255)
                                                                                            cAAGTACTTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGTGATGCCAAGGTTTACACGCa                                                      >  1:1432746/1‑150 (MQ=255)
                                                                                            cAAGTACTTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGCGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGTGATGCCAAGGTTTACACGCa                                                      >  1:183029/1‑150 (MQ=255)
                                                                                               gTACTTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGTGATGCCAAGGTTTACACGCAGgc                                                   >  1:4370068/1‑150 (MQ=255)
                                                                                                   ttCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGTGATGCCAAGGTTTACACGCAGGCGCTg                                               >  2:2083148/1‑150 (MQ=255)
                                                                                                   ttCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGTGATGCCAAGGTTTACACGCAGGCGCTg                                               >  1:3698194/1‑150 (MQ=255)
                                                                                                   ttCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGTGATGCCAAGGTTTACACGCAGGCGCTg                                               <  2:4354258/150‑1 (MQ=255)
                                                                                                    tCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGTGATGCCAAGGTTTACACGCAGGCGCTgg                                              >  2:652291/1‑150 (MQ=255)
                                                                                                     ccAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGTGATGCCAAGGTTTACACGCAGGCGCTGga                                             >  2:3569854/1‑150 (MQ=255)
                                                                                                     ccAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGTGATGCCAAGGTTTACACGCAGGCGCTGga                                             >  1:2618422/1‑150 (MQ=255)
                                                                                                      cAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGTGATGCCAAGGTTTACACGCAGGCGCTGgac                                            <  2:3538772/150‑1 (MQ=255)
                                                                                                      cAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGTGATGCCAAGGTTTACACGCAGGCGCTGgac                                            >  2:1820520/1‑150 (MQ=255)
                                                                                                       aGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGTGATGCCAAGCTTTACACGCAGGCGCTGgacg                                           >  2:1041527/1‑150 (MQ=255)
                                                                                                         aTCGACCTCAATGCCGATGAACACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGTGATGCCAAGGTTTACACGCAGGCGCTGgacgac                                         <  1:625956/150‑1 (MQ=255)
                                                                                                             aCTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGTGATGCCAAGGTTTACACGCAGGCGCTGGACGACGccc                                     >  1:4601849/1‑150 (MQ=255)
                                                                                                              cTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGTGATGCCAAGGTTTACACGCAGGCGCTGGACGACGCCCg                                    <  1:1156081/150‑1 (MQ=255)
                                                                                                              cTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGTGATGCCAAGGTTTACACGCAGGCGCTGGACGACGCCCg                                    <  2:4454605/150‑1 (MQ=255)
                                                                                                               ttCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGTGATGCCAAGGTTTACACGCAGGCGCTGGACGACGCCCGc                                   >  2:3853174/1‑150 (MQ=255)
                                                                                                               ttCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGTGATGCCAAGGTTTACACGCAGGCGCTGGACGACGCCCGc                                   >  1:3398631/1‑150 (MQ=255)
                                                                                                               ttCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGTGATGCCAAGGTTTACACGCAGGCGCTGGACGACGCCCGc                                   >  1:33707/1‑150 (MQ=255)
                                                                                                                    tGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGTGATGCCAAGGTTTACACGCAGGCGCTGGACGACGCCCGCAGCGt                              >  2:4121125/1‑150 (MQ=255)
                                                                                                                    tGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGTGATGCCAAGGTTTACACGCAGGCGCTGGACGACGCCCGCAGCGt                              >  2:4121187/1‑150 (MQ=255)
                                                                                                                        gATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGTGATGCCAAGGTTTACACGCAGGCGCTGGACGACGCCCGCAGCGTActg                          >  1:2518821/1‑150 (MQ=255)
                                                                                                                         aTGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGTGATGCCAAGGTTTACACGCAGGCGCTGGACGACGCCCGCAGCGTActgc                         >  2:3053223/1‑150 (MQ=255)
                                                                                                                           gAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGTGATGCCAAGGTTTACACGCAGGCGCTGGACGACGCCCGCAGCGTActgctg                       <  2:2502521/150‑1 (MQ=255)
                                                                                                                                    cGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGTGATGCCAAGGTTTACACGCAGGCGCTGGACGACGCCCGCAGCGTACTGCTGGCCAACTTc              <  2:57584/150‑1 (MQ=255)
                                                                                                                                     ttCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGTGATGCCAAGGTTTACACGCAGGCGCTGGACGACGCCCGCAGCGTACTGCTGGCCAACTTCa             >  2:4611372/2‑150 (MQ=255)
                                                                                                                                      tCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATg                                                                            >  1:1548644/1‑86 (MQ=255)
                                                                                                                                      tCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATg                                                                            <  2:1548644/86‑1 (MQ=255)
                                                                                                                                      tCCGTGGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGTGATGCCAAGGTTTACACGCAGGCGCTGGACGACGCCCGCAGCGTACTGCTGGCCAACTTCaa            <  1:2693751/150‑1 (MQ=255)
                                                                                                                                       ccGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGTGATGCCAAGGTTTACACGCAGGCGCTGGACGACGCCCGCAGCGTACTGCTGGCCAACTTCAAt           <  1:3839326/150‑1 (MQ=255)
                                                                                                                                       ccGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGTGATGCCAAGGTTTACACGCAGGCGCTGGACGACGCCCGCAGCGTACTGCTGGCCAACTTCAAt           <  1:4148950/150‑1 (MQ=255)
                                                                                                                                       ccGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGTGATGCCAAGGTTTACACGCAGGCGCTGGACGACGCCCGCAGCGTACTGCTGGCCAACTTCAAt           <  1:1065926/150‑1 (MQ=255)
                                                                                                                                       ccGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACAATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGTGATGCCAAGGTTTACACGCAGGCGCTGGACGACGCCCGCAGCGTACTGCTGGCCAACTTCAAt           <  2:5093803/150‑1 (MQ=255)
                                                                                                                                        cGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGTGATGCCAAGGTTTACACGCAGGCGCTGGACGACGCCCGCAGCGTACTGCTGGCCAACTTCAATg          <  2:4661429/150‑1 (MQ=255)
                                                                                                                                         gTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGTGATGCCAAGGTTTACACGCAGGCGCTGGACGACGCCCGCAGCGTACTGCTGGCCAACTTCAATGc         <  2:1923844/150‑1 (MQ=255)
                                                                                                                                         gTTGCTGCCCGGCCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGTGATGCCAAGGTTTACACGCAGGCGCTGGACGACGCCCGCAGCGTACTGCTGGCCAACTTCAATGc         >  2:3708919/1‑150 (MQ=255)
                                                                                                                                            gCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGTGATGCCAAGGTTTACACGCAGGCGCTGGACGACGCCCGCAGCGTACTGCTGGCCAACGTCAATGCAGa      <  2:524048/150‑1 (MQ=255)
                                                                                                                                               ggCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGTGATGCCAAGGTTTACACGCAGGCGCTGGACGACGCCCGCAGCGTACTGCTGGCCAACTTCAATGCAGACaa   >  2:219907/1‑150 (MQ=255)
                                                                                                                                               ggCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGTGATGCCAAGGTTTACACGCAGGCGCTGGACGACGCCCGCAGCGTACTGCTGGCCAACTTCAATGCAGACaa   <  1:547855/150‑1 (MQ=255)
                                                                                                                                                gCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGTGATGCCAAGGTTTACACGCAGGCGCTGGACGACGCCCGCAGCGTACTGCTGGCCAACTTCAATGCAGACAAc  <  1:3345516/150‑1 (MQ=255)
                                                                                                                                                gCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGTGATGCCAAGGTTTACACGCAGGCGCTGGACGACGCCCGCAGCGTACTGCTGGCCAACTTCAATGCAGACAAc  >  1:564791/1‑150 (MQ=255)
                                                                                                                                                gCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGTGATGCCAAGGTTTACACGCAGGCGCTGGACGACGCCCGCAGCGTACTGCTGGCCAACTTCAATGCAGACAAc  <  2:2618422/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
CCCGAGTTCGACAGTAACGCCGATGCGCTTGGCGCCCTGACCGCCGACGGTGATGGTGCCAGCAGGCTGTCGCAGTGGTGGGCAGAAATATCCAAGTACTTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGTGATGCCAAGGTTTACACGCAGGCGCTGGACGACGCCCGCAGCGTACTGCTGGCCAACTTCAATGCAGACAAC  >  NC_002947/240387‑240680

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: