Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,716,168 C→A 6.4% A41S (GCG→TCG)  paaI ← hydroxyphenylacetyl‑CoA thioesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,716,1680CA 6.4% 193.2 / 3.6 63A41S (GCG→TCG) paaIhydroxyphenylacetyl‑CoA thioesterase
Reads supporting (aligned to +/- strand):  ref base C (32/27);  new base A (2/2);  total (34/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.85e-01

CCAGGTAGTCGATGGTGCAGCCCAGCGCCACCGTGGCCTGGTCGTAGCTATTGCAGGCGAAGGCGAACGCCGAATCGGCCAACGCGAACAGGAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGATACCCAGGCCCTGGGTGGCCGGGTCGCGGGCGTACATGGCGTCGGCACAGGCTTGTGCCAGTTGCTCTTCAGTCATGTAGGGTGCCTCCTTTG  >  NC_002947/3716027‑3716305
                                                                                                                                             |                                                                                                                                         
ccAGGTAGTCGATGGTGCAGCCCAGCGCCACCGTGGCCTGGTCGTAGCTATTGCAGGCGAAGGCGAACGCCGAATCGGCCAACTCGAACAGGCAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATgcgc                                                                                                                                    <  5:250778/149‑1 (MQ=255)
ccAGGTAGTCGATGGTGCAGCCCAGCGCCACCGTGGCCTGGTCGTAGCTATTGCAGGCGAAGGCGAACGCCGAATCGGCCAACGCGAACAGGAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATgcg                                                                                                                                     >  2:70092/1‑148 (MQ=255)
  aGGTAGTCGATGGTGCAGCCCAGCGCCACCGTGGCCTGGTCGTAGCTATTGCAGGCGAAGGCGAACGCCGAATCGGCCAACGCGAACAGGAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAg                                                                                                                                  >  5:71999/1‑149 (MQ=255)
      aGTCGATGGTGCAGCCCAGCGCCACCGTGGCCTGGTCGTAGCTATTGCAGGCGAAGGCGAACGCCGAATCGGCCAACGCGAACAGGAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTg                                                                                                                              >  5:135592/1‑149 (MQ=255)
       gTCGATGGTGCAGCCCAGCGCCACCGTGGCCTGGTCGTAGCTATTGCAGGCGAAGGCGAACGCCGAATCGGCCAACGCGAACAGGAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTg                                                                                                                              <  1:12993/148‑1 (MQ=255)
        tCGATGGTGCAGCCCAGCGCCACCGTGGCCTGGTCGTAGCTATTGCAGGCGAAGGCGAACGCCGAATCGGCCAACGCGAACAGGAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTggc                                                                                                                            >  3:38388/1‑149 (MQ=255)
             ggTGCAGCCCAGCGCCACCGTGGCCTGGTCGTAGCTATTGCAGGCGAAGGCGAACGCCGAATCGGCCAACGCGAACAGGAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTggcccggc                                                                                                                       >  5:222998/1‑149 (MQ=255)
                gCAGCCCAGCGCCACCGTGGCCTGGTCGTAGCTATTGCAGGCGAAGGCGAACGCCGAATCGGCCAACGCGAACAGGAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGgcccggcccg                                                                                                                    >  3:120925/1‑149 (MQ=255)
                   gCCCAGCGCCACCGTGGCCTGGTCGTAGCTATTGCAGGCGAAGGCGAACGCCGAATCGGCCAACGCGAACAGGAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGc                                                                                                                 >  5:214962/1‑149 (MQ=255)
                     ccAGCGCCACCGTGGCCTGGTCGTAGCTATTGCAGGCGAAGGCGAACGCCGAATCGGCCAACGCGAACAGGAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCg                                                                                                               <  8:214075/149‑1 (MQ=255)
                     ccAGCGCCACCGTGGCCTGGTCGTAGCTATTGCAGGCGAAGGCGAACGCCGAATCGGCCAACGCGAACAGGAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCg                                                                                                               <  3:187213/149‑1 (MQ=255)
                          gcCACCGTGGCCTGGTCGTAGCTATTGCAGGCGAAGGCGAACGCCGAATCGGCCAACGCGAACAGGAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCATc                                                                                                          <  6:214962/149‑1 (MQ=255)
                                 tGGCCTGGTCGTAGCTATTGCAGGCGAAGGCGAACGCCGAATCGGCCAACGCGAACAGGAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAgg                                                                                                   <  4:120925/149‑1 (MQ=255)
                                   gCCTGGTCGTAGCTATTGCAGGCGAAGGCGAACGCCGAATCGGCCAACGCGAACAGGAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCt                                                                                                 <  1:70092/149‑1 (MQ=255)
                                   gCCTGGTCGTAGCTATTGCAGGCGAAGGCGAACGCCGAATCGGCCAACGCGAACAGGAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCt                                                                                                 <  4:38388/149‑1 (MQ=255)
                                        gTCGTAGCTATTGCAGGCGAAGGCGAACGCCGAATCGGCCAACGCGAACAGGAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGATa                                                                                             <  1:139125/148‑1 (MQ=255)
                                           gTAGCTATTGCAGGCGAAGGCGAACGCCGAATCGGCCAACGCGAACAGGAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGATAcc                                                                                           >  3:277710/1‑147 (MQ=255)
                                                    gCATGCGAAGGCGAACGCCGAATCGGCGAACGAGAACAGGAATCCGCAATGGCAGGTGACATGCCCCTGGATCATGTCGTCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGTGCCGGCATCCAGCAGGCTTATACCCAGGCCCTgg                                                                                 <  7:166975/148‑1 (MQ=255)
                                                     cAGGCGAAGGCGAACGCCGAATCGGCCAACGCGAACAGGAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGATACCCAGGCCCTGGGt                                                                               >  2:293245/1‑149 (MQ=255)
                                                     cAGGCGAAGGCGAACGCCGAATCGGCCAACGCGAACAGGAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGATACCCAGGCCCTGGGt                                                                               >  3:29472/1‑149 (MQ=255)
                                                      aggcgaaggcgaACGCCGAATCGGCCAACGCGAACAGGAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGATACCCAGGCCCTGGGTg                                                                              <  6:71999/149‑1 (MQ=255)
                                                        gcgaaggcgaaCGCCGAATCGGCCAACGCGAACAGGAAACCGCCATGGCAGGTGCCATCCCCCTGGATCATCTCGCCACGTACCGCCATGCCCAGGCTGGCCCGGCCCCGGCCGGCATCCAGCAGGCTGATACCCAGGCCCAGGGTGGc                                                                            >  5:29439/1‑149 (MQ=255)
                                                          gaaggcgaaCGCCGAATCGGCCAACGCGAACAGGATACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGATACCCAGGCCCTGGGTGGCCg                                                                          >  1:54558/1‑149 (MQ=255)
                                                           aaggcgaaCGCCGAATCGGCCAACGCGAACAGGAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGATACCCAGGCCCTGGGTGGCCgg                                                                         >  6:157154/1‑149 (MQ=255)
                                                           aaggcgaaCGCCGAATCGGCCAACGCGAACAGGAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGATACCCAGGCCCTGGGTGGCCgg                                                                         <  2:54558/149‑1 (MQ=255)
                                                             ggcgaaCGCCGAATCGGCCAACGCGAACAGGAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGATACCCAGGCCCTGGGTGGCCGGGt                                                                       <  6:222998/149‑1 (MQ=255)
                                                                gaaCGCCGAATCGGCCAACGCGAACAGGAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCg                                                                                                               <  2:320304/106‑1 (MQ=255)
                                                                gaaCGCCGAATCGGCCAACGCGAACAGGAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCg                                                                                                               >  1:320304/1‑106 (MQ=255)
                                                                gaaCGCCGAATCGGCCAACGCGAACAGGAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGATACCCAGGCCCTGGGTGGCCGGGTcgc                                                                    >  4:221735/1‑149 (MQ=255)
                                                                             gCCAACGCGAACAGGAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGATACCCCGGCCCTGGGTGGCCGGGTCGCGGGCGTACATGGc                                                       >  2:171986/1‑149 (MQ=255)
                                                                             gCCAACGCGAACAGGAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGATACCCAGGCCCTGGGTGGCCGGGTCGCGGGCGTACATGGc                                                       >  6:211022/1‑149 (MQ=255)
                                                                               cAACGCGAACAGGAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGATACCCAGGCCCTGGGTGGCCGGGTCGCGGGCGTACATGGCGt                                                     <  5:211022/149‑1 (MQ=255)
                                                                                   gcgAACAGGAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGATACCCAGGCCCTGGGTGGCCGGGTCGCGGGCGTACATGGCGTCGGc                                                 >  7:234664/1‑149 (MQ=255)
                                                                                   gcgAACAGGAAACCGCCATGGCAGGTGCCAAGCCCCTGGATCATGTCGGCACGTACCGCCAGGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGATACCCAGGCCCTGGGTGGCCGGGTCGCGGGCGTACATGGCGTCgg                                                  >  1:30689/1‑148 (MQ=255)
                                                                                    cgAACAGGAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGATACCCAGGCCCTGGGTGGCCGGGTCGCGGGCGTACATg                                                         >  5:96246/1‑140 (MQ=255)
                                                                                    cgAACAGGAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGATACCCAGGCCCTGGGTGGCCGGGTCGCGGGCGTACATg                                                         <  6:96246/140‑1 (MQ=255)
                                                                                     gAACAGGAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGATACCCAGGCCCTGGGTGGCCGGGTCGCGGGCGTACATGGCGTCGGca                                                >  4:227644/1‑148 (MQ=255)
                                                                                         aGGAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGATACCCAGGCCCTGGGTGGCCGGGTCGCGGGCGTACATGGCGTCGGCACAGGc                                           >  7:65783/1‑149 (MQ=255)
                                                                                         aGGAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGATACCCAGGCCCTGGGTGGCCGGGTCGCGGGCGTACATGGCGTCGGCACAGGc                                           <  1:293245/149‑1 (MQ=255)
                                                                                           gAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCTGGCATCCAGCAGGCTGATACCCAGGCCCTGGGTGGCCGGGTCGCGGGCGTACAt                                                          >  2:301864/1‑132 (MQ=255)
                                                                                                cGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGATACCCAGGCCCTGGGTGGCCGGGTCGCGGGCGTACATGGCGTCGGCACAGGCTTGTGc                                     >  2:74922/1‑148 (MQ=255)
                                                                                                  ccATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGATACCCAGGCCCTGGGTGGCCGGGTCGCGGGCGTACATGGCGTCGGCACAGGCTTGTGCCAg                                  <  8:65783/149‑1 (MQ=255)
                                                                                                    aTGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGATACCCAGGCCCTGGGTGGCCGGGTCGCGGGCGTACATGGCGTCGGCACAGGCTTGTGCCAGtt                                >  4:100927/1‑149 (MQ=255)
                                                                                                     tGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGATACCCAGGCCCTGGGTGGCCGGGTCGCGGGCGTACATGGCGTCGGCACAGGCTTGTGCCAGTTg                               <  4:277710/149‑1 (MQ=255)
                                                                                                       gCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGATACCCAGGCCCTGGGTGGCCGGGTCGCGGGCGTACATGGCGTCGGCACAGGCTTGTGCCAGTTGct                             <  6:135592/149‑1 (MQ=255)
                                                                                                       gCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGATACCCAGGCCCTGGGTGGCCGGGTCGCGGGCGTACATGGCGTCGGCACAAGCCTGTGccgacgcca                             <  1:171986/149‑8 (MQ=255)
                                                                                                                  gCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGATACCCAGGCCCTGGGTGGCCGGGTCGCGGGCGTACATGGCGTCGGCACAGGCTTGTGCCAGTTGCTCTTCAGTCATg                  >  7:324392/1‑149 (MQ=255)
                                                                                                                    cccTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGATACCCAGGCCCTGGGTGGCCGGGTCGCGGGCGTACATGGCGTCGGCACAGGCTTGTGCCAGTTGCTCtt                          >  5:16209/1‑139 (MQ=255)
                                                                                                                    cccTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGATACCCAGGCCCTGGGTGGCCGGGTCGCGGGCGTACATGGCGTCGGCACAGGCTTGTGCCAGTTGCTCtt                          <  6:16209/139‑1 (MQ=255)
                                                                                                                     ccTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGATACCCAGGCCCTGGGTGGCCGGGTCGCGGGCGTACATGGCGTCGGCACAGGCTTGTGCCAGTTGCTCTTCAGTCATGTAg               >  5:125124/1‑149 (MQ=255)
                                                                                                                     ccTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGATACCCAGGCCCTGGGTGGCCGGGTCGCGGGCGTACATGGCGTCGGCACAGGCTTGTGCCAGTTGCTCTTCAGTCATGTAg               <  3:221735/149‑1 (MQ=255)
                                                                                                                      cTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGATACCCAGGCCCTGGGTGGCCGGGTCGCGGGCGTACATGGCGTCGGCACAAGCCTGTGccgac                                 >  6:295337/1‑127 (MQ=255)
                                                                                                                           gcATGGCGGTACGTGCCGACATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGATACCCAGGCCCTGGGTGGCCGGGTCGCGGGCGTACATGGCGTCGGCACAGGCTTGTGCCAGTTGctc                            <  7:86893/129‑1 (MQ=255)
                                                                                                                           gcATGGCGGTACGTGCCGACATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGATACCCAGGCCCTGGGTGGCCGGGTCGCGGGCGTACATGGCGTCGGCACAGGCTTGTGCCAGTTGctc                            <  7:91704/129‑1 (MQ=255)
                                                                                                                           gcATGGCGGTACGTGCCGACATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGATACCCAGGCCCTGGGTGGCCGGGTCGCGGGCGTACATGGCGTCGGCACAGGCTTGTGCCAGTTGctc                            >  8:86893/2‑130 (MQ=255)
                                                                                                                           gcATGGCGGTACGTGCCGACATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGATACCCAGGCCCTGGGTGGCCGGGTCGCGGGCGTACATGGCGTCGGCACAGGCTTGTGCCAGTTGctc                            >  8:91704/2‑130 (MQ=255)
                                                                                                                             aTGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGATACCCAGGCCCTGGGTGGCCGGGTCGCGGGCGTACATGGCGTCGGCACAGGCTTGTGCCAGTTGCTCTTCAGTCATGTAGGGTGcctc       >  7:34626/1‑149 (MQ=255)
                                                                                                                                tCGGTACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGATCCCCAGGCCCTGGGTGGCCTTGTCGCGGGCGTACATGGCGTCGGCAGAGGCTTGTGCCAGTTGCTCTTCGGTCATGTAGGGTGCCTCCtt    <  6:29439/149‑1 (MQ=255)
                                                                                                                                tCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGATACCCAGGCCCTGGGTGGCCGGGTCGCGGGCGTACATGGCGTCGGCACAGGCTTGTGCCAGTTGCTCTTCAGTCATGTAGGGTGCCTCCtt    >  1:288363/1‑149 (MQ=255)
                                                                                                                                  ggCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGATACCCAGGCCCTGGGTGGCCGGGTCGCGGGCGTACATGGCGTCGGCACAGGCTTGTGCCAGTTGCTCTTCAGTCATGTAGGGTGCCTCCttt   <  5:157154/148‑1 (MQ=255)
                                                                                                                                  ggCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGATACCCAGGCCCTGGGTGGCCGGGTCGCGGGCGTACATGGCGTCGGCACAGGCTTGTGCCAGTTGCTCTTCAGTCATGTAGGGTGCCTCCTTTg  <  8:234664/149‑1 (MQ=255)
                                                                                                                                         aCCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGATACCCAGGCCCTGGGTGGCCGGGTCGCGGGCGTACATGGCGTCGGCACAGGCTTGTGCCAGTTGCTCTTCAGTCAt                   <  4:66043/125‑1 (MQ=255)
                                                                                                                                         aCCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGATACCCAGGCCCTGGGTGGCCGGGTCGCGGGCGTACATGGCGTCGGCACAGGCTTGTGCCAGTTGCTCTTCAGTCAt                   >  3:66043/1‑125 (MQ=255)
                                                                                                                                             |                                                                                                                                         
CCAGGTAGTCGATGGTGCAGCCCAGCGCCACCGTGGCCTGGTCGTAGCTATTGCAGGCGAAGGCGAACGCCGAATCGGCCAACGCGAACAGGAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGATACCCAGGCCCTGGGTGGCCGGGTCGCGGGCGTACATGGCGTCGGCACAGGCTTGTGCCAGTTGCTCTTCAGTCATGTAGGGTGCCTCCTTTG  >  NC_002947/3716027‑3716305

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑