Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 54,103 G→T 8.9% R200R (CGC→CGA surA ← peptidyl‑prolyl cis‑trans isomerase (PPIase)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00927354,1030GT 8.9% 946.8 / 43.2 258R200R (CGC→CGAsurApeptidyl‑prolyl cis‑trans isomerase (PPIase)
Reads supporting (aligned to +/- strand):  ref base G (118/117);  new base T (17/6);  total (135/123)
Fisher's exact test for biased strand distribution p-value = 4.70e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 275 total aligned reads displayed.

CGGTGCTTAATGCCTGGGCGAAGATCCCGGGCAACTCCTGAATACGGCCCCAGCCCATCTGGCCGCCGTTCAGCGCCTGCTGGTCGGCAGAATGAGCAATCGCCAGCTTACCGAAATCAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACCTGCTGCGCCAGGGATTCGACTTCCTGCGGCAGGATGGT  >  CP009273/53955‑54252
                                                                                                                                                    |                                                                                                                                                     
cgGTGCTTAATGCCTGGGCGAAGATCCCGGGCAACTCCTGAATACGGCCCCAGCCCATCTGGCCGCCGTTCAGCGCCTGCTGGTCGGCAGAATGAGCAATCGCCAGCTTACCGAAATCAGCGCCGTTACGCGCCTGATCGACAATGgcgc                                                                                                                                                      <  2:3732472/150‑1 (MQ=255)
cgGTGCTTAATGCCTGGGCGAAGATCCCGGGCAACTCCTGAATACGGCCCCAGCCCATCTGGCCGCCGTTCAGCGCCTGCTGGTCGGCAGAATGAGCAATCGCCAGCTTACCGAAATCAGCGCCGTTACGCGCCTGATCGACAATGgcgc                                                                                                                                                      >  1:1517689/1‑150 (MQ=255)
   tGCTTAATGCCTGGGCGAAGATCCCGGGCAACTCCTGAATACGGCCCCAGCCCATCTGGCCGCCGTTCAGCGCCTGCTGGTCGGCAGAATGAGCAATCGCCAGCTTACCGAAATCAGCGCCGTTACGCGCCTGATCGACAATGgctcgcg                                                                                                                                                   <  2:310478/150‑1 (MQ=255)
        aaTGCCTGGGCGAAGATCCCGGGCAACTCCTGAATACGGCCCCAGCCCATCTGGCCGCCGTTCAGCGCCTGCTGGTCGGCAGAATGAGCAATCGCCAGCTTACCGAAATCAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTgg                                                                                                                                              <  1:2446784/150‑1 (MQ=255)
        aaTGCCTGGGCGAAGATCCCGGGCAACTCCTGAATACGGCCCCAGCCCATCTGGCCGCCGTTCAGCGCCTGCTGGTCGGCAGAATGAGCAATCGCCAGCTTACCGAAATCAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTgg                                                                                                                                              <  2:1755941/150‑1 (MQ=255)
        aaTGCCTGGGCGAAGATCCCGGGCAACTCCTGAATACGGCCCCAGCCCATCTGGCCGCCGTTCAGCGCCTGCTGGTCGGCAGAATGAGCAATCGCCAGCTTACCGAAATCAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTgg                                                                                                                                              <  2:2488025/150‑1 (MQ=255)
         aTGCCTGGGCGAAGATCCCGGGCAACTCCTGAATACGGCCCCAGCCCATCTGGCCGCCGTTCAGCGCCTGCTGGTCGGCAGAATGAGCAATCGCCAGCTTACCGAAATCAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGc                                                                                                                                             <  1:1119273/150‑1 (MQ=255)
         aTGCCTGGGCGAAGATCCCGGGCAACTCCTGAATACGGCCCCAGCCCATCTGGCCGCCGTTCAGCGCCTGCTGGTCGGCAGAATGAGCAATCGCCAGCTTACCGAAATCAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGc                                                                                                                                             <  1:2326639/150‑1 (MQ=255)
           gCCTGGGCGAAGATCCCGGGCAACTCCTGAATACGGCCCCAGCCCATCTGGCCGCCGTTCAGCGCCTGCTGGTCGGCAGAATGAGCAATCGCCAGCTTACCGAAATCAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCtt                                                                                                                                           >  1:1470912/1‑150 (MQ=255)
                 gCGAAGATCCCGGGCAACTCCTGAATACGGCCCCAGCCCATCTGGCCGCCGTTCAGCGCCTGCTGGTCGGCAGAATGAGCAATCGCCAGCTTACCGAAATCAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCt                                                                                                                                     <  1:162457/150‑1 (MQ=255)
                 gCGAAGATCCCGGGCAACTCCTGAATACGGCCCCAGCCCATCTGGCCGCCGTTCAGCGCCTGCTGGTCGGCAGAATGAGCAATCGCCAGCTTACCGAAATCAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCt                                                                                                                                     <  2:2512959/150‑1 (MQ=255)
                  cGAAGATCCCGGGCAACTCCTGAATACGGCCCCAGCCCATCTGGCCGCCGTTCAGCGCCTGCTGGTCGGCAGAATGAGCAATCGCCAGCTTACCGAAATCAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCtt                                                                                                                                    <  1:1929486/150‑1 (MQ=255)
                  cGAAGATCCCGGGCAACTCCTGAATACGGCCCCAGCCCATCTGGCCGCCGTTCAGCGCCTGCTGGTCGGCAGAATGAGCAATCGCCAGCTTACCGAAATCAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCtt                                                                                                                                    <  1:3188391/150‑1 (MQ=255)
                   gAAGATCCCGGGCAACTCCTGAATACGGCCCCAGCCCATCTGGCCGCCGTTCAGCGCCTGCTGGTCGGCAGAATGAGCAATCGCCAGCTTACCGAAATCAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTc                                                                                                                                   >  2:2367537/1‑150 (MQ=255)
                      gATCCCGGGCAACTCCTGAATACGGCCCCAGCCCATCTGGGCGCCGTTCAGCGCCTGCTGGTCGGCAGAATGAGCAATCGCCAGCTTACCGAAATCAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGtt                                                                                                                                <  2:119410/150‑1 (MQ=255)
                        tCCCGGGCAACTCCTGAATACGGCCCCAGCCCATCTGGCCGCCGTTCAGCGCCTGCTGGTCGGCAGAATGAGCAATCGCCAGCTTACCGAAATCAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCa                                                                                                                              >  2:1205951/1‑150 (MQ=255)
                           cGGGCAACTCCTTAATACGGCCCCCGCCCATCTGGCCGCCGTTCAGCGCCTGCTGGTCGGCAGAATGAGCAATCGCCAGCTTACCGAAATCAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCt                                                                                                                           <  2:2915861/150‑1 (MQ=255)
                           cGGGCAACTCCTGAATACGGCCCCAGCCCATCTGGCCGCCGTTCAGCGCCTGCTGGTCGGCAGAATGAGCAATCGCCAGCTTACCGAAATCAGCGCCGTTACGCGCCTGATCAACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCt                                                                                                                           >  2:2182434/1‑150 (MQ=255)
                            gggCAACTCCTGAATACGGCCCCAGCCCATCTGGCCGCCGTTCAGCGCCTGCTGGTCGGCAGAATGAGCAATCGCCAGCTTACCGAAATCAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTg                                                                                                                          >  2:4253618/1‑150 (MQ=255)
                             ggCAACTCCTGAATACGGCCCCAGCCCATCTGGCCGCCGTTCAGCGCCTGCTGGTCGGCAGAATGAGCAATCGCCAGCTTACCGAAATCAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGa                                                                                                                         <  2:2309042/150‑1 (MQ=255)
                              gCAACTCCTGACTACGGCCCCAGCCCATCTGGCCGCCGTTCAGCGCCTGCTGGTCGGCAGAATGAGCAATCGCCAGCTTACCGAAATCAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGAt                                                                                                                        <  2:3228960/150‑1 (MQ=255)
                               cAACTCCTGAATACGGCCCCAGCCCATCTGGCCGCCGTTCAGCGCCTGCTGGTCGGCAGAATGAGCAATCGCCAGCTTACCGAAATCAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATc                                                                                                                       >  2:958141/1‑150 (MQ=255)
                                aaCTCCTGAATACGGCCCCAGCCCATCTGGCCGCCGTTCAGCGCCTGCTGGTCGGCAGAATGAGCAATCGCCAGCTTACCGAAATCAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCa                                                                                                                      <  2:2764853/150‑1 (MQ=255)
                                  cTCCTGAATACGGCCCCAGCCCATCTGGCCGCCGTTCAGCGCCTGCTGGTCGGCAGAATGAGCAATCGCCAGCTTACCGAAATCAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCaga                                                                                                                    <  2:4525100/150‑1 (MQ=255)
                                   tCCTGAATACGGCCCCAGCCCATCTGGCCGCCGTTCAGCGCCTGCTGGTCGGCAGAATGAGCAATCGCCAGCTTACCGAAATCAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCagag                                                                                                                   <  1:2080318/150‑1 (MQ=255)
                                    ccTGAATACGGCCCCAGCCCATCTGGCCGCCGTTCAGCGCCTGCTGGTCGGCAGAATGAGCAATCGCCAGCTTACCGAAATCAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAgg                                                                                                                  >  1:846907/1‑150 (MQ=255)
                                      tGAATACGGCCCCAGCCCATCTGGCCGCCGTTCAGCGCCTGCTGGTCGGCAGAATGAGCAATCGCCAGCTTACCGAAATCAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTc                                                                                                                >  1:4291216/1‑150 (MQ=255)
                                        aaTACGGCCCCAGCCCATCTGGCCGCCGTTCAGCGCCTGCTGGTCGGCAGAATGAGCAATCGCCAGCTTACCGAAATCAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCgg                                                                                                              <  2:2364930/150‑1 (MQ=255)
                                         aTACGGCCCCAGCCCATCTGGCCGCCGTTCAGCGCCTGCTGGTCGGCAGAATGAGCAATCGCCAGCTTACCGAAATCAGCGCCGTTACGCGCCTGATCGACAATGGCTCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCggg                                                                                                             >  1:694453/1‑150 (MQ=255)
                                            cGGCCCCAGCCCATCTGGCCGCCGTTCAGCGCCTGCTGGTCGGCAGAATGAGCAATCGCCAGCTTACCGAAATCAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGttt                                                                                                          >  2:2965958/1‑150 (MQ=255)
                                            cGGCCCCAGCCCATCTGGCCGCCGTTCAGCGCCTGCTGGTCGGCAGAATGAGCAATCGCCAGCTTACCGAAATCAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGttt                                                                                                          >  1:852777/1‑150 (MQ=255)
                                            cGGCCCCAGCCCATCTGGCCGCCGTTCAGCGCCTGCTGGTCGGCAGAATGAGCAATCGCCAGCTTACCGAAATCAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGttt                                                                                                          >  1:1184583/1‑150 (MQ=255)
                                               ccccAGCCCATCTGGCCGCCGTTCAGCGCCTGCTGGTCGGCAGAATGAGCAATCGCCAGCTTACCGAAATCAGCGCCGTTACGCGCCTGTTCGACAATTGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTcc                                                                                                       <  1:1016180/150‑1 (MQ=255)
                                               ccccAGCCCATCTGGCCGCCGTTCAGCGCCTGCTGGTCGGCAGAATGAGCAATCGCCAGCTTACCGAAATCAGCGCCGTTACGCGCCTGATCGACAATGGCTCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTcc                                                                                                       >  2:910066/1‑150 (MQ=255)
                                                cccAGCCCATCTGGCCGCCGTTCAGCGCCTGCTGGTCGGCAGAATGAGCAATCGCCAGCTTACCGAAATCAGCGCCGTTACGCGCCTGATCGACAATGGCTCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCg                                                                                                      >  1:4370931/1‑150 (MQ=255)
                                                 ccAGCCCATCTGGCCGCCGTTCAGCGCCTGCTGGTCGGCAGAATGAGCAATCGCCAGCTTACCGAAATCAGCGCCGTTACGCGCCTGATCGACAATGGCTCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCgg                                                                                                     >  1:324585/1‑150 (MQ=255)
                                                  cAGCCCATCTGGCCGCCGTTCAGCGCCTGCTGGTCGGCAGAATGAGCAATCGCCAGCTTACCGAAATCAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGc                                                                                                    >  2:1482661/1‑150 (MQ=255)
                                                  cAGCCCATCTGGCCGCCGTTCAGCGCCTGCTGGTCGGCAGAATGAGCAATCGCCAGCTTACCGAAATCAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTGTTCCGGc                                                                                                    >  2:3319009/1‑150 (MQ=255)
                                                      ccATCTGGCCGCCGTTCAGCGCCTGCTGGTCGGCAGAATGAGCAATCGCCAGCTTACCGACATCAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCg                                                                                                >  2:1835275/1‑150 (MQ=255)
                                                      ccATCTGGCCGCCGTTCAGCGCCTGCTGGTCGGCAGAATGAGCAATCGCCAGCTTACCGAAATCAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCg                                                                                                >  2:3542116/1‑150 (MQ=255)
                                                       cATCTGGCCGCCGTTCAGCGCCTGCTGGTCGGCAGAATGAGCAATCGCCAGCTTACCGAAATCAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCgg                                                                                               <  1:4130994/150‑1 (MQ=255)
                                                               cgccgTTCAGCGCCTGCTGGTCGGCAGAATGAGCAATCGCCAGCTTACCGAAATCAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCGGGATCAGGa                                                                                       >  1:573570/1‑150 (MQ=255)
                                                                  cgTTCAGCGCCTGCTGGTCGGCAGAATGAGCAATCGCCAGCTTACCGAAATCAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCGGGATCAGGAtgt                                                                                    >  2:2584602/1‑150 (MQ=255)
                                                                  cgTTCAGCGCCTGCTGGTCGGCAGAATGAGCAATCGCCAGCTTACCGAAATCAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCAGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCGGGATCAGGAtgt                                                                                    >  2:3136547/1‑150 (MQ=255)
                                                                    ttCAGCGCCTGCTGGTCGGCAGAATGAGCAATCGCCAGCTTACCGAAATCAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCGGGATCAGGATGTgg                                                                                  <  1:377704/150‑1 (MQ=255)
                                                                          gcCTGCTGGTCGGCAGAATGAGCAATCGCCAGCTTACCGAAATCAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCGGGATCAGGATGTGGCTCAgg                                                                            >  2:1444375/1‑150 (MQ=255)
                                                                          gcCTGCTGGTCGGCAGAATGAGCAATCGCCAGCTTACCGAAATCAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCGGGATCAGGATGTGGCTCAgg                                                                            <  1:1609140/150‑1 (MQ=255)
                                                                            ctgctgGTCGGCAGAATGAGCAATCGCCAGCTTACCGAAATCAGCGCCGTTACGCGCCTGATCGACAATGGCTCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCGGGATCAGGATGTGGCTCAGGtt                                                                          >  1:2880743/1‑150 (MQ=255)
                                                                            ctgctgGTCGGCAGAATGAGCAATCGCCAGCTTACCGAAATCAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCGGGATCAGGATGTGGCTCAGGtt                                                                          >  1:4206570/1‑150 (MQ=255)
                                                                            ctgctgGTCGGCAGAATGAGCAATCGCCAGCTTACCGAAATCAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCGGGATCAGGATGTGGCTCAGGtt                                                                          <  2:860087/150‑1 (MQ=255)
                                                                              gctgGTCGGCAGAATGAGCAATCGCCAGCTTACCGAAATCAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCGGGATCAGGATGTGGCTCAGGTTCa                                                                        <  2:559854/150‑1 (MQ=255)
                                                                               ctgGTCGGCAGAATGAGCAATCGCCAGCTTACCGAAATCAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCGGGATCAGGATGTGGCTCAGGTTCAg                                                                       >  1:1058780/1‑150 (MQ=255)
                                                                                  gTCGGCAGAATGAGCAATCGCCAGCTTACCGAAATCAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCGGGATCAGGATGTGGCTCAGGTTCAGCTc                                                                    >  1:1561572/1‑150 (MQ=255)
                                                                                  gTCGGCAGAATGAGCAATCGCCAGCTTACCGAAATCAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCGGGATCAGGATGTGGCTCAGGTTCAGCTc                                                                    >  1:3551922/1‑150 (MQ=255)
                                                                                    cGGCAGAATGAGCAATCGCCAGCTTACCGAAATCAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCAGGATCAGGATGTGGCTCAGGTTCAGCTCAg                                                                  >  1:2004931/1‑150 (MQ=255)
                                                                                       cAGAATGAGCAATCGCCAGCTTACCGAAATCAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGc                                                               >  2:4064211/1‑150 (MQ=255)
                                                                                       cAGAATGAGCAATCGCCAGCTTACCGAAATCAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGc                                                               >  2:1598008/1‑150 (MQ=255)
                                                                                        aGAATGAGCAATCGCCAGCTTACCGAAATCAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCt                                                              >  1:3330073/1‑150 (MQ=255)
                                                                                           aTGAGCAATCGCCAGCTTACCGAAATCAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGc                                                           <  1:2421385/150‑1 (MQ=255)
                                                                                             gAGCAATCGCCAGCTTACCGAAATCAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGcgt                                                         >  1:2104583/1‑150 (MQ=255)
                                                                                             gAGCAATCGCCAGCTTACCGAAATCAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGcgt                                                         >  1:4349245/1‑150 (MQ=255)
                                                                                              aGCAATCGCCAGCTTACCGAAATCAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGcgtc                                                        >  2:371713/1‑150 (MQ=255)
                                                                                              aGCAATCGCCAGCTTACCGAAATCAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGcgtc                                                        <  1:2367529/150‑1 (MQ=255)
                                                                                                cAATCGCCAGCTTACCGAAATCAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGcgtcgt                                                      >  2:1346121/1‑150 (MQ=255)
                                                                                                cAATCGCCAGCTTACCGAAATCAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGcgtcgt                                                      >  1:221535/1‑150 (MQ=255)
                                                                                                 aaTCGCCAGCTTACCGAAATCAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGtt                                                     <  1:4311636/150‑1 (MQ=255)
                                                                                                  aTCGCCAGCTTACCGAAATCAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGttt                                                    >  1:265859/1‑150 (MQ=255)
                                                                                                    cGCCAGCTTACCGAAATCAGCGCCGTTACGCGCCTGATCGACAATGGCTCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTg                                                  >  2:776830/1‑150 (MQ=255)
                                                                                                    cGCCAGCTTACCGAAATCAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTg                                                  <  2:1967053/150‑1 (MQ=255)
                                                                                                    cGCCAGCTTACCGAAATCAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTg                                                  <  2:3055829/150‑1 (MQ=255)
                                                                                                     gCCAGCTTACCGAAATCAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTgg                                                 <  2:3780610/150‑1 (MQ=255)
                                                                                                       cAGCTTACCGAAATCAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGtt                                               <  2:846907/150‑1 (MQ=255)
                                                                                                       cAGCTTACCGAAATCAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGtt                                               <  2:1339508/150‑1 (MQ=255)
                                                                                                          cTTACCGAAATCAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTAcc                                            <  2:34597/150‑1 (MQ=255)
                                                                                                          cTTACCGAAATCAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTAcc                                            >  1:375957/1‑150 (MQ=255)
                                                                                                              ccGAAATCAGCGCCGTTACGCGCCTGATCGACAATGGCTCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCAcc                                        <  2:4322659/150‑1 (MQ=255)
                                                                                                                 aaaTCAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACctgc                                     >  2:3827147/1‑150 (MQ=255)
                                                                                                                  aaTCAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACctgct                                    >  2:2142319/1‑150 (MQ=255)
                                                                                                                     cAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACCTGCTgcg                                 <  2:4495973/150‑1 (MQ=255)
                                                                                                                     cAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACCTGCTgcg                                 <  1:1717424/150‑1 (MQ=255)
                                                                                                                      aGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACCTGCTgcgc                                >  2:1048597/1‑150 (MQ=255)
                                                                                                                        cgcCGTTACGCGCCTGAGCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACCTGCTGCGCCa                              <  2:2916972/150‑1 (MQ=255)
                                                                                                                             ttACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACCTGCTGCGCCAGGGAt                         >  1:3301122/1‑150 (MQ=255)
                                                                                                                             ttACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACCTGCTGCGCCAGGGAt                         <  1:2984154/150‑1 (MQ=255)
                                                                                                                              tACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACCTGCTGCGCCAGGGAtt                        >  2:4080302/1‑150 (MQ=255)
                                                                                                                                   gcCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACCTGCTGCGCCAGGGATTCGACt                   >  2:657480/1‑150 (MQ=255)
                                                                                                                                   gcCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACCTGCTGCGCCAGGGATTCGACt                   <  1:1835176/150‑1 (MQ=255)
                                                                                                                                   gcCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACCTGCTGCGCCAGGGATTCGACt                   <  2:3299966/150‑1 (MQ=255)
                                                                                                                                    ccTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACCTGCTGCGCCAGGGATTCGACtt                  <  2:33262/150‑1 (MQ=255)
                                                                                                                                       gATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACCTGCTGCGCCAGGGATTCGACTTCCt               <  1:3713901/150‑1 (MQ=255)
                                                                                                                                        aTCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACCTGCTGCGCCAGGGATTCGACTTCCTg              <  2:3690062/150‑1 (MQ=255)
                                                                                                                                        aTCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACCTGCTGCGCCAGGGATTCGACTTCCTg              <  2:3122774/150‑1 (MQ=255)
                                                                                                                                         tCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACCTGCTGCGCCAGGGATTCGACTTCCTGc             >  1:2517086/1‑150 (MQ=255)
                                                                                                                                         tCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACCTGCTGCGCCAGGGATTCGACTTCCTGc             >  1:2803620/1‑150 (MQ=255)
                                                                                                                                          cGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACCTGCTGCGCCAGGGATTCGACTTCCTGCg            >  1:2351102/1‑150 (MQ=255)
                                                                                                                                          cGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACCTGCTGCGCCAGGGATTCGACTTCCTGCg            <  1:762081/150‑1 (MQ=255)
                                                                                                                                             cAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACCTGCTGCGCCAGGGATTCGACTTCCTGCGGCa         <  1:2965958/150‑1 (MQ=255)
                                                                                                                                                tGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACCTGCTGCGCCAGGGATTCGACTTCCTGCGGCAGGa      >  2:1192944/1‑150 (MQ=255)
                                                                                                                                                  gcgcgcgcCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACCTGCTGCGCCAGGGATTCGACTTCCTGCGGCAGGATg    <  2:3249282/150‑1 (MQ=255)
                                                                                                                                                   cgcgcgcCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACCTGCTGCGCCAGGGATTCGACTTCCTGCGGCAGGATgg   <  1:2363507/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CGGTGCTTAATGCCTGGGCGAAGATCCCGGGCAACTCCTGAATACGGCCCCAGCCCATCTGGCCGCCGTTCAGCGCCTGCTGGTCGGCAGAATGAGCAATCGCCAGCTTACCGAAATCAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACCTGCTGCGCCAGGGATTCGACTTCCTGCGGCAGGATGGT  >  CP009273/53955‑54252

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: