Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,716,768 T→A 6.3% D344V (GAC→GTC)  paaH ← 3‑hydroxyacyl‑CoA dehydrogenase PaaC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,716,7680TA 6.3% 193.8 / 2.3 63D344V (GAC→GTC) paaH3‑hydroxyacyl‑CoA dehydrogenase PaaC
Reads supporting (aligned to +/- strand):  ref base T (31/28);  new base A (2/2);  total (33/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

TGTCGCTGGTATCCGCTGACCAGCTGATGGCAATGCGCGTGGCGCTGCTGTAGTCCAGCGCAAGGTCGAGCAGCACCAGGTTGCGCAGCCCGTCCTCACGGGCACGTTGGCTGGCCAGGCGGCCGTCGGACAGGGCCAGGGTGGCGTCGCCCACCTGGATCAGGCCACTACCGGCGCGCTGGGTCACGGCAATGCCGCTTTGGCGAAGGCGCTCGACCAGCGGCTGCATCACCCCCAGGTTGCCCTCGACAACGCAGGCCTCTGCGCAGGCGGAGCTGTGCAGTTCGGCCGGCTG  >  NC_002947/3716622‑3716916
                                                                                                                                                  |                                                                                                                                                    
tGTCGCTGGTATCCGCTGACCAGCTGATGGCAATGCGCGTGGCGCTGCTGTAGTCCAGCGCAAGGTCGAGCAGCACCAGGTTGCGCAGCCCGTCCTCACGGGCACGTTGGCTGGCCAGGCGGCCGTCGGACAGGGCCAGGGTGGCGTCg                                                                                                                                                    >  1:275580/1‑149 (MQ=255)
    gCTGGTATCCGCTGACCAGCTGATGGCAATGCGCGTGGCGCTGCTGTAGTCCAGCGCAAGGTCGAGCAGCACCAGGTTGCGCAGCCCGTCCTCACGGGCACGTTGGCTGGCCAGGCGGCCGTCGGACAGGGCCAGGGTGGCGTCGCCCa                                                                                                                                                >  5:206544/1‑149 (MQ=255)
    gCTGGTATCCGCTGACCAGCTGATGGCAATGCGCGTGGCGCTGCTGTAGTCCAGCGCAAGGTCGAGCAGCACCAGGTTGCGCAGCCCGTCCTCACGGGCACGTTGGCTGGCCAGGCGGCCGTCGGACAGGGCCAGGGTGGCGTCGCCCa                                                                                                                                                >  2:198085/1‑149 (MQ=255)
    gCTGGTATCCGCTGACCAGCTGATGGCAATGCGCGTGGCGCTGCTGTAGTCCAGCGCAAGGTCGAGCAGCACCAGGTTGCGCAGCCCGTCCTCACGGGCACGTTGGCTGGCCAGGCGGCCGTCGGACAGGGCCAGGGTGGCGTCGCCCa                                                                                                                                                >  2:321536/1‑149 (MQ=255)
     cTGGTATCCGCTGACCAGCTGATGGCAATGCGCGTGGCGCTGCTGTAGTCCAGCGCAAGGTCGAGCAGCACCAGGTTGCGCAGCCCGTCCTCACGGGCACGTTGGCTGGCCAGGCGGCCGTCGGACAGGGCCAGGGTGGCGTCGCCCAc                                                                                                                                               >  7:125882/1‑149 (MQ=255)
        gTATCCGCTGACCAGCTGATGGCAATGCGCGTGGCGCTGCTGTAGTCCAGCGCAAGGTCGAGCAGCACCAGGTTGCGCAGCCCGTCCTCACGGGCACGTTGGCTGGCCAGGCGGCCGTCGGACAGGGCCAGGGTGGCGTCGCCCACCTg                                                                                                                                            <  3:205933/149‑1 (MQ=255)
             cGCTGACCAGCTGATGGCAATGCGCGTGGCGCTGCTGTAGTCCAGCGCAAGGTCGAGCAGCACCAGGTTGCGCAGCCCGTCCTCACGGGCACGTTGGCTGGCCAGGCGGCCGTCGGACAGGGCCAGGGTGGCGTCGCCCACCTGGATc                                                                                                                                        <  4:73134/148‑1 (MQ=255)
             cGCTGACCAGCTGATGGCAATGCGCGTGGCGCTGCTGTAGTCCAGCGCAAGGTCGAGCAGCACCAGGTTGCGCAGCCCGTCCTCACGGGCACGTTGGCTGGCCAGGCGGCCGTCGGACAGGGCCAGGGTGGCGTCGCCCACCTGGATc                                                                                                                                        <  2:324929/148‑1 (MQ=255)
              gCTGACCAGCTGATGGCAATGCGCGTGGCGCTGCTGTAGTCCAGCGCAAGGTCGAGCAGCACCAGGTTGCGCAGCCCGTCCTCACGGGCACGTTGGCTGGCCAGGCGGCCGTCGGACAGGGCCAGGGTGGCGTCGCCCACCTGGATCAg                                                                                                                                      >  6:64151/1‑149 (MQ=255)
                         gATGGCAATGCGCGTGGCGCTGCTGTAGTCCAGCGCAAGGTCGAGCAGCACCAGGTTGCGCAGCCCGTCCTCACGGGCACGTTGGCTGGCCAGGCGGCCGTCGGACAGGGCCAGGGTGGCGTCGCCCACCTGGATCAGGCCACTACCgg                                                                                                                           >  8:124357/1‑149 (MQ=255)
                          aTGGCAATGCGCGTGGCGCTGCTGTAGTCCAGCGCAAGGTCGAGCAGCACCAGGCTGCGCAGCCCGTCCTCACGGGCACGTTGGCTGGCCAGGCGGCCGTCGGACAGTGCCAGGGTGGCGTCGCCCACCTGGATCAGGCCACTACCGgc                                                                                                                          <  1:159694/149‑1 (MQ=255)
                             gCAATGCGCGTGGCGCTGCTGTAGTCCAGCGCAAGGTCGAGCAGCACCAGGTTGCGCAGCCCGTCCTCACGGGCACGTTGGCTGGCCAGGCGGCCGTCGGACAGGGCCAGGGTGGCGTCGCCCACCTGGATCAGGCCACTACCGgcgcg                                                                                                                       >  8:90951/1‑149 (MQ=255)
                                 tGCGCGTGGCGCTGCTGTAGTCCAGCGCAAGGTCGAGCAGCACGAGGTTGCGCAGCCCGTCCTCACGGGCACGTTGGCTGGCCAGGCGGCCGTCGGACAGGGCCAGGGTGGCGTCGCCCACCTGGATCAGGCCACTACCGGCGCGCTgg                                                                                                                   <  7:84863/149‑1 (MQ=255)
                                 tGCGCGTGGCGCTGCTGGAGTCCAGCGCAAGGTCGAGCAGCACCAGGTTGCGCAGCCCGTCCTCACGGGCACGTTGGCTGGCCAGGCGGCCGTCGGACAGGGCCAGGGTGGCGTCGCCCACCTGGATCAGGCCACTACCGGCGCGCTgg                                                                                                                   <  1:123242/149‑1 (MQ=255)
                                   cgcgTGGCGCTGCTGTAGTCCAGCGCAAGGTCGAGCAGCACCAGGTTGCGCAGCCCGTCCTCACGGGCACGTTGGCTGGCCAGGCGGCCGTCGGACAGGGCCAGGGTGGCGTCGCCCACCTGGATCAGGCCACTACCGGCGCGCTGGGt                                                                                                                 <  2:275580/149‑1 (MQ=255)
                                       tGGCGCTGCTGTAGTCCAGCGCAAGGTCGAGCAGCACCAGGTTGCGCAGCCCGTCCTCACGGGCACGTTGGCTGGCCAGGCGGCCGTCGGACAGGGCCAGGGTGGCGTCGCCCACCTGGATCAGGCCACTACCGGCGCGCTGGGTCACg                                                                                                             <  7:268809/149‑1 (MQ=255)
                                                tgTAGTCCAGCGCAAGGTCGAGCAGCACCAGGTTGCGCAGCCCGTCCTCACGGGCACGTTGGCTGGCCAGGCGGCCGTCGGACAGGGCCAGGGTGGCGTCGCCCACCTGGATc                                                                                                                                        <  2:291368/113‑1 (MQ=255)
                                                tgTAGTCCAGCGCAAGGTCGAGCAGCACCAGGTTGCGCAGCCCGTCCTCACGGGCACGTTGGCTGGCCAGGCGGCCGTCGGACAGGGCCAGGGTGGCGTCGCCCACCTGGATc                                                                                                                                        >  1:291368/1‑113 (MQ=255)
                                                  tAGTCCAGCGCAAGGTCGAGCAGCACCAGGTTGCGCAGCCCGTCCTCACGGGCACGTTGGCTGGCCAGGCGGCCGTCGGACAGGGCCAGGGTGGCGTCGCCCACCTGGATCAGGCCACTACCGGCGCGCTGGGTCACGGCAATGCCGCt                                                                                                  >  3:78555/1‑149 (MQ=255)
                                                  tAGTCCAGCGCAAGGTCGAGCAGCACCAGGTTGCGCAGCCCGTCCTCACGGGCACGTTGGCTGGCCAGGCGGCCGTCGGACAGGGCCAGGGTGGCGTCGCCCACCTGGATCAGGCCACTACCGGCGCGCTGGGTCACGGCAATGCCGCt                                                                                                  >  3:83398/1‑149 (MQ=255)
                                                    gTCCAGCGCAAGGTCGAGCAGCACCAGGTTGCGCAGCCCGTCCTCACGGGCACGTTGGCTGGCCAGGCGGCCGTCGGACAGGGCCAGGGTGGCGTCGCCCACCTGGATCAGGCCACTACCGGCGCGCTGGGTCACGGCAATGTCGCttt                                                                                                >  7:100772/1‑149 (MQ=255)
                                                     tCCAGCGCAAGGTCGAGCAGCACCAGGTTGCGCAGCCCGTCCTCACGGGCACGTTGGCTGGCCAGGCGGCCGTCGGACAGGGCCAGGGTGGCGTCGCCCACCTGGATCAGGCCACTACCGGCGCGCTGGGTCACGGCAATGCCGCttt                                                                                                <  3:30640/148‑1 (MQ=255)
                                                     tCCAGCGCAAGGTCGAGCAGCACCAGGTTGCGCAGCCCGTCCTCACGGGCACGTTGGCTGGCCAGGCGGCCGTCGGACAGGGCCAGGGTGGCGTCGCCCACCTGGATCAGGCCACTACCGGCGCGCTGGGTCACGGCAATGCCGCTTTg                                                                                               <  5:294061/149‑1 (MQ=255)
                                                     tCCAGCGCAAGGTCGAGCAGCACCAGGTTGCGCAGCCCGTCCTCACGGGCACGTTGGCTGGCCAGGCGGCCGTCGGACAGGGCCAGGGTGGCGTCGCCCACCTGGATCAGGCCACTACCGGCGCGCTGGGTCACGGCAATGCCGCTTTg                                                                                               >  5:235530/1‑149 (MQ=255)
                                                     tCCAGCGCAAGGTCGAGCAGCACCAGGTTGCGCAGCCCGTCCTCACGGGCACGTTGGCTGGCCAGGCGGCCGTCGGACAGGGCCAGGGTGGCGTCGCCCACCTGGATCAGGCCACTACCGGCGCGCTGGGTCACGGCAATGCCGCTTTg                                                                                               <  7:90951/149‑1 (MQ=255)
                                                       cAGCGCAAGGTCGAGCAGCACCAGGTTGCGCAGCCCGTCCTCACGGGCACGTTGGCTGGCCAGGCGGCCGTCGGACAGGGCCAGGGTGGCGTCGCCCACCTGGATCAGGCCACTACCGGCGCGCTGGGTCACTGCAATGCAGCTTTgg                                                                                              <  6:201793/148‑1 (MQ=255)
                                                       cAGCGCAAGGTCGAGCAGCACCAGGTTGCGCAGCCCGTCCTCACGGGCACGTTGGCTGGCCAGGCGGCCGTCGGACAGGGCCAGGGTGGCGTCGCCCACCTGGATCAGGCCACTACCGGCGCGCTGGGTCACGGCAATGCCGCTTTgg                                                                                              <  3:187705/148‑1 (MQ=255)
                                                       cAGCGCAAGGTCGAGCAGCACCAGGTTGCGCAGCCCGTCCTCACGGGCACGTTGGCTGGCCAGGCGGCCGTCGGACAGGGCCAGGGTGGCGTCGCCCACCTGGATCAGGCCACTACCGGCGCGCTGGGTCACGGCAATGCCGCTTTGGc                                                                                             >  2:65489/1‑149 (MQ=255)
                                                        aGCGCAAGGTCGAGCAGCACCAGGTTGCGCAGCCCGTCCTCACGGGCACGTTGGCTGGCCAGGCGGCCGTCGGACAGGGCCAGGGTGGCGTCGCCCACCTGGATCAGGCCACTACCGGCGCGCTGGGTCACGGCAATGCCGCTTTGGCg                                                                                            <  4:70068/149‑1 (MQ=255)
                                                            cAAGGTCGAGCAGCACCAGGTTGCGCAGCCCGTCCTCACGGGCACGTTGGCTGGCCAGGCGGCCGTCGGACAGGGCCAGGGTGGCGTCGCCCACCTGGATCAGGCCACTACCGGCGCGCTGGGTCACGGCAATGCCGCTTTGGCGAAgg                                                                                        >  3:189196/1‑149 (MQ=255)
                                                                 tCGAGCAGCACCAGGTTGCGCAGCCCGTCCTCACGGGCACGTTGGCTGGCCAGGCGGCCGTCGGACAGGGCCAGGGTGGCGTCGCCCACCTGGATCAGGCCACTACCGGCGCGCTGGGTCACGGCAATGCCGCTTTGGCGAAGGCGCTc                                                                                   <  5:64151/149‑1 (MQ=255)
                                                                    agcagcACCAGGTTGCGCAGCCCGTCCTCACGGGCACGTTGGCTGGCCAGGCGGCCGTCGGACAGGGCCAGGGTGGCGTCGCCCACCTGGATCAGGCCACTACCGGCGCGCTGGGTCACGGCAATGTCGCTTTGGCGAAGGCGCTCGAc                                                                                <  8:100772/149‑1 (MQ=255)
                                                                    agcagcACCAGGTTGCGCAGCCCGTCCTCACGGGCACGTTGGCTGGCCAGGCGGCCGTCGGACAGGGCCAGGGTGGCGTCGCCCACCTGGATCAGGCCACTACCGGCGCGCTGGGTCACGGCAATGCCGCTTTGGCGAAGGCGCTCGAc                                                                                <  6:206544/149‑1 (MQ=255)
                                                                       agcaCCAGGTTGCGCAGCCCGTCCTCACGGGCACGTTGGCTGGCCAGGCGGCCGTCGGACAGGGCCAGGGTGGCGTCGCCCACCTGGATCAGGCCACTACCGGCGCGCTGGGTCACGGCAATGCCGCTTTGGCGAAGGCGCTCGACCAg                                                                             <  2:198315/149‑1 (MQ=255)
                                                                              ggTTGCGCAGCCCGTCCTCACGGGCACGTTGGCTGGCCAGGCGGCCGTCGGACAGGGCCAGGGTGGCGTCGCCCACCTGGATCAGGCCACTACCGGCGCGCTGGGTCACGGCAATGCCGCTTTGGCGAAGGCGCTCGACCAGCGGCTGc                                                                      <  1:65489/149‑1 (MQ=255)
                                                                              ggTTGCGCAGCCCGTCCTCACGGGCACGTTGGCTGGCCAGGCGGCCGTCGGACAGGGCCAGGGTGGCGTCGCCCACCTGGATCAGGCCACTACCGGCGCGCTGGGTCACGGCAATGCCGCTTTGGCGAAGGCGCTCGACCAGCGGCTGc                                                                      <  6:235530/149‑1 (MQ=255)
                                                                                   cgcAGCCCGTCCTCACGGGCACGTTGGCTGGCCAGGCGGCCGTCGGACAGGGCCAGGGTGGCGTCGCCCACCTGGATCAGGCCACTACCGGCGCGCTGGGTCACGGCAATGCCGCTTTGGCGAAGGCGCTCGACCAGCGGCTGCATCAc                                                                 <  2:43381/149‑1 (MQ=255)
                                                                                   cgcAGCCCGTCCTCACGGGCACGTTGGCTGGCCAGGCGGCCGTCGGACAGGGCCAGGGTGGCGTCGCCCACCTGGATCAGGCCACTACCGGCGCGCTGGGTCACGGCAATGCCGCTTTGGCGAAGGCGCTCGACCAGCGGCTGCATCAc                                                                 <  1:321536/149‑1 (MQ=255)
                                                                                   cgcAGCCCGTCCTCACGGGCACGTTGGCTGGCCAGGCGGCCGTCGGACAGGGCCAGGGTGGCGTCGCCCACCTGGATCAGGCCACTACCGGCGCGCTGGGTCACGGCAATGCCGCTTTGGCGAAGGCGCTCGACCAGCGGCTGCATCAc                                                                 >  3:143207/1‑149 (MQ=255)
                                                                                    gcAGCCCGTCCTCACGGGCACGTTGGCTGGCCAGGCGGCCGTCGGACAGGGCCAGGGTGGCGTCGCCCACCTGGATCAGGCCACTACCGGCGCGCTGGGTCACGGCAATGCCGCTTTGGCGAAGGCGCTCGACCAGCGGCTGCATCAcc                                                                >  7:124156/1‑149 (MQ=255)
                                                                                    gcAGCCCGTCCTCACGGGCACGTTGGCTGGCCAGGCGGCCGTCGGACAGGGCCAGGGTGGCGTCGCCCACCTGGATCAGGCCACTACCGGCGCGCTGGGTCACGGCAATGCCGCTTTGGCGAAGGCGCTCGACCAGCGGCTGCATCAcc                                                                >  5:172840/1‑149 (MQ=255)
                                                                                    gcAGCCCGTCCTCACAGGCACGTTGGCTGGCCAGGCGGCCGTCGGACAGGGCCAGGGTGGCGTCGCCCACCTGGATCAGGCCACTACCGGCGCGCTgg                                                                                                                   >  6:54816/1‑98 (MQ=255)
                                                                                    gcAGCCCGTCCTCACAGGCACGTTGGCTGGCCAGGCGGCCGTCGGACAGGGCCAGGGTGGCGTCGCCCACCTGGATCAGGCCACTACCGGCGCGCTgg                                                                                                                   <  5:54816/98‑1 (MQ=255)
                                                                                              cTCCCGGGCACGTTGGCTGGCCAGGCGGCCGTAGGACAGGGCCAGGGTGGCGTCGACCACCTGGATCAGGCCACTCCCGGCGCGCTGGGTCACGGCAATGCCGCTTTGGCGAAGGCGCTCGGCAAGCGGCGGCATCACCCCCAGGTTGc                                                      >  4:104092/1‑149 (MQ=255)
                                                                                               tCACGGGCACGTTGGCTGGCCAGGCGGCCGTCGGACAGGGCCAGGGTGGCGTCGCCCACCTGGATCAGGCCACTACCGGCGCGCTGGGTCACGGCAATGCCGCTTTGGCGAAGGCGCTCGACCAGCGGCTGCATCACCCCCAGGTTGcc                                                     >  1:199325/1‑149 (MQ=255)
                                                                                                  cGGGCACGTTGGCTGGCCAGGCGGCCGTCGGACAGGGCCAGGGTGGCGTCGCCCACCTGGATCAGGCCACTACCGGCGCGCTGGGTCACGGCAATGCCGCTTTGGCGAAGGCGCTCGACCAGCGGCTGCATCACCCCCAGGTTGCCCTc                                                  <  2:85347/149‑1 (MQ=255)
                                                                                                   gggCACGTTGGCTGGCCAGGCGGCCGTCGGACAGGGCCAGGGTGGCGTCGCCCACCTGGATCAGGCCACTACCGGCGCGCTGGGTCACGGCAATGCCGCTTTGGCGAAGGCGCTCGACCAGCGGCTGCATCACCCCCAGGTTGCCCTc                                                  >  5:238637/1‑148 (MQ=255)
                                                                                                   gggCACGTTGGCTGGCCAGGCGGCCGTCGGACAGGGCCAGGGTGGCGTCGCCCACCTGGATCAGGCCACTACCGGCGCGCTGGGTCACGGCAATGCCGCTTTGGCGAAGGCGCTCGACCAGCGGCTGCATCACCCCCAGGTTGCCCTc                                                  >  5:243847/1‑148 (MQ=255)
                                                                                                   gggCACGTTGGCTGGCCAGGCGGCCGTCGGACAGGGCCAGGGTGGCGTCGCCCACCTGGATCAGGCCACTACCGGCGCGCTGGGTCACGGCAATGCCGCTTTGGCGAAGGCGCTCGACCAGCGGCTGCATCACCCCCAGGTTGCCCTc                                                  >  4:113987/1‑148 (MQ=255)
                                                                                                                 gcCAGGCGGCCGTCGGACAGGGCCAGGGTGGCGTCGCCCACCTGGATCAGGCCACTACCGGCGCGCTGGGTCACGGCAATGCCGCTTTGGCGAAGGCGCTCGACCAGCGGCTGCATCACCCCCAGGTTGCCCTCGACAACGCAGGCctc                                   >  1:302117/1‑149 (MQ=255)
                                                                                                                               ggACAGGGCCAGGGTGGCGTCGCCCACCTGGATCAGGCCACTACCGGCGCGCTGGGTCACGGCAATGCCGCTTTGGCGAAGGCGCTCGACCAGCGGCTGCATCACCCCCAGGTTGCCCTCGACAACGCAGGCCTCTGCGCAGGCGGAGc                     <  7:124357/149‑1 (MQ=255)
                                                                                                                                gACAGGGCCAGGGTGGCGTCGCCCACCTGGATCAGGCCACTACCGGCGCGCTGGGTCACGGCAATGCCGCTTTGGCGAAGGCGCTCGACCAGCGGCTGCATCCCCCCCAGGTTGCCCTCGCCAACGCAGGCCTCTGCGCAGGCGGAGCt                    >  8:320210/1‑149 (MQ=255)
                                                                                                                                gACAGGGCCAGGGTGGCGTCGCCCACCTGGATCAGGCCACTACCGGCGCGCTGGGTCACGGCAATGCCGCTTTGGCGAAGGCGCTCGACCAGCGGCTGCATCACCCCCAGGTTGCCCTCGACAACGCAGGc                                      <  4:85819/131‑1 (MQ=255)
                                                                                                                                gACAGGGCCAGGGTGGCGTCGCCCACCTGGATCAGGCCACTACCGGCGCGCTGGGTCACGGCAATGCCGCTTTGGCGAAGGCGCTCGACCAGCGGCTGCATCACCCCCAGGTTGCCCTCGACAACGCAGGc                                      >  3:85819/1‑131 (MQ=255)
                                                                                                                                 cctgatccAGGTGGGCGACGCCACCCTGGATCAGGCCCCTACCGGCGCGCTGGGTCACGGCAATGCCGCTTTGGCGAAGGCGCTCGACCAGCGGCTGCATCAACCCCAGGGTGCCCGCGCCAACGCAGGCCTCTGCGCAGGCGGAGCtg                   >  8:135319/7‑149 (MQ=255)
                                                                                                                                 aCAGGGCCAGGGTGGCGTCGCCCACCTGGATCAGGCCACTACCGGCGCGCTGGGTCACGGCAATGCCGCTTTGGCGAAGGCGCTCGACCAGCGGCTGCATCACCCCCAGGTTGCCCTCGACAACGCAGGCCTCTGCGCAGGCGGAGCtg                   >  8:117568/1‑149 (MQ=255)
                                                                                                                                   tgatccAGGTGGGCGACGCCACCCTGGATCAGGCCACTACCGGCGCGCTGGGTCACGGCAATGCCGCTTTGGCGAAGGCGCTCGACCAGCGGCTGCATCACCCCCAg                                                           >  1:349193/5‑107 (MQ=38)
                                                                                                                                   tgatccAGGTGGGCGACGCCACCCTGGATCAGGCCACTACCGGCGCGCTGGGTCACGGCAATGCCGCTTTGGCGAAGGCGCTCGACCAGCGGCTGCATCACCCCCAg                                                           <  2:349193/103‑1 (MQ=38)
                                                                                                                                     ggCCAGGGTGGCGTCGCCCACCTGGATCAGGCCACTACCGGCGCGCTGGGTCACGGCAATGCCGCTTTGGCGAAGGCGCTCGACCAGCGGCTGCATCACCCCCAGGTTGCCCTCGACAACGCAGGCCTCTGCGCAGGCGGAGCtgtg                 >  5:157337/1‑147 (MQ=255)
                                                                                                                                     ggCCAGGGTGGCGTCGCCCACCTGGATCAGGCCACTACCGGCGCGCTGGGTCACGGCAATGCCGCTTTGGCGAAGGCGCTCGACCAGCGGCTGCATCACCCCCAGGTTGCCCTCGACAACGCAGGCCTCTGCGCAGGCGGAGCtgt                  >  3:327273/1‑146 (MQ=255)
                                                                                                                                     ggCCAGGGTGGCGTCGCCCACCTGGATCAGGCCACTACCGGCGCGCTGGGTCACGGCAATGCCGCTTTGGCGAAGGCGCTCGACCAGCGGCTGCATCACCCCCAGGTTGCCCTCGACAACGCAGGCCTCTGCGCAGGCGGAGCtgt                  <  4:327273/146‑1 (MQ=255)
                                                                                                                                       ccAGGGTGGCGTCGCCCACCTGGATCAGGCCACTACCGGCGCGCTGGGTCACGGCAATGCCGCTTTGGCGAAGGCGCTCGACCAGCGGCTGCATCACCCCCAGGTTGCCCTCGACAACGCAGGCCTCTGCGCAGGCGGAGCTGTGCAGt             >  4:176424/1‑149 (MQ=255)
                                                                                                                                             gggCGACGCCACCCTGGATCAGGCCACTACCGGCGCGCTGGGTCACGGCAATGCCGCTTTGGCGAAGGCGCTCGACCAGCGGCTGCATCAACCCCAGGTTGCCCTCGACAACGCAGGCCTCTGCGCAGGCGGAGCTGTGCAGTTcggcc       <  7:135319/148‑1 (MQ=255)
                                                                                                                                                  tCGCCCACCTGGATCAGGCCACTACCGGCGCGCTGGGTCACGGCAATGCCGCTTTGGCGAAGGCGCTCGACCAGCGGCTGCATCACCCCCAGGTTGCCCTCGACAACGCAGGCCTCTGCGCAGGCGGAGCTGTGCAGTTCGGCCGGCTg  >  1:312962/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TGTCGCTGGTATCCGCTGACCAGCTGATGGCAATGCGCGTGGCGCTGCTGTAGTCCAGCGCAAGGTCGAGCAGCACCAGGTTGCGCAGCCCGTCCTCACGGGCACGTTGGCTGGCCAGGCGGCCGTCGGACAGGGCCAGGGTGGCGTCGCCCACCTGGATCAGGCCACTACCGGCGCGCTGGGTCACGGCAATGCCGCTTTGGCGAAGGCGCTCGACCAGCGGCTGCATCACCCCCAGGTTGCCCTCGACAACGCAGGCCTCTGCGCAGGCGGAGCTGTGCAGTTCGGCCGGCTG  >  NC_002947/3716622‑3716916

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