Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,419,176 G→A 5.3% R160C (CGC→TGC)  gudD ← bifunctional D‑glucarate dehydratase/L‑idarate epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,419,1760GA 5.3% 473.4 / 10.5 150R160C (CGC→TGC) gudDbifunctional D‑glucarate dehydratase/L‑idarate epimerase
Reads supporting (aligned to +/- strand):  ref base G (65/77);  new base A (4/4);  total (69/81)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 151 total aligned reads displayed.

GAGCTTGAAGTCGCTGAAGCCATAGCGCGCCTTGGCGGCCTCGGCCAGCCGCACCACTGCTTCCGGGGTCAGGGCCTTTTCGTGGCGCAGGCGAAACCAGTCGTCATCGCTGTCGGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCGGTCGCCGATGTAGAACAGATAGCCGAGCATCTTCACCGCATCGCGTTGCTGGCCCTCGCCCAGCAATGCTGCCATCGGCACATTCAGGTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAAT  >  NC_002947/5419027‑5419323
                                                                                                                                                     |                                                                                                                                                   
gAGCTTGAAGTCGCTGAAGCCATAGCGCGCCTTGGCGGCCTCGGCCAGCCGCACCACTGCTTCCGGGGTCAGGGCCTTTTCGTGGCGCAGGCGAAACCAGTCGTCATCGCTGTCAGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCg                                                                                                                                                     <  1:1252108/150‑1 (MQ=255)
    ttGAAGTCGCTGAAGCCATAGCGCGCCTTGGCGGCCTCGGCCAGCCGCACCACTGCTTCCGGGGTCAGGGCCTTTTCGTGGCGCAGGCGAAACCAGTCGTCATCGCTGTCGGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCGGTCg                                                                                                                                                 <  1:2410735/150‑1 (MQ=255)
    ttGAAGTCGCTGAAGCCATAGCGCGCCTTGGCGGCCTCGGCCAGCCGCACCACTGCTTCCGGGGTCAGGGCCTTTTCGTGGCGCAGGCGAAACCAGTCGTCATCGCTGTCGGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCGGTCg                                                                                                                                                 >  1:1701070/1‑150 (MQ=255)
    ttGAAGTCGCTGAAGCCATAGCGCGCCTTGGCGGCCTCGGCCAGCCGCACCACTGCTTCCGGGGTCAGGGCCTTTTCGTGGCGCAGGCGAAACCAGTCGTCATCGCTGTCGGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCCGGCGGTCg                                                                                                                                                 >  2:4272281/1‑150 (MQ=255)
     tGAAGTCGCTGAAGCCATAGCGCGCCTTGGCGGCCTCGGCCAGCCGCACCACTGCTTCCGGGGTCAGGGCCTTTTCGTGGCGCAGGCGAAACCAGTCGTCATCGCTGTCGGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCGGTCGc                                                                                                                                                <  2:4685790/150‑1 (MQ=255)
      gAAGTCGCTGAAGCCATAGCGCGCCTTGGCGGCCTCGGCCAGCCGCACCACTGCTTCCGGGGTCAGGGCCTTTTCGTGGCGCAGGCGAAACCAGTCGTCATCGCTGTCGGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCGGTCGcc                                                                                                                                               <  2:1458580/150‑1 (MQ=255)
        aGTCGCTGAAGCCATAGCGCGCCTTGGCGGCCTCGGCCAGCCGCACCACTGCTTCCGGGGTCAGGGCCTTTTCGTGGCGCAGGCGAAACCAGTCGTCATCGCTGTCGGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCGGTCGCCGa                                                                                                                                             <  1:2280534/150‑1 (MQ=255)
        aGTCGCTGAAGCCATAGCGCGCCTTGGCGGCCTCGGCCAGCCGCACCACTGCTTCCGGGGTCAGGGCCTTTTCGTGGCGCAGGCGAAACCAGTCGTCATCGCTGTCGGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCGGTCGCCGa                                                                                                                                             >  2:1775906/1‑150 (MQ=255)
         gTCGCTGAAGCCATAGCGCGCCTTGGCGGCCTCGGCCAGCCGCACCACTGCTTCCGGGGTCAGGGCCTTTTCGTGGCGCAGGCGAAACCAGTCGTCATCGCTGTCGGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCGGTCGCCGAt                                                                                                                                            <  1:2454791/150‑1 (MQ=255)
         gTCGCTGAAGCCATAGCGCGCCTTGGCGGCCTCGGCCAGCCGCACCACTGCTTCCGGGGTCAGGGCCTTTTCGTGGCGCAGGCGAAACCAGTCGTCATCGCTGTCGGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCGGTCGCCGAt                                                                                                                                            <  2:3626369/150‑1 (MQ=255)
         gTCGCTGAAGCCATAGCGCGCCTTGGCGGCCTCGGCCAGCCGCACCACTGCTTCCGGGGTCAGGGCCTTTTCGTGGCGCAGGCGAAACCAGTCGTCATCGCTGTCGGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCGGTCGCCGAt                                                                                                                                            <  2:3617743/150‑1 (MQ=255)
           cGCTGAAGCCATAGCGCGCCTTGGCGGCCTCGGCCAGCCGCACCACTGCTTCCGGGGTCAGGGCCTTTTCGTGGCGCAGGCGAAACCAGTCGTCATCGCTGTCGGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCGGTCGCCGATGt                                                                                                                                          <  1:821809/150‑1 (MQ=255)
            gCTGAAGCCATAGCGCGCCTTGGCGGCCTCGGCCAGCCGCACCACTGCTTCCGGGGTCAGGGCCTTTTCGTGGCGCAGGCGAAACCAGTCGTCATCGCTGTCGGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCGGTCGCCGATGTa                                                                                                                                         <  2:2846389/150‑1 (MQ=255)
            gCTGAAGCCATAGCGCGCCTTGGCGGCCTCGGCCAGCCGCACCACTGCTTCCGGGGTCAGGGCCTTTTCGTGGCGCAGGCGAAACCAGTCGTCATCGCTGTCGGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCGGTCGCCGATGTa                                                                                                                                         >  2:1523718/1‑150 (MQ=255)
            gCTGAAGCCATAGCGCGCCTTGGCGGCCTCGGCCAGCCGCACCACTGCTTCCGGGGTCAGGGCCTTTTCGTGGCGCAGGCGAAACCAGTCGTCATCGCTGTCGGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCGGTCGCCGATGTa                                                                                                                                         >  2:4970082/1‑150 (MQ=255)
            gCTGAAGCCATAGCGCGCCTTGGCGGCCTCGGCCAGCCGCACCACTGCTTCCGGGGTCAGGGCCTTTTCGTGGCGCAGGCGAAACCAGTCGTCATCGCTGTCGGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCGGTCGCCGATGTa                                                                                                                                         <  2:3149561/150‑1 (MQ=255)
            gCTGAAGCCAGAGCGCGCCTTGGCGGCCTCGGCCAGCCGCACCACTGCTTCCGGGGTCAGGGCCTTTTCGTGGCGCAGGCGAAACCAGTCGTCATCGCTGTCGGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCGGTCGCCGATGTa                                                                                                                                         <  2:1323812/150‑1 (MQ=255)
              tGAAGCCATAGCGCGCCTTGGCGGCCTCGGCCAGCCGCACCACTGCTTCCGGGGTCAGGGCCTTTTCGTGGCGCAGGCGAAACCAGTCGTCATCGCTGTCGGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCGGTCGCCGATGTAGa                                                                                                                                       >  1:1350096/1‑150 (MQ=255)
                aaGCCATAGCGCGCCTTGGCGGCCTCGGCCAGCCGCACCACTGCTTCCGGGGTCAGGGCCTTTTCGTGGCGCAGGCGAAACCAGTCGTCATCGCTGTCGGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCGGTCGCCGATGTAGAAc                                                                                                                                     <  2:4302385/150‑1 (MQ=255)
                 aGCCATGGCGCGCCTTGGCGGCCTCGGCCAGCCGCACCACTGCTTCCGGGGTCAGGGCCTTTTCGTGGCGCAGGCGAAACCAGTCGTCATCGCTGTCGGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCGGTCGCCGATGTAGAACa                                                                                                                                    <  2:1194219/150‑1 (MQ=255)
                 aGCCATAGCGCGCCTTGGCGGCCTCGGCCAGCCGCACCACTGCTTCCGGGGTCAGGGCCTTTTCGTGGCGCAGGCGAAACCAGTCGTCATCGCTGTCGGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCGGTCGCCGATGTAGAACa                                                                                                                                    <  2:4256312/150‑1 (MQ=255)
                   ccATAGCGCGCCTTGGCGGCCTCGGCCAGCCGCACCACTGCTTCCGGGGTCAGGGCCTTTTCGTGGCGCAGGCGAAACCAGTCGTCATCGCTGTCGGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCGGTCGCCGATGTAGAACAGa                                                                                                                                  >  2:2306306/1‑150 (MQ=255)
                   ccATAGCGCGCCTTGGCGGCCTCGGCCAGCCGCACCACTGCTTCCGGGGTCAGGGCCTTTTCGTGGCGCAGGCGAAACCAGTCGTCATCGCTGTCGGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCGGTCGCCGATGTAGAACAGa                                                                                                                                  >  1:2833395/1‑150 (MQ=255)
                    cATAGCGCGCCTTGGCGGCCTCGGCCAGCCGCACCACTGCTTCCGGGGTCAGGGCCTTTTCGTGGCGCAGGCGAAACCAGTCGTCATCGCTGTCGGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCGGTCGCCGATGTAGAACAGAt                                                                                                                                 >  2:4643022/1‑150 (MQ=255)
                    cATAGCGCGCCTTGGCGGCCTCGGCCAGCCGCACCACTGCTTCCGGGGTCAGGGCCTTTTCGTGGCGCAGGCGAAACCAGTCGTCATCGCTGTCGGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCGGTCGCCGATGTAGAACAGAt                                                                                                                                 >  2:4643042/1‑150 (MQ=255)
                    cATAGCGCGCCTTGGCGGCCTCGGCCAGCCGCACCACTGCTTCCGGGGTCAGGGCCTTTTCGTGGCGCAGGCGAAACCAGTCGTCATCGCTGTCGGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCGGTCGCCAATGTAGAACAGAt                                                                                                                                 <  2:1324602/150‑1 (MQ=255)
                     aTAGCGCGCCTTGGCGGCCTCGGCCAGCCGCACCACTGCTTCCGGGGTCAGGGCCTTTTCGTGGCGCAGGCGAAACCAGTCGTCATCGCTGTCGGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCGGTCGCCGATGTAGAACAGATa                                                                                                                                <  2:813770/150‑1 (MQ=255)
                     aTAGCGCGCCTTGGCGGCCTCGGCCAGCCGCACCACTGCTTCCGGGGTCAGGGCCTTTTCGTGGCGCAGGCGAAACCAGTCGTCATCGCTGTCGGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCGGTCGCCGATGTAGAACAGATa                                                                                                                                <  2:813968/150‑1 (MQ=255)
                     aTAGCGCGCCTTGGCGGCCTCGGCCAGCCGCACCACTGCTTCCGGGGTCAGGGCCTTTTCGTGGCGCAGGCGAAACCAGTCGTCATCGCTGTCGGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCGGTCGCCGATGTAGAACAGATa                                                                                                                                >  1:3942557/1‑150 (MQ=255)
                     aTAGCGCGCCTTGGCGGCCTCGGCCAGCCGCACCACTGCTTCCGGGGTCAGGGCCTTTTCGTGGCGCAGGCGAAACCAGTCGTCATCGCTGTCGGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCGGTCGCCGATGTAGAACAGATa                                                                                                                                >  1:4040393/1‑150 (MQ=255)
                     aTAGCGCGCCTTGGCGGCCTCGGCCAGCCGCACCACTGCTTCCGGGGTCAGGGCCTTTTCGTGGCGCAGGCGAAACCAGTCGTCATCGCTGTCGGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCGGTCGCCGATGTAGAACAGATa                                                                                                                                >  2:236285/1‑150 (MQ=255)
                        gcgcgcCTTGGCGGCCTCGGCCAGCCGCACCACTGCTTCCGGGGTCAGGGCCTTTTCGTGGCGCAGGCGAAACCAGTCGTCATCGCTGTCGGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCGGTCGCCGATGTAGAACAGATAGcc                                                                                                                             >  2:199241/1‑150 (MQ=255)
                             ccTTGGCGGCCTCGGCCAGCCGCACCACTGCTTCCGGGGTCAGGGCCTTTTCGTGGCGCAGGCGAAACCAGTCGTCATCGCTGTCGGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCGGTCGCCGATGTAGAACAGATAGCCGAGCa                                                                                                                        >  1:200451/1‑150 (MQ=255)
                             ccTTGGCGGCCTCGGCCAGCCGCACCACTGCTTCCGGGGTCAGGGCCTTTTCGTGGCGCAGGCGAAACCAGTCGTCATCGCTGTCGGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCGGTCGCCGATGTAGAACAGATAGCCGAGCa                                                                                                                        >  1:1354218/1‑150 (MQ=255)
                             ccTTGGCGGCCTCGGCCAGCCGCACCACTGCTTCCGGGGTCAGGGCCTTTTCGTGGCGCAGGCGAAACCAGTCGTCATCGCTGTCGGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCGGTCGCCGATGTAGAACAGATAGCCGAGCa                                                                                                                        >  1:1860111/1‑150 (MQ=255)
                              cTTGGCGGCCTCGGCCAGCCGCACCACTGCTTCCGGGGTCAGGGCCTTTTCGTGGCGCAGGCGAAACCAGTCGTCATCGCTGTCGGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCGGTCGCCGATGTAGAACAGATAGCCGAGCAt                                                                                                                       >  2:4141165/1‑150 (MQ=255)
                                 ggcggcCTCGGCCAGCCGCACCACTGCTTCCGGGGTCAGGGCCTTTTCGTGGCGCAGGCGAAACCAGTCGTCATCGCTGTCGGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCGGTCGCCGATGTAGAACAGATAGCCGAGCATCtt                                                                                                                    <  1:2394785/150‑1 (MQ=255)
                                   cggcCTCGGCCAGCCGCACCACTGCTTCCGGGGTCAGGGCCTTTTCGTGGCGCAGGCGAAACCAGTCGTCATCGCTGTCGGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCGGTCGCCGATGTAGAACAGATAGCCGAGCATCTTCa                                                                                                                  <  1:1504736/150‑1 (MQ=255)
                                    ggcCTCGGCCAGCCGCACCACTGCTTCCGGGGTCAGGGCCTTTTCGTGGCGCAGGCGAAACCAGTCGTCATCGCTGTCGGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCGGTCGCCGATGTAGAACAGATAGCCGAGCATCTTCAc                                                                                                                 <  2:1938235/150‑1 (MQ=255)
                                      ccTCGGCCAGCCGCACCACTGCTTCCGGGGTCAGGGCCTTTTCGTGGCGCAGGCGAAACCAGTCGTCATCGCTGTCGGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCGGTCGCCGATGTAGAACAGATAGCCGAGCATCTTCACCg                                                                                                               >  1:1573158/1‑150 (MQ=255)
                                        tCGGCCAGCCGCACCACTGCTTCCGGGGTCAGGGCCTTTTCGTGGCGCAGGCGAAACCAGTCGTCATCGCTGTCGGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCGGTCGCCGATGTAGAACAGATAGCCGAGCATCTTCACCGCa                                                                                                             >  2:2104354/1‑150 (MQ=255)
                                        tCGGCCAGCCGCACCACTGCTTCCGGGGTCAGGGCCTTTTCGTGGCGCAGGCGAAACCAGTCGTCATCGCTGTCGGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCGGTCGCCGATGTAGAACAGATAGCCGAGCATCTTCACCGCa                                                                                                             >  1:4616595/1‑150 (MQ=255)
                                          ggCCAGCCGCACCACTGCTTCCGGGGTCAGGGCCTTTTCGTGGCGCAGGCGAAACCAGTCGTCATCGCTGTCGGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCGGTCGCCGATGTAGAACAGATAGCCGAGCATCTTCACCGCATc                                                                                                           <  1:184808/150‑1 (MQ=255)
                                                 cGCACCACTGCTTCCGGGGTCAGGGCCTTTTCGTGGCGCAGGCGAAACCAGTCGTCATCGCTGTCGGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCGGTCGCCGATGTAGAACAGATAGCCGAGCATCTTCACCGCATCGCGTTGc                                                                                                    <  1:1919782/150‑1 (MQ=255)
                                                  gCACCACTGCTTCCGGGGTCAGGGCCTTTTCGTGGCGCAGGCGAAACCAGTCGTCATCGCTGTCGGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCGGTCGCCGATGTAGAACAGATAGCCGAGCATCTTCACCGCATCGCGTTGCt                                                                                                   >  2:915449/1‑150 (MQ=255)
                                                   caccacTGCTTCCGGGGTCAGGGCCTTTTCGTGGCGCAGGCGAAACCAGTCGTCATCGCTGTCGGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCGGTCGCCGATGTAGAACAGATAGCCGAGCATCTTCACCGCATCGCGTTGCTg                                                                                                  <  2:4296868/150‑1 (MQ=255)
                                                   caccacTGCTTCCGGGGTCAGGGCCTTTTCGTGGCGCAGGCGAAACCAGTCGTCATCGCTGTCGGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCGGTCGCCGATGTAGAACAGATAGCCGAGCATCTTCACCGCATCGCGTTGCTg                                                                                                  <  1:3607256/150‑1 (MQ=255)
                                                   caccacTGCTTCCGGGGTCAGGGCCTTTTCGTGGCGCAGGCGAAACCAGTCGTCATCGCTGTCGGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCGGTCGCCGATGTAGAACAGATAGCCGAGCATCTTCACCGCATCGCGTTGCTg                                                                                                  <  2:3644235/150‑1 (MQ=255)
                                                     ccacTGCTTCCGGGGTCAGGGCCTTTTCGTGGCGCAGGCGAAACCAGTCGTCATCGCTGTCGGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCGGTCGCCGATGTAGAACAGATAGCCGAGCATCTTCACCGCATCGCGTTGCTGGc                                                                                                >  2:1059579/1‑150 (MQ=255)
                                                      cacTGCTTCCGGGGTCAGGGCCTTTTCGTGGCGCAGGCGAAACCAGTCGTCATCGCTGTCGGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCGGTCGCCGATGTAGAACAGATAGCCGAGCATCTTCACCGCATCGCGTTGCTGGcc                                                                                               >  2:3069396/1‑150 (MQ=255)
                                                      cacTGCTTCCGGGGTCAGGGCCTTTTCGTGGCGCAGGCGAAACCAGTCGTCATCGCTGTCGGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCGGTCGCCGATGTAGAACAGATAGCCGAGCATCTTCACCGCATCGCGTTGCTGGcc                                                                                               >  2:3446200/1‑150 (MQ=255)
                                                         tGCTTCCGGGGTCAGGGCCTTTTCGTGGCGCATGCGAAACCAGTCGTCATCGCTGTCGCCGTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCGGTCGCCGATGTAGAACAGATAGCCGAGCATCTTCACCGCATCGCGTTGCTGGCCCTc                                                                                            >  2:2419786/1‑150 (MQ=255)
                                                           cTTCCGGGGTCAGGGCCTTTTCGTGGCGCAGGCGAAACCAGTCGTCATCGCTGTCGGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCGGTCGCCGATGTAGAACAGATAGCCGAGCATCTTCACCGCATCGCGTTGCTGGCCCTCGc                                                                                          >  1:4738069/1‑150 (MQ=255)
                                                            ttCCGGGGTCAGGGCCTTTTCGTGGCGCAGGCGAAACCAGTCGTCATCGCTGTCGGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCAGTCGCCGATGTAGAACAGATAGCCGAGCATCTTCACCGCATCGCGTTGCTGGCCCTCGcc                                                                                         >  2:4926746/1‑150 (MQ=255)
                                                            ttCCGGGGTCAGGGCCTTTTCGTGGCGCAGGCGAAACCAGTCGTCATCGCTGTCGGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCAGTCGCCGATGTAGAACAGATAGCCGAGCATCTTCACCGCATCGCGTTGCTGGCCCTCGcc                                                                                         >  2:4926497/1‑150 (MQ=255)
                                                            ttCCGGGGTCAGGGCCTTTTCGTGGCGCAGGCGAAACCAGTCGTCATCGCTGTCGGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCAGTCGCCGATGTAGAACAGATAGCCGAGCATCTTCACCGCATCGCGTTGCTGGCCCTCGcc                                                                                         >  2:4926461/1‑150 (MQ=255)
                                                                        ggCCTTTTCGTGGCGCAGGCGAAACCAGTCGTCATCGCTGTCGGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCGGTCGCCGATGTAGAACAGATAGCCGAGCATCTTCACCGCATCGCGTTGCTGGCCCTCGCCCAGCAAtgctgc                                                                             >  1:3151222/1‑150 (MQ=255)
                                                                         gCCTTTTCGTGGCGCAGGCGAAACCAGTCGTCATCGCTGTCGGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCGGTCGCCGATGTAGAACAGATAGCCGAGCATCTTCACCGCATCGCGTTGCTGGCCCTCGCCCAGCAATGCTGcc                                                                            <  1:2353374/150‑1 (MQ=255)
                                                                         gCCTTTTCGTGGCGCAGGCGAAACCAGTCGTCATCGCTGTCGGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCGGTCGCCGATGTAGAACAGATAGCCGAGCATCTTCACCGCATCGCGTTGCTGGCCCTCGCCCAGCAATGCTGcc                                                                            >  1:3462607/1‑150 (MQ=255)
                                                                           cTTTTCGTGGCGCAGGCGAAACCAGTCGTCATCGCTGTCGGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCGGTCGCCGATGTAGAACAGATAGCCGAGCATCTTCACCGCATCGCGTTGCTGGCCCTCGCCCAGCAATGCTGCCAt                                                                          >  1:3137508/1‑150 (MQ=255)
                                                                           cTTTTCGTGGCGCAGGCGAAACCAGTCGTCATCGCTGTCGGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCGGTCGCCGATGTAGAACAGATAGCCGAGCATCTTCACCGCATCGCGTTGCTGGCCCTCGCCCAGCAATGCTGCCAt                                                                          <  1:4413365/150‑1 (MQ=255)
                                                                           cTTTTCGTGGCGCAGGCGAAACCAGTCGTCATCGCTGTCGGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCGGTCGCCGATGTAGAACAGATAGCCGAGCATCTTCACCGCATCGCGTTGCTGGCCCTCGCCCAGCAATGCTGCCAt                                                                          <  2:2254481/150‑1 (MQ=255)
                                                                           cTTTTCGTGGCGCAGGCGAAACCAGTCGTCATCGCTGTCGGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCAGTCGCCGATGTAGAACAGATAGCCGAGCATCTTCACCGCATCGCGTTGCTGGCCCTCGCCCAGCAATGCTGCCAt                                                                          <  1:4926746/150‑1 (MQ=255)
                                                                           cTTTTCGTGGCGCAGGCGAAACCAGTCGTCATCGCTGTCGGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCAGTCGCCGATGTAGAACAGATAGCCGAGCATCTTCACCGCATCGCGTTGCTGGCCCTCGCCCAGCAATGCTGCCAt                                                                          <  1:4926497/150‑1 (MQ=255)
                                                                           cTTTTCGTGGCGCAGGCGAAACCAGTCGTCATCGCTGTCGGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCAGTCGCCGATGTAGAACAGATAGCCGAGCATCTTCACCGCATCGCGTTGCTGGCCCTCGCCCAGCAATGCTGCCAt                                                                          <  1:4926461/150‑1 (MQ=255)
                                                                                cGTGGCGCAGGCGAAACCAGTCGTCATCGCTGTCGGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCGGTCGCCGATGTAGAACAGATAGCCGAGCATCTTCACCGCATCGCGTTGCTGGCCCTCGCCCAGCAATGCTGCCATCGGca                                                                     <  2:1967255/150‑1 (MQ=255)
                                                                                       cAGGCGAAACCAGTCGTCATCGCTGTCGGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCGGTCGCCGATGTAGAACAGATAGCCGAGCATCTTCACCGCATCGCGTTGCTGGCCCTCGCCCAGCAATGCTGCCATCGGCACATTCAg                                                              <  2:5089023/150‑1 (MQ=255)
                                                                                            gAAACCAGTCGTCATCGCTGTCGGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCGGTCGCCGATGTAGAACAGATAGCCGAGCATCTTCACCGCATCGCGTTGCTGGCCCTCGCCCAGCAATGCTGCCATCGGCACATTCAGGTGCt                                                         >  1:3820206/1‑150 (MQ=255)
                                                                                                  aGTCGTCATCGCTGTCGGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCGGTCGCCGATGTAGAACAGATAGCCGAGCATCTTCACCGCATCGCGTTGCTGGCCCTCGCCCAGCAATGCTGCCATCGGCACATTCAGGTGCTGGCCGa                                                   >  2:2846947/1‑150 (MQ=255)
                                                                                                    tcgtcATCGCTGTCGGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCGGTCGCCGATGTAGAACAGATAGCCGAGCATCTTCACCGCATCGCGTTGCTGGCCCTCGCCCAGCAATGCTGCCATCGGCACATTCAGGTGCTGGCCGAGc                                                 <  1:507985/150‑1 (MQ=255)
                                                                                                    tcgtcATCGCTGTCGGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCGGTCGCCGATGTAGAACAGATAGCCGAGCATCTTCACCGCATCGCGTTGCTGGCCCTCGCCCAGCAATGCTGCCATCGGCACATTCAGGTGCTGGCCGAGc                                                 >  2:3215994/1‑150 (MQ=255)
                                                                                                     cgtcATCGCTGTCGGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCGGTCGCCGATGTAGAACAGATAGCCGAGCATCTTCACCGCATCGCGTTGCTGGCCCTCGCCCAGCAATGCTGCCATCGGCACATTCAGGTGCTGGCCGAGCa                                                <  1:2971331/150‑1 (MQ=255)
                                                                                                     cgtcATCGCTGTCGGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCGGTCGCCGATGTAGAACAGATAGCCGAGCATCTTCACCGCATCGCGTTGCTGGCCCTCGCCCAGCAATGCTGCCATCGGCACATTCAGGTGCTGGCCGAGCa                                                <  2:141862/150‑1 (MQ=255)
                                                                                                      gtcATCGCTGTCGGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCGGTCGCCGATGTAGAACAGATAGCCGAGCATCTTCACCGCATCGCGTTGCTGGCCCTCGCCCAGCAATGCTGCCATCGGCACATTCAGGTGCTGGCCGAGCAg                                               <  1:3207107/150‑1 (MQ=255)
                                                                                                        cATCGCTGTCGGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCGGTCGCCGATGTAGAACAGATAGCCGAGCATCTTCACCGCATCGCGTTGCTGGCCCTCGCCCAGCAATGCTGCCATCGGCACATTCAGGTGCTGGCCGAGCAGGt                                             <  2:1701070/150‑1 (MQ=255)
                                                                                                         aTCGCTGTCGGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCGGTCGCCGATGTAGAACAGATAGCCGAGCATCTTCACCGCATCGCGTTGCTGGCCCTCGCCCAGCAATGCTGCCATCGGCACATTCAGGTGCTGGCCGAGCAGGTc                                            >  1:1995179/1‑150 (MQ=255)
                                                                                                           cGCTGTCGGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCGGTCGCCGATGTAGCACAGATAGCCGAGCATCTTCACCGCATCGCGTTGCTGGCCCTCGCCCAGCAATGCTGCCATCGGCACATTCAGGTGCTGGCCGAGCAGGTCGa                                          <  1:4590030/150‑1 (MQ=255)
                                                                                                           cGCTGTCGGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCGGTCGCCGATGTAGAACAGATAGCCGAGCATCTTCACCGCATCGCGTTGCTGGCCCTCGCCCAGCAATGCTGCCATCGGCACATTCAGGTGCTGGCCGAGCAGGTCGa                                          <  1:4590018/150‑1 (MQ=255)
                                                                                                           cGCTGTCGGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCGGTCGCCGATGTAGAACAGATAGCCGAGCATCTTCACCGCATCGCGTTGCTGGCCCTCGCCCAGCAATGCTGCCATCGGCACATTCAGGTGCTGGCCGAGCAGGTCGa                                          <  1:2207397/150‑1 (MQ=255)
                                                                                                            gCTGTCGGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCGGTCGCCGATGTAGAACAGATAGCCGAGCATCTTCACCGCATCGCGTTGCTGGCCCTCGCCCAGCAATGCTGCCATCGGCACATTCAGGTGCTGGCCGAGCAGGTCGAg                                         >  2:3715182/1‑150 (MQ=255)
                                                                                                                tCGGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCGGTCGCCGATGTAGAACAGATAGCCGAGCATCTTCACCGCATCGCGTTGCTGGCCCTCGCCCAGCAATGCTGCCATCGGCACATTCAGGTGCTGGCCGAGCAGGTCGAGCAAc                                     >  2:1682249/1‑150 (MQ=255)
                                                                                                                tCGGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCGGTCGCCGATGTAGAACAGATAGCCGAGCATCTTCACCGCATCGCGTTGCTGGCCCTCGCCCAGCAATGCTGCCATCGGCACATTCAGGTGCTGGCCGAGCAGGTCGAGCAAc                                     >  2:1738998/1‑150 (MQ=255)
                                                                                                                tCGGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCGGTCGCCGATGTAGAACAGATAGCCGAGCATCTTCACCGCATCGCGTTGCTGGCCCTCGCCCAGCAATGCTGCCATCGGCACATTCAGGTGCTGGCCGAGCAGGTCGAGCAAc                                     >  1:4264571/1‑150 (MQ=255)
                                                                                                                  ggCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCGGTCGCCGATGTAGAACAGATAGCCGAGCATCTTCACCGCATCGCGTTGCTGGCCCTCGCCCAGCAATGCTGCCATCGGCACATTCAGGTGCTGGCCGAGCAGGTCGAGCAACGc                                   <  2:2520824/150‑1 (MQ=255)
                                                                                                                   gcgtCGTTGCGGTAGGCCAGGGCGGTCTGCGGGCGGGCGCCGATGTAGAACAGATAGCCGAGCATCTTCACCGCCTCGCGTTGCTGGCCCTCGCCCAGCAATGCTGCCAACGGCACATGCAGGGGCTGGCCGAGCAGGTCGAGCAACGcc                                  <  1:2250620/147‑1 (MQ=255)
                                                                                                                     gtcgtcgCGGTAGGCCAGGTCGGTCTGCTGGCGGTCGCCGATGTAGAACAGATAGCCGAGCATCTTCACCGCATCGCGTTGCTGGCCCTCGCCCAGCAATGCTGCCATCGGCACATTCAGGTGCTGGCCGAGCAGGTCGAGCAACGCCGa                                <  1:4315920/149‑1 (MQ=255)
                                                                                                                               tAGGCCAGGTCGGTCTGCTGGCGGTCGCCGATGTAGAACAGATAGCCGAGCATCTTCACCGCATCGCGTTGCTGGCCCTCGCCCAGCAATGCTGCCATCGGCACATTCAGGTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGATGGcc                      <  2:5110681/150‑1 (MQ=255)
                                                                                                                                  gCCAGGTCGGTCTGCTGGCGGTCGCCGATGTAGAACAGATAGCCGAGCATCTTCACCGCATCGCGTTGCTGGCCCTCGCCCAGCAATGCTGCCATCGGCACATTCAGGTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGATGGCCGtt                   >  1:311280/1‑150 (MQ=255)
                                                                                                                                    cAGGTCGGTCTGCTGGCGGTCGCCGATGTAGAACAGATAGCCGAGCATCTTCACCGCATCGCGTTGCTGGCCCTCGCCCAGCAATGCTGCCATCGGCACATTCAGGTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTAc                 >  1:2015615/1‑150 (MQ=255)
                                                                                                                                    cAGGTCGGTCTGCTGGCGGTCGCCGAAGTAGAACAGATAGCCGAGCATCTTCACCGCATCGCGTTGCTGGCCCTCGCCAAGCAATGCTGCCATAGGCACATTCAGGTGCTGGCCGAGCAGGTCGAGCAACTCCGACTCGATGGCCGTTAc                 >  1:469147/1‑150 (MQ=255)
                                                                                                                                    aaGGTCGGTCTGCTGGCGGTCGCCGATGTAGAACAGATAGCCGAGCATCTTCACCGCATCGCGTTGCTGGCCCTCGCCCAGCAATGCTGCCATCGGCACATTCAGGTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTAc                 >  2:2351068/2‑150 (MQ=255)
                                                                                                                                      ggtcggtcTGCTGGCGGTCGCCGATGTAGAACAGATAGCCGAGCTTCTTCACCGCATCGCGTTGCTGGCGCTCGCCCAGCAATGCTGCCATGGGCACAGTCAGGTGCTGGCCGAGCAGGTGGAGCAACGCCGACTCGAGGGCCGGGACCg               >  2:3761060/1‑150 (MQ=255)
                                                                                                                                      ggtcggtcTGCTGGCGGTCGCCGATGTAGAACAGATAGCCGAGCATCTTCACCGCATCGCGTTGCTGGCCCTCGCCCAGCAATGCTGCCATCGGCACATTCAGGTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTACCg               >  1:3570770/1‑150 (MQ=255)
                                                                                                                                       gtcggtcTGCTGGCGGTCGCCGATGTAGAACAGATAGCCGAGCATCTTCACCGCATCGCGTTGCTGGCCCTCGCCCTGCAATGCTGCCATCGGCACATTCAGGTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTACCGc              <  1:2207590/150‑1 (MQ=255)
                                                                                                                                         cggtcTGCTGGCGGTCGCCGATGTAGAACAGATAGCCGAGCATCTTCACCGCATCGCGTTGCTGGCCCTCGCCCAGCAATGCTGCCATCGGCACATTCAGGTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTACCGCAt            <  2:458582/150‑1 (MQ=255)
                                                                                                                                             ctgctgGCGGTCGCCGATGTAGAACAGATAGCCGAGCATCTTCACCGCATCGCGTTGCTGGCCCTCGCCCAGCAATGCTGCCATCGGCACATTCAGGTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTACCGCATGTAc        <  1:1244428/150‑1 (MQ=255)
                                                                                                                                              tgctgGCGGTCGCCGATGTAGAACAGATAGCCGAGCATCTTCACCGCATCGCGTTGCTGGCCCTCGCCCAGCAATGCTGCCATCGGCACATTCAGGTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTACCGCATGTACg       >  1:4592558/1‑150 (MQ=255)
                                                                                                                                              tgctgGCGGTCGCCGATGTAGAACAGATAGCCGAGCATCTTCACCGCATCGCGTTGCTGGCCCTCGCCCAGCAATGCTGCCATCGGCACATTCAGGTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTACCGCATGTACg       <  2:1461396/150‑1 (MQ=255)
                                                                                                                                                ctgGCGGTCGCCGATGTAGAACAGATAGCCGAGCATCTTCACCGCATCGCGTTGCTGGCCCTCGCCCAGCAATGCTGCCATCGGCACATTCAGGTGCTGGCCGAGCAGGTCGAGCAACGc                                   <  1:1730110/120‑1 (MQ=255)
                                                                                                                                                ctgGCGGTCGCCGATGTAGAACAGATAGCCGAGCATCTTCACCGCATCGCGTTGCTGGCCCTCGCCCAGCAATGCTGCCATCGGCACATTCAGGTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGt     >  1:3303351/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GAGCTTGAAGTCGCTGAAGCCATAGCGCGCCTTGGCGGCCTCGGCCAGCCGCACCACTGCTTCCGGGGTCAGGGCCTTTTCGTGGCGCAGGCGAAACCAGTCGTCATCGCTGTCGGCTTCGTTGCGGTAGGCCAGGTCGGTCTGCTGGCGGTCGCCGATGTAGAACAGATAGCCGAGCATCTTCACCGCATCGCGTTGCTGGCCCTCGCCCAGCAATGCTGCCATCGGCACATTCAGGTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAAT  >  NC_002947/5419027‑5419323

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