Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 239,949 A→G 6.7% T82A (ACC→GCC)  PP_0188 → uroporphyrin‑III C‑methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947239,9490AG 6.7% 256.6 / 7.0 89T82A (ACC→GCC) PP_0188uroporphyrin‑III C‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (42/41);  new base G (4/2);  total (46/43)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.22e-01

CCCACCTGCCAAGCGCTCCGGCAGTGGCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGCGACCAGCAGCGTCTCAGCCAGCGCCTGG  >  NC_002947/239801‑240093
                                                                                                                                                    |                                                                                                                                                 
cccaccTGCCAAGCGCTCCGGCAGTGGCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGACAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGa                                                                                                                                                   >  7:155864/1‑149 (MQ=255)
 ccaccTGCCAAGCGCTCCGGCAGTGGCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGACAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGAc                                                                                                                                                  >  1:133254/1‑149 (MQ=255)
 ccaccTGCCAAGCGCTCCGGCAGTGGCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGACAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGAc                                                                                                                                                  <  6:105019/149‑1 (MQ=255)
  caccTGCCAAGCGCTCCGGCAGTGGCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGACAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGAcc                                                                                                                                                 <  5:89893/149‑1 (MQ=255)
    ccTGCCAAGCGCTCCGGCAGTGGCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGACAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGAcc                                                                                                                                                 >  4:82310/1‑147 (MQ=255)
     cTGCCAAGCGCTCCGGCAGTGGCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGACAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAg                                                                                                                                              >  1:110118/1‑149 (MQ=255)
      tGCCAAGCGCTCCGGCAGTGGCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGACAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAgc                                                                                                                                            >  6:265056/1‑149 (MQ=255)
        ccAAGCGCTCCGGCAGTGGCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGACAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAgcgc                                                                                                                                          <  2:133254/149‑1 (MQ=255)
        ccAAGCGCTCCGGCAGTGGCCTGGCAGCGCTGGCCTTGCTGCTTGGGGCGGCCGGGGTGGCGGTAGGTGGCTGGGGGGTCTGACAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAgcgc                                                                                                                                          >  2:237885/1‑149 (MQ=255)
          aaGCGCTCCGGCAGTGGCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGACAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAgcgcgc                                                                                                                                        >  1:62979/1‑149 (MQ=255)
          aaGCGCTCCGGCAGTGGCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGACAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAgcgcgc                                                                                                                                        <  3:82310/149‑1 (MQ=255)
          aaGCGCTCCGGCAGTGGCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGACAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAgcgcgc                                                                                                                                        >  2:265074/1‑149 (MQ=255)
          aaGCGCTCCGGCAGGGGCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGACAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGGGTCAgcgcct                                                                                                                                        <  3:276879/149‑3 (MQ=255)
               cTCCGGCAGTGGCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGACAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAgg                                                                                                                                   <  3:185045/149‑1 (MQ=255)
                tCCGGCAGTGGCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGACAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGgc                                                                                                                                  <  2:244085/149‑1 (MQ=255)
                tCCGGCAGTGGCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGACAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGgc                                                                                                                                  >  5:331623/1‑149 (MQ=255)
                  cGGCAGTGGCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGgcgc                                                                                                                                >  7:28859/1‑149 (MQ=255)
                      aGTGGCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGACAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTgca                                                                                                                            <  8:30907/149‑1 (MQ=255)
                      aGTGGCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGACAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTgca                                                                                                                            <  2:110118/149‑1 (MQ=255)
                      aGTGGCCTGGCAGCGCTGGCCTTGCTGCTGGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGACAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTg                                                                                                                              >  1:28222/1‑147 (MQ=255)
                               gCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGACAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGt                                                                                                                     <  2:62979/147‑1 (MQ=255)
                               gCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGACAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTg                                                                                                                    <  6:331623/148‑1 (MQ=255)
                                   cgcTGGCCTTGCTGCTTGGTGCGGCCGCGGTGGCGGTAGGTGGCTGGGGTGTCTGACAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAAcag                                                                                                               >  4:222615/1‑149 (MQ=255)
                                     cTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGACAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCTAGCACATCGAGGCGCTGACACAGCGCGCCGAGGCGCTGCAGCAACGTGAAcagca                                                                                                             >  8:335813/1‑149 (MQ=255)
                                     cTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGACAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCCGCGCCGAGGCGCTGCAGCAACGTGAAcagc                                                                                                              <  1:40174/149‑1 (MQ=255)
                                        gCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGACAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGAt                                                                                                          <  5:265056/149‑1 (MQ=255)
                                          cTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGACAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGAt                                                                                                          >  1:161613/1‑147 (MQ=255)
                                          cTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGACAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCa                                                                                                        >  3:211272/1‑149 (MQ=255)
                                                gctTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGACAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTgcg                                                                                                   >  7:43524/1‑148 (MQ=255)
                                                 ctTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGACAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTAACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTgcg                                                                                                   <  2:155884/147‑1 (MQ=255)
                                                     gTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGACAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACt                                                                                              <  6:318161/148‑1 (MQ=255)
                                                     gTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGACAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTg                                                                                             <  1:237885/149‑1 (MQ=255)
                                                        cGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGACAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGc                                                                                           <  2:28222/148‑1 (MQ=255)
                                                        cGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGACAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGc                                                                                           >  3:162248/1‑148 (MQ=255)
                                                        cGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGACAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGc                                                                                           <  3:116559/148‑1 (MQ=255)
                                                        cGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGACAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGc                                                                                           <  8:155864/148‑1 (MQ=255)
                                                                      gTAGGTGGCTGGGGTGTCTGACAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAgccgc                                                                            >  8:221268/1‑149 (MQ=255)
                                                                      gTAGGTGGCTGGGGTGTCTGACAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAgccgc                                                                            <  4:162248/149‑1 (MQ=255)
                                                                        aGGTGGCTGGGGTGTCTGACAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAgccgc                                                                            >  1:106900/1‑147 (MQ=255)
                                                                            ggCTGGGGTGTCTGACAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACAACGAGGCGCAGACCAAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGa                                                                      >  8:355246/1‑149 (MQ=255)
                                                                             gCTGGGGTGTCTGACAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAg                                                                     <  4:211272/149‑1 (MQ=255)
                                                                                ggggTGTCTGACAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTa                                                                  <  4:56459/149‑1 (MQ=255)
                                                                                   gtgtCTGACAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTagaa                                                               <  8:69507/149‑1 (MQ=255)
                                                                                     gtCTGACAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGcc                                                                                  <  1:68598/128‑1 (MQ=255)
                                                                                     gtCTGACAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGcc                                                                                  >  2:68598/1‑128 (MQ=255)
                                                                                      tCTGACAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTgcagca                                                                                                                         >  4:163334/1‑88 (MQ=255)
                                                                                      tCTGACAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTgcagca                                                                                                                         <  3:163334/88‑1 (MQ=255)
                                                                                      tCTGACAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAAcag                                                                                                               >  8:287515/1‑98 (MQ=255)
                                                                                      tCTGACAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAAcag                                                                                                               <  7:287515/98‑1 (MQ=255)
                                                                                       cTGACAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGAcc                                                           >  6:90683/1‑149 (MQ=255)
                                                                                        tgacAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGAccg                                                          >  3:356091/4‑149 (MQ=255)
                                                                                             ggTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGc                                                     >  2:9315/1‑149 (MQ=255)
                                                                                                  gTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCa                                                                                         >  1:250588/1‑108 (MQ=255)
                                                                                                  gTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCa                                                                                         <  2:250588/108‑1 (MQ=255)
                                                                                                   tCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTgg                                               >  5:325951/1‑149 (MQ=255)
                                                                                                   tCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTgg                                               <  2:161613/149‑1 (MQ=255)
                                                                                                     aaCTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCTCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCATCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTgtt                                              <  8:28859/148‑3 (MQ=255)
                                                                                                     aaCTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTa                                             <  8:317487/149‑1 (MQ=255)
                                                                                                          cAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCa                                        <  4:1660/149‑1 (MQ=255)
                                                                                                            gggCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGt                                      >  3:102639/1‑149 (MQ=255)
                                                                                                                aGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCa                                  <  5:90683/149‑1 (MQ=255)
                                                                                                                 gTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAg                                 >  6:63390/1‑149 (MQ=255)
                                                                                                                  tGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAgg                                <  7:221268/149‑1 (MQ=255)
                                                                                                                    aGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAggg                               <  4:102639/148‑1 (MQ=255)
                                                                                                                    aGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGc                              <  6:325951/149‑1 (MQ=255)
                                                                                                                     gTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGc                              <  8:302819/148‑1 (MQ=255)
                                                                                                                      ttCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGCGa                            >  8:145287/1‑149 (MQ=255)
                                                                                                                      ttCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGCGa                            >  7:280308/1‑149 (MQ=255)
                                                                                                                      ttCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGGCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGCGa                            <  3:222615/149‑1 (MQ=255)
                                                                                                                         gcgctgggTCAGCGCCTCGATGTGCTGGCCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGCGAc                           >  8:310179/6‑147 (MQ=255)
                                                                                                                           gctgggTCAGCGCCTCGATGTGCTGGCCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGCGACcag                        >  6:194617/4‑148 (MQ=255)
                                                                                                                           ccAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCa                                                                                 <  6:72030/91‑1 (MQ=255)
                                                                                                                           ccAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCa                                                                                 >  5:72030/1‑91 (MQ=255)
                                                                                                                           ccAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTa                                             <  6:53687/127‑1 (MQ=255)
                                                                                                                           ccAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTa                                             >  5:53687/1‑127 (MQ=255)
                                                                                                                             tgggTCAGCGCCTCGATGTGCTGGCCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTGGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGCGACcagcag                     >  6:233283/2‑149 (MQ=255)
                                                                                                                                    gCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGCGACCAGCAGCGTCTca              >  2:160793/1‑149 (MQ=255)
                                                                                                                                     cacaTCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGCGACCAGCAGCGTCTcag             >  7:328610/1‑149 (MQ=255)
                                                                                                                                        ctCGATGTGCTGGCCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGCGACCAGCAGCGTCTcagc            >  7:34613/2‑147 (MQ=255)
                                                                                                                                        aTCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGCGACCAGCAGCGTCTcagc            >  7:260774/1‑147 (MQ=255)
                                                                                                                                            aGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGCGAcc                          >  3:277654/1‑129 (MQ=255)
                                                                                                                                            aGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGCGAcc                          <  4:277654/129‑1 (MQ=255)
                                                                                                                                            aGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGCGACCAGCAGCGTCTCAGCCAgcgc      >  3:28956/1‑149 (MQ=255)
                                                                                                                                              gcgcTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGCGACCAGCAGCGTCTCAGCCAGCGCCt    <  2:106900/149‑1 (MQ=255)
                                                                                                                                               tgcTGGCCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGCGACCAGCAGCGTCTCAGCCAGCGCCt    <  5:194617/147‑1 (MQ=255)
                                                                                                                                                gcTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGTCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGc                                           <  2:65024/108‑1 (MQ=255)
                                                                                                                                                gcTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGCGACCAGCAGCGTCTCAGCCAGCGCCTg   >  5:182946/1‑148 (MQ=255)
                                                                                                                                                gcTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGCGACCAGCAGCGTCTCAGCCAGCGCCTg   >  5:188327/1‑148 (MQ=255)
                                                                                                                                                gcTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGACGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGc                                           >  1:65024/1‑108 (MQ=255)
                                                                                                                                                 cTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGCGACCAGCAGCGTCTCAGCCAGCGCCTgg  >  8:101073/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
CCCACCTGCCAAGCGCTCCGGCAGTGGCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGCGACCAGCAGCGTCTCAGCCAGCGCCTGG  >  NC_002947/239801‑240093

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: