Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,443,122 A→T 10.4% intergenic (‑53/+123) ghrB ← / ← PP_1262 2‑ketoaldonate reductase/hydroxypyruvate/glyoxylate reductase/LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,443,1220AT10.4% 225.8 / 15.1 77intergenic (‑53/+123)ghrB/PP_12622‑ketoaldonate reductase/hydroxypyruvate/glyoxylate reductase/LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (29/40);  new base T (3/5);  total (32/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.03e-01

GGGGTTGGGCAGGATCACGTTGAAGTCTTGCTGCAGGCGTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGGTGCCGGAACAGGCGAAATCAATGGTTAGCAACCTACCATTGTCCGCCCCCGCATTGCAGCAGCTTTTTAGTTGGCGATCAGCACCTGATG  >  NC_002947/1442980‑1443268
                                                                                                                                               |                                                                                                                                                  
ggggTTGGGCAGGATCACGTTGAAGTCTTGCTGCAGGCGTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGc                                                                                                                                             >  2:814248/1‑150 (MQ=255)
  ggTTGGGCAGGATCACGTTGAAGTCTTGCTGCAGGCGTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTc                                                                                                                                           >  1:140692/1‑150 (MQ=255)
   gTTGGGCAGGATCACGTTGAAGTCTTGCTGCAGGCGTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCTTGCCCGcgaa                                                                                                                                          <  1:632228/150‑4 (MQ=255)
   gTTGGGCAGGATCACGTTGAAGTCTTGCTGCAGGCGTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCTTGCCCGcgaa                                                                                                                                          <  1:172369/150‑4 (MQ=255)
        gCAGGATCACGTTGAAGTCTTGCTGCAGGCGTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACgg                                                                                                                                     >  2:287072/1‑150 (MQ=255)
              tCACGTTGAAGTCTTGCTGCAGGCGTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGgc                                                                                                                               >  2:734081/1‑150 (MQ=255)
              tCACGTTGAAGTCTTGCTGCAGGCGTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGgc                                                                                                                               >  1:1398314/1‑150 (MQ=255)
              tCACGTTGAAGTCTTGCTGCAGGCGTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGGGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGgc                                                                                                                               >  2:337195/1‑150 (MQ=255)
                   ttGAAGTCTTGCTGCAGGCGTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGt                                                                                                                          <  1:814248/150‑1 (MQ=255)
                        gTCTTGCTGCAGGCGTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACggg                                                                                                                                    <  1:142929/135‑1 (MQ=255)
                        gTCTTGCTGCAGGCGTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACggg                                                                                                                                    >  2:142929/1‑135 (MQ=255)
                         tCTTGCTGCAGGCGTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCTTGCCCGcgaa                                                                                                                                          >  1:1232654/1‑125 (MQ=255)
                         tCTTGCTGCAGGCGTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCTTGCCCGcgaa                                                                                                                                          <  2:1232654/128‑4 (MQ=255)
                          cTTGCTGCAGGCGTTCGGCCATGGCCGGGGTGATACGGCTGAAGGACAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGt                                                                                                                   <  2:670049/150‑1 (MQ=255)
                           ttGCTGCAGGCGTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGTa                                                                                                                  >  2:5453/1‑150 (MQ=255)
                               tgcAGGCGTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCACTCCTACGGGGGGGCGCGGTGCCtt                                                                                                                     <  1:1159750/143‑1 (MQ=255)
                               tgcAGGCGTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCACTCCTACGGGGGGGCGCGGTGCCtt                                                                                                                     >  2:1159750/1‑143 (MQ=255)
                                gcAGGCGTTCGGCCATGGCCGGGGTGATACGGCGGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGTAGGAGc                                                                                                             <  1:1425475/150‑1 (MQ=255)
                                  aGGCGTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCCACTTCTCTGCAATGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGTAGGAGCgg                                                                                                           <  1:1240858/150‑1 (MQ=255)
                                  aGGCGTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGTAGGAGCgg                                                                                                           <  2:763226/150‑1 (MQ=255)
                                      gTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCTTGCCCGcgaag                                                                                                                                         >  2:568324/1‑112 (MQ=255)
                                      gTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCTTGCCCGcgaag                                                                                                                                         <  1:568324/116‑5 (MQ=255)
                                      gTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGTAGGAGCGGGCtt                                                                                                       >  2:1099292/1‑150 (MQ=255)
                                       ttCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCTTGCCCGcgaa                                                                                                                                          >  1:730141/1‑111 (MQ=255)
                                       ttCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCTTGCCCGcgaa                                                                                                                                          <  2:730141/114‑4 (MQ=255)
                                        tCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGTAGGAGCGGGCTTGc                                                                                                     <  1:709480/150‑1 (MQ=255)
                                          ggCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGTAGGAGCGGGCTTGccc                                                                                                   <  1:734081/150‑1 (MQ=255)
                                               tGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGTAGGAGCGGGCTTGCCCGCGaa                                                                                              <  2:140692/150‑1 (MQ=255)
                                                 gCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGCCCCTTCGCGGGCAAGCCCGCTCCTACAAGGCACCGCGCCCCCCCGTAGGAGCGGGCTTGCCCGCGAAgg                                                                                            <  2:615565/150‑1 (MQ=18)
                                                   cGGGGTGATACGGCTGAAGGCCAGGGCGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGGTCCTACGGGGGGGCGCGGTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGGTg                                                                                          <  1:889186/150‑1 (MQ=255)
                                                   cGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGGTg                                                                                          >  2:1392030/1‑150 (MQ=255)
                                                     gggTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGGTGcc                                                                                        <  1:655620/150‑1 (MQ=255)
                                                      ggTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACggg                                                                                                                                    <  1:1201664/105‑1 (MQ=255)
                                                      ggTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACggg                                                                                                                                    >  2:1201664/1‑105 (MQ=255)
                                                      ggTGATACGGCGGAGGGCAAGGAGGGTTTTTTTCAGCGGGCAACTTCTCTGTAGGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCCCCTACGGGGGGGCGCGGTGCCTTGTAGGAGCGGGCTGCCCCGCGAAGGTGCCg                                                                                       <  2:1398314/150‑1 (MQ=255)
                                                             cGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGGTGCCGGAACAgg                                                                                >  1:52388/1‑150 (MQ=255)
                                                              ggCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGGTGCCGGAACAGGc                                                                               <  2:511914/150‑1 (MQ=255)
                                                                   agggCCGGGAGGTTTTTTTTCTTCGGCCACTTTCTCGGAAGGGGGGGGGGGGGCTGGCGGGGCCCTTTCGCGGGCAAGCCCGCTCTTCGGGGGGGGCGCGGTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGGTGCCGGAACAGGCGAAAt                                                                          <  2:1023715/148‑1 (MQ=255)
                                                                      gCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGGTGCCGGAACAGGCGAAATCaa                                                                       <  1:172343/150‑1 (MQ=255)
                                                                       ccAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGGTGCCGGAACAGGCGAAATCAAt                                                                      <  1:1042639/150‑1 (MQ=255)
                                                                        cAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGGTGCCGGAACa                                                                                  <  1:894908/137‑1 (MQ=255)
                                                                        cAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGGTGCCGGAACa                                                                                  >  2:894908/1‑137 (MQ=255)
                                                                          ggACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGGTGCCGGAACAGGCGAAATCAAt                                                                      <  1:1277002/147‑1 (MQ=255)
                                                                          ggACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGGTGCCGGAACAGGCGAAATCAAt                                                                      >  2:1277002/1‑147 (MQ=255)
                                                                          ggACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGATATTGGGGTCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGGTGCCGGAACAGGCGAAATCAATGGt                                                                   <  1:1156114/150‑1 (MQ=255)
                                                                           gACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGTAGGAGCGGGCTTGCCCGCGAAg                                                                                             <  2:704669/123‑1 (MQ=255)
                                                                           gACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGTAGGAGCGGGCTTGCCCGCGAAg                                                                                             >  1:704669/1‑123 (MQ=255)
                                                                           gACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGGTGCCGGAACAGGCGAAATCAATGGtt                                                                  >  1:336480/1‑150 (MQ=255)
                                                                              ggTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCtt                                                                                                                     >  2:614234/1‑96 (MQ=255)
                                                                              ggTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCtt                                                                                                                     <  1:614234/96‑1 (MQ=255)
                                                                                 cTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGTAGGAGCGGGCTTGCCAGCGAAGGTGCCGGAACAGGCGAAATCAATGGTTAGCAAc                                                            <  1:337195/150‑1 (MQ=255)
                                                                                  tttttCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCt                                                                                                                                            >  1:375965/1‑69 (MQ=255)
                                                                                  tttttCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCt                                                                                                                                            <  2:375965/69‑1 (MQ=255)
                                                                                            ggCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGGTGCCGGAACAGGCGAAATCAATGGTTAGCAACCTACCATTGTc                                                 <  1:287072/150‑1 (MQ=255)
                                                                                                  ttCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGGTGCCGGAACAGGCGAAATCAATGGTTAGCAACCTACCATTGTCCGcccc                                           >  1:648290/1‑150 (MQ=255)
                                                                                                  ttCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGGTGCCGGAACAGGCGAAATCAATGGTTAGCAACCGACCATTGTCCGcccc                                           >  2:994295/1‑150 (MQ=255)
                                                                                                   tctcTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGGTGCCGGAACAGGCGAAATCAAt                                                                      <  2:266086/122‑1 (MQ=255)
                                                                                                   tctcTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGGTGCCGGAACAGGCGAAATCAAt                                                                      >  1:266086/1‑122 (MQ=255)
                                                                                                           aGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGGTGCCGGAACAGGCGAAATCAATGGTTAGCAACCTACCATTGTCCGCCCCCGCATTgca                                  <  2:496974/150‑1 (MQ=255)
                                                                                                           aGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGGTGCCGGAACAGGCGAAATCAATGGTTAGCAACCTACCATTGTCCGCCCCCGCATTgca                                  <  2:336480/150‑1 (MQ=255)
                                                                                                                   tGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGGTGCCGGAACAGGCGAAATCAATGGTTAGCAACCTACCATTGTCCGCCCCCGCATTGCAGCAGCttt                          <  1:994295/150‑1 (MQ=255)
                                                                                                                      ggCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGGTGCCGGAACAGGCGAAATCAATGGTTAGCAACCTACCATTGTCCGCCCCCGCATTGCAGCAGCTTTTTa                       <  1:1099292/150‑1 (MQ=255)
                                                                                                                             cTGGTCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGGTGCCGGAACAGGCGAAATCAAt                                                                      >  2:571982/1‑96 (MQ=255)
                                                                                                                             cTGGTCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGGTGCCGGAACAGGCGAAATCAAt                                                                      <  1:571982/96‑1 (MQ=255)
                                                                                                                              tGGTCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGcgcg                                                                                                                            >  2:1357364/1‑41 (MQ=255)
                                                                                                                              tGGTCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGcgcg                                                                                                                            <  1:1357364/41‑1 (MQ=255)
                                                                                                                                gTCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTg                                                                                                                    <  2:165302/47‑1 (MQ=255)
                                                                                                                                gTCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTg                                                                                                                    >  1:165302/1‑47 (MQ=255)
                                                                                                                                 tCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGGTGCCGGAACAGGCGAAATCAATGGTTAGCAACCTACCATTGTCCGCCCCCGCATTg                                    >  2:492554/1‑126 (MQ=255)
                                                                                                                                 tCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGGTGCCGGAACAGGCGAAATCAATGGTTAGCAACCTACCATTGTCCGCCCCCGCATTg                                    <  1:492554/126‑1 (MQ=255)
                                                                                                                                 tCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGGTGCCGGAACAGGCGAAATCAATGGTTAGCAACCTACCATTGTCCGCCCCCGCATTGCAGCAGCTTTTTAGTTGGCGATCa            <  1:865947/150‑1 (MQ=255)
                                                                                                                                    cTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGGTGCCGGAACCGGCGACAGCAAGGGTTAGCAACGTACGAGGTTCCTCCCCCGCATTGCAGCAGCTTTTTAGTTGGCGATCAGCa         >  2:870743/1‑150 (MQ=38)
                                                                                                                                    cTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGGTGCCGGAACAGGCGAAATCAATGGTTAGCAACCTACCATTGTCCGCCCCCGCATTGCAGCAGCTTTTTAGTTGGCGATCAGCa         <  2:648290/150‑1 (MQ=255)
                                                                                                                                       tCGCTGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGGTGCCGGAACAGGCGAAATCAATGGTTAGCAACCTACCATTGTCCGCCCCCGCATTGCAGCAGCTTTTTAGTTGGCGATCAGCAcc       >  2:1007377/1‑150 (MQ=255)
                                                                                                                                           gggCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGGTGCCGGAACAGGCGAAATCAATGGTTAGCAACCTACCATTGTCCGCCCCCGCATTGCAGCAGCTTTTTAGTTGGCGATCAGCACCTGa    >  2:1071546/1‑149 (MQ=255)
                                                                                                                                           gggCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGGTGCCGGAACAGGCGAAATCAATGGTTAGCAACCTACCATTGTCCGCCCCCGCATTGCAGCAGCTTTTTAGTTGGCGATCAGCACCTGa    <  1:1071546/149‑1 (MQ=255)
                                                                                                                                            ggCAAGCACCCTCCTACGGGGGGGCGCGGTGCCTTGTAGGAGCGGGCTTGCCAGCGAATGAGCTGCGACATGCTAAAGCAATGGTTAGCAACCTACCATTGTCCGCCCCCGCATTGCAGCAGCTTTTTAGTTGGCGATCAGCACCTGATg  <  1:1390944/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
GGGGTTGGGCAGGATCACGTTGAAGTCTTGCTGCAGGCGTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGGTGCCGGAACAGGCGAAATCAATGGTTAGCAACCTACCATTGTCCGCCCCCGCATTGCAGCAGCTTTTTAGTTGGCGATCAGCACCTGATG  >  NC_002947/1442980‑1443268

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