Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,417,437 A→G 10.5% S99P (TCG→CCG)  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,417,4370AG10.5% 309.7 / 19.6 114S99P (TCG→CCG) gudDbifunctional D‑glucarate dehydratase/L‑idarate epimerase
Reads supporting (aligned to +/- strand):  ref base A (57/45);  new base G (9/3);  total (66/48)
Fisher's exact test for biased strand distribution p-value = 2.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.28e-01

Only 100 of 117 total aligned reads displayed.

CACCGCATCGCGTTGCTGGCCCTCGCCCAGCAATGCTGCCATCGGCACATTCAGGTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCGGGCGTCCTCCAGGGTCTGGCGAATCTTCTCGCCGCCGGGCACTTCGCCCAGGCCGGTATTGCC  >  NC_002947/5417288‑5417581
                                                                                                                                                      |                                                                                                                                                
ccccGCACCGCGTTGCTGGCCCTCCCACAGCAATGCTGCCATCGGCACATTCAGGTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGa                                                                                                                                                  <  1:1136261/148‑1 (MQ=255)
 aCCGCATCGCGTTGCTGGCCCTCGCCCAGCAATGCTGCCATCGGCACATTCAGGTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAg                                                                                                                                                 >  2:374737/1‑150 (MQ=255)
 aCCGCATCGCGTTGCTGGCCCTCGCCCAGCAATGCTGCCATCGGCACATTCAGGTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAg                                                                                                                                                 >  2:1640014/1‑150 (MQ=255)
   cGCATCGCGTTGCTGGCCCTCGCCCAGCAATGCTGCCATCGGCACATTCAGTTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTc                                                                                                                                               <  1:146928/150‑1 (MQ=255)
   cGCATCGCGTTGCTGGCCCTCGCCCAGCAATGCTGCCATCGGCACATTCAGGTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTc                                                                                                                                               >  2:663865/1‑150 (MQ=255)
    gCATCGCGTTGCTGGCCCTCGCCCAGCAATGCTGCCATCGGCACATTCAGGTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTcg                                                                                                                                              >  2:2231485/1‑150 (MQ=255)
     cATCGCGTTGCTGGCCCTCGCCCAGCAATGCTGCCATCGGCACATTCAGGTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACTCGGCGGGTcgc                                                                                                                                             <  2:1542220/150‑1 (MQ=255)
         GCGTTGCTGGCCCTCGCCCAGCAATGCTGCCATCGGCACATTCAGGTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACTCGGCGGGTCgcggcc                                                                                                                                         >  1:1120110‑M1/1‑148 (MQ=255)
                tGGCCCTCGCCCAGCAATGCTGCCATCGGCACATTCAGGTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGtt                                                                                                                                         <  2:821841/143‑1 (MQ=255)
                tGGCCCTCGCCCAGCAATGCTGCCATCGGCACATTCAGGTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGtt                                                                                                                                         >  1:821841/1‑143 (MQ=255)
                    ccTCGCCCAGCAATGCTGCCATCGGCACATTCAGGTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGaa                                                                                                                                   >  1:2191649/1‑145 (MQ=255)
                    ccTCGCCCAGCAATGCTGCCATCGGCACATTCAGGTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGaa                                                                                                                                   <  2:2191649/145‑1 (MQ=255)
                    ccTCGCCCAGCAATGCTGCCATCGGCACATTCAGGTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCt                                                                                                                              >  2:553512/1‑150 (MQ=255)
                      tCGCCCAGCAATGCTGCCATCGGCACATTCAGGTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTgg                                                                                                                            <  2:359441/150‑1 (MQ=255)
                          ccAGCAATTCTGCCATCGGCACATTCAGGTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCa                                                                                                                        <  1:1915456/150‑1 (MQ=255)
                           cAGCAATGCTGCCATCGGCACATTCAGGTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCAt                                                                                                                       <  1:1640014/150‑1 (MQ=255)
                                 tgctgcCATCGGCACATTCAGGTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCAtt                                                                                                                 <  1:979246/150‑1 (MQ=255)
                                   cTGCCATCGGCACATTCAGGTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGa                                                                                                               <  2:1277763‑M2/149‑1 (MQ=255)
                                    tgcCATCGGCACATTCAGGTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGgc                                                                                                                           <  1:1282248/137‑1 (MQ=255)
                                    tgcCATCGGCACATTCAGGTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGgc                                                                                                                           >  2:1282248/1‑137 (MQ=255)
                                    tgcCATCGGCACATTCAGGTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGa                                                                                                               >  2:1992712/1‑149 (MQ=255)
                                    tgcCATCGGCACATTCAGGTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGa                                                                                                               <  1:1992712/149‑1 (MQ=255)
                                    tgcCATCGGCACATTCAGGTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAg                                                                                                              >  2:346547/1‑150 (MQ=255)
                                    tgcCATCGGCACATTCAGGTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAg                                                                                                              >  1:879073/1‑150 (MQ=255)
                                    tgcCATCGGCACATTCAGGTGCTGGCCAAGCAGGTGAAGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAg                                                                                                              <  2:192047/150‑1 (MQ=255)
                                      ccATCGGCACATTCAGGTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCa                                                                                                            >  2:1853119/1‑150 (MQ=255)
                                      ccATCGGCACATTCAGGTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCa                                                                                                            <  2:879073/150‑1 (MQ=255)
                                      ccATCGGCACATTCAGGTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCa                                                                                                            >  1:2024611/1‑150 (MQ=255)
                                        aTCGGCACATTCAGGTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGc                                                                                                                                      <  1:2375786/122‑1 (MQ=255)
                                        aTCGGCACATTCAGGTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGc                                                                                                                                      >  2:2375786/1‑122 (MQ=255)
                                           ggCACATTCAGGTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACgcgc                                                                                                       >  2:2342133/1‑150 (MQ=255)
                                            gCACATTCAGGTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCt                                                                                                      >  1:956281/1‑150 (MQ=255)
                                              acaTTCAGGTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTACCGCGGCGACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTgg                                                                                                    <  1:1835730/150‑1 (MQ=255)
                                              acaTTCAGGTGCTGGCCGAGCAGGTAGACCTACGCCGACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTgg                                                                                                    >  1:1643007/1‑150 (MQ=255)
                                               caTTCAGGTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGt                                                                                                   <  1:374737/150‑1 (MQ=255)
                                                 ttCAGGTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTaa                                                                                                 <  2:1965203/150‑1 (MQ=255)
                                                 ttCAGGTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTaa                                                                                                 >  1:728599/1‑150 (MQ=255)
                                                    aGGTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATgg                                                                                              <  2:2400397/150‑1 (MQ=255)
                                                    aGGTGCTGGCCGAGCACGTCGAGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATgg                                                                                              >  2:2021301/1‑150 (MQ=255)
                                                      gTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGcc                                                                                            <  2:2024611/150‑1 (MQ=255)
                                                            accGAGCACGTCGAGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAATGCGCACCTCCCAGGTCTGCCGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAAGGCTAGCGCATGGCATTGAGCACGCGCTCGTAATGGCCGAACgc                                                                                      >  2:585891/2‑149 (MQ=255)
                                                              cGAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCt                                                                                    <  1:1765298/150‑1 (MQ=255)
                                                                  cAGGTCGAGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGcc                                                                                >  2:410733/1‑150 (MQ=255)
                                                                      tCGAGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGacc                                                                            >  2:1668826/1‑150 (MQ=255)
                                                                        gAGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGaccac                                                                          >  2:2338705/1‑150 (MQ=255)
                                                                           cAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAg                                                                       >  2:1276619/1‑150 (MQ=255)
                                                                                  gACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCgg                                                                <  1:410733/150‑1 (MQ=255)
                                                                                  gACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCgg                                                                >  1:1123874/1‑150 (MQ=255)
                                                                                  gACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCgg                                                                >  2:2003852/1‑150 (MQ=255)
                                                                                     tCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACgcgc                                                                                                       <  2:1747170/108‑1 (MQ=255)
                                                                                     tCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACgcgc                                                                                                       >  1:1747170/1‑108 (MQ=255)
                                                                                           gCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCGGGCGtcctcc                                                       >  1:934614/1‑150 (MQ=255)
                                                                                              gTTACCGCATGTACGGTAATGCGCAGGTCGAGGGTCTGCGGGCCGCGCCCCGCCCAGTCGCGGTTGGCCCTTGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCGGGCGTCCTCCAgg                                                    <  1:1853119/150‑1 (MQ=255)
                                                                                              gTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCGGGCGTCCTCCAgg                                                    >  1:1073608/1‑150 (MQ=255)
                                                                                              gTTACCGCATGTACGGTAATGCGCAGGTCGAAAGTCTGCCGGCGGGCCTCGACCGAGGCGTGGTTGGCGAAGGTATGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCGGGCGTCCTCCAgg                                                    <  1:2342133/150‑1 (MQ=255)
                                                                                                  ccGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCGGGCGTCCTCCAGGGTCt                                                >  1:868498/1‑150 (MQ=255)
                                                                                                      aTGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGaccac                                                                          >  1:1676000/1‑120 (MQ=255)
                                                                                                      aTGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGaccac                                                                          <  2:1676000/120‑1 (MQ=255)
                                                                                                      aTGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCGGGCGTCCTCCAGGGTCTGGCg                                            >  2:184788/1‑150 (MQ=255)
                                                                                                       tGTACGTTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCGGGCGTCCTCCAGGGTCTGGCGa                                           <  2:868498/150‑1 (MQ=255)
                                                                                                        gTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCGGGCGTCCTCCAGGGTCTGGCGaa                                          <  1:553512/150‑1 (MQ=255)
                                                                                                         tACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGc                                                                                             >  1:447588/1‑98 (MQ=255)
                                                                                                            ggTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCGGGCGTCCTCCAGGGTCTGGCGAAtctt                                      <  1:1276619/150‑1 (MQ=255)
                                                                                                            ggTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCGGGCGTCCTCCAGGGTCTGGCGAAtctt                                      >  1:1189454/1‑150 (MQ=255)
                                                                                                              tAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCGGGCGTCCTCCAGGGTCTGGCGAAtcttct                                    <  1:2388893/150‑1 (MQ=255)
                                                                                                               aaTGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCTAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCGGGCGTCCTCCAGGGTCTGGCGAATCTtctc                                   <  2:1643007/150‑1 (MQ=255)
                                                                                                               aaTGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCGGGCGTCCTCCAGGGTCTGGCGAATCTtctc                                   <  1:1668826/150‑1 (MQ=255)
                                                                                                                   cgcAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCGGGCGTCCTCCAGGGTCTGGCGAATCTTCTcgccg                               >  1:2274376/1‑150 (MQ=255)
                                                                                                                      aGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGaccacc                                                                         <  1:1725313/105‑1 (MQ=255)
                                                                                                                      aGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGaccacc                                                                         >  2:1725313/1‑105 (MQ=255)
                                                                                                                      aGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCGGGCGTCCTCCAGGGTCTGGCGAATCTTCTCGccgccg                            >  2:2208017/1‑150 (MQ=255)
                                                                                                                       ggTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCGGGCGTCCTCCAGGGTCTGGCGAATCTTCTCGCCGCCgg                           <  2:2212155/150‑1 (MQ=255)
                                                                                                                        gTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCGGGCGTCCTCCAGGGTCTGGCGAATCTTCTCGCCGCCggg                          >  2:47580/1‑150 (MQ=255)
                                                                                                                           gAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCGGGCGTCCTCCAGGGTCTGGCGAATCTTCTCGCCGCCGGGCAc                       >  1:2259283/1‑150 (MQ=255)
                                                                                                                           gAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCGGGCGTCCTCCAGGGTCTGGCGAATCTTCTCGCCGCCGGGCAc                       >  1:31192/1‑150 (MQ=255)
                                                                                                                           gAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCGGGCGTCCTCCAGGGTCTGGCGAATCTTCTCGCCGCCGGGCAc                       >  2:2189747/1‑150 (MQ=255)
                                                                                                                            aaGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCGGGCGTCCTCCAGGGTCTGGCGAATCTTCTCGCCGCCGGGCACt                      >  2:1586770/1‑150 (MQ=255)
                                                                                                                                    gatccctcgACCCTCCGAGTCGCGGTTGGCGAAGATCTTGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCTGCTGGCCGACCACCAGGCTGCGGGCGTCCTCCAGGGTCTGGCGAATCTTCTCGCCGCCGGGCACTTCGCCCAg              <  1:1138034/143‑1 (MQ=255)
                                                                                                                                    ccaaccGCGACTCGGCGGGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGAACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCGGGCGTCCTCCAGGGTCTGGCGACTCTTCTCGCCGCCGGGCACTTCGCCCAg              >  2:655650/5‑150 (MQ=255)
                                                                                                                                     caaccGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGAccacca                                                                        >  2:1920734/4‑91 (MQ=255)
                                                                                                                                     caaccGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGAccacca                                                                        <  1:1920734/88‑1 (MQ=255)
                                                                                                                                     aGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCGGGCGTCCTCCAGGGTCTGGCGAATCTTCTCGCCGCCGGGCACTTCGCCCAgg             <  1:75853/150‑1 (MQ=255)
                                                                                                                                     aGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATGGAGCACGCGGTGCTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCGGGCGTCCTCCAGGGTCTGGCGAATCTTCTCGCCGCCGGGCACTTCGCCCAgg             >  2:727458/1‑150 (MQ=255)
                                                                                                                                      ggCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCGGGCGTCCTCCAGGGTCTGGCGAATCTTCTCGCCGCCGGGCACTTCGCCCAGGc            <  1:426986/150‑1 (MQ=255)
                                                                                                                                      aaccGCGACTCGGCGGGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCGGGCGTCCTCCAGGGTCTGGCGAATCTTCTCGCCGCCGGGCACTTCGCCCAGGc            >  1:2125296/3‑150 (MQ=255)
                                                                                                                                       accGCGACTCGGCGGGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCGGGCGTCCTCCAGGGTCTGGCGAATCTTCTCGCCGCCGGGCACTTCGCCCAGGcc           >  2:1080500/2‑150 (MQ=255)
                                                                                                                                        ccGCGACTCGGCGGGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCGGGCGTCCTCCAGGGTCTGGCGAATCTTCTCGCCGCCGGGCACTTCGCCCAGGCCg          >  2:1573769/1‑150 (MQ=255)
                                                                                                                                        ccGCGACTCGGCGGGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCGGGCGTCCTCCAGGGTCTGGCGAATCTTCTCGCCGCCGGGCACTTCGCCCAGGCCg          >  1:1060327/1‑150 (MQ=255)
                                                                                                                                         cgcgACTCGGCGGGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGa                                                                                          >  1:1055919/1‑69 (MQ=38)
                                                                                                                                         cgcgACTCGGCGGGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGa                                                                                          <  2:1055919/69‑1 (MQ=38)
                                                                                                                                         cgcgACTCGGCGGGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCGGGCGTCCTCCAGGGTCt                                                <  2:1914549/111‑1 (MQ=255)
                                                                                                                                         cgcgACTCGGCGGGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCGGGCGTCCTCCAGGGTCt                                                >  1:1914549/1‑111 (MQ=255)
                                                                                                                                         cgcgACTCGGCGGGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCGGGCGTCCTCCAGGGTCTGGCGAATCTTCTCGCCGCCGGGCACTTCGCCCAGGCCgg         >  2:2406015/1‑150 (MQ=255)
                                                                                                                                           cgACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCGGGCGTCCTCCAGGGTCTGGCGAATCTTCTCGCCGCCGGGCACTTCGCCCAGGCCGGTa       >  1:2210760/1‑150 (MQ=255)
                                                                                                                                             aCCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCGGGCGTCCTCCAGGGTCTGGCGAATCTTCTCGCCGCCGGGCACTTCGCCCAGGCCGGTAtt     <  2:728599/150‑1 (MQ=255)
                                                                                                                                               ccgctgAGTCGCGGGTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCGGGCGTCCTCCAGGGTCTGGCGAATCTTCTCGCCGCCGGGCACTTCGCCCAGGCCGGTATTGc   <  2:515750/150‑1 (MQ=255)
                                                                                                                                               ccgccgAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCGGGCGTCCTCCAGGGTCTGGCGAATCTTCTCGCCGCCGGGCACTTCGCCCAGGCCGGTATTGc   >  1:765383/1‑150 (MQ=255)
                                                                                                                                               ccgccgAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCGGGCGTCCTCCAGGGTCTGGCGAATCTTCTCGCCGCCGGGCACTTCGCCCAGGCCGGTATTGc   >  1:2271328/1‑150 (MQ=255)
                                                                                                                                                cgccgAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCGGGCGTCCTCCAGGGTCTGGCGAATCTTCTCGCCGCCGGGCACTTCGCCCAGGCCGGTATTGcc  >  1:1077788/1‑150 (MQ=255)
                                                                                                                                                      aGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCt                                                                                    <  1:2201673/63‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                
CACCGCATCGCGTTGCTGGCCCTCGCCCAGCAATGCTGCCATCGGCACATTCAGGTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCGGGCGTCCTCCAGGGTCTGGCGAATCTTCTCGCCGCCGGGCACTTCGCCCAGGCCGGTATTGCC  >  NC_002947/5417288‑5417581

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