Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,417,434 C→G 10.8% A100P (GCG→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,4340CG10.8% 335.5 / 27.7 111A100P (GCG→CCG) gudDbifunctional D‑glucarate dehydratase/L‑idarate epimerase
Reads supporting (aligned to +/- strand):  ref base C (54/45);  new base G (9/3);  total (63/48)
Fisher's exact test for biased strand distribution p-value = 2.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.38e-01

Only 100 of 119 total aligned reads displayed.

CTTCACCGCATCGCGTTGCTGGCCCTCGCCCAGCAATGCTGCCATCGGCACATTCAGGTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCGGGCGTCCTCCAGGGTCTGGCGAATCTTCTCGCCGCCGGGCACTTCGCCCAGGCCGGTATTGCC  >  NC_002947/5417285‑5417581
                                                                                                                                                      |                                                                                                                                                   
cTTCACCGCATCGCGTTGCTGGCCCTCGCCCAGCAATGCTGCCATCGGCACATTCAGGTGCTGTCCGGGCGGGGCGAGCAACGCCGACTCGATGGCCATTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACgcgc                                                                                                                                                     >  2:75853/1‑150 (MQ=255)
 ttCACCGCATCGCGTTGCTGGCCCTCGCCCAGCAATGCTGCCATCGGCACATTCAGGTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACccgcc                                                                                                                                                    <  1:1085740/150‑1 (MQ=255)
 ttCACCGCATCACGTTGCTGGCCCTCGCCCAGCAATGCTGCCATCGGCACATTCAGGTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACccgcc                                                                                                                                                    >  2:1956883/1‑150 (MQ=255)
   ccccGCACCGCGTTGCTGGCCCTCCCACAGCAATGCTGCCATCGGCACATTCAGGTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGa                                                                                                                                                  <  1:1136261/148‑1 (MQ=255)
    aCCGCATCGCGTTGCTGGCCCTCGCCCAGCAATGCTGCCATCGGCACATTCAGGTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAg                                                                                                                                                 >  2:374737/1‑150 (MQ=255)
    aCCGCATCGCGTTGCTGGCCCTCGCCCAGCAATGCTGCCATCGGCACATTCAGGTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAg                                                                                                                                                 >  2:1640014/1‑150 (MQ=255)
     ccGCATCGCGTTGCTGGCCCTCGCCCAGCAATGCTGCCATCGGCACATTCAGGTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGt                                                                                                                                                <  1:1758696/150‑1 (MQ=255)
      cGCATCGCGTTGCTGGCCCTCGCCCAGCAATGCTGCCATCGGCACATTCAGTTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTc                                                                                                                                               <  1:146928/150‑1 (MQ=255)
      cGCATCGCGTTGCTGGCCCTCGCCCAGCAATGCTGCCATCGGCACATTCAGGTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTc                                                                                                                                               >  2:663865/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)
                   tGGCCCTCGCCCAGCAATGCTGCCATCGGCACATTCAGGTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAg                                                                                                                                  <  1:302601/150‑1 (MQ=255)
                       ccTCGCCCAGCAATGCTGCCATCGGCACATTCAGGTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGaa                                                                                                                                   <  2:2191649/145‑1 (MQ=255)
                       ccTCGCCCAGCAATGCTGCCATCGGCACATTCAGGTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGaa                                                                                                                                   >  1:2191649/1‑145 (MQ=255)
                       ccTCGCCCAGCAATGCTGCCATCGGCACATTCAGGTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCt                                                                                                                              >  2:553512/1‑150 (MQ=255)
                       ccTCGCCCAGCAATGCTGCCATCGGCACATTCAGGTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCt                                                                                                                              >  2:1765298/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)
                                  aaTGCTGCCATCGGCACATTCAGGTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTCCAGTGTGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCa                                                                                                                   <  1:1956883/150‑1 (MQ=255)
                                    tgctgcCATCGGCACATTCAGGTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCAtt                                                                                                                 <  1:979246/150‑1 (MQ=255)
                                      cTGCCATCGGCACATTCAGGTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGa                                                                                                               <  2:1277763‑M2/149‑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                                                                                                              >  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)
                                                       aGGTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATgg                                                                                              <  2:2400397/150‑1 (MQ=255)
                                                        ggTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGc                                                                                             <  1:2231485/150‑1 (MQ=255)
                                                         gTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGcc                                                                                            <  2:2024611/150‑1 (MQ=255)
                                                               gCCGAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCgg                                                                                      >  2:1588802/1‑150 (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                                                                >  2:2003852/1‑150 (MQ=255)
                                                                                     gACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCgg                                                                <  1:410733/150‑1 (MQ=255)
                                                                                     gACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCgg                                                                >  1:1125517/1‑150 (MQ=255)
                                                                                        tCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACgcgc                                                                                                       >  1:1747170/1‑108 (MQ=255)
                                                                                          gATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCGGGCGtc                                                           >  2:1971192/1‑150 (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                                                                          <  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                                                                                             <  2:447588/98‑1 (MQ=255)
                                                                                                            tACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGc                                                                                             >  1:447588/1‑98 (MQ=255)
                                                                                                             aCGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCGGGCGTCCTCCAGGGTCTGGCGAAtc                                        >  1:1599303/1‑150 (MQ=255)
                                                                                                               ggTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCGGGCGTCCTCCAGGGTCTGGCGAAtctt                                      >  1:1189454/1‑150 (MQ=255)
                                                                                                               ggTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCGGGCGTCCTCCAGGGTCTGGCGAAtctt                                      <  1:1276619/150‑1 (MQ=255)
                                                                                                                 tAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCGGGCGTCCTCCAGGGTCTGGCGAAtcttct                                    <  1:2388893/150‑1 (MQ=255)
                                                                                                                 tAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCGGGCGTCCTCCAGGGTCTGGCGAAtcttct                                    >  2:2380903/1‑150 (MQ=255)
                                                                                                                  aaTGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCTAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCGGGCGTCCTCCAGGGTCTGGCGAATCTtctc                                   <  2:1643007/150‑1 (MQ=255)
                                                                                                                  aaTGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCGGGCGTCCTCCAGGGTCTGGCGAATCTtctc                                   <  1:1668826/150‑1 (MQ=255)
                                                                                                                         aGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGaccacc                                                                         >  2:1725313/1‑105 (MQ=255)
                                                                                                                         aGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGaccacc                                                                         <  1:1725313/105‑1 (MQ=255)
                                                                                                                         aGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCGGGCGTCCTCCAGGGTCTGGCGAATCTTCTCGccgccg                            >  2:2208017/1‑150 (MQ=255)
                                                                                                                          ggTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCGGGCGTCCTCCAGGGTCTGGCGAATCTTCTCGCCGCCgg                           <  2:2212155/150‑1 (MQ=255)
                                                                                                                           gTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCGGGCGTCCTCCAGGGTCTGGCGAATCTTCTCGCCGCCggg                          >  2:47580/1‑150 (MQ=255)
                                                                                                                              gAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCGGGCGTCCTCCAGGGTCTGGCGAATCTTCTCGCCGCCGGGCAc                       >  2:2189747/1‑150 (MQ=255)
                                                                                                                              gAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCGGGCGTCCTCCAGGGTCTGGCGAATCTTCTCGCCGCCGGGCAc                       >  1:2259283/1‑150 (MQ=255)
                                                                                                                              gAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCGGGCGTCCTCCAGGGTCTGGCGAATCTTCTCGCCGCCGGGCAc                       >  1:31192/1‑150 (MQ=255)
                                                                                                                               aaGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCGGGCGTCCTCCAgg                                                    <  2:671979/120‑1 (MQ=255)
                                                                                                                               aaGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCGGGCGTCCTCCAgg                                                    >  1:671979/1‑120 (MQ=255)
                                                                                                                                       gatccctcgACCCTCCGAGTCGCGGTTGGCGAAGATCTTGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCTGCTGGCCGACCACCAGGCTGCGGGCGTCCTCCAGGGTCTGGCGAATCTTCTCGCCGCCGGGCACTTCGCCCAg              <  1:1138034/143‑1 (MQ=255)
                                                                                                                                       ccaaccGCGACTCGGCGGGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGAACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCGGGCGTCCTCCAGGGTCTGGCGACTCTTCTCGCCGCCGGGCACTTCGCCCAg              >  2:655650/5‑150 (MQ=255)
                                                                                                                                        caaccGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGAccacca                                                                        <  1:1920734/88‑1 (MQ=255)
                                                                                                                                        aGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCGGGCGTCCTCCAGGGTCTGGCGAATCTTCTCGCCGCCGGGCACTTCGCCCAgg             <  1:75853/150‑1 (MQ=255)
                                                                                                                                         ggCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCGGGCGTCCTCCAGGGTCTGGCGAATCTTCTCGCCGCCGGGCACTTCGCCCAGGc            <  1:426986/150‑1 (MQ=255)
                                                                                                                                          accGCGACTCGGCGGGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCGGGCGTCCTCCAGGGTCTGGCGAATCTTCTCGCCGCCGGGCACTTCGCCCAGGcc           >  2:1080500/2‑150 (MQ=255)
                                                                                                                                           ccGCGACTCGGCGGGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCGGGCGTCCTCCAGGGTCTGGCGAATCTTCTCGCCGCCGGGCACTTCGCCCAGGCCg          >  1:1060327/1‑150 (MQ=255)
                                                                                                                                           ccGCGACTCGGCGGGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCGGGCGTCCTCCAGGGTCTGGCGAATCTTCTCGCCGCCGGGCACTTCGCCCAGGCCg          >  2:1573769/1‑150 (MQ=255)
                                                                                                                                            cgcgACTCGGCGGGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGa                                                                                          <  2:1055919/69‑1 (MQ=38)
                                                                                                                                            cgcgACTCGGCGGGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGa                                                                                          >  1:1055919/1‑69 (MQ=38)
                                                                                                                                            cgcgACTCGGCGGGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCGGGCGTCCTCCAGGGTCt                                                <  2:1914549/111‑1 (MQ=255)
                                                                                                                                            cgcgACTCGGCGGGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCGGGCGTCCTCCAGGGTCTGGCGAATCTTCTCGCCGCCGGGCACTTCGCCCAGGCCgg         >  2:2406015/1‑150 (MQ=255)
                                                                                                                                              cgACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCGGGCGTCCTCCAGGGTCTGGCGAATCTTCTCGCCGCCGGGCACTTCGCCCAGGCCGGTa       >  1:2210760/1‑150 (MQ=255)
                                                                                                                                                aCCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCGGGCGTCCTCCAGGGTCTGGCGAATCTTCTCGCCGCCGGGCACTTCGCCCAGGCCGGTAtt     <  2:728599/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)
                                                                                                                                                      |                                                                                                                                                   
CTTCACCGCATCGCGTTGCTGGCCCTCGCCCAGCAATGCTGCCATCGGCACATTCAGGTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCGGGCGTCCTCCAGGGTCTGGCGAATCTTCTCGCCGCCGGGCACTTCGCCCAGGCCGGTATTGCC  >  NC_002947/5417285‑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: