Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,442,663 G→C 6.0% P136R (CCG→CGG)  ghrB ← 2‑ketoaldonate reductase/hydroxypyruvate/glyoxylate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,442,6630GC 6.0% 203.4 / 5.0 67P136R (CCG→CGG) ghrB2‑ketoaldonate reductase/hydroxypyruvate/glyoxylate reductase
Reads supporting (aligned to +/- strand):  ref base G (33/30);  new base C (2/2);  total (35/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.16e-01

CTTCTCCAGCGCGGTCTTGCGGCTGTTGCCGCTGTACAGGATCGGCATGTTGAAGCCGAAACGGCCACGGCGTGCCACTGCTGCACCAATGTTGCCCATGCCGACGATGCCCAGGGTCTTGCCATGCACATCGCTGCCGAAGTGGGCCGGGCCCACGGTGGCCTGCCAGTTGCCGGCCTTGGTCCAGGCATCCAGCTCGGCGGTGCGCCGGGCGCAGCCCATGATCAGCGAGAAGCCCAGGTCGGCGGTGCTTTCGGTCAGCACGTCGGGGGTGTTGGTCAGGGCGATGCCGCGC  >  NC_002947/1442515‑1442809
                                                                                                                                                    |                                                                                                                                                  
cTTCTCCAGCGCGGTCTTGCGGCTGTTGCCGCTGTACAGGATCGGCATGTTGAAGCCGAAACGGCCACGGCGTGCCACTGCTGCACCAATGTTGCCCATGCCGACGATGCCCAGGGTCTTGCCATGCACATCGCTGCCGAAGTgggccgg                                                                                                                                                   <  2:1647397/150‑1 (MQ=255)
cTTCTCCAGCGCGGTCTTGCGGCTGTTGCCGCTGTACAGGATCGGCATGTTGAAGCCGAAACGGCCACGGCGTGCCACTGCTGCACCAATGTTGCCCATGCCGACGATGCCCAGGGTCTTGCCATGCACATCGCTGCCGAAGTgggccgg                                                                                                                                                   >  1:254004/1‑150 (MQ=255)
 ttCTCCAGCGCGGTCTTGCGGCTGTTGCCGCTGTACAGGATCGGCATGTTGAAGCCGAAACGGCCACGGCGTGCCACTGCTGCACCAATGTTGCCCATGCCGACGATGCCCAGGGTCTTGCCATGCACATCGCTGCCGAAGTgggccggg                                                                                                                                                  >  1:633773/1‑150 (MQ=255)
  tctcCAGCGCGGTCTTGCGGCTGTTGCCGCTGTACAGGATCGGCATGTTGAAGCCGAAACGGCCACGGCGTGCCACTGCTGCACCAATGTTGCCCATGCCGACGATGCCCAGGGTCTTGCCATGCACATCGCTGCCGAAGTgggccgggc                                                                                                                                                 <  2:891179/150‑1 (MQ=255)
      cAGCGCGGTCTTGCGGCTGTTGCCGCTGTACAGGATCGGCATGTTGAAGCCGAAACGGCCACGGCGTGCCACTGCTGCACCAATGTTGCCCATGCCGACGATGCCCAGGGTCTTGCCATGCACATCGCTGCCGAAGTGGGCCGGGCCCAc                                                                                                                                             >  2:300623/1‑150 (MQ=255)
       aGCGCGGTCTTGCGGCTGTTGCCGCTGTACAGGATCGGCATGTTGAAGCCGAAACGGCCACGGCGTGCCACTGCTGCACCAATGTTGCCCATGCCGACGATGCCCAGGGTCTTGCCATGCACATCGCTGCCGAAGTGGGCCGGGCCCACg                                                                                                                                            <  1:300623/150‑1 (MQ=255)
       aGCGCGGTCTTGCGGCTGTTGCCGCTGTACAGGATCGGCATGTTGAAGCCGAAACGGCCACGGCGTGCCACTGCTGCACCAATGTTGCCCATGCCGACGATGCCCAGGGTCTTGCCATGCACATCGCTGCCGAAGTGGGCCGGGCCCACg                                                                                                                                            <  1:1617094/150‑1 (MQ=255)
        gcgcGGTCTTGCGGCTGTTGCCGCTGTACAGGATCGGCATGTTGAAGCCGAAACGGCCACGGCGTGCCACTGCTGCACCAATGTTGCCCCAGGCGACTATGCACCGGGCCCTGCCATGCACATCGCTGCCGAAGTGGGCCGGGCCCACgg                                                                                                                                           <  2:67356/150‑1 (MQ=255)
         cgcgGTCTTGCGGCTGTTGACGCTGTACAGGATCGGCATGTTGAAGCCGAAACGGCCACGGCGTGCCACTGCTGCACCAATGTTGCCCATGCCGACGATGCCCAGGGTCTTGCCATGCACATCGCTGCCGAAGTGGGCCCGGCCCActtc                                                                                                                                          <  1:1177185/150‑4 (MQ=255)
            ggTCTTGCGGCTGTGGCCGCTGTACAGGATCGGCATGTTGAAGCCGAAACGGCCACGGCGTGCCACTGCTGCACCAATGTTGCCCATGCCGACGATGCGCAGGGTCTTGCCATGCACAGCGCTGCCGAAGTGGGCCGGGCCCACGGTGGc                                                                                                                                       <  2:421239/150‑1 (MQ=255)
                ttGCGGCTGTTGCCGCTGTACAGGATCGGCATGTTGAAGCCGAAACGGCCACGGCGTGCCACTGCTGCACCAATGTTGCCCATGCCGACGATGCCCAGGGTCTTGCCATGCACATCGCTGCCGAAGTGGGCCGGGCCCACGGTGGCCTGc                                                                                                                                   >  2:495987/1‑150 (MQ=255)
                   cGGCTGTTGCCGCTGTACAGGATCGGCATGTTGAAGCCGAAACGGCCACGGCGTGCCACTGCTGCACCAATGTTGCCCATGCCGACGATGCCCAGGGTCTTGCCATGCACATCGCTGCCGAAGTGGGCCGGGCCCACGGTGGCCTGCCAg                                                                                                                                >  1:1388482/1‑150 (MQ=255)
                        gTTGCCGCTGTACAGGATCGGCATGTTGAAGCCGAAACGGCCACGGCGTGCCACTGCTGCACCAATGTTGCCCATGCCGACGATGCCCAGGGTCTTGCCATGCACATCGCTGCCGAAGTGGGCCGGGCCCACGGTGGCCTGCCAGTTGcc                                                                                                                           <  1:684274/150‑1 (MQ=255)
                         ttGCCGCTGTACAGGATCGGCATGTTGAAGCCGAAACGGCCACGGCGTGCCACTGCTGCACCAATGTTGCCCATGCCGACGATGCCCAGGGTCTTGCCATGCACATCGCTGCCGAAGTGGGCCGGGCCCACGGTGGCCTGCCAGTTGCCg                                                                                                                          >  2:759500/1‑150 (MQ=255)
                            ccGCTGTACAGGATCGGCCTGTTGAAGCCGAAACGGCCACGGCGTGCCACTGCTGCACCAGTGTTGCCCCTGCCGACGAGGCCCAGGGTCTTGCCCTGCACATCGCTGCCGAAGTTGGCCGGGCCCGCGGTGGCCTGCCAGATgggaaga                                                                                                                       >  2:1065692/1‑144 (MQ=255)
                             cGCTGTACAGGATCGGCATGTTGAAGCCGAAACGGCCACGGCGTGCCACTGCTGCACCAATGTTGCCCATGCCGACGATGCCCAGGGTCTTGCCATGCACATCGCTGCCGAAGTGGGCCGGGCCCACGGTGGCCTGCCAGTTGCCGGCCt                                                                                                                      >  2:164075/1‑150 (MQ=255)
                             cGCTGTACAGGATCGGCATGTTGAAGCCGAAACGGCCACGGCGTGCCACTGCTGCACCAATGTTGCCCATGCCGACGATGCCCAGGGTCTTGCCATGCACATCGCTGCCGAAGTGGGCCGGGCCCACGGTGGCCTGCCAGTTGCCGGCCt                                                                                                                      <  2:99961/150‑1 (MQ=255)
                              gCTGTACAGGATCGGCATGTTGAAGCCGAAACGGCCACGGCGTGCCACTGCTGCACCAATGTTGCCCATGCCGACGATGCCCAGGGTCTTGCCATGCACATCGCTGCCGAAGTGGGCCGGGCCCACGGTGGCCTGCCAGTTGCCGGCCtt                                                                                                                     >  1:1284776/1‑150 (MQ=255)
                               cTGTACAGGATCGGCATGTTGAAGCCGAAACGGCCACGGCGTGCCACTGCTGCACCAATGTTGCCCATGCCGACGATGCCCAGGGTCTTGCCATGCACATCGCTGCCGAAGTGGGCCGGGCCCACGGTGGCCTGCCAGTTGCCGGCCTTg                                                                                                                    <  1:773944/150‑1 (MQ=255)
                                    cAGGATCGGCATGTTGAAGCCGAAACGGACAAGGCGTGCTACTGGGTCACCAAGGTGGCCCATGCCGACGATGAGCAGGGTCTTGCCATGCACATCGCTGCCGAAGTGGGCCGGGCCCACGGTGGCCTGCCAGTTGCCGGCCTTGGTCCa                                                                                                               <  1:1884974/150‑1 (MQ=255)
                                    cAGGAGCGGCATGTTGAAGCCGAAACGGCCACGGCGTGCCACTGCTGCACCAATGTTGCCCATGCCGACGATGCCCAGGGTCTTGCCATGCACATCGCTGCCGAAGTGGGCCGGGCCCACGGTGGCCTGCCAGTTGCCGGCCTTGGTCCa                                                                                                               <  1:1031870/150‑1 (MQ=255)
                                       gATCGGCATGTTGAAGCCGAAACGGCCACGGCGTGCCACTGCTGCACCAATGTTGCCCATGCCGACGATGCCCAGGGTCTTGCCATGCACATCGCTGCCGAAGTgggccggg                                                                                                                                                  >  1:1932308/1‑112 (MQ=255)
                                       gATCGGCATGTTGAAGCCGAAACGGCCACGGCGTGCCACTGCTGCACCAATGTTGCCCATGCCGACGATGCCCAGGGTCTTGCCATGCACATCGCTGCCGAAGTgggccggg                                                                                                                                                  <  2:1932308/112‑1 (MQ=255)
                                       gATCGGCATGTTGAAGCCGAAACGGCCACGGCGTGCCACTGCTGCACCAATGTTGCCCATGCCGACGATGCCCAGGGTCTTGCCATGCACATCGCTGCCGAAGTGGGCCGGGCCCACGGTGGCCTGCCAGTTGCCGGCCTTGGTCCAg                                                                                                              <  2:149637/148‑1 (MQ=255)
                                       gATCGGCATGTTGAAGCCGAAACGGCCACGGCGTGCCACTGCTGCACCAATGTTGCCCATGCCGACGATGCCCAGGGTCTTGCCATGCACATCGCTGCCGAAGTGGGCCGGGCCCACGGTGGCCTGCCAGTTGCCGGCCTTGGTCCAg                                                                                                              >  1:149637/1‑148 (MQ=255)
                                       gATCGGCATGTTGAAGCCGAAACGGCCACGGCGTGCCACTGCTGCACCAATGTTGCCCATGCCGACGATGCCCAGGGTCTTGCCATGCACATCGCTGCCGAAGTGGGCCGGCCCCACGGTTGCCTGCCAGTTGCCGGCCTTGGTCCAGGc                                                                                                            >  1:1655874/1‑150 (MQ=255)
                                           ggCATGTTGAAGCCGAAACGGCCACGGCGTGCCACTGCTGCACCAATGTTGCCCATGCCGACGATGCCCAGGGTCTTGCCATGCACATCGCTGCCGAAGTGGGCCGGGCCCACGGTGGCCTGCCAGTTGCCGGCCTTGGTCCAGGCATcc                                                                                                        >  1:1193939/1‑150 (MQ=255)
                                            ggaTGTTGACGTCGAAGTTGCCACGGCGTACCACTGCTGCACCAATGGTGTCCATGCCGACGATGCCCAGGGTCTTGCCATGCACATCGCTGCCGAAGTGGGCCGGGCCCACGGTGGCCTGCCAGTTGCCGGCCTTGGTCCAGGCATCCa                                                                                                       >  2:1044352/3‑150 (MQ=255)
                                            gCATGTTGAAGCCGAAACGGCCACGGCGTGCCACTGCTGCACCAATGTTGCCCATGCCGACGATGCCCAGGGTCTTGCCATGCACATCGCTGCCGAAGTGGGCCGGGCCCACGGTGGc                                                                                                                                       >  1:280784/1‑118 (MQ=255)
                                            gCATGTTGAAGCCGAAACGGCCACGGCGTGCCACTGCTGCACCAATGTTGCCCATGCCGACGATGCCCAGGGTCTTGCCATGCACATCGCTGCCGAAGTGGGCCGGGCCCACGGTGGc                                                                                                                                       <  2:280784/118‑1 (MQ=255)
                                                  tGAAGCCGAAACGGCCACGGCGTGCCACTGCTGCACCAATGTTGCCCATGCCGACGATGCCCAGGGTCTTGCCATGCACATCGCTGCCGAAGTGGGCCGGGCCCACGGTGGCCTGCCAGTTGCCGGCCTTGGTCCAGGCATCCAGCTcgg                                                                                                 <  1:164075/150‑1 (MQ=255)
                                                   gAAGCCGAAACGGCCACGGCGTGCCACTGCTGCACCAATGTTGCCCATGCCGACGATGCCCAGGGTCTTGCCATGCACATGCCGGCCGAAGTGGGCCGGGCCCACGGTGGCCTGCCAGTTGCCGGCCTTGGTCCAGGCATCCAGCTcggc                                                                                                >  1:1393045/1‑150 (MQ=255)
                                                   gAAGCCGAAACGGCCACGGCGTGCCACTGCTGCACCAATGTTGCCCATGCCGACGATGCCCAGGGTCTTGCCATGCACATCGCTGCCGAAGTGGGCCGGGCCCACGGTGGCCTGCCAGTTGCCGGCCTTGGTCCAGGCATCCAGCTcggc                                                                                                >  1:162226/1‑150 (MQ=255)
                                                     aGCCGAAACGGCCACGGCGTGCCACTGCTGCACCAATGTTGCCCATGCCGACGATGCCCAGGGTCTTGCCATGCACATCGCTGCCAAAGTGGGCCGGGCCCACGGTGGCCTGCCAGTTGCCGGCCTTGGTCCCGGCATCCAGCTcggccg                                                                                              >  1:342885/1‑148 (MQ=255)
                                                       ccGAAACGGCCACGGCGTGCCACTGCTGCACCAATGTTGCCCATGCCGACGATGCCCAGGGTCTTGCCATGCACATCGCTGCCGAAGTGGGCCGGGCCCACGGTGGCCTGCCAGTTGCCGGCCTTGGTCCa                                                                                                               >  1:442335/1‑131 (MQ=255)
                                                       ccGAAACGGCCACGGCGTGCCACTGCTGCACCAATGTTGCCCATGCCGACGATGCCCAGGGTCTTGCCATGCACATCGCTGCCGAAGTGGGCCGGGCCCACGGTGGCCTGCCAGTTGCCGGCCTTGGTCCa                                                                                                               <  2:442335/131‑1 (MQ=255)
                                                       ccGAAACGGCCACGGCGTGCCACTGCTGCACCAATGTTGCCCATGCCGACGATGCCCAGGGTCTTGCCATGCACATCGCTGCCGAAGTGGGCCGGGCCCACGGTGGCCTGCCAGTTGCCGGCCTTGGTCCAGGCATCCAGCTCGGCGGTg                                                                                            <  2:564391/150‑1 (MQ=255)
                                                             cGGCCACGGCGTGCCACTGCTGCACCAATGTTGCCCATGCCGACGATGCCCAGGGTCTTGCCATGCACATCGCTGCCGAAGTGGGCCGGGCCCACGGTGGCCTGCCAGTTGCCGGCCTTGGTCCAGGCATCCAGCTCGGCGGTGCGCCgg                                                                                      <  2:1284776/150‑1 (MQ=255)
                                                                ccACGGCGTGCCACTGCTGCACCAATGTTGCCCATGCCGACGATGCCCAGGGTCTTGCCATGCACATCGCTGCCGAAGTGGGCCGGGCCCACGGTGGCCTGCCAGTTGCCGGCCTTGGTCCAGGCATCCAGCTCGGCGGTGCGCCGGgcg                                                                                   <  1:461378/150‑1 (MQ=255)
                                                                   cGGCGTGCCACTGCTGCACCAATGTTGCCCATGCCGACGATGCCCAGCGTCTTGCCATGCACATCGCTGCCGAAGTGGGCCGGGCCCACGGTGGCCTGCCCGTTGCCGGCCTTGGTCCAGGCATCCATCTCGGCGGTGGCCCGGCCGCAg                                                                                <  2:565512/150‑1 (MQ=255)
                                                                       gTGCCACTGCTGCACCAATGTTGCCCATGCCGACGATGCCCAGGGTCTTGCCATGCACATCGCTGCCGAAGTGGGCCGGGCCCACGGTGGCCTGCCAGTTGCCGGCCTTGGTCCAGGCATCCAGCTTGGCGGTGCGCCGGGGGCCGcatc                                                                            >  2:1288043/1‑147 (MQ=255)
                                                                            aCTGCTGCACCAATGTTGCCCATGCCGACGATGCCCAGGGTCTTGCCATGCACATCGCTGCCGAAGTGGGCCGGGCCCACGGTGGCCTGCCAGTTGCCGGCCTTGGTCCAGGCATCCAGCTCGGCGGTGCGCCGGGCGCAGCCCATGATc                                                                       >  1:461407/1‑150 (MQ=255)
                                                                               gctgcACCAATGTTGCCCATGCCGACGATGCCCAGGGTCTTGCCATGCACATCGCTGCCGAAGTGGGCCGGGCCCACGGTGGCCTGCCAGTTGCCGGCCTTGGTCCAGGCATCCAGCTCGGCGGTGCGCCGGGCGCAGCCCATGATCAGc                                                                    <  1:759500/150‑1 (MQ=255)
                                                                                      cAATGTTGCCGATGCCGACGATGCCCAGGGTCTTGCCATGCACATCGCTGCCGAAGTGGGCCGGGCCCACGGTGGCCTGCCAGTTGCCGGCCTTGGTCCAGGCATCCAGCTCGGCGGTGCGCCGGGCGCAGCCCATGATCAGCGAGAAGc                                                             >  2:519032/1‑150 (MQ=255)
                                                                                      cAATGTTGCCCATGCCGACGATGCCCAGGGTCTTGCCATGCACATCGCTGCCGGAGAGGGCCGGGCCCCCGGTGGCCTGCCAGTTGCCGGCCTTGGTCCAGGCATCCAGCTCGGCGGTGCGCCGGGCGCAGCCCATGATCCGAGAGAAGc                                                             >  2:1665238/1‑150 (MQ=255)
                                                                                       aaTGTTGCCCATGCCGACGATGCCCAGGGTCTTGCCATGCACATCGCTGCCGAAGTGGGCCGGGCCCACGGTGGCCTGCCAGTTGCCGGCCTTGGTCCAGGCATCCAGCTCGGCGGTGCGCCGGGCGCAGCCCATGATCAGCGAGAAGcc                                                            <  1:41380/150‑1 (MQ=255)
                                                                                       aaTGTTGCCCATGCCGACGATGCCCAGGGTCTTGCCATGCACATCGCTGCCGAAGTGGGCCGGGCCCACGGTGGACTTTCAGTTGCCGGAATAGGTCCACGTAACGTGGTTTGCGGTTCGCCGGGCGCAGCCCATGATCAGCGAGAAGcc                                                            >  1:918499/1‑150 (MQ=255)
                                                                                        aTGTTGCCCATGCCGACGATGCCCAGGGTCTTGCCATGCACATCGCTGCCGAAGTGGGCCGGGCCCACGGTGGCCTGCCAGTTGCCGGCCTTGGTCCAGGCATCCAGCTCGGCGGTGCGCCGGGCGCAGCCCATGCTCAGCGAGAAGccc                                                           <  2:1393045/150‑1 (MQ=255)
                                                                                               ccATGCCGACGATGCCCAGGGTCTTGCCATGCACATCGCTGCCGAAGTGGGCCGGGCCCACGGTGGCCTGCCAGTTGCCGGCCTTGGTCCAGGCATCCAGCTCGGCGGTGCGCCGGGCGCAGCCCATGATCAGCGAGAAGCCCAGGTcga                                                    >  1:424252/1‑149 (MQ=255)
                                                                                                  tGCCGACGATGCCCAGGGTCTTGCCATGCACATCGCTGCCGAAGTGGGCCGGGCCCACGGTGGCCTGCCAGTTGCCGGCCTTGGTCCAGGCATCCAGCTCGGCGGTGCGCCGGGCGCAGCCCATGATCAGCGAGAAGCCCAGGTcggcgg                                                 <  2:1193939/150‑1 (MQ=255)
                                                                                                    ccGACGATGCCCAGGGTCTTGCCATGCACATCGCTGCCGAAGTGGGCCGGGCCCACGGTGGCCTGCCAGTTGCCGGCCTTGGTCCAGGCATCCAGCTCGGCGGTGCGCCGGGCGCAGCCCATGATCAGCGAGAAGCCCAGGTCGGCGGTg                                               <  2:1618300/150‑1 (MQ=255)
                                                                                                     cgacgaTGCCCAGGGTCTTGCCATGCACATCGCTGCCGAAGTGGGCCGGGCCCACGGTGGCCTGCCAGTTGCCGGCCTTGGTCCAGGCATCCAGCTCGGCGGTGCGCCGGGCGCAGCCCATGATCAGCGAGAAGCCCAGGTcggcg                                                  >  2:1490860/1‑146 (MQ=255)
                                                                                                     cgacgaTGCCCAGGGTCTTGCCATGCACATCGCTGCCGAAGTGGGCCGGGCCCACGGTGGCCTGCCAGTTGCCGGCCTTGGTCCAGGCATCCAGCTCGGCGGTGCGCCGGGCGCAGCCCATGATCAGCGAGAAGCCCAGGTcggcg                                                  <  1:1490860/146‑1 (MQ=255)
                                                                                                     cgacgaTGCCCAGGGTCTTGCCATGCACATCGCTGCCGAAGTGGGCCGGGCCCACGGTGGCCTGCCAGTTGCCGGCCTTGGTCCAGGCATCCAGCTCGGCGGTGCGCCGGGCGCAGCCCATGATCAGCGAGAAGCCCAGGTCGGCGGTGc                                              >  2:583183/1‑150 (MQ=255)
                                                                                                          aTGCCCAGGGTCTTGCCATGCACATCGCTGCCGAAGTGGGCCGGGCCCACGGTGGCCTGCCAGTTGCCGGCCTTGGTCCAGGCATCCAGCTCGGCGGTGCGCCGGGCGCAGCCCATGATCAGCGAGAAGCCCAGGTCGGCGGTGCTTTCg                                         <  1:479837/150‑1 (MQ=255)
                                                                                                              ccAGGGTCTTGCCATGCACATCGCTGCCGAAGTGGGCCGGGCCCACGGTGGCTTGCCAGTTGCCGGCCTTGGTCCAGGCATCCAGCTCGGCGGTGCGCCGGGCGCAGCCCATGATCAGCGAGAAGCCCAg                                                         >  1:721817/1‑130 (MQ=255)
                                                                                                              ccAGGGTCTTGCCATGCACATCGCTGCCGAAGTGGGCCGGGCCCACGGTGGCCTGCCAGTTGCCGGCCTTGGTCCAGGCATCCAGCTCGGCGGTGCGCCGGGCGCAGCCCATGATCAGCGAGAAGCCCAg                                                         <  2:721817/130‑1 (MQ=255)
                                                                                                              ccAGGGTCTTGCCATGCACATCGCTGCCGAAGTGGGCCGGGCCCACGGTGGCCTGCCAGTTGCCGGCCTTGGTCCAGGCATCCAGCTCGGCGGTGCGCCGGGCGCAGCCCATGATCAGCGAGAAGCCCAGGTCGGCGGTGCTTTCGGTCa                                     >  1:694466/1‑150 (MQ=255)
                                                                                                               cAGGGTCTTGCCATGCACATCGCTGCCGAAGTGGGCCGGGCCCACGGTGGCCTGCCAGTTGCCGGCCTTGGTCCAGGCATCCAGCTCGGCGGTGCGCCGGGCGCAGCCCATGATCAGCGAGAAGCCCAGGTCGGCGGTGCTTTCGGTCAg                                    >  2:1809157/1‑150 (MQ=255)
                                                                                                               cAGGGTCTTGCCATGCACATCGCTGCCGAAGTGGGCCGGGCCCACGGTGGCCTGCCAGTTGCCGGCCTTGGTCCAGGCATCCAGCTCGGCGGTGCGCCGGGCGCAGCCCATGATCAGCGAGAAGCCCAGGTCGGCGGTGCTTTCGGTCAg                                    >  2:1095230/1‑150 (MQ=255)
                                                                                                                      ttGCCATGCACATCGCTGCCGAAGTGGGCCGGGCCCACGGTGGCCTGCCAGTTGCCGGCCTTGGTCCAGGCATCCAGCTCGGCGGTGCGCCGGGCGCAGCCCATGATCAGCGAGAAg                                                              <  1:134648/117‑1 (MQ=255)
                                                                                                                      ttGCCATGCACATCGCTGCCGAAGTGGGCCGGGCCCACGGTGGCCTGCCAGTTGCCGGCCTTGGTCCAGGCATCCAGCTCGGCGGTGCGCCGGGCGCAGCCCATGATCAGCGAGAAg                                                              >  2:134648/1‑117 (MQ=255)
                                                                                                                         ccATGCACATCGCTGCCGAAGTGGGCCGGGCCCACGGTGGCCTGCCAGTTGCCGGCCTTGGTCCAGGCATCCAGCTCGGCGGTGCGCCGGGCGCAGCCCATGATCAGCGAGAAGCCCAGGTCGGCGGTGCTTTCGGTCAGCACGTCgggg                          <  2:1388482/150‑1 (MQ=255)
                                                                                                                           aTGCACATCGCTGCCGAAGTGGGCCGGGCCCACGGTGGCCTGCCAGTTGCCGGCCTTGGTCCAGGCATCCAGCTCGGCGGTGCGCCGGGCGCAGCCCATGATCAGCGAGAAGCCCAGGTCGGCGGTGCTTTCGGTCAGCACGTCGGGGgt                        >  1:531622/1‑150 (MQ=255)
                                                                                                                                    gCTGCCGAAGTGGGCCGGGCCCACGGTGGCCTGCCAGTTGCCGGCCTTGGTCCAGGCATCCAGCTCGGCGGTGCGCCGGGCGCAGCCCATGATCAGCGAGaa                                                               <  2:416858/102‑1 (MQ=255)
                                                                                                                                    gCTGCCGAAGTGGGCCGGGCCCACGGTGGCCTGCCAGTTGCCGGCCTTGGTCCAGGCATCCAGCTCGGCGGTGCGCCGGGCGCAGCCCATGATCAGCGAGaa                                                               >  1:416858/1‑102 (MQ=255)
                                                                                                                                      tGCCGAAGTGGGCCGGGCCCACGGTGGCCTGCCAGTTGCCGGCCTTGGTCCAGGCATCCAGCTCGGCGGTGCGCCGGGCGCAGCCCATGATCAGCGAGAAGCCCAGGTCGGCGGTGCtt                                            >  1:271030/1‑119 (MQ=255)
                                                                                                                                      tGCCGAAGTGGGCCGGGCCCACGGTGGCCTGCCAGTTGCCGGCCTTGGTCCAGGCATCCAGCTCGGCGGTGCGCCGGGCGCAGCCCATGATCAGCGAGAAGCCCAGGTCGGCGGTGCtt                                            <  2:271030/119‑1 (MQ=255)
                                                                                                                                      tGCCGAAGTGGGCCGGGCCCACGGTGGCCTGCCAGTTGCCGGCCTTGGTCCAGGCATCCAGCTCGGCGGTGCGCCGGGCGCAGCCCATGATCAGCGAGAAGCCCAGGTCGGCGGTGCTTTCGGTCAGCACGTCGGGGGTGTTGGTCAggg             >  1:371551/1‑150 (MQ=255)
                                                                                                                                      tGCCGAAGTGGGCCGGGCCCACGGTGGCCTGCCAGTTGCCGGCCTTGGTCCAGGCATCCAGCTCGGCGGTGCGCCGGGCGCAGCCCATGATCAGCGAGAAGCCCAGGTCGGCGGTGCTTTAGGTCAGCACGTCGGGGGTGTTGGTCAggg             <  1:1095230/150‑1 (MQ=255)
                                                                                                                                         caccgTGGGCCCGGCCCACGGTGGCCTGCCAGTTGCCGGCCTTGGTCCAGGCATCCAGCTCGGCGGTGCGCCGGGCGCAGCCCATGATCAGCGAGAAGCCCAGGTCGGCGGTGCTTTCGGTCa                                     <  2:1341661/119‑1 (MQ=255)
                                                                                                                                         caccgTGGGCCCGGCCCACGGTGGCCTGCCAGTTGCCGGCCTTGGTCCAGGCATCCAGCTCGGCGGTGCGCCGGGCGCAGCCCATGATCAGCGAGAAGCCCAGGTCGGCGGTGCTTTCGGTCa                                     >  1:1341661/5‑123 (MQ=255)
                                                                                                                                         caccgTGGGCCCGGCCCACGGTGGCCTGCCAGTTGCCGGCCTTGGTCCAGGCATCCAGCTCGGCGGTGCGCCGGGCGCAGCCCATGATCAGCGAGAAGCCCAGGTCGGCGGTGCTTTCGGTCAGCACGTCGGGGGTGTTGGTCAGGGCGa          >  2:743993/5‑150 (MQ=255)
                                                                                                                                                 gccgggccCACGGTGGCCTGCCAGTTGCCGGCCTTGGTCCAGGCATCCAGCTCGGCGGTGCGCCGGGCGCAGCCCATGATCAGCGAGAAGCCCAGGTCGGCGGTGCTTTCGGTCAGCACGTCGGGGGTGTTGGTCAGGGCGATGCCgcgc  >  2:1706701/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CTTCTCCAGCGCGGTCTTGCGGCTGTTGCCGCTGTACAGGATCGGCATGTTGAAGCCGAAACGGCCACGGCGTGCCACTGCTGCACCAATGTTGCCCATGCCGACGATGCCCAGGGTCTTGCCATGCACATCGCTGCCGAAGTGGGCCGGGCCCACGGTGGCCTGCCAGTTGCCGGCCTTGGTCCAGGCATCCAGCTCGGCGGTGCGCCGGGCGCAGCCCATGATCAGCGAGAAGCCCAGGTCGGCGGTGCTTTCGGTCAGCACGTCGGGGGTGTTGGTCAGGGCGATGCCGCGC  >  NC_002947/1442515‑1442809

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