Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,716,159 C→T 8.0% A44T (GCC→ACC)  paaI ← hydroxyphenylacetyl‑CoA thioesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,716,1590CT 8.0% 180.5 / 3.2 62A44T (GCC→ACC) paaIhydroxyphenylacetyl‑CoA thioesterase
Reads supporting (aligned to +/- strand):  ref base C (31/26);  new base T (2/3);  total (33/29)
Fisher's exact test for biased strand distribution p-value = 6.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGCCGGGGCCAGGTAGTCGATGGTGCAGCCCAGCGCCACCGTGGCCTGGTCGTAGCTATTGCAGGCGAAGGCGAACGCCGAATCGGCCAACGCGAACAGGAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGATACCCAGGCCCTGGGTGGCCGGGTCGCGGGCGTACATGGCGTCGGCACAGGCTTGTGCCAGTTGCTCTTCAGTCATGTAGGGTGCCTCCTTTG  >  NC_002947/3716018‑3716305
                                                                                                                                             |                                                                                                                                                  
gcgcCGGGGCCAGGTAGTCGATGGTGCAGCCCAGCGCCACCGTGGCCTGGTCGTAGCTATTGCAGGCGAAGGCGAACGCCGAATCGGCCAACGCGAACAGGAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTAcc                                                                                                                                             <  4:10247/149‑1 (MQ=255)
         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                                                                                                 <  4:38388/149‑1 (MQ=255)
                                            gCCTGGTCGTAGCTATTGCAGGCGAAGGCGAACGCCGAATCGGCCAACGCGAACAGGAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCt                                                                                                 <  1:70092/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                                                                               >  3:29472/1‑149 (MQ=255)
                                                              cAGGCGAAGGCGAACGCCGAATCGGCCAACGCGAACAGGAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGATACCCAGGCCCTGGGt                                                                               >  2:293245/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                                                                                                               >  1:320304/1‑106 (MQ=255)
                                                                         gaaCGCCGAATCGGCCAACGCGAACAGGAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCg                                                                                                               <  2:320304/106‑1 (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                                                         <  6:96246/140‑1 (MQ=255)
                                                                                             cgAACAGGAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGATACCCAGGCCCTGGGTGGCCGGGTCGCGGGCGTACATg                                                         >  5:96246/1‑140 (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                          <  6:16209/139‑1 (MQ=255)
                                                                                                                             cccTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGATACCCAGGCCCTGGGTGGCCGGGTCGCGGGCGTACATGGCGTCGGCACAGGCTTGTGCCAGTTGCTCtt                          >  5:16209/1‑139 (MQ=255)
                                                                                                                              ccTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGATACCCAGGCCCTGGGTGGCCGGGTCGCGGGCGTACATGGCGTCGGCACAGGCTTGTGCCAGTTGCTCTTCAGTCATGTAg               <  3:221735/149‑1 (MQ=255)
                                                                                                                              ccTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGATACCCAGGCCCTGGGTGGCCGGGTCGCGGGCGTACATGGCGTCGGCACAGGCTTGTGCCAGTTGCTCTTCAGTCATGTAg               >  5:125124/1‑149 (MQ=255)
                                                                                                                               cTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGATACCCAGGCCCTGGGTGGCCGGGTCGCGGGCGTACATGGCGTCGGCACAAGCCTGTGccgac                                 >  6:295337/1‑127 (MQ=255)
                                                                                                                                    gcATGGCGGTACGTGCCGACATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGATACCCAGGCCCTGGGTGGCCGGGTCGCGGGCGTACATGGCGTCGGCACAGGCTTGTGCCAGTTGctc                            >  8:86893/2‑130 (MQ=255)
                                                                                                                                    gcATGGCGGTACGTGCCGACATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGATACCCAGGCCCTGGGTGGCCGGGTCGCGGGCGTACATGGCGTCGGCACAGGCTTGTGCCAGTTGctc                            >  8:91704/2‑130 (MQ=255)
                                                                                                                                    gcATGGCGGTACGTGCCGACATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGATACCCAGGCCCTGGGTGGCCGGGTCGCGGGCGTACATGGCGTCGGCACAGGCTTGTGCCAGTTGctc                            <  7:91704/129‑1 (MQ=255)
                                                                                                                                    gcATGGCGGTACGTGCCGACATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGATACCCAGGCCCTGGGTGGCCGGGTCGCGGGCGTACATGGCGTCGGCACAGGCTTGTGCCAGTTGctc                            <  7:86893/129‑1 (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)
                                                                                                                                             |                                                                                                                                                  
GCGCCGGGGCCAGGTAGTCGATGGTGCAGCCCAGCGCCACCGTGGCCTGGTCGTAGCTATTGCAGGCGAAGGCGAACGCCGAATCGGCCAACGCGAACAGGAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGATACCCAGGCCCTGGGTGGCCGGGTCGCGGGCGTACATGGCGTCGGCACAGGCTTGTGCCAGTTGCTCTTCAGTCATGTAGGGTGCCTCCTTTG  >  NC_002947/3716018‑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: ‑