Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,786,381 C→T 6.9% G291S (GGC→AGC)  mhpT ← 3‑(3‑hydroxy‑phenyl)propionate transporter MhpT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,786,3810CT 6.9% 336.5 / 11.4 115G291S (GGC→AGC) mhpT3‑(3‑hydroxy‑phenyl)propionate transporter MhpT
Reads supporting (aligned to +/- strand):  ref base C (56/51);  new base T (3/5);  total (59/56)
Fisher's exact test for biased strand distribution p-value = 4.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.75e-01

Only 100 of 115 total aligned reads displayed.

GCGCACCTGGGTCGAATACAGCGTCGGCGCCAGCGCGTACAGCACCAACTGGCCACCAATGGCGCAATACCCTGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACTGCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAACGATCCAGCCGCCCCACCCAGGTTGAACAGAATCTGCACGGCGCCCGCCTGCGGCCGGCTGAAGCCTTGGCCGATCAA  >  NC_002947/3786235‑3786522
                                                                                                                                                  |                                                                                                                                             
gcgcACCTGGGTCGAATACAGCGTCGGCGCCAGCGCGTTCAGCACCAACTGGCCACCAATGGCGCAATACCCTGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACTGCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGc                                                                                                                                            >  2:1838059/1‑150 (MQ=255)
gcgcACCTGGGTCGAATACAGCGTCGGCGCCAGCGCGTACAGCACCAACTGGCCACCAATGGCGCAATACCCTGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACTGCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGc                                                                                                                                            >  1:2356516/1‑150 (MQ=255)
  gcACCTGGGTCGAATACAGCGTCGGCGCCAGCGCGTACAGCACCAACTGTCCACCAATGGCGCAATACCCTGCAGTGAACCCGGCCAGCAACATCAGCCAGAAACTGCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGt                                                                                                                                          >  2:58727/1‑150 (MQ=255)
   cACCTGGGTCGAATACAGCGTCGGCGCCAGCGCGTACAGCACCAACTGGCCACCAATGGCGCAATACCCTGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACTGCTCGACAGCCCCAGCCCCGCCAACGACGCCGGCATGCTGGCGtc                                                                                                                                         <  2:2061819/150‑2 (MQ=255)
   cACCTGGGTCGAATACAGCGTCGGCGCCAGCGCGTACAGCACCAACTGGCCACCAATGGCGCAATACCCTGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACTGCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTa                                                                                                                                         <  1:1105932/150‑1 (MQ=255)
    aCCTGGGTCGAATACAGCGTCGGCGCCAGCGCGTACAGCACCAACTGGCCACCAATGGCGCAATACCCTGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACTGCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAg                                                                                                                                        <  1:1468289/150‑1 (MQ=255)
    aCCTGGGTCGAATACAGCGTCGGCGCCAGCGCGTACAGCACCAACTGGCCACCAATGGCGCAATACCCTGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACTGCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAg                                                                                                                                        <  1:2266659/150‑1 (MQ=255)
    aCCTGGGTCGAATACAGCGTCGGCGCCAGCGCGTACAGCACCAACTGGCCACCAATGGCGCAATACCCTGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACTGCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAg                                                                                                                                        <  1:2230212/150‑1 (MQ=255)
    aCCTGGGTCGAATACAGCGTCGGAGCGGCAATGTACAGCACCAACTGGCCACCAATGGCGCAATACCCTGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACTGCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAg                                                                                                                                        >  1:663990/1‑150 (MQ=255)
     ccTGGGTCGAATACAGCGTCGGCGCCAGCGCGTACAGCACCAACTGGCCACCAATGGCGCAATACCCTGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACTGCTCGACAGCCCCAGCCCCGCCAACGACGCCGGCATGCTGGCGtcgt                                                                                                                                       <  2:560312/150‑4 (MQ=255)
        gggTCGAATACAGCGTCGGCGCCAGCGCGTACAGCACCAACTGGCCACCAATGGCGCAATACCCTGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACTGCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGa                                                                                                                                    >  2:868279/1‑150 (MQ=255)
        ggctggAAGACAGCGTCGGCGCAAGCGCGTACAGCACCAACTGGCCACCAATGGCGCAATACCCGGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACTGCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGa                                                                                                                                    <  2:544552/145‑1 (MQ=255)
           tCGAATACAGCGTCGGCGCCAGCGCGTACAGCACCAACTGGCCACCAATGGCGCAATACCCTGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACTGCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCa                                                                                                                                 >  2:1000661/1‑150 (MQ=255)
           tCGAATACAGCGTCGGCGCCAGCGCGTACAGCACCAACTGGCCACCAATGGCGCAATACCCTGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACTGCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCa                                                                                                                                 >  1:2135606/1‑150 (MQ=255)
             gAATACAGCGTCGGCGCCAGCGCGTACAGCACCAACTGGCCACCAATGGCGCAATATCCATCAGTGAACCCGGCCAGTAACATCAGCCCGAAACTGCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCTTAGTCGATCAAt                                                                                                                               <  1:982863/150‑1 (MQ=255)
              aaTACAGCGTCGGCGCCAGCGCGTACAGCACCAACTGGCCACCAATGGCGCAATACCCTGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACTGCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATa                                                                                                                              >  2:375024/1‑150 (MQ=255)
               aTACAGCGTCGGCGCCAGCGCGTACAGCACCAACTGGCCACCAATGGCGCAATACCCTGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACTGCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATAc                                                                                                                             >  1:1436538/1‑150 (MQ=255)
                 aCAGCGTCGGCGCCAGCGCGTACAGCACCAACTGGCCACCAATGGCGCAATACCCTGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACTGCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATAccg                                                                                                                           <  2:1893021/150‑1 (MQ=255)
                 aCAGCGTCGGCGCCAGCGCGTACAGCACCAACTGGCCACCAATGGCGCAATACCCTGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACTGCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATAccg                                                                                                                           >  2:2440542/1‑150 (MQ=255)
                  cAGCGTCGGCGCCAGCGCGTACAGCACCAACTGGCCACCAATGGCGCAATACCCTGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACTGCTCGACAGCCCCAGCCCCGCCCACTTCGCCAGCATGCCGGCGTAGGCGATCAATAccgc                                                                                                                          <  1:1773361/150‑1 (MQ=255)
                  cAGCGTCGGCGCCAGCGCGTACAGCACCAACTGGCCACCAATGGCGCAATACCCTGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACTGCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATAccgc                                                                                                                          <  2:2225694/150‑1 (MQ=255)
                  cAGCGTCGGCGCCAGCGCGTACAGCACCAACTGGCCACCAATGGCGCAATACCCTGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACTGCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATAccgc                                                                                                                          <  2:1682159/150‑1 (MQ=255)
                             ccAGCGCGTACAGCACCAACTGGCCACCAATGGCGCAATACCCTGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACTGCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCaa                                                                                                               >  1:761006/1‑150 (MQ=255)
                              cAGCGCGTACAGCACCAACTGGCCACCAATGGCGCAATACCCTGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACTGCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCaaa                                                                                                              >  1:1768346/1‑150 (MQ=255)
                               aGCGCGTACAGCACCAACTGGCCACCAATGGCGCAATACCCTGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACTGCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAAc                                                                                                             <  1:868279/150‑1 (MQ=255)
                                  gcgTACAGCACCCACTGGCCACCAATGGCGCAATACCCTGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACTGCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGgcgggcggtattgatcgcctacgccggcatgctgg                                                                                                          <  1:343929‑M1/150‑36 (MQ=255)
                                    gTACAGCACCAACTGGCCACCAATGGCGCAATACCCTGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACTGCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCcgcg                                                                                                        >  2:441103/1‑150 (MQ=255)
                                        aGCACCAACTGGCCACCAATGGCGCAATACCCTGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACTGCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTcc                                                                                                    <  1:2329318/150‑1 (MQ=255)
                                         gCACCAACTGGCCACCAATGGCGCAATACCCTGCAGTGAACCCGGCCAGCAACATCAGCCCGAACCTGCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCca                                                                                                   <  1:834483/150‑1 (MQ=255)
                                           aCCAACTGGCCACCAATGGCGCAATACCCTGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACTGCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCcatc                                                                                                 >  2:656734/1‑150 (MQ=255)
                                             cAACTGGCCACCAAGGGCGCAATACCCTGCATGGAACCCGGCCAGCAACATCAGCCCGAAACTGCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCcatcat                                                                                               <  2:137367/150‑1 (MQ=255)
                                                 tGGCCACCAATGGCGCAATACCCTGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACTGCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGa                                                                                           <  1:1584198/150‑1 (MQ=255)
                                                     caccaATGGCGCAATACCCTGCAGTGAACCCGGCCAGCATCATCAGCCCGAAACTGCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGgcgggcggtattgatcgcctacgccggcatgctggcgtc                                                                                                      >  1:187802‑M1/1‑96 (MQ=255)
                                                     caccaATGGCGCAATACCCTGCAGTGAACCCGGCCAGCATCATCAGCCCGAAACTGCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGgcgggcggtattgatcgcctacgccggcatgctggcgtc                                                                                                      <  2:187802‑M1/135‑40 (MQ=255)
                                                     caccaATGGCGCAATACCCTGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACTGCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAAt                                                                                                                               <  1:363007/110‑1 (MQ=255)
                                                     caccaATGGCGCAATACCCTGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACTGCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAAt                                                                                                                               >  2:363007/1‑110 (MQ=255)
                                                     caccaATGGCGCAATACCCTGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACTGCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCcatca                                                                                                <  1:1812958/141‑1 (MQ=255)
                                                      accaATGGCGCAATACCCTGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACTGCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGt                                                                                      >  1:1581508/1‑150 (MQ=255)
                                                         aaTGGCGCAATACCCTGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACTGCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCcc                                                                                                                    >  1:1818878‑M2/1‑115 (MQ=255)
                                                         aaTGGCGCAATACCCTGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACTGCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCcc                                                                                                                    <  2:1818878‑M2/117‑3 (MQ=255)
                                                         aaTGGCGCAATACCCTGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACTGCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAg                                                                                   >  1:2009559/1‑150 (MQ=255)
                                                            ggCGCAATACCCTGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACTGCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGaa                                                                                <  2:2009559/150‑1 (MQ=255)
                                                               gcAATACCCTGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACTGCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAACGa                                                                             <  2:761006/150‑1 (MQ=255)
                                                                cAATACCCTGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACTGCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAACGAt                                                                            >  1:2319409/1‑150 (MQ=255)
                                                                       cTGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACTGCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAACGATCCAgccg                                                                     >  2:1576591/1‑150 (MQ=255)
                                                                       cTGCAGTGAACCCGGCAAGCAACATCAGCCCGAAACTGCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCcg                                                                                                          >  2:627346/1‑113 (MQ=255)
                                                                         gCAGTGAACCCGGCCAGCAACATCAGCCCGAAACTGCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCCGTCAGGAACGATCCAgccgcc                                                                   <  1:55782/150‑1 (MQ=255)
                                                                            gTGAACCCGGCCAGCAACATCAGCCCGAAACTGCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGTTCCATCATCCGACCGGTCAg                                                                                   <  1:1310990/131‑1 (MQ=255)
                                                                            gTGAACCCGGCCAGCAACATCAGCCCGAAACTGCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGAc                                                                                          >  1:1624599/1‑124 (MQ=255)
                                                                            gTGAACCCGGCCAGCAACATCAGCCCGAAACTGCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGAc                                                                                          <  2:1624599/124‑1 (MQ=255)
                                                                                  ccGGCCAGCAACATCAGCCCGAAACTGCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAACGATCCAGCCGCCCCACCCAgg                                                          <  1:1838059/150‑1 (MQ=255)
                                                                                  ccGGCCAGCAACATCAGCCCGAAACTGCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAACGATCCAGCCGCCCCACCCAgg                                                          >  1:1034985/1‑150 (MQ=255)
                                                                                  ccGGCCAGCAACATCAGCCCGAAACTGCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAACGATCCAGCCGCCCCACCCAgg                                                          <  1:2355312/150‑1 (MQ=255)
                                                                                       cAGCAACATCAGCCCGAAACTGCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAACGATCCAGCCGCCCCACCCAGGTTGaa                                                     <  2:902321/150‑1 (MQ=255)
                                                                                        aGCAACATCAGCCCGAAACTGCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAACGATCCAGCCGCCCCACCCAGGTTGaa                                                     <  1:181911/149‑1 (MQ=255)
                                                                                        aGCAACATCAGCCCGAAACTGCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAACGATCCAGCCGCCCCACCCAGGTTGaa                                                     >  2:181911/1‑149 (MQ=255)
                                                                                        aGCAACATCAGCCCGAAACTGCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAACGATCCAGCCGCCCCACCCAGGTTGAAc                                                    >  1:1079615/1‑150 (MQ=255)
                                                                                         gCAACATCAGCCCGAAACTGCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAACGATCCAGCCGCCCCACCCAGGTTGAACa                                                   >  2:1673328/1‑150 (MQ=255)
                                                                                          cAACATCAGCCCGAAACTGCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAACGATCCAGCCGCCCCACCCAGGTTGAACAg                                                  >  1:1835460/1‑150 (MQ=255)
                                                                                          cAACATCAGCCCGAAACTGCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAACGATCCAGCCGCCCCACCCAGGTTGAACAg                                                  >  2:534928/1‑150 (MQ=255)
                                                                                           aaCATCAGCCCGAAACTGCTCGACAGCCCCAGCCCCGCCAACGACGCCGGCATGCTGGCGtcgt                                                                                                                                       >  1:726165/1‑61 (MQ=37)
                                                                                           aaCATCAGCCCGAAACTGCTCGACAGCCCCAGCCCCGCCAACGACGCCGGCATGCTGGCGtcgt                                                                                                                                       <  2:726165/64‑4 (MQ=37)
                                                                                           aaCATCAGCCCGAAACTGCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAACGATCCAGCCGCCCCACCCAGGTTGAACAGa                                                 <  1:12774/150‑1 (MQ=255)
                                                                                            aCATCAGCCCGAAACTGCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAACGATCCAGCCGCCCCACCCAGGTTGAACAGaa                                                >  2:1557797/1‑150 (MQ=255)
                                                                                            aCATCAGACCGAAACTCCGCGACAGCCCCAGCCCCGCCAACGACGCCATCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGGGGTCCATCATCCGACCGGTCAGGAACGATCCAGCCGCCCCACCCAGGTTGAACAGaa                                                <  1:1951366/150‑1 (MQ=255)
                                                                                             cATCAGCCCGAAACTGCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCaa                                                                                                                                <  2:1710230/69‑1 (MQ=255)
                                                                                             cATCAGCCCGAAACTGCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCaa                                                                                                                                >  1:1710230/1‑69 (MQ=255)
                                                                                               tCAGCCCGAAACTGCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAACGATCCAGCCGCCCCACCCAGGTTGAACAGAATCt                                             >  2:739219/1‑150 (MQ=255)
                                                                                                cAGCCCGAAACTGCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAACGATCCAGCCGCCCCACCCAGGTTGAACAGAATCTg                                            >  1:1409309/1‑150 (MQ=255)
                                                                                                  gCCCGAAACTGCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAACGATCCAGCCGCCCCACCCAGGTTGAACAGAATCTGCa                                          <  2:2319409/150‑1 (MQ=255)
                                                                                                   cccGAAACTGCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAACGATcc                                                                          <  2:1201644/117‑1 (MQ=255)
                                                                                                   cccGAAACTGCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAACGATcc                                                                          >  1:1201644/1‑117 (MQ=255)
                                                                                                   cccGAAACTGCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAACGATCCAGCCGCCCCACCCAGGTTGAACAGAATCTGCAc                                         <  1:375024/150‑1 (MQ=255)
                                                                                                   cccGAAACTGCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAACGATCCAGCCGCCCCACCCAGGTTGAACAGAATCTGCAc                                         <  2:2135606/150‑1 (MQ=255)
                                                                                                     cGAAACTGCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAACGATCCAGCCGCCCCACCCAGGTTGAACAGAATCTGCACgg                                       >  2:864259/1‑150 (MQ=255)
                                                                                                        aaCTGCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAACGATCCAGCCGCCCCACCCAGGTTGAACAGAATCTGCACGgcgc                                    <  1:441103/150‑1 (MQ=255)
                                                                                                           tGCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAACGATCCAGCCGCCCCACCCAGGTTGAACAGAATCTGCACGGcgcccg                                 >  1:2129528/1‑150 (MQ=255)
                                                                                                           tGCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAACGATCCAGCCGCCCCACCCAGGTTGAACAGAATCTGCACGGcgcccg                                 >  1:666934/1‑150 (MQ=255)
                                                                                                           tGCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAACGATCCAGCCGCCCCACCCAGGTTGAACAGAATCTGCACGGcgcccg                                 >  2:396733/1‑150 (MQ=255)
                                                                                                                gACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGAGCTCGCGTCTTCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAACGATCCAGCCGCCCCACCCAGGTTGAACAGAATCTGCACGGCGCCCGCCTGc                            <  2:663990/150‑1 (MQ=255)
                                                                                                                  cagccccagcccCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAAGAACGATCCAGCCGCCCCACCCAGGTTGAACAGAATCTGCACGGCGCCCGCCTGcgg                          >  1:1345769/1‑150 (MQ=255)
                                                                                                                       ccagccccGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAACGATCCAGCCGCCCCACCCAGGTTGAACAGAATCTGCACGGCGCCCGCCTGcggccggc                     >  2:234211/1‑150 (MQ=255)
                                                                                                                        cagccccGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAACGATCCAGCCGCCCCACCCAGGTTGAACAGAATCTGCACGGCGCCCGCCTGCGGCCGGCt                    <  1:656734/150‑1 (MQ=255)
                                                                                                                            cccGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCTGCCCTCGTACCGGTCAGGAACGATCCAGCCGCCCCACCCAGGTTGAACAGAATCTGCACGGCGCCCGCCTGCGGCCGGCTGAAg                <  1:1000661/150‑1 (MQ=255)
                                                                                                                              cGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAACGATCCAGCCGCCCCACCCAGGTTGAACAGAATCTGCACGGCGCCCGCCTGCGGCCGGCTGAAGcc              >  1:2165091/1‑150 (MQ=255)
                                                                                                                              cGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAACGATCCAGCCGCCCCACCCAGGTTGAACAGAATCTGCACGGCGCCCGCCTGCGGCCGGCTGAAGcc              >  1:450496/1‑150 (MQ=255)
                                                                                                                                ccAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAACGATCCAGCCGCCCCACCCAGGTTGAACAGAATCTGCACGGCGCCCGCCTGCGGCCGGCTGAAGCCtt            <  2:1034985/150‑1 (MQ=255)
                                                                                                                                    ctacgCCGGCATGCTGGCGTCGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAACGATCCAGCCGCCCCACCCAGGTTGAACAGAATCTGCACGGCGCCCGCCTGCGGCCGGCTGAAGc               >  2:1056756/3‑143 (MQ=255)
                                                                                                                                    ctacgCCGGCATGCTGGCGTCGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAACGATCCAGCCGCCCCACCCAGGTTGAACAGAATCTGCACGGCGCCCGCCTGCGGCCGGCTGAAGc               <  1:1056756/141‑1 (MQ=255)
                                                                                                                                     gacgCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCccggcatgctggcgtcgttggcggggctggggctgtcgagcagtttcgggctgatgatgctggccgggttcactgcagggtattgcgccattggtg                      >  1:187802‑M2/1‑39 (MQ=255)
                                                                                                                                     gacgCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCccggcatgctggcgtcgttggcggggctggggctgtcgagcagtttcgggctgatgatgctggccgggttcactgcagggtattgcgccattggtg                      <  2:187802‑M2/135‑97 (MQ=255)
                                                                                                                                        gCCAGCATGCCGGCGTAGGCGGTCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAACGATCCAGCCGCCCCACCCAGGTTGAACAGAATCTGCACGGCGCCCGCCTGCGGCCGGCTGAAGCCTTGGCCGATc    >  1:1573326/1‑150 (MQ=255)
                                                                                                                                        gCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAACGATCCAGCCGCCCCACCCAGGTTGAACAGaa                                                >  1:91274/1‑106 (MQ=255)
                                                                                                                                        gCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAACGATCCAGCCGCCCCACCCAGGTTGAACAGAATCTGCACGGCGCCCGCCTGCGGCCGGCTGAAGCCTTGGCCGATc    <  2:2129528/150‑1 (MQ=255)
                                                                                                                                        gCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAACGATCCAGCCGCCCCACCCAGGTTGAACAGAATCTGCACGGCGCCCGCCTGCGGCCGGCTGAAGCCTTGGCCGATc    <  1:2135755/150‑1 (MQ=255)
                                                                                                                                        gCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAACGATCCAGCCGCCCCACCCAGGTTGAACAGAATCTGCACGGCGCCCGCCTGCGGCCGGCTGAAGCCTTGGCCGATc    <  2:1835460/150‑1 (MQ=255)
                                                                                                                                         ccAGCATGCCGGCGTAGGCGATCAATACCGCCCGCccggcatgctggcgtcgttggcggggctggggctgtcgagcagtttcgggctgatgttgctggccgggttcactgcagggtattgcgccattggtggccagtgggtgctgtacgc   <  1:343929‑M2/150‑116 (MQ=255)
                                                                                                                                          cAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAACGATCCAGCCGCCCCACCCAGGTTGAACAGAATCTGCACGGCGCCCGCCTGCGGCCGGCTGAAGCCTTGGCCGATCaa  <  1:1557797/150‑1 (MQ=255)
                                                                                                                                           aGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAACGATCCAGCCGCCCCACCCAGGTTGAACAGAATCTGCACGGCGCCCGCCTGCGGCCGGCTGAAGCCtt            <  2:2057228/139‑1 (MQ=255)
                                                                                                                                           aGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAACGATCCAGCCGCCCCACCCAGGTTGAACAGAATCTGCACGGCGCCCGCCTGCGGCCGGCTGAAGCCtt            >  1:2057228/1‑139 (MQ=255)
                                                                                                                                                  |                                                                                                                                             
GCGCACCTGGGTCGAATACAGCGTCGGCGCCAGCGCGTACAGCACCAACTGGCCACCAATGGCGCAATACCCTGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACTGCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAACGATCCAGCCGCCCCACCCAGGTTGAACAGAATCTGCACGGCGCCCGCCTGCGGCCGGCTGAAGCCTTGGCCGATCAA  >  NC_002947/3786235‑3786522

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