Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,572,931 C→A 6.6% A337A (GCC→GCA pcaB → 3‑carboxy‑cis,cis‑muconate cycloisomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,572,9310CA 6.6% 747.2 / 27.7 225A337A (GCC→GCApcaB3‑carboxy‑cis,cis‑muconate cycloisomerase
Reads supporting (aligned to +/- strand):  ref base C (96/114);  new base A (5/10);  total (101/124)
Fisher's exact test for biased strand distribution p-value = 4.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 227 total aligned reads displayed.

GCCGCCGTGTTGATCGGTGCCGCGACCCGCGTGCCGGGCCTGCTGTCGACGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTGGTGCTGGCCGAAGCGGTGAGCATCGTCCTCGCCCAGCGTCTGGGTCGCGACCGTGCCCACCA  >  NC_002947/1572782‑1573077
                                                                                                                                                     |                                                                                                                                                  
         ttGATCGGTGCCGCGACCCGCGTGCCGGGCCTGCTGTCGACGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAg                                                                                                                                           <  2:5226940/150‑1 (MQ=255)
         ttGATCGGTGCCGCGACCCGCGTGCCGGGCCTGCTGTCGACGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCAGGGCGCCAg                                                                                                                                           <  2:1748112/150‑1 (MQ=255)
          tGATCGGTGCCGCGACCCGCGTGCCGGGCCTGCTGTCGACGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAgg                                                                                                                                          <  1:894248/150‑1 (MQ=255)
          tGATCGGTGCCGCGACCCGCGTGCCGGGCCTGCTGTCGACGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCAGGGCGCCAga                                                                                                                                          <  2:697245/150‑2 (MQ=255)
            aTCGGTGCCGCGACCCGCGTGCCGGGCCTGCTGTCGACGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCAGGGCGCCAgaga                                                                                                                                        <  1:3507427/150‑4 (MQ=255)
            aTCGGTGCCGCGACCCGCGTGCCGGGCCTGCTGTCGACGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCAGGGCGCCAgaga                                                                                                                                        <  1:772776/150‑4 (MQ=255)
            aTCGGTGCCGCGACCCGCGTGCCGGGCCTGCTGTCGACGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCAGGGCGCCAgaga                                                                                                                                        <  1:2232815/150‑4 (MQ=255)
               ggTGCCGCGACCCGCGTGCCGGGCCTGCTGTCGACGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCCCTGGCGCAGGGCGCCAgagacca                                                                                                                                     <  2:1172908/150‑7 (MQ=255)
                gTGCCGCGACCCGCGTGCCGGGCCTGCTGTCGACGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAg                                                                                                                                    >  2:3020821/1‑150 (MQ=255)
                   ccGCGACCCGCGTGCCGGGCCTGCTGTCGACGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGtca                                                                                                                                 >  2:5191206/1‑148 (MQ=255)
                    cgcgACCCGCGTGCCGGGCCTGCTGTCGACGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCGGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGAt                                                                                                                                <  2:5022729/150‑1 (MQ=255)
                     gcgACCCGCGTGCCGGGCCTGCTGTCGACGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGAtt                                                                                                                               >  1:3770342/1‑150 (MQ=255)
                       gACCCGCGTGCCGGGCCTGCTGTCGACGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGc                                                                                                                             >  1:3012381/1‑150 (MQ=255)
                          ccGCGTGCCGGGCCTGCTGTCGACGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGa                                                                                                                          <  2:3747182/150‑1 (MQ=255)
                           cgcgTGCCGGGCCTGCTGTCGACGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAg                                                                                                                         >  2:2075246/1‑150 (MQ=255)
                           cgcgTGCCGGGCCTGCTGTCGACGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAg                                                                                                                         <  1:4174616/150‑1 (MQ=255)
                             cgTGCCGGGCCTGCTGTCGACGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAggg                                                                                                                       <  2:5381289/150‑1 (MQ=255)
                              gTGCCGGGCCTGCTGTCGACGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCATCCTGGGCCTATGGCATGCCGAGTGGGTAACCCTGCCGGATTTCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGc                                                                                                                      <  2:4501681/150‑1 (MQ=255)
                                  cGGGCCTGCTGTCGACGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATgg                                                                                                                  <  1:4923470/150‑1 (MQ=255)
                                  cGGGCCTGCTGTCGACGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATgg                                                                                                                  >  1:1208693/1‑150 (MQ=255)
                                   gggCCTGCTGTCGACGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGa                                                                                                                 >  1:3365339/1‑150 (MQ=255)
                                        tgctgTCGACGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTgg                                                                                                            <  2:1208693/150‑1 (MQ=255)
                                        tgctgTCGACGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTgg                                                                                                            <  2:1811098/150‑1 (MQ=255)
                                        tgctgTCGACGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTgg                                                                                                            >  2:4097934/1‑150 (MQ=255)
                                         gctgTCGACGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGa                                                                                                           <  1:1676453/150‑1 (MQ=255)
                                          ctgTCGACGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGAt                                                                                                          <  1:1471034/150‑1 (MQ=255)
                                          ctgTCGACGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGAt                                                                                                          <  1:4677896/150‑1 (MQ=255)
                                          ctgTCGACCCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGAt                                                                                                          <  1:3617380/150‑1 (MQ=255)
                                           tgTCGACGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCAGGGCGCCAgagac                                                                                                                                       >  2:1409140/1‑116 (MQ=255)
                                            gTCGACGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGc                                                                                                        <  2:1714108/150‑1 (MQ=255)
                                             tCGACGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGcc                                                                                                       >  1:4099037/1‑150 (MQ=255)
                                               gACGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCcgc                                                                                                     >  1:3270189/1‑150 (MQ=255)
                                                aCGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCAGGGCGCCAgaga                                                                                                                                        >  2:3846582/1‑111 (MQ=255)
                                                aCGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCAGGGCGCCAgaga                                                                                                                                        <  1:3846582/114‑4 (MQ=255)
                                                 cGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCcgcgc                                                                                                   <  2:4279971/150‑1 (MQ=255)
                                                 cGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCcgcgc                                                                                                   <  1:1795371/150‑1 (MQ=255)
                                                 cGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATATGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCcgcgc                                                                                                   <  1:109142/150‑1 (MQ=255)
                                                    tGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCa                                                                                                >  1:1385244/1‑150 (MQ=255)
                                                    tGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCa                                                                                                >  1:162067/1‑150 (MQ=255)
                                                    tGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCa                                                                                                >  1:3320292/1‑150 (MQ=255)
                                                      ttCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATg                                                                                              >  2:2722176/1‑150 (MQ=255)
                                                      ttCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATg                                                                                              >  1:1579356/1‑150 (MQ=255)
                                                         gCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATgcgc                                                                                           <  2:1579365/150‑1 (MQ=255)
                                                         gCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATgcgc                                                                                           >  1:3217943/1‑150 (MQ=255)
                                                              cATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTaa                                                                                      <  2:3368905/150‑1 (MQ=255)
                                                               aTGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAAc                                                                                     <  2:318681/150‑1 (MQ=255)
                                                                  ccTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTc                                                                                  <  1:3404964/150‑1 (MQ=255)
                                                                    tCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGa                                                                                >  1:3607070/1‑150 (MQ=255)
                                                                    tCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGa                                                                                >  1:3191125/1‑150 (MQ=255)
                                                                    tCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGa                                                                                <  1:1840741/150‑1 (MQ=255)
                                                                    tCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGa                                                                                <  1:2870851/150‑1 (MQ=255)
                                                                    tCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGa                                                                                >  1:2545567/1‑150 (MQ=255)
                                                                     cAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGAc                                                                               <  2:1156909/150‑1 (MQ=255)
                                                                     cAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGAc                                                                               <  2:3928076/150‑1 (MQ=255)
                                                                     cAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGAc                                                                               <  1:434263/150‑1 (MQ=255)
                                                                      aGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGAcc                                                                              <  2:1040970/150‑1 (MQ=255)
                                                                          gCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGAc                                                                          <  2:2837097/150‑1 (MQ=255)
                                                                           cACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCCGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCAGCGCGCATGCGCCGTAACCTCGACCTGAcc                                                                         >  2:3389314/1‑150 (MQ=255)
                                                                             cGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCa                                                                       <  2:3321536/150‑1 (MQ=255)
                                                                              gAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAg                                                                                                                >  2:1923680/1‑108 (MQ=255)
                                                                              gAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAg                                                                                                                <  1:1923680/108‑1 (MQ=255)
                                                                              gAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCaa                                                                      >  1:1181206/1‑150 (MQ=255)
                                                                              gAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCaa                                                                      <  2:5198392/150‑1 (MQ=255)
                                                                              gAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCaa                                                                      <  1:4770185/150‑1 (MQ=255)
                                                                                   cAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCt                                                                 <  2:3472265/150‑1 (MQ=255)
                                                                                    aGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTg                                                                <  2:4212193/150‑1 (MQ=255)
                                                                                     gCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTgg                                                               <  2:4827330/150‑1 (MQ=255)
                                                                                       cTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTGGTg                                                             >  1:1421812/1‑150 (MQ=255)
                                                                                        tGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAg                                                                                                                                           >  2:4259875/1‑71 (MQ=255)
                                                                                        tGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAg                                                                                                                                           <  1:4259875/71‑1 (MQ=255)
                                                                                            ccTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTGGTGCTGGc                                                        <  1:4734363/150‑1 (MQ=255)
                                                                                              tATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTGGTGCTGGCCg                                                      >  2:2798264/1‑150 (MQ=255)
                                                                                              tATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTGGTGCTGGCCg                                                      >  2:2232583/1‑150 (MQ=255)
                                                                                                tGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTGGTGCTGGCCGaa                                                    <  1:5151550/150‑1 (MQ=255)
                                                                                                   cATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTGGTg                                                             <  2:2948562/138‑1 (MQ=255)
                                                                                                   cATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTGGTg                                                             >  1:2948562/1‑138 (MQ=255)
                                                                                                     tGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTGGTGCTGGCCGAAGCGGt                                               <  2:202883/150‑1 (MQ=255)
                                                                                                       ccGGGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTGGTGCTGGCCGAAGCGGTGa                                             >  2:4544186/1‑150 (MQ=255)
                                                                                                       ccGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTGGTGCTGGCCGAAGCGGTGa                                             <  1:1140210/150‑1 (MQ=255)
                                                                                                       ccGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTGGTGCTGGCCGAAGCGGTGa                                             <  1:1659554/150‑1 (MQ=255)
                                                                                                       ccGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTGGTGCTGGCCGAAGCGGTGa                                             <  1:2670311/150‑1 (MQ=255)
                                                                                                            tGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTGGTGCTGGCCGAAGCGGTGAGCATc                                        <  1:970380/150‑1 (MQ=255)
                                                                                                             gggAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTGGTGCTGGCCGAAGCGGTGAGCATCg                                       <  2:5065903/150‑1 (MQ=255)
                                                                                                                           gATATCTGCTGCCTGTTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTGGTGCTGGCCGAAGCGGTGAGCATCGTCCTCGCCCAGCGt                         >  2:582116/1‑150 (MQ=255)
                                                                                                                            atatCTGCTGCCTGGTATCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTGGTGCTGGCCGAAGCGGTGAGCATCGTCCTCGCCCAGCGTc                        >  2:2491009/1‑150 (MQ=255)
                                                                                                                               tCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTGGTGCTGGCCGAAGCGGTGAGCATCGTCCTCGCCCAGCGTCTgg                     >  2:2561522/1‑150 (MQ=255)
                                                                                                                                  gctgcCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTGGTGCTGGCCGAAGCGGTGAGCATCGTCCTCGCCCAGCGTCTGGGTc                  >  1:5331295/1‑150 (MQ=255)
                                                                                                                                   ctgcCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTGGTGCTGGCCGAAGCGGTGAGCATCGTCCTCGCCCAGCGTCTGGGTcg                 >  2:589761/1‑150 (MQ=255)
                                                                                                                                    tgcCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTGGTGCTGGCCGAAGCGGTGAGCATCGTCCTCGCCCAGCGTCTGGGTcgc                <  2:1421876/150‑1 (MQ=255)
                                                                                                                                      ccTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTGGTGCTGGCCGAAGCGGTGAGCATCGTCCTCGCCCAGCGTCTGGGTCGCGa              >  1:1998928/1‑150 (MQ=255)
                                                                                                                                       cTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTGGTGCTGGCCGAAGCGGTGAGCATCGTCCTCGCCCAGCGTCTGGGTCGCGAc             >  1:4962936/1‑150 (MQ=255)
                                                                                                                                        tGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTGGTGCTGGCCGAAGCGGTGAGCATCGTCCTCGCCCAGCGTCTGGGTCGCGAcc            >  1:635718/1‑150 (MQ=255)
                                                                                                                                         ggTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGCTGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTGGTGCTGGCCGAAGCGGTGAGCATCGTCCTCGCCCAGCGTCTGGGTCGCGACCg           >  2:2390416/1‑150 (MQ=255)
                                                                                                                                           tctcTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTGGTGCTGGCCGAAGCGGTGAGCATCGTCCTCGCCCAGCGTCTGGGTCGCGACCGTg         <  2:4665590/150‑1 (MQ=255)
                                                                                                                                           tctcTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTGGTGCTGGCCGAAGCGGTGAGCATCGTCCTCGCCCAGCGTCTGGGTCGCGACCGTg         <  1:2075246/150‑1 (MQ=255)
                                                                                                                                             tctGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTGGTGCTGGCCGAAGCGGTGAGCATCGTCCTCGCCCAGCGTCTGGGTCGCGACCGTGcc       <  2:2108019/150‑1 (MQ=255)
                                                                                                                                             cctGGCGCAGGGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTGGTGCTGGCCGAAGCGGTGAGCAGCGTCCTCGCCCCGCGTCTGGGTCGCGACCGTGcc       >  2:4456446/2‑150 (MQ=255)
                                                                                                                                               tGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTGGTGCTGGCCGAAGCGGTGAGCATCGTCCTCGCCCAGCGTCTGGGTCGCGACCGTGCcca     >  2:439788/1‑150 (MQ=255)
                                                                                                                                                ggCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTGGTGCTGGCCGAAGCGGTGAGCATCGTCCTCGCCCAGCGTCTGGGTCGCGACCGTGCccac    <  1:3020821/150‑1 (MQ=255)
                                                                                                                                                 gcgcCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTGGTGCTGGCCGAAGCGGTGAGCATCGTCCTCGCCCAGCGTCTGGGTCGCGACCGTGCccacc   <  2:3566932/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GCCGCCGTGTTGATCGGTGCCGCGACCCGCGTGCCGGGCCTGCTGTCGACGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTGGTGCTGGCCGAAGCGGTGAGCATCGTCCTCGCCCAGCGTCTGGGTCGCGACCGTGCCCACCA  >  NC_002947/1572782‑1573077

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: