Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 239,466 C→T 10.5% L183F (CTC→TTC)  hemD → uroporphyrinogen‑III synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947239,4660CT10.5% 205.8 / 11.5 76L183F (CTC→TTC) hemDuroporphyrinogen‑III synthase
Reads supporting (aligned to +/- strand):  ref base C (34/34);  new base T (3/5);  total (37/39)
Fisher's exact test for biased strand distribution p-value = 7.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.75e-01

GATCGTTCGCGGCGTCGGAGGTCGCGAACTGCTGGCAGAGCGTCTTGCAGAGCAAGGTGCTAGTGTCGAATATCTGGAACTGTATCGTCGCTGCCT‑‑GCCGGCCTACCCGGCGGG‑‑TA‑‑CGCTGATGCGCCGCATTGAAGCGGAACGCCTCAATGGCCTGGTGGTCAGCAGTGGGCAGGGTTTTGAACATTTGCAGCAGATGGCCGGTGCCGATTGGCTGCGGTTGGCCCAGTTGCCGCTGTTCGTGCCCAGCCCGCGGGTCGCCGAGCAGGCCAAGGCCGCCGGGGCCCAACAGGT  >  NC_002947/239321‑239614
                                                                                                                                                        |                                                                                                                                                    
gATCGTTTGCGGCGTCGGAGGTCGCGAACTGCTGGCAGAGCGTCTTGCAGAGCAAGGTGCTAGTGTCGAATATCTGGAACTGTATCGTCGCAGCCT‑‑GCCGGCCTACCCGGCGGG‑‑TA‑‑CGCTGATGCGCCGCATTGAAGCGGAACGCCTc                                                                                                                                                    >  2:246248/1‑148 (MQ=255)
    gTTCGCGGCGTCGGAGGTCGCGAACTGCTGGCAGAGCGTCTTGCAGAGCAAGGTGCTAGTGTCGAATATCTGGAACTGTATCGTCGCTGCCT‑‑GCCGGCCTACCCGGCGGG‑‑TA‑‑CGCTGATGCGCCGCATTGAGGCGTTCCGCTTCAATgc                                                                                                                                               <  3:38354/149‑2 (MQ=255)
     ttCGCGGCGTCGGAGGTCGCGAACTGCTGGCAGAGCGTCTTGCAGAGCAAGGTGCTAGTGTCGAATATCTGGAACTGTATCGTCGCTGCCT‑‑GCCGGCCTACCCGGCGGG‑‑TA‑‑CGCTGATGCGCCGCATTGAGGCGTTCCGCTTCAATgcg                                                                                                                                              <  1:167882/149‑3 (MQ=255)
        gcggcgTCGGAGGTCGCGAACTGCTGGCAGAGCGTCTTGCAGAGCAAGGTGCTAGTGTCGAATATCTGGAACTGTATCGTCGCTGCCT‑‑GCCGGCCTACCCGGCGGG‑‑TA‑‑CGCTGATGCGCCGCATTGAGGCGTTCCGCTTCAATgcggc                                                                                                                                            <  3:50371/148‑5 (MQ=255)
        gcggcgTCGGAGGTCGCGAACTGCTGGCAGAGCGTCTTGCAGAGCAAGGTGCTAGTGTCGAATATCTGGAACTGTATCGTCGCTGCCT‑‑GCCGGCCTACCCGGCGGG‑‑TA‑‑CGCTGATGCGCCGCATTGAAGCGGAACGCCTCAATGGCCt                                                                                                                                            <  6:63695/148‑1 (MQ=255)
         cggcgTCGGAGGTCGCGAACTGCTGGCAGAGCGTCTTGCAGAGCAAGGTGCTAGTGTCGAATATCTGGAACTGTATCGTCGCTGCCT‑‑GCCGGCCTACCCGGCGGG‑‑TA‑‑CGCTGATGCGCCGCATTGAAGCGGAACGCCTCAATGGCCtgg                                                                                                                                          <  4:23795/149‑1 (MQ=255)
             gTCGGAGGTCGCGAACTGCTGGCAGAGCGTCTTGCAGAGCAAGGTGCTAGTGTCGAATATCTGGAACTGTATCGTCGCTGCCT‑‑GCCGGCCTACCCGGCGGG‑‑TA‑‑CGCTGATGCGCCGCATTGAAGCGGAACGCCTCAATGGCCTggtggt                                                                                                                                      <  8:80647/149‑1 (MQ=255)
              tCGGAGGTCGCGAACTGCTGGCAGAGCGTCTTGCAGAGCAAGGTGCTAGTGTCGAATATCTGGAACTGTATCGTCGCTGCCT‑‑GCCGGCCTACCCGGCGGG‑‑TA‑‑CGCTGATGCGCCGCATTGAAGCGGAACGCCTCAATGGCCTGGTGGTc                                                                                                                                     <  5:54333/149‑1 (MQ=255)
                ggAGGTCGCGAACTGCTGGCAGAGCGTCTTGCAGAGCAAGGTGCTAGTGTCGAATATCTGGAACTGTATCGTCGCTGCCT‑‑GCCGGCCTACCCGGCGGG‑‑TA‑‑CGCTGATGCGCCGCATTGAAGCGGAACGCCTCAATGGCCTGGTGGTcag                                                                                                                                   <  5:231538/149‑1 (MQ=255)
                    gTCGCGAACTGCTTGCAGATCGTTTTGCAGATCAATGTGCTTGTGTCGTATATCTTGAACTGTATTGTCGCTTCCCTCAGCCGGCCTACCCGGCGGG‑‑TA‑‑CGCTGATGCGCCGCATTGAAGCGGAACGCCTCAATGGCCTGGTGGTcag                                                                                                                                   <  1:299893/148‑1 (MQ=255)
                    gTCGCGAACTGCTGGCAGAGCGTCTTGCAGAGCAAGGTGCTAGTGTCGAATATCTGGAACTGTATCGTCGCTGCCT‑‑GCCGGCCTACCCGGCGGG‑‑TA‑‑CGCTGATGCGCCGCATTGAAGCGGAACGCCTCAATGGCCTGGTGGTCAGCAGt                                                                                                                               >  7:283128/1‑149 (MQ=255)
                           aCTGCTGGCAGAGCGTCTTGCAGAGCAAGGTGCTAGTGTCGAATATCTGGAACTGTATCGTCGCTGCCT‑‑GCCGGCCTACCCGGCGGG‑‑TA‑‑CGCTGATGCGCCGCATTGAAGCGGAACGCCTCAATGGCCTGGTGGTCAGCAGTGGGCAgg                                                                                                                        <  2:252102/149‑1 (MQ=255)
                           aCTGCTGGCAGAGCGTCTTGCAGAGCAAGGTGCTAGTGTCGAATATCTGGAACTGTATCGTCGCTGCCT‑‑GCCGGCCTACCCGGCGGG‑‑TA‑‑CGCTGATGCGCCGCATTGAAGCGGAACGCCTCAATGGCCTGGTGGTCAGCAGTGGGCAgg                                                                                                                        <  5:43021/149‑1 (MQ=255)
                               ctgGCAGAGCGTCTTGCAGAGCAAGGTGCTAGTGTCGAATATCTGGAACTGTATCGTCGCTGCCT‑‑GCCGGCCTACCCGGCGGG‑‑TA‑‑CGCTGATGCGCCGCATTGAAGCGGAACGCCTCaa                                                                                                                                                  >  7:113121/1‑119 (MQ=255)
                               ctgGCAGAGCGTCTTGCAGAGCAAGGTGCTAGTGTCGAATATCTGGAACTGTATCGTCGCTGCCT‑‑GCCGGCCTACCCGGCGGG‑‑TA‑‑CGCTGATGCGCCGCATTGAAGCGGAACGCCTCaa                                                                                                                                                  <  8:113121/119‑1 (MQ=255)
                                tgGCAGAGCGTCTTGCAGAGCAAGGTGCTAGTGTCGAATATCTGGAACTGTATCGTCGCTGCCT‑‑GCCGGCCTACCCGGCGGG‑‑TA‑‑CGCTGATGCGCCGCATTGAAGCGGAACGCCTCAATGGCCTGGTGGTCAGCAGTGGGCAGGGtttt                                                                                                                   >  7:183088/1‑149 (MQ=255)
                                tgGCAGAGCGTCTTGCAGAGCAAGGTGCTAGTGTCGAATATCTGGAACTGTATCGTCGCTGCCT‑‑GCCGGCCTACCCGGCGGG‑‑TA‑‑CGCTGATGCGCCGCATTGAAGCGGAACGCCTCAATGGCCTGGTGGTCAGCAGTGGGCAGGGtttt                                                                                                                   >  7:241701/1‑149 (MQ=255)
                                 ggCAGAGCGTCTTGCAGAGCAAGGTGCTAGTGTCGAATATCTGGAACTGTATCGTCGCTGCCT‑‑GCCGGCCTACCCGGCGGG‑‑TA‑‑CGCTGATGCGCCGCATTGAGGCGTTCCGCTTCAATg                                                                                                                                                >  7:329065/1‑119 (MQ=255)
                                 ggCAGAGCGTCTTGCAGAGCAAGGTGCTAGTGTCGAATATCTGGAACTGTATCGTCGCTGCCT‑‑GCCGGCCTACCCGGCGGG‑‑TA‑‑CGCTGATGCGCCGCATTGAGGCGTTCCGCTTCAATg                                                                                                                                                <  8:329065/119‑1 (MQ=255)
                                  gCAGAGCGTCTTGCAGAGCAAGGTGCTAGTGTCGAATATCTGGAACTGTATCGTCGCTGCCT‑‑GCCGGCCTACCCGGCGGG‑‑TA‑‑CGCTGATGCGCCGCATTGAAGCGGAACGCCTCAATGGCCTGGTGGTCAGCAGTGGGCAGGGTTTTg                                                                                                                  >  3:21854/1‑148 (MQ=255)
                                     gagCGTCTTGCAGAGCAAGGTGCTAGTGTCGAATATCTGGAACTGTATCGTCGCTGCCT‑‑GCCGGCCTACCCGGCGGG‑‑TA‑‑CGCTGATGCGCCGCATTGAAGCGGAACGCCTCAATGGCCTGGTGGTCAGCAGTGGGCAGGGTTTTGAACa                                                                                                              >  7:206917/1‑149 (MQ=255)
                                          tCTTGCAGAGCAAGGTGCTAGTGTCGAATATCTGGAACTGTATCGTCGCTGCCT‑‑GCCGGCCTACCCGGCGGG‑‑TA‑‑CGCTGATGCGCCGCATTGAAGCGGAACGCCTCAATGGCCTGGTGGTCAGCAGTGGGCAGGGTTTTGAACATTTgc                                                                                                         <  8:183088/149‑1 (MQ=255)
                                          tCTTGCAGAGCAAGGTGCTAGTGTCGAATATCTGGAACTGTATCGTCGCGGCCT‑‑GCCGGCCTACCCGGCGGG‑‑TA‑‑CGCTGATGCGCCGCATTGAAGCGGAACGCCTCAATGGCCTGGTGGTCAGCAGTGGGCAGGGTTTTGAACATTTgc                                                                                                         <  1:246248/149‑1 (MQ=255)
                                            ttGCAGAGCAAGGTGCTAGTGTCGAATATCTGGAACTGTATCGTCGCTGCCT‑‑GCCGGCCTACCCGGCGGG‑‑TA‑‑CGCTGATGCGCCGCATTGAAGCGGAACGCCTCAATGGCCTGGTGGTCAGCAGTGGGCAGGGTTTTGAACATTTgcag                                                                                                       >  4:152089/1‑149 (MQ=255)
                                              gCAGAGCAAGGTGCTAGTGTCGAATATCTGGAACTGTATCGTCGCTGCCT‑‑GCCGGCCTACCCGGCGGGTTTAGGCACACATGCGCCGCATTCAAGCGGAACGCCACAATGCCCTGGTGGTCAGCAGTGGGCAGGTTTTTGAACAtttc                                                                                                          >  1:114948/1‑147 (MQ=255)
                                              gCAGAGCAAGGTGCTAGTGTCGAATATCTGGAACTGTATCGTCGCTGCCT‑‑GCCGGCCTACCCGGCGGG‑‑TA‑‑CGCTGATGCGCCGCATTGAGGCGTTCCGCTTCAATg                                                                                                                                                <  2:164758/106‑1 (MQ=255)
                                              gCAGAGCAAGGTGCTAGTGTCGAATATCTGGAACTGTATCGTCGCTGCCT‑‑GCCGGCCTACCCGGCGGG‑‑TA‑‑CGCTGATGCGCCGCATTGAGGCGTTCCGCTTCAATg                                                                                                                                                >  1:164758/1‑106 (MQ=255)
                                              gCAGAGCAAGGTGCTAGTGTCGAATATCTGGAACTGTATCGTCGCTGCCT‑‑GCCGGCCTACCCGGCGGG‑‑TA‑‑CGCTGATGCGCCGCATTGAAGCGGAACGCCTCAATGGCCTGGTGGTCAGCAGTGGGCAGGGTTTTGAACATTTgcag                                                                                                       >  6:254881/1‑147 (MQ=255)
                                              gCAGAGCAAGGTGCTAGTGTCGAATATCTGGAACTGTATCGTCGCTGCCT‑‑GCCGGCCTACCCGGCGGG‑‑TA‑‑CGCTGATGCGCCGCATTGAAGCGGAACGCCTCAATGGCCTGGTGGTCAGCAGTGGGCAGGGTTTTGAACATTTgcag                                                                                                       >  1:28525/1‑147 (MQ=255)
                                              gCAGAGCAAGGTGCTAGTGTCGAATATCTGGAACTGTATCGTCGCTGCCT‑‑GCCGGCCTACCCGGCGGG‑‑TA‑‑CGCTGATGCGCCGCATTGAAGCGGAACGCCTCAATGGCCTGGTGGTCAGCAGTGGGCAGGGTTTTGAACATTTgcag                                                                                                       <  8:51450/147‑1 (MQ=255)
                                                agagCAAGGTGCTAGTGTCGAATATCTGGAACTGTATCGTCGCTGCCT‑‑GCCGGCCTACCAGGCGGG‑‑TA‑‑CGCTGATGCGCCGCATTGAAGCGGAACGCCTCAATGGCCTGGTGGTCAGCAGTCGGAAGGGTTTTGAACATTTGaagcag                                                                                                    >  8:263266/1‑148 (MQ=255)
                                                  agCAAGGTGCTAGTGTCGAATATCTGGAACTGTATCGTCGCTGCCT‑‑GCCGGCCTACCCGGCGGG‑‑TA‑‑CGCTGATGCGCCGCATTGAAGCGGAACGCCTCAATGGCCTGGTGGTCAGCAGTGGGCAGGGTTTTGAACATTTg                                                                                                          <  7:216602/140‑1 (MQ=255)
                                                  agCAAGGTGCTAGTGTCGAATATCTGGAACTGTATCGTCGCTGCCT‑‑GCCGGCCTACCCGGCGGG‑‑TA‑‑CGCTGATGCGCCGCATTGAAGCGGAACGCCTCAATGGCCTGGTGGTCAGCAGTGGGCAGGGTTTTGAACATTTg                                                                                                          >  8:216602/1‑140 (MQ=255)
                                                             aGTGTCGAATATCTGGAACTGTATCGTCGCTGCCT‑‑GCCGGCCTACCCGGCGGG‑‑TA‑‑CGCTGATGCGCCGCATTGAAGCGGAACGCCTCAATGGCCTGGTGGTCAGCAGTGGGCAGGGTTTTGaa                                                                                                                >  2:217332/1‑123 (MQ=255)
                                                             aGTGTCGAATATCTGGAACTGTATCGTCGCTGCCT‑‑GCCGGCCTACCCGGCGGG‑‑TA‑‑CGCTGATGCGCCGCATTGAAGCGGAACGCCTCAATGGCCTGGTGGTCAGCAGTGGGCAGGGTTTTGaa                                                                                                                <  1:217332/123‑1 (MQ=255)
                                                              gtgtCGAATATCTGGAACTGTATCGTCGCTGCCT‑‑GCCGGCCTACCCGGCGGG‑‑TA‑‑CGCTGATGCGCCGCATTGAAGCGGAACGCCTCAATGGCCTGGTGGTCAGCAGTGGGCAGGGTTTTGAACATTTGCAGCAGATGGCCGGTGCCGAt                                                                                     <  7:123475/149‑1 (MQ=255)
                                                               tgtCGAATATCTGGAACTGTATCGTCGCTGCCT‑‑GCCGGCCTACCCGGCGGG‑‑TA‑‑CGCTGATGCGCCGCATTGAAGCGGAACGCCTCAATGGCCTGGTGGTCAGCAGTGGGCAGGGTTTTGAACATTTGCAGCAGATGGCCGGTGCCGAtt                                                                                    <  1:125304/149‑1 (MQ=255)
                                                               tgtCGAATATCTGGAACTGTATCGTCGCTGCCT‑‑GCCGGCCTACCCGGCGGG‑‑TA‑‑CGCTGATGCGCCGCATTGAAGCGGAACGCCTCAATGGCCTGGTGGTCAGCAGTGGGCAGGGTTTTGAACATTTGCAGCAGATGGCCGGTGCCGAt                                                                                     <  5:254881/148‑1 (MQ=255)
                                                                gtCGAATATCTGGAACTGTATCGTCGCTGCCT‑‑GCCGGCCTACCCGGCGGG‑‑TA‑‑CGCTGATGCGCCTCATTGAAGCGGAACGCCTCAATGGCCTGGTGGTCAGCAGTGGGCAGGGTTTTGAACATTTGCGGCAGATGGCCGGTGCCGAtt                                                                                    <  4:205313/148‑1 (MQ=255)
                                                                gtCGAATATCTGGAACTGTATCGTCGCTGCCT‑‑GCCGGCCTACCCGGCGGG‑‑TA‑‑CGCTGATGCGCCGCATTGAAGCGGAACGCCTCAATGGCCTGGTGGTCAGCAGTGGGCAGGGTTTTGAACATTTGCGGCAGATGGCCGGTGCCGAtt                                                                                    >  3:205313/1‑148 (MQ=255)
                                                                 tCGAATATCTGGAACTGTATCGTCGCTGCCT‑‑GCCGGCCTACCCGGCGGG‑‑TA‑‑CGCTGATGCGCCGCATTGAAGCGGAACGCCTCAATGGCCTGGTGGTCAGCAGTGGGCAGGGTTTTGAACATTTGCAGCAGATGGCCGGTGCCGATTgg                                                                                  <  5:40007/149‑1 (MQ=255)
                                                                 tCGAATATCTGGAACTGTATCGTCGCTGCCT‑‑GCCGGCCTACCCGGCGGG‑‑TA‑‑CGCTGATGCGCCGCATTGAAGCGGAACGCCTCAATGGCCTGGTGGTCAGCAGTGGGCAGGGTTTTGAACATTTGCAGCAGATGGCCGGTGCCGATTgg                                                                                  >  5:43084/1‑149 (MQ=255)
                                                                    aaTATCTGGAACTGTATCGTCGCTGCCT‑‑GCCGGCCTACCCGGCGGG‑‑TA‑‑CGCTGATGCGCCGCATTGAAGCGGAACGCCTCAATGGCCTGGTGGTCAGCAGTGGGCAGGGTTTTGAACATTTGCAGCAGATGGCCGGTGCCGATTGGCTg                                                                               >  6:156047/1‑149 (MQ=255)
                                                                     atatCTGGAACTGTATCGTCGCTGCCT‑‑GCCGGCCTACCCGGCGGG‑‑TA‑‑CGCTGATGCGCCGCATTGAAGCGGAACGCCTCAATGGCCTGGTGGTCAGCAGTGGGCAGGGTTTTGAACATTTGCAGCAGATGGCCGGTGCCGATTGGCTGc                                                                              >  3:31050/1‑149 (MQ=255)
                                                                         cTGGAACTGTATCGTCGCTGCCT‑‑GCCGGCCTACCCGGCGGG‑‑TA‑‑CGCTGATGCGCCGCATTGAAGCGGAACGCCTCAATGGCCTGGTGGTCAGCAGTGGGCAGGGTTTTGAACATTTGCAGCAGATGGCCGGTGCCGATTGGCTGCGGtt                                                                          >  1:134818/1‑149 (MQ=255)
                                                                           ggAACTGTATCGTCGCTGCCT‑‑GCCGGCCTACCCGGCGGG‑‑TA‑‑CGCTGATGCGCCGCATTGAAGCGGAACGCCTCAATGGCCTGGTGGTCAGCAGTGGGCAGGGTTTTGAACATTTGCAGCAGATGGCCGGt                                                                                           >  4:14236/1‑130 (MQ=255)
                                                                           ggAACTGTATCGTCGCTGCCT‑‑GCCGGCCTACCCGGCGGG‑‑TA‑‑CGCTGATGCGCCGCATTGAAGCGGAACGCCTCAATGGCCTGGTGGTCAGCAGTGGGCAGGGTTTTGAACATTTGCAGCAGATGGCCGGt                                                                                           <  3:14236/130‑1 (MQ=255)
                                                                            gAACTGTATCGTCGCTGCCT‑‑GCCGGCCTACCCGGCGGG‑‑TA‑‑CGCTGATGCGCCGCATTGAAGCGGAACGCCTCAATGGCCTGGTGGTCAGCAGTGGGCAGGGTTTTGAACATTTGCAGCAGATGGCCGGTGCCGATTGGCTGCGGTTg                                                                         <  5:156047/147‑1 (MQ=255)
                                                                                tGTATCGTCGCTGCCT‑‑GCCGGCCTACCCGGCGGG‑‑TA‑‑CGCTGATGCGCCGCATTGAAGCGGAACGCCTCAATGGCCTGGTGGTCAGCAGTGGGCAGGGTTTTGAACATTTGCAGCAGATGGCCGGTGCCGATTGGCTGCGGTTGGCCCAg                                                                   >  8:75571/1‑149 (MQ=255)
                                                                                tGTATCGTCGCTGCCT‑‑GCCGGCCTACCCGGCGGG‑‑TA‑‑CGCTGATGCGCCGCATTGAAGCGGAACGCCTCAATGGCCTGGTGGTCAGCAGTGGGCAGGGTTTTGAACATTTGCAGCAGATGGCCGGTGCCGATTGGCTGCGGTTGGCCCAg                                                                   >  3:249863/1‑149 (MQ=255)
                                                                                  tATCGTCGCTGCCT‑‑GCCGGCCTACCCGGCGGG‑‑TA‑‑CGCTGATGCGCCGCATTGAAGCGGAACGCCTCAATGGCCTGGTGGTCAGCAGTGGGCAGGGTTTTGAACATTTGCAGCAGATGGCCGGTGCCGATTGGCTGCGGTTGGCCCAg                                                                   <  3:152089/147‑1 (MQ=255)
                                                                                  tATCGTCGCTGCCT‑‑GCCGGCCTACCCGGCGGG‑‑TA‑‑CGCTGATGCGCCGCATTGAAGCGGAACGCCTCAATGGCCTGGTGGTCAGCAGTGGGCAGGGTTTTGAACATTTGCAGCAGATGGCCGGTGCCGATTGGCTGCGGTTGGCCCAGtt                                                                 >  3:258086/1‑149 (MQ=255)
                                                                                  tATCGTCGCTGCCT‑‑GCCGGCCTACCCGGCGGG‑‑TA‑‑CGCTGATGCGCCGCATTGAAGCGGAACGCCTCAATGGCCTGGTGGTCAGCAGTGGGCAGGGTTTTGAACATTTGCAGCAGATGGCCGGTGCCGATTGGCTGCGGTTGGCCCAGtt                                                                 <  2:134818/149‑1 (MQ=255)
                                                                                        cgCTGCCT‑‑CCCGGCCTACCCGGGGGG‑‑TA‑‑CGCTGTTCCGCCGCATTGCAGGGGATCGCCTGAATGGCCTGGAGTTCTGCAGTGGGCAGGGTTTTGAACATTTGCTGGAGATGTCCGGTGCCGATTGGCTTCGGTTGGCCCAGTTGCCGCt                                                           <  4:288598/149‑1 (MQ=255)
                                                                                         gCTGCCT‑‑GCCGGCCTACCCGGCGGG‑‑TA‑‑CGCTGATGCGCCGCATTGAAGCGGAACGCCTCAATGGCCTGGTGGTCAGCAGTGGGCAGGGTTTTGAACATTTGCAGCAGATGGCCGGTGCCGATTGGCTGCGGTTGGCCCAGTTGCCGCTg                                                          <  6:43084/149‑1 (MQ=255)
                                                                                         gCTGCCT‑‑GCCGGCCTACCCGGCGGG‑‑TA‑‑CGCTGATGCGCCGCATTGAAGCGGAACGCCTCAATGGCCTGGTGGTCAGCAGTGGGCAGGGTTTTGAACATTTGCAGCAGATGGCCGGTGCCGATTGGCTGCGGTTGGCCCAGTTGCCGCTg                                                          <  8:241701/149‑1 (MQ=255)
                                                                                              cct‑‑gccGGCCTACCCGGCGGG‑‑TA‑‑CGCTGATGCGCCGCATTGAAGCGGAACGCCTCAATGGCCTGGTGGTCAGCAGTGGGCAGGGTTTTGAACATTTGCAGCAGATGGCCGGTGCCGATTGGCTGCGGTTGGCCCAGTTGCCGCTGTTc                                                       <  4:31050/148‑1 (MQ=255)
                                                                                                t‑‑gccGGCCTACCCGGCGGG‑‑TA‑‑CGCTGATGCGCCGCATTGAAGCGGAACGCCTCAATGGCCTGGTGGTCAGCAGTGGGCAGGGTTTTGCACATTTGCAGCAGATGGCCGGTGCCGATTGGCTGCGGTTGGCCCAGTTGCCGCTGTTCGTg                                                    >  1:61050/1‑149 (MQ=255)
                                                                                                t‑‑gccGGCCTACCCGGCGGG‑‑TA‑‑CGCTGATGCGCCGCATTGAAGCGGAACGCCTCAATGGCCTGGTGGTCAGCAGTGGGCAGGGTTTTGAACATTTgc                                                                                                         >  2:153008/1‑96 (MQ=255)
                                                                                                t‑‑gccGGCCTACCCGGCGGG‑‑TA‑‑CGCTGATGCGCCGCATTGAAGCGGAACGCCTCAATGGCCTGGTGGTCAGCAGTGGGCAGGGTTTTGAACATTTgc                                                                                                         <  1:153008/96‑1 (MQ=255)
                                                                                                t‑‑gccGGCCTACCCGGCGGG‑‑TA‑‑CGCTGATGCGCCGCATTGAAGCGGAACGCCTCAATGGCCTGGTGGTCAGCAGTGGGCAGGGTTTTGAACATTTGCAGCAGATGGCCGGTGCCGATTGGCTGCGGTTGGCCCAGTTGCCGCTGTTCGTg                                                    >  4:198724/1‑149 (MQ=255)
                                                                                                            cccGGCGGG‑‑TA‑‑CGCTGATGCGCCGCATTGAAGCGGAACGCCTCAATGGCCTGGGGGTCAGCAGTGGGCAGGGTTTTGAACATTTGCAGCAGATGGCCGGTGCCGATTGGCTGCGGTTGGCCCAGTTGCCGCTGTTCGTGCCCAGCCcgc                                          <  2:114948/149‑1 (MQ=255)
                                                                                                             ccGGCGGG‑‑TA‑‑CGCTGATGCGCCGCATTGAAGCGGAACGCCTCAATGGCCTGGTGGTCAGCAGTGGGCAGGGTTTTGAACATTTGCAGCAGATGGCCGGTGCCGATTGGCTGCGGTTGGCCCAGTTGCCGCTGTTCGTGCCCAGCCcgcg                                         >  8:111724/1‑149 (MQ=255)
                                                                                                              cggcggG‑‑TA‑‑CGCTGATGCGCCGCATTGAAGCGGAACGCCTCAATGTCCTGGTGGTCAGCAGTGGGCAGGGTTTTGAACATTTGCAGGAGATGGCCGGTGCCGATTGGCTGCGGTTGGCCCAGTTGCCGCTGTTCGTGCCCAGCCCGCgg                                        <  7:263266/149‑1 (MQ=255)
                                                                                                              cggcggG‑‑TA‑‑CGCTGATGCGCCGCATTGAAGCGGAACGCCTCAATGGCCTGGTGGTCAGCAGt                                                                                                                               >  6:379900/1‑62 (MQ=255)
                                                                                                              cggcggG‑‑TA‑‑CGCTGATGCGCCGCATTGAAGCGGAACGCCTCAATGGCCTGGTGGTCAGCAGt                                                                                                                               <  5:379900/62‑1 (MQ=255)
                                                                                                                           cGCTGATGCGCCGCATTGAAGCGGAACGCCTCAATGGCCTGGTGGTCAGCAGTGGGCAGGGTTTTGAACATTTGCAGCAGATGGCCGGTGCCGATTGGCTGCGGTTGGCCCAGTTGCCGCTGTTCGTGCCCAGCCCGCGGGTCGCCGa                                <  2:28525/148‑1 (MQ=255)
                                                                                                                                  gcgcCGCATTGAAGCGGAACGCCTCAATGGCCTGGTGGTCAGCAGTGGGCAGGGTTTTGAACATTTGCAGCAGATGGCCGGTGCCGATTGGCTGCGGTTGGCCCAGTTGCCGCTGTTCGTGCCCAGCCCGCGGGTCGCCGAGCAGGcc                         >  1:231569/1‑148 (MQ=255)
                                                                                                                                    gccgcATTGAAGCGGAACGCCTCAATGGCCTGGTGGTCAGCAGTGGGCAGGGTTTTGAACATTTGCAGCAGATGGCCGGTGCCGATTGGCTGCGGTTGGCCCAGTTGCCGCTGTTCGTGCCCAGCCCGCGGGTCGCCGAGCAGGCCAAg                      >  6:246650/1‑149 (MQ=255)
                                                                                                                                        cATTGAAGCGGAACGCCTCAATGGCCTGGTGGTCAGCAGTGGGCAGGGTTTTGAACATTTGCAGCAGATGGCCGGTGCCGATTGGCTGCGGTTGGCCCAGTTGCCGCTGTTCGTGCCCAGCCCGCGGGTCGCCGAGCAGGCCAAGgccg                  >  2:11096/1‑149 (MQ=255)
                                                                                                                                         aTTGAAGCGGAACGCCTCAATGGCCTGGTGGTCAGCAGTGGGCAGGGTTTTGAACATTTGCAGCAGATGGCCGGTGCCGATTGGCTGCGGTTGGCCCAGTTGCCGCTGTTCGTGCCCAGCCCGCGGGTCGCCGAGCAGGCCAAGgccgc                 >  2:252254/1‑149 (MQ=255)
                                                                                                                                            gAAGCGGAACGCCTCAATGGCCTGGTGGTCAGCAGTGGGCAGGGTTTTGAACATTTGCAGCAGATGGCCGGTGCCGATTGGCTGCGGTTGGCCCAGTTGCCGCTGTTCGTg                                                    <  2:205179/111‑1 (MQ=255)
                                                                                                                                            gAAGCGGAACGCCTCAATGGCCTGGTGGTCAGCAGTGGGCAGGGTTTTGAACATTTGCAGCAGATGGCCGGTGCCGATTGGCTGCGGTTGGCCCAGTTGCCGCTGTTCGTg                                                    >  1:205179/1‑111 (MQ=255)
                                                                                                                                             aaGCGGAACGCCTCAATGGCCTGGTGGTCAGCAGTGGGCAGGGTTTTGAACATTTGCAGCAGATGGCCGGTGCCGATTGGCTGCGGTTGGCCCAGTTGCCGCTGTTCGTGCCCAGCCCGCGGGTCGCCGAGCAGGCCAAGGCCGCCggg             >  6:132360/1‑149 (MQ=255)
                                                                                                                                             aaGCGGAACGCCTCAATGGCCTGGTGGTCAGCAGTGGGCAGGGTTTTGAACATTTGCAGCAGATGGCCGGTGCCGATTGGCTGCGGTTGGCCCAGTTGCCGCTGTTCGTGCCCAGCCCGCGGGTCGCCGAGCAGGCCAAGGCCGCCggg             <  1:61537/149‑1 (MQ=255)
                                                                                                                                               gcgttccGCTTCAATGGCCTGGTGGTCAGCAGTGGGCAGGGTTTTGAACATTTGCAGCAGATGGCCGGTGCCGATTGGCTGCGGTTGGCCCAGTTGCCGCTGTTCGTGCCCAGCCCGCGGGTCGCCGAGCAGGCCAAGGCCGCCGGGGc           >  5:348034/7‑149 (MQ=255)
                                                                                                                                                       ccTCAATGGCCTGGTGGTCAGCAGTGGGCAGGGTTTTGAACATTTGCAGCAGATGGCCGGTGCCGATTGGCTGCGGTTGGCCCAGTTGCCGCTGTTCGTGCCCAGCCCGCGGGTCGCCGAGCAGGCCAAGGCCGCCGGGGCCCAACAgg   >  4:343711/1‑149 (MQ=255)
                                                                                                                                                        cTCAATGGCCTGGTGGTCAGCAGTGGGCAGGGTTTTGAACATTTGCAGCAGATGGCCGGTGCCGATTGGCTGCGGTTGGCCCAGTTGCCGCTGTTCGTGCCCa                                                <  5:34407/103‑1 (MQ=255)
                                                                                                                                                        cTCAATGGCCTGGTGGTCAGCAGTGGGCAGGGTTTTGAACATTTGCAGCAGATGGCCGGTGCCGATTGGCTGCGGTTGGCCCAGTTGCCGCTGTTCGTGCCCa                                                >  6:34407/1‑103 (MQ=255)
                                                                                                                                                        cTCAATGGCCTGGTGGTCAGCAGTGGGCAGGGTTTTGAACATTTGCAGCAGATGGCCGGTGCCGATTGGCTGCGGTTGGCCCAGTTGCCGCTGTTCGTGCCCAGCCCGCGGGTCGCCGAGCAGGCCAAGGCCGCCGGGGCCCAACAGGt  >  5:350745/1‑149 (MQ=255)
                                                                                                                                                        cTCAATGGCCTGGTGGTCAGCAGTGGGCAGGGTTTTGAACATTTGCAGCAGATGGCCGGTGCCGATTGGCTGCGGTTGGCCCAGTTGCCGCTGTTCGTGCCCAGCCCGCGGGTCGCCGAGCAGGCCAAGGCCGCCGGGGCCCAACAGGt  >  5:55969/1‑149 (MQ=255)
                                                                                                                                                        |                                                                                                                                                    
GATCGTTCGCGGCGTCGGAGGTCGCGAACTGCTGGCAGAGCGTCTTGCAGAGCAAGGTGCTAGTGTCGAATATCTGGAACTGTATCGTCGCTGCCT‑‑GCCGGCCTACCCGGCGGG‑‑TA‑‑CGCTGATGCGCCGCATTGAAGCGGAACGCCTCAATGGCCTGGTGGTCAGCAGTGGGCAGGGTTTTGAACATTTGCAGCAGATGGCCGGTGCCGATTGGCTGCGGTTGGCCCAGTTGCCGCTGTTCGTGCCCAGCCCGCGGGTCGCCGAGCAGGCCAAGGCCGCCGGGGCCCAACAGGT  >  NC_002947/239321‑239614

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: