Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 239,460 G→T 9.2% E181* (GAA→TAA) ‡ hemD → uroporphyrinogen‑III synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947239,4600GT 9.2% 230.8 / 11.9 76E181* (GAA→TAA) ‡hemDuroporphyrinogen‑III synthase
Reads supporting (aligned to +/- strand):  ref base G (35/34);  new base T (2/5);  total (37/39)
Fisher's exact test for biased strand distribution p-value = 4.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.75e-01

TCTTGATCGTTCGCGGCGTCGGAGGTCGCGAACTGCTGGCAGAGCGTCTTGCAGAGCAAGGTGCTAGTGTCGAATATCTGGAACTGTATCGTCGCTGCCT‑‑GCCGGCCTACCCGGCGGG‑‑TA‑‑CGCTGATGCGCCGCATTGAAGCGGAACGCCTCAATGGCCTGGTGGTCAGCAGTGGGCAGGGTTTTGAACATTTGCAGCAGATGGCCGGTGCCGATTGGCTGCGGTTGGCCCAGTTGCCGCTGTTCGTGCCCAGCCCGCGGGTCGCCGAGCAGGCCAAGGCCGCCGGG  >  NC_002947/239317‑239603
                                                                                                                                                      |                                                                                                                                               
tCTTGATCGTTCGCGGCGTCGGAGGTCGCGAACTGCTGGCAGAGCGTCTTGCAGAGCAAGGTGCTAGTGTCGAATATCTGGAACTGTATCGTCGCTGCCT‑‑GCCGGCCTACCCGGCGGG‑‑TA‑‑CGCTGATGCGCCGCATTGAAGCGGAACGc                                                                                                                                            >  6:40007/1‑149 (MQ=255)
    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                                                                                                             <  5:43021/149‑1 (MQ=255)
                               aCTGCTGGCAGAGCGTCTTGCAGAGCAAGGTGCTAGTGTCGAATATCTGGAACTGTATCGTCGCTGCCT‑‑GCCGGCCTACCCGGCGGG‑‑TA‑‑CGCTGATGCGCCGCATTGAAGCGGAACGCCTCAATGGCCTGGTGGTCAGCAGTGGGCAgg                                                                                                             <  2:252102/149‑1 (MQ=255)
                                   ctgGCAGAGCGTCTTGCAGAGCAAGGTGCTAGTGTCGAATATCTGGAACTGTATCGTCGCTGCCT‑‑GCCGGCCTACCCGGCGGG‑‑TA‑‑CGCTGATGCGCCGCATTGAAGCGGAACGCCTCaa                                                                                                                                       <  8:113121/119‑1 (MQ=255)
                                   ctgGCAGAGCGTCTTGCAGAGCAAGGTGCTAGTGTCGAATATCTGGAACTGTATCGTCGCTGCCT‑‑GCCGGCCTACCCGGCGGG‑‑TA‑‑CGCTGATGCGCCGCATTGAAGCGGAACGCCTCaa                                                                                                                                       >  7:113121/1‑119 (MQ=255)
                                    tgGCAGAGCGTCTTGCAGAGCAAGGTGCTAGTGTCGAATATCTGGAACTGTATCGTCGCTGCCT‑‑GCCGGCCTACCCGGCGGG‑‑TA‑‑CGCTGATGCGCCGCATTGAAGCGGAACGCCTCAATGGCCTGGTGGTCAGCAGTGGGCAGGGtttt                                                                                                        >  7:241701/1‑149 (MQ=255)
                                    tgGCAGAGCGTCTTGCAGAGCAAGGTGCTAGTGTCGAATATCTGGAACTGTATCGTCGCTGCCT‑‑GCCGGCCTACCCGGCGGG‑‑TA‑‑CGCTGATGCGCCGCATTGAAGCGGAACGCCTCAATGGCCTGGTGGTCAGCAGTGGGCAGGGtttt                                                                                                        >  7:183088/1‑149 (MQ=255)
                                     ggCAGAGCGTCTTGCAGAGCAAGGTGCTAGTGTCGAATATCTGGAACTGTATCGTCGCTGCCT‑‑GCCGGCCTACCCGGCGGG‑‑TA‑‑CGCTGATGCGCCGCATTGAGGCGTTCCGCTTCAATg                                                                                                                                     <  8:329065/119‑1 (MQ=255)
                                     ggCAGAGCGTCTTGCAGAGCAAGGTGCTAGTGTCGAATATCTGGAACTGTATCGTCGCTGCCT‑‑GCCGGCCTACCCGGCGGG‑‑TA‑‑CGCTGATGCGCCGCATTGAGGCGTTCCGCTTCAATg                                                                                                                                     >  7:329065/1‑119 (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                                                                                            >  1:28525/1‑147 (MQ=255)
                                                  gCAGAGCAAGGTGCTAGTGTCGAATATCTGGAACTGTATCGTCGCTGCCT‑‑GCCGGCCTACCCGGCGGG‑‑TA‑‑CGCTGATGCGCCGCATTGAAGCGGAACGCCTCAATGGCCTGGTGGTCAGCAGTGGGCAGGGTTTTGAACATTTgcag                                                                                            <  8:51450/147‑1 (MQ=255)
                                                  gCAGAGCAAGGTGCTAGTGTCGAATATCTGGAACTGTATCGTCGCTGCCT‑‑GCCGGCCTACCCGGCGGG‑‑TA‑‑CGCTGATGCGCCGCATTGAAGCGGAACGCCTCAATGGCCTGGTGGTCAGCAGTGGGCAGGGTTTTGAACATTTgcag                                                                                            >  6:254881/1‑147 (MQ=255)
                                                    agagCAAGGTGCTAGTGTCGAATATCTGGAACTGTATCGTCGCTGCCT‑‑GCCGGCCTACCAGGCGGG‑‑TA‑‑CGCTGATGCGCCGCATTGAAGCGGAACGCCTCAATGGCCTGGTGGTCAGCAGTCGGAAGGGTTTTGAACATTTGaagcag                                                                                         >  8:263266/1‑148 (MQ=255)
                                                      agCAAGGTGCTAGTGTCGAATATCTGGAACTGTATCGTCGCTGCCT‑‑GCCGGCCTACCCGGCGGG‑‑TA‑‑CGCTGATGCGCCGCATTGAAGCGGAACGCCTCAATGGCCTGGTGGTCAGCAGTGGGCAGGGTTTTGAACATTTg                                                                                               >  8:216602/1‑140 (MQ=255)
                                                      agCAAGGTGCTAGTGTCGAATATCTGGAACTGTATCGTCGCTGCCT‑‑GCCGGCCTACCCGGCGGG‑‑TA‑‑CGCTGATGCGCCGCATTGAAGCGGAACGCCTCAATGGCCTGGTGGTCAGCAGTGGGCAGGGTTTTGAACATTTg                                                                                               <  7:216602/140‑1 (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:43084/1‑149 (MQ=255)
                                                                     tCGAATATCTGGAACTGTATCGTCGCTGCCT‑‑GCCGGCCTACCCGGCGGG‑‑TA‑‑CGCTGATGCGCCGCATTGAAGCGGAACGCCTCAATGGCCTGGTGGTCAGCAGTGGGCAGGGTTTTGAACATTTGCAGCAGATGGCCGGTGCCGATTgg                                                                       <  5:40007/149‑1 (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                                                                                <  3:14236/130‑1 (MQ=255)
                                                                               ggAACTGTATCGTCGCTGCCT‑‑GCCGGCCTACCCGGCGGG‑‑TA‑‑CGCTGATGCGCCGCATTGAAGCGGAACGCCTCAATGGCCTGGTGGTCAGCAGTGGGCAGGGTTTTGAACATTTGCAGCAGATGGCCGGt                                                                                >  4:14236/1‑130 (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                                                      <  2:134818/149‑1 (MQ=255)
                                                                                      tATCGTCGCTGCCT‑‑GCCGGCCTACCCGGCGGG‑‑TA‑‑CGCTGATGCGCCGCATTGAAGCGGAACGCCTCAATGGCCTGGTGGTCAGCAGTGGGCAGGGTTTTGAACATTTGCAGCAGATGGCCGGTGCCGATTGGCTGCGGTTGGCCCAGtt                                                      >  3:258086/1‑149 (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                                                                                                                    <  5:379900/62‑1 (MQ=255)
                                                                                                                  cggcggG‑‑TA‑‑CGCTGATGCGCCGCATTGAAGCGGAACGCCTCAATGGCCTGGTGGTCAGCAGt                                                                                                                    >  6:379900/1‑62 (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                                         >  1:205179/1‑111 (MQ=255)
                                                                                                                                                gAAGCGGAACGCCTCAATGGCCTGGTGGTCAGCAGTGGGCAGGGTTTTGAACATTTGCAGCAGATGGCCGGTGCCGATTGGCTGCGGTTGGCCCAGTTGCCGCTGTTCGTg                                         <  2:205179/111‑1 (MQ=255)
                                                                                                                                                 aaGCGGAACGCCTCAATGGCCTGGTGGTCAGCAGTGGGCAGGGTTTTGAACATTTGCAGCAGATGGCCGGTGCCGATTGGCTGCGGTTGGCCCAGTTGCCGCTGTTCGTGCCCAGCCCGCGGGTCGCCGAGCAGGCCAAGGCCGCCggg  <  1:61537/149‑1 (MQ=255)
                                                                                                                                                 aaGCGGAACGCCTCAATGGCCTGGTGGTCAGCAGTGGGCAGGGTTTTGAACATTTGCAGCAGATGGCCGGTGCCGATTGGCTGCGGTTGGCCCAGTTGCCGCTGTTCGTGCCCAGCCCGCGGGTCGCCGAGCAGGCCAAGGCCGCCggg  >  6:132360/1‑149 (MQ=255)
                                                                                                                                                      |                                                                                                                                               
TCTTGATCGTTCGCGGCGTCGGAGGTCGCGAACTGCTGGCAGAGCGTCTTGCAGAGCAAGGTGCTAGTGTCGAATATCTGGAACTGTATCGTCGCTGCCT‑‑GCCGGCCTACCCGGCGGG‑‑TA‑‑CGCTGATGCGCCGCATTGAAGCGGAACGCCTCAATGGCCTGGTGGTCAGCAGTGGGCAGGGTTTTGAACATTTGCAGCAGATGGCCGGTGCCGATTGGCTGCGGTTGGCCCAGTTGCCGCTGTTCGTGCCCAGCCCGCGGGTCGCCGAGCAGGCCAAGGCCGCCGGG  >  NC_002947/239317‑239603

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: