Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,560,063 T→G 5.0% D162A (GAC→GCC)  hisZ ← ATP phosphoribosyltransferase regulatory subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,560,0630TG 5.0% 316.9 / 6.0 109D162A (GAC→GCC) hisZATP phosphoribosyltransferase regulatory subunit
Reads supporting (aligned to +/- strand):  ref base T (38/65);  new base G (4/2);  total (42/67)
Fisher's exact test for biased strand distribution p-value = 2.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.15e-01

Only 100 of 110 total aligned reads displayed.

TGCTGCCCAGGTCTTTCGGCAGGTCGGCCGTCAGCGCCTGCACTTCATCGACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAGGCTGATGACTTCGACATCGCTGGTAGGGCTGGCATCGCCGTACAGTTCGGCACCCAGCTGGATCGGGCTGCGCGAGGTGGACAAGGCACGTGGCTGA  >  NC_002947/5559914‑5560209
                                                                                                                                                     |                                                                                                                                                  
tgctgcCCAGGTCTTTCGGCAGGTCGGCCGTCAGCGCCTGCACTTCATCGACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGt                                                                                                                                                    >  2:328803/1‑150 (MQ=255)
   tgcCCAGGTCTTTCGGCAGGTCGGCCGTCAGCGCCTGCACTTCATCGACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCa                                                                                                                                                 >  1:1272838/1‑150 (MQ=255)
     cccAGGTCTTTCGGCAGGTCGGCCGTCAGCGCCTGCACTTCATCGACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACGGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCAtg                                                                                                                                               >  2:1229742/1‑150 (MQ=255)
     cccAGGTCTTTCGGCAGGTCGGCCGTCAGCGCCTGCACTTCATCGACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGTCGACATGGCCCAGGTCCAtg                                                                                                                                               <  1:676386/150‑1 (MQ=255)
     cccAGGTCTTTCGGCAGGTCGGCCGTCAGCGCCTGCACTTCATCGACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCAtg                                                                                                                                               <  1:1229742/150‑1 (MQ=255)
       cAGGTCTTTCGGCAGGTCGGCCGTCAGCGCCTGCACTTCATCGACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCAtgtg                                                                                                                                             <  1:1237588/150‑1 (MQ=255)
       cAGGTCTTTCGGCAGGTCGGCCGTCAGCGCCTGCACTTCATCGACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCAtgtg                                                                                                                                             <  1:1507892/150‑1 (MQ=255)
         ggTCTTTCGGCAGGTCGGCCGTCAGCGCCTGCACTTCATCGACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGaa                                                                                                                                           <  1:2140475/150‑1 (MQ=255)
           tCTTTCGGCAGGTCGGCCGTCAGCGCCTGCACTTCATCGACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACa                                                                                                                                         >  2:2059049/1‑150 (MQ=255)
            cTTTCGGCAGGTCGGCCGTCAGCGCCTGCACTTCATCGACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACAt                                                                                                                                        >  1:699452/1‑150 (MQ=255)
             tttCGGCAGGTCGGCCGTCAGCGCCTGCACTTCATCGACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATc                                                                                                                                       >  2:2117382/1‑150 (MQ=255)
              ttCGGCAGGTCGGCCGTCAGCGCCTGCACTTCATCGACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATcc                                                                                                                                      <  1:617084/150‑1 (MQ=255)
              ttCGGCAGGTCGGCCGTCAGCGCCTGCACTTCATCGACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGCAGGTAAGGGCCAGGCCGAGGTAGATGCCTACATGGCCCAGGTCCATGTGAACATcc                                                                                                                                      <  1:298658/150‑1 (MQ=255)
                cGGCAGGTCGGCCGTCAGCGCCTGCACTTCATCGACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCgg                                                                                                                                    <  2:2290295/150‑1 (MQ=255)
                cGGCAGGTCGGCCGTCAGCGCCTGCACTTCATCGACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCgg                                                                                                                                    <  2:173199/150‑1 (MQ=255)
                cGGCAGGTCGGCCGTCAGCGCCTGCACTTCATCGACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCgg                                                                                                                                    <  2:358718/150‑1 (MQ=255)
                  gCAGGTCGGCCGTCAGCGCCTGCACTTCATCGACAGCCTTGCGTTGCATGGCGTCGAAGAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGca                                                                                                                                  >  2:346773/1‑150 (MQ=255)
                  gCAGGTCGGCCGTCAGCGCCTGCACTTCATCGACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGca                                                                                                                                  >  2:1611064/1‑150 (MQ=255)
                      gTCGGCCGTCAGCGCCTGCACTTCATCGACGGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATc                                                                                                                              <  1:1584962/150‑1 (MQ=255)
                      gTCGGCCGTCAGCGCCTGCACTTCATCGACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATc                                                                                                                              <  1:1960516/150‑1 (MQ=255)
                         ggCCGTCAGCGCCTGCACTTCATCGACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGc                                                                                                                           <  2:2363795/150‑1 (MQ=255)
                         ggCCGTCAGCCCCTGCACTTCATCGACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTCCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGc                                                                                                                           <  2:734524/150‑1 (MQ=255)
                           ccGTCAGCGCCTGCACTTCATCGACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCa                                                                                                                         <  2:1989160/150‑1 (MQ=255)
                           ccGTCAGCGCCTGCACTTCATCGACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCa                                                                                                                         <  2:737442/150‑1 (MQ=255)
                            cGTCAGCGCCTGCACTTCATCGACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAg                                                                                                                        <  2:1817022/150‑1 (MQ=255)
                             gTCAGCGCCTGCACTTCATCGACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGt                                                                                                                       >  2:1249099/1‑150 (MQ=255)
                             gTCAGCGCCTGCACTTCATCGACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGt                                                                                                                       <  1:2292406/150‑1 (MQ=255)
                                 gcgcCTGCACTTCATCGACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCAtgt                                                                                                                                              <  2:2017279/123‑1 (MQ=255)
                                      tGCACTTCATCGACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTg                                                                                                              <  1:1338601/150‑1 (MQ=255)
                                       gCACTTCATCGACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTgg                                                                                                             <  1:582244/150‑1 (MQ=255)
                                          cTTCATCGACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGtt                                                                                                                      >  1:1937456/1‑138 (MQ=255)
                                          cTTCATCGACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGtt                                                                                                                      <  2:1937456/138‑1 (MQ=255)
                                           ttCATCGACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAg                                                                                                         <  2:2215232/150‑1 (MQ=255)
                                           ttCATCGACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAg                                                                                                         >  1:1146279/1‑150 (MQ=255)
                                           ttCATCGACAGACTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAg                                                                                                         >  2:699862/1‑150 (MQ=255)
                                             cATCGACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCa                                                                                                       >  1:2418912/1‑150 (MQ=255)
                                             cATCGACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCa                                                                                                       <  2:2340298/150‑1 (MQ=255)
                                               tCGACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATc                                                                                                     >  1:232796/1‑150 (MQ=255)
                                                cGACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCACGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCa                                                                                                    <  1:931787/150‑1 (MQ=255)
                                                 gACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCAtgt                                                                                                                                              >  2:1456568/1‑107 (MQ=255)
                                                   cAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTATCGAGCATCAGGc                                                                                                 <  1:575160/150‑1 (MQ=255)
                                                      ccTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGCCCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAGGCtga                                                                                              <  1:1541165/150‑1 (MQ=255)
                                                      ccTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAGGCtga                                                                                              >  1:2061226/1‑150 (MQ=255)
                                                      ccTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAGGCtga                                                                                              >  1:2322276/1‑150 (MQ=255)
                                                       cTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAGGCtgat                                                                                             <  1:2059049/150‑1 (MQ=255)
                                                            gTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGTAGCGTGGCGAGCATCAGGCTGATGACtt                                                                                        <  2:1249768/150‑1 (MQ=255)
                                                            gTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAGGCTGATGACtt                                                                                        >  1:1426812/1‑150 (MQ=255)
                                                             ttGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAGGCTGATGACTTc                                                                                       <  2:2329291/150‑1 (MQ=255)
                                                             ttGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAGGCTGATGACTTc                                                                                       <  2:1393355/150‑1 (MQ=255)
                                                                     ggTCGAACAGCAGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAGGCTGATGACTTCGACATCGc                                                                               >  1:184059/2‑150 (MQ=255)
                                                                      gTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAGGCTGATGACTTCGACATCGCt                                                                              <  2:2124881/150‑1 (MQ=255)
                                                                      gTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAGGATGATGACTTCGACATCGCt                                                                              <  1:682700/150‑1 (MQ=255)
                                                                      gTCGAACAGCAGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAGGCTGATGACTTCGACATCGCt                                                                              >  2:2445435/1‑150 (MQ=255)
                                                                              gctgctGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAGGCTGATGACTTCGACATCGCTGGTAGGGc                                                                      <  2:2061226/150‑1 (MQ=255)
                                                                               ctgctgTTCGACCGCGCCGGACAGACCGGCGGCAGGGGTCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAGGCTGATGACTTCGACATCGCTGGTAGNGCt                                                                     <  2:1505929/150‑1 (MQ=255)
                                                                               cagctgTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAGGCTGATGACTTCGACATCGCTGGTAGGGCt                                                                     >  1:612231/3‑150 (MQ=255)
                                                                                     ttCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAGGCTGATGACtt                                                                                        >  2:812829/1‑125 (MQ=255)
                                                                                     ttCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAGGCTGATGACtt                                                                                        <  1:812829/125‑1 (MQ=255)
                                                                                     ttCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAGGCTGATGACTTCGACATCGCTGGTAGGGCTGGCATc                                                               <  2:803453/150‑1 (MQ=255)
                                                                                         aCCNCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAGGCTGATGACTTCGTCATCGCTGGTAGGGCTGGCATCGCCg                                                           >  2:1495773/1‑150 (MQ=255)
                                                                                         aCCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAGGCTGATGACTTCGACATCGCTGGTAGGGCTGGCATCGCCg                                                           <  2:232796/150‑1 (MQ=255)
                                                                                          ccGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAGGCTGATGACTTCGACATCGCTg                                                                             >  2:1922066/1‑131 (MQ=255)
                                                                                          ccGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAGGCTGATGACTTCGACATCGCTg                                                                             <  1:1922066/131‑1 (MQ=255)
                                                                                          ccGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAGGCTGATGACTTCGACATCGCTGGTAGGGCTGGCATCGCCGt                                                          >  2:338599/1‑150 (MQ=255)
                                                                                           cgcgCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAGGCTGATGACTTCGACATCGCTGGTAGGGCTGGCATCGCCGTa                                                         <  1:1524628/150‑1 (MQ=255)
                                                                                           cgcgCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAGGCTGATGACTTCGACATCGCTGGTAGGGCTGGCATCGCCGTa                                                         <  2:1426812/150‑1 (MQ=255)
                                                                                              gcCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAGGCTGATGACTTCGACATCGCTGGTAGGGCTGGCATCGCCGTACAg                                                      >  1:1412579/1‑150 (MQ=255)
                                                                                              gcCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAGGCTGATGACTTCGACATCGCTGGTAGGGCTGGCATCGCCGTACAg                                                      >  2:365255/1‑150 (MQ=255)
                                                                                                cGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAGGCTGATGACTTCGACATCGCTGGTAGGGCTGGCATCGCCGTACAGtt                                                    <  1:1674539/150‑1 (MQ=255)
                                                                                                  gACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAGGCTGATGACTTCGACATCGCTGGTAGGGCTGGCATCGCCGTACAGTTc                                                   >  1:574498/1‑149 (MQ=255)
                                                                                                        ccGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCAt                                                                                                                                                <  1:575938/50‑1 (MQ=255)
                                                                                                        ccGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCAt                                                                                                                                                >  2:575938/1‑50 (MQ=255)
                                                                                                        ccGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAGGCTGATGACTTCGACATCGCTGGTAGGGCTGGCATCGCCGTACAGTTCGGCAccc                                            >  1:1265930/1‑150 (MQ=255)
                                                                                                          ggcggcACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTTAACTGCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAGGCTGATGACTTCGACATCGCTGGTAGGGCTGGCATCGCCGTACAGTTCGGCACCCAg                                          <  2:2276814/150‑1 (MQ=255)
                                                                                                                aCGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAGGCTGATGACTTCGACATCGCTGGTAGGGCTGGCATCGCCGTACAGTTCGGCACCCAGCTGGAt                                    <  2:295808/150‑1 (MQ=255)
                                                                                                                 cGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAGGCTGATGACTTCGACATCGCTGGTAGGGCTGGCATCGCCGTACAGTTCGGCACCCAGCTGGATc                                   <  1:1611064/150‑1 (MQ=255)
                                                                                                                    gCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAGGCTGATGACTTCGACATCGCTGGTAGGGCTGGCATCGCCGTACAGTTCGGCACCCAGCTGGATCggg                                >  1:39781/1‑150 (MQ=255)
                                                                                                                    gCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAGGCTGATGACTTCGACATCGCTGGTAGGGCTGGCATCGCCGTACAGTTCGGCACCCAGCTGGATCggg                                >  2:2105647/1‑150 (MQ=255)
                                                                                                                    gCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAGGCTGATGACTTCGACATCGCTGGTAGGGCTGGCATCGCCGTACAGTTCGGCACCCAGCTGGATCggg                                <  2:39781/150‑1 (MQ=255)
                                                                                                                      cAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAGGCTGATGACTTCGACATCGCTGGTAGGGCTGGCATCGCCGTACAGTTCGGCACCCAGCTGGATCGGGCt                              <  2:1439270/150‑1 (MQ=255)
                                                                                                                       aGGCCGCGGTAGATGTCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAGGCTGATGACTTCGACATCGCTGGTAGGGCTGGCATCGCCGTACAGTTCGGCACCCAGCTGGATCGGGCTg                             >  2:1748065/1‑150 (MQ=255)
                                                                                                                        ggCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAGGCTGATGACTTCGACATCGCTGGTAGGGCTGGCATCGCCGTACAGTTCGGCACCCAGCTGGATCGGGCTgc                            <  2:612231/150‑1 (MQ=255)
                                                                                                                        ggCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAGGCTGATGACTTCGACATCGCTGGTAGGGCTGGCATCGCCGTACAGTTCGGCACCCAGCTGGATCGGGCTgc                            <  1:699862/150‑1 (MQ=255)
                                                                                                                         gCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAGGCTGATGACTTCGACATCGCTGGTAGGGCTGGCATCGCCGTACAGTTCGGCACCCAGCTGGATCGGGCTgcg                           >  2:74660/1‑150 (MQ=255)
                                                                                                                              ggTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAGGCTGATGACTTCGACATCGCTGGTAGGGCTGGCAt                                                                <  2:415181/108‑1 (MQ=255)
                                                                                                                              ggTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAGGCTGATGACTTCGACATCGCTGGTAGGGCTGGCAt                                                                >  1:415181/1‑108 (MQ=255)
                                                                                                                              ggTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAGGCTGATGACTTCGACATCGCTGGTAGGGCTGGCATCGCCGTACAGTTCGGCACCCAGCTGGATCGGGCTGCGCGAgg                      <  1:1249099/150‑1 (MQ=255)
                                                                                                                              ggTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAGGCTGATGACTTCGACATCGCTGGTAGGGCTGGCATCGCCGTACAGTTCGGCACCCAGCTGGATCGGGCTGCGCGAgg                      <  1:2445435/150‑1 (MQ=255)
                                                                                                                                tAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAGGCTGATGACTTCGACATCGCTGGTAGGGCTGGCATCGCCGTACAGTTCGGCACCCAGCTGGATCGGGCTGCGCGAGGTg                    >  2:2215318/1‑150 (MQ=255)
                                                                                                                                     gttcaCATGGACCTGGGCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAGGCTGATGACTTCGACATCGCTGGTAGGGCTGGCAt                                                                <  1:931979/97‑1 (MQ=255)
                                                                                                                                     gttcaCATGGACCTGGGCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAGGCTGATGACTTCGACATCGCTGGTAGGGCTGGCAt                                                                >  2:931979/5‑101 (MQ=255)
                                                                                                                                       cGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAGGCTGATGACTTCGACATCGCTGGTAGGGCTGGCATCGCCGTACAGTTCGGCACCCAGCTGGATCGGGCTGCGCGAGGTGGACAAgg             <  2:20167/150‑1 (MQ=255)
                                                                                                                                        gACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAGGCTGATGACTTCGACGTCTCTGGTAGCGCTCGCATCGCCGTACAGTTCGGCACCCAGCTGGATCGGGCTGCGCGAGGTGGACAAGGc            >  2:566285/1‑150 (MQ=255)
                                                                                                                                        gACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAGGCTGATGACTTCGACATCGCTGGTAGGGCTGGCATCGCCGTACAGTTCGGCACCCAGCTGGATCGGGCTGCGCGAGGTGGACAAGGc            <  2:77259/150‑1 (MQ=255)
                                                                                                                                        caCATGGACCTGGGCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAGGCTGATGACTTCGACATCGCTGGTAGGGCTGGCATCGCCGTACAGTTCGGCACCCAGCTGGATCGGGCTGCGCGa                        <  1:2352265/137‑1 (MQ=255)
                                                                                                                                        caCATGGACCTGGGCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAGGCTGATGACTTCGACATCGCTGGTAGGGCTGGCATCGCCGTACAGTTCGGCACCCAGCTGGATCGGGCTGCGCGa                        >  2:2352265/2‑138 (MQ=255)
                                                                                                                                        caCATGGACCTGGGCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAGGCTGATGACTTCGACATCGCTGGTAGGGCTGGCATCGCCGTACAGTTCGGCACCCAGCTGGATCGGGCTGCGCGAGGTGGACAAGGc            >  1:866789/2‑150 (MQ=255)
                                                                                                                                              gCCCAGGGCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAGGCTGATGACTTCGACATCGCTGGTAGGGCTGGCATCGCCGTACAGTTCGGCACCCAGCTGGATCGGGCTGCGCGAGGTGGACACGGCACGTgg      >  1:927174/1‑150 (MQ=255)
                                                                                                                                                ccAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAGGCTGATGACTTCGACATCGCTGGTAGGGCTGGCATCGCCGTACAGTTCGGCACCCAGCTGGATCGGGCTGCGCGAGGTGGACAAGGCACGTGGCt    >  2:2222476/1‑150 (MQ=255)
                                                                                                                                                  aGGTCCATGTGAACATCCGGCACATCGGCAAGTTGCAGACTGGCGAGCATGTGTTTGTTGACTTCGACATCGCTGCTGGGGCTGGCATCGCCGTACAGTTCGGCACCCAGCTGGATCGGGCTGCGCGAGGTGGACAAGGCACGTGGCTGa  <  1:1345402/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
TGCTGCCCAGGTCTTTCGGCAGGTCGGCCGTCAGCGCCTGCACTTCATCGACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAGGCTGATGACTTCGACATCGCTGGTAGGGCTGGCATCGCCGTACAGTTCGGCACCCAGCTGGATCGGGCTGCGCGAGGTGGACAAGGCACGTGGCTGA  >  NC_002947/5559914‑5560209

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: