Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,716,155 T→G 6.6% D45A (GAC→GCC)  paaI ← hydroxyphenylacetyl‑CoA thioesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,716,1550TG 6.6% 180.6 / 2.6 60D45A (GAC→GCC) paaIhydroxyphenylacetyl‑CoA thioesterase
Reads supporting (aligned to +/- strand):  ref base T (31/25);  new base G (2/2);  total (33/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.90e-01

TCATCGCGCAGCGCCGGGGCCAGGTAGTCGATGGTGCAGCCCAGCGCCACCGTGGCCTGGTCGTAGCTATTGCAGGCGAAGGCGAACGCCGAATCGGCCAACGCGAACAGGAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGATACCCAGGCCCTGGGTGGCCGGGTCGCGGGCGTACATGGCGTCGGCACAGGCTTGTGCCAGTTGCTCTTCAGTCATGTAGGGTGCCTCCTT  >  NC_002947/3716008‑3716303
                                                                                                                                                   |                                                                                                                                                    
tCATCGCGCAGCGCCGGGGCCAGGTAGTCGATGGTGCAGCCCAGCGCCACCGTGGCCTGGTCGTAGCTATTGCAGGCTAAGGCGAACGCCGAATCGGCCAACGCGAACAGGAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTc                                                                                                                                                     <  8:224695/149‑1 (MQ=255)
tCATCGCGCAGCGCCGGGGCCAGGTAGTCGATGGTGCAGCCCAGCGCCACCGTGGCCTGGTCGTAGCTATTGCAGGCGAAGGCGAACGCCGAATCGGCCAACGCGAACAGGAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTc                                                                                                                                                     >  3:99610/1‑149 (MQ=255)
          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                                                                                               <  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                                                                             >  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                                                                       <  2:54558/149‑1 (MQ=255)
                                                                              aaggcgaaCGCCGAATCGGCCAACGCGAACAGGAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGATACCCAGGCCCTGGGTGGCCgg                                                                       >  6:157154/1‑149 (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                                         <  1:293245/149‑1 (MQ=255)
                                                                                                            aGGAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGATACCCAGGCCCTGGGTGGCCGGGTCGCGGGCGTACATGGCGTCGGCACAGGc                                         >  7:65783/1‑149 (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             >  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)
                                                                                                                                                   |                                                                                                                                                    
TCATCGCGCAGCGCCGGGGCCAGGTAGTCGATGGTGCAGCCCAGCGCCACCGTGGCCTGGTCGTAGCTATTGCAGGCGAAGGCGAACGCCGAATCGGCCAACGCGAACAGGAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGATACCCAGGCCCTGGGTGGCCGGGTCGCGGGCGTACATGGCGTCGGCACAGGCTTGTGCCAGTTGCTCTTCAGTCATGTAGGGTGCCTCCTT  >  NC_002947/3716008‑3716303

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: ‑