Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,549,778 C→A 12.7% M293I (ATG→ATT eutC ← ethanolamine ammonia‑lyase, small subunit (light chain)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,549,7780CA12.7% 968.4 / 70.1 277M293I (ATG→ATTeutCethanolamine ammonia‑lyase, small subunit (light chain)
Reads supporting (aligned to +/- strand):  ref base C (132/109);  new base A (21/15);  total (153/124)
Fisher's exact test for biased strand distribution p-value = 7.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 280 total aligned reads displayed.

GAGTCGGCAGAAATCAGACCGAGGCTACGAATATGCGGCGGCAATTTCAGTTCACGCGCAAAATCGGCGTTAACGGAGGCAATCACCCGCATAGCGGTGACAGACGGTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTGCGATCGGCCTCAACGGTGGTCGCCATACGCGGCGAGTA  >  CP009273/2549629‑2549927
                                                                                                                                                     |                                                                                                                                                     
  gTCGGCAGAAATCAGACCGAGGCTACGAATATGCGGCGGCAATTTAAGTTCACGCGCAAAATCGGCGTTAACGGAGGCAATCACCCGCATAGCGGTGACAGACGGTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTCAt                                                                                                                                                     <  2:418562/150‑1 (MQ=255)
    cGGCAGAAATCAGACCGAGGCTACGAATATGCGGCGGCAATTTAAGTTCACGCGCAAAATCGGCGTTAACGGAGGCAATCACCCGCATAGCGGTGACAGACGGTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGt                                                                                                                                                   >  1:1948365/1‑150 (MQ=255)
    cGGCAGAAATCAGACCGAGGCTACGAATATGCGGCGGCAATTTAAGTTCACGCGCAAAATCGGCGTTAACGGAGGCAATCACCCGCATAGCGGTGACAGACGGTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGt                                                                                                                                                   >  1:2707784/1‑150 (MQ=255)
     ggCAGAAATCAGACCGAGGCTACGAATATGCGGCGGCAATTTAAGTTCACGCGCAAAATCGGCGTAAACGGAGGCAATCACCCGCATAGCGGTGACAGACGGTCGAATCAAATCTAAAGCTGGCATGATGCATCCTTAGCGGGTCATGtt                                                                                                                                                  >  1:3310453/1‑150 (MQ=255)
        aGAAATCAGACCGAGGCTACGAATATGCGGCGGCAATTTAAGTTCACGCGCAAAATCGGCGTTAACGGAGGCAATCACCCGCATAGCGGTGACAGACGGTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGAt                                                                                                                                               <  1:1994275/150‑1 (MQ=255)
         gAAATCAGACCGAGGCTACGAATATGCGGCGGCAATTTCAGTTCACGCGCAAAATCGGCGTTAACGGAGGCAATCACCCGCATAGCGGTGACAGACGGTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATg                                                                                                                                              >  2:4397960/1‑150 (MQ=255)
           aaTCAGACCGAGGCTACGAATATGCGGCGGCAATTTAAGTTCACGCGCAAAATCGGCGTTAACGGAGGCAATCACCCGCATAGCGGTGACAGACGGTCGAATCAAATCTAAAGCTGGCACGATGCCTCCTTATCGGGTCATGTTGATGcc                                                                                                                                            <  1:2380571/150‑1 (MQ=255)
               aGACCGAGGCTACGAATATGCGGCGGCAATTTAAGTTCACGCGCAAAATCGGCGTTAACGGAGGCAATCACCCGCATAGCGGTGACAGACGGTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGAt                                                                                                                                        >  2:1154627/1‑150 (MQ=255)
                   cGAGGCTACGAATATGCGGCGGCAATTTAAGTTCACGCGCAAAATCGGCGTTAACGGAGGCAATCACCCGCATAGCGGTGACAGACGGTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCtt                                                                                                                                    <  2:4243732/150‑1 (MQ=255)
                    gAGGCTACGAATATGCGGCGGCAATTTAAGTTCACGCGCAAAATCGGCGTTAACGGAGGCAATCACCCGCATAGCGGTGACAGACGGTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCttt                                                                                                                                   <  2:1306740/150‑1 (MQ=255)
                    gAGGCTACGAATATGCGGCGGCAATTTAAGTTCACGCGCAAAATCGGCGTTAACGGAGGCAATCACCCGCATAGCGGTGACAGACGGTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCttt                                                                                                                                   >  2:2653746/1‑150 (MQ=255)
                     aGGCTACGAATATGCGGCGGCAATTTAAGTTCACGCGCAAAATCGGCGTTAACGGAGGCAATCACCCGCATAGCGGTGACAGACGGTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTAATGTTGATGCCGGATGCTTTc                                                                                                                                  <  2:2519623/150‑1 (MQ=255)
                       gCTACGAATATGCGGCGGCAATTTAAGTTCACGCGCAAAATCGGCGTTAACGGAGGCAATCACCCGCATAGCGGTGACAGACGGTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGcc                                                                                                                                            <  1:1023190/138‑1 (MQ=255)
                       gCTACGAATATGCGGCGGCAATTTAAGTTCACGCGCAAAATCGGCGTTAACGGAGGCAATCACCCGCATAGCGGTGACAGACGGTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTg                                                                                                                                >  2:3485188/1‑150 (MQ=255)
                         tACGAATATGCGGCGGCAATTTAAGTTCACGCGCAAAATCGGCGTTAACGGAGGCAATCACCCGCATAGCGGTGACAGACGGTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCt                                                                                                                              <  1:2617949/150‑1 (MQ=255)
                         tACGAATATGCGGCGGCAATTTAAGTTCACGCGCAAAATCGGCGTTAACGGAGGCAATCACCCGCATAGCGGTGACAGACGGTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTAATGTTGATGCCGGATGCTTTCTGCt                                                                                                                              >  2:1024591/1‑150 (MQ=255)
                               tatGCGGCGGCAATTTAAGTTCACGCGCAAAATCGGCGTTAACGGAGGCAATCACCCGCATAGCGGTGACAGACGGTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCa                                                                                                                        >  2:360694/1‑150 (MQ=255)
                                 tGCGGCGGCAATTTAAGTTCACGCGCAAAATCGGCGTTAACGGAGGCAATCACCCGCATAGCGGTGACAGACGGTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATa                                                                                                                      <  2:4643957/150‑1 (MQ=255)
                                  gcggcgGCAATTTAAGTTCACGCGCAAAATCGGCGTTAACGGAGGCAATCACCCGCATAGCGGTGACAGACGGTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATAc                                                                                                                     <  1:3578259/150‑1 (MQ=255)
                                    ggcggcAATTTAAGTTCACGCGCAAAATCGGCGTTAACGGAGGCAATCACCCGCATAGCGGTGACAGACGGTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGt                                                                                                                   >  1:3703741/1‑150 (MQ=255)
                                    ggcggcAATTTAAGTTCACGCGCAAAATCGGCGTTAACGGAGGCAATCACCCGCATAGCGGTGACAGACGGTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGt                                                                                                                   >  1:3703717/1‑150 (MQ=255)
                                    ggcggcAATTTAAGTTCACGCGCAAAATCGGCGTTAACGGAGGCAATCACCCGCATAGCGGTGACAGACGGTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGt                                                                                                                   >  2:1334322/1‑150 (MQ=255)
                                     gcggcAATTTAAGTTCACGCGCAAAATCGGCGTTAACGGAGGCAATCACCCGCATAGCGGTGACAGACGGTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTAATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGtt                                                                                                                  <  2:845092/150‑1 (MQ=255)
                                         cAATTTAAGTTCACGCGCAAAATCGGCGTTAACGGAGGCAATCACCCGCATAGCGGTGACAGACGGTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGc                                                                                                              <  2:3012027/150‑1 (MQ=255)
                                          aaTTTAAGTTCACGCGCAAAATCGGCGTTAACGGAGGCAATCACCCGCATAGCGGTGACAGACGGTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGcc                                                                                                             >  2:3848373/1‑150 (MQ=255)
                                          aaTTTAAGTTCACGCGCAAAATCGGCGTTAACGGAGGCAATCACCCGCATAGCGGTGACAGACGGTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGcc                                                                                                             >  1:2188714/1‑150 (MQ=255)
                                          aaTTTAAGTTCACGCGCAAAATCGGCGTTAACGGAGGCAATCACCCGCATAGCGGTGACAGACGGTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGcc                                                                                                             >  1:2947478/1‑150 (MQ=255)
                                          aaTTTAAGTTCACGCGCAAAATCGGCGTTAACGGAGGCAATCACCCGCATAGCGGTGACAGACGGTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTAATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGcc                                                                                                             >  1:2258140/1‑150 (MQ=255)
                                          aaTTTAAGTTCACGCGCAAAATCGGCGTTAACGGAGGCAATCACCCGCATAGCGGTGACAGACGGTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTAATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGcc                                                                                                             >  1:207968/1‑150 (MQ=255)
                                           aTTTAAGTTCACGCGCAAAATCGGCGTTAACGGAGGCAATCACCCGCATAGCGGTGACAGACGGTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCa                                                                                                            <  1:4003127/150‑1 (MQ=255)
                                              taaGTTCACGCGCAAAATCGGCGTTAACGGAGGCAATCACCCGCATAGCGGTGACAGACGGTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTAATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAAt                                                                                                         >  1:1730458/3‑150 (MQ=255)
                                              taaGTTCACGCGCAAAATCGGCGTTAACGGAGGCAATCACCCGCATAGCGGTGACAGACGGTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTAATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAAt                                                                                                         >  2:4142600/3‑150 (MQ=255)
                                                aGTTCACGCGCAAAATCGGCGTTAACGGAGGCAATCACCCGCATAGCGGTGACAGACGGTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATcc                                                                                                       <  2:2166346/150‑1 (MQ=255)
                                                aGTTCACGCGCAAAATCGGCGTTAACGGAGGCAATCACCCGCATAGCGGTGACAGACGGTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTAATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATcc                                                                                                       >  1:1233949/1‑150 (MQ=255)
                                                 gTTCACGCGCAAAATCGGCGTTAACGGAGGCAATCACCCGCATAGCGGTGACAGACGGTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTAATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCca                                                                                                      >  2:22788/1‑150 (MQ=255)
                                                           aaaaTCGGCGTTAACGGAGGCAATCACCCGCATAGCGGTGACAGACGGTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGc                                                                                            >  2:1412608/1‑150 (MQ=255)
                                                           aaaaTCGGCGTTAACGGAGGCAATCACCCGCATAGCGGTGACAGACGGTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGc                                                                                            <  2:101502/150‑1 (MQ=255)
                                                            aaaTCGGCGTTAACGGAGGCAATCACCCGCATAGCGGTGACAGACGGTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCt                                                                                           <  1:555919/150‑1 (MQ=255)
                                                            aaaTCGGCGTTAACGGAGGCAATCACCCGCATAGCGGTGACAGACGGTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCt                                                                                           <  2:3730803/150‑1 (MQ=255)
                                                             aaTCGGCGTTAACGGAGGCAATCACCCGCATAGCGGTGACAGACGGTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTAATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTg                                                                                          <  2:2364952/150‑1 (MQ=255)
                                                              aTCGGCGTTAACGGAGGCAATCACCCGCATAGCGGTGACAGACGGTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGc                                                                                         >  1:3100837/1‑150 (MQ=255)
                                                              aTCGGCGTTAACGGAGGCAATCACCCGCATAGCGGTGACAGACGGTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTAATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGc                                                                                         <  1:263607/150‑1 (MQ=255)
                                                               tCGGCGTTAACGGAGGCAATCACCCGCATAGCGGTGACAGACGGTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCg                                                                                        <  1:4054168/150‑1 (MQ=255)
                                                                    gTTAACGGAGGCAATCACCCGCATAGCGGTGACAGACGGTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTAATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTc                                                                                   >  2:4602933/1‑150 (MQ=255)
                                                                      tAACGGAGGCAATCACCCGCATAGCGGTGACAGACGGTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTAATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCaa                                                                                 >  2:3944137/1‑150 (MQ=255)
                                                                        aCGGAGGCAATCACCCGCATAGCGGTGACAGACGGTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAAcc                                                                               <  2:2449517/150‑1 (MQ=255)
                                                                          ggAGGCAATCACCCGCATAGCGGTGACAGACGGTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACcgg                                                                             >  2:3242320/1‑150 (MQ=255)
                                                                             ggCAATCACCCGCATAGCGGTGACAGACGGTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACcggcgg                                                                          >  1:538160/1‑150 (MQ=255)
                                                                              gCAATCACCCGCATAGCGGTGACAGACGGTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCggcggc                                                                         <  1:1082105/150‑1 (MQ=255)
                                                                               cAATCACCCGCATAGCGGTGACAGACGGTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTAATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGgcggcg                                                                        <  2:4110396/150‑1 (MQ=255)
                                                                                aaTCACCCGCATAGCGGTGACAGACGGTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGt                                                                       >  2:3433033/1‑150 (MQ=255)
                                                                                  tCACCCGCATAGCGGTGACAGACGGTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTgc                                                                     <  1:2004903/150‑1 (MQ=255)
                                                                                  tCACCCGCATAGCGGTGACAGACGGTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTgc                                                                     >  1:2366532/1‑150 (MQ=255)
                                                                                        gCATAGCGGTGACAGACGGTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCt                                                               <  1:2786755/150‑1 (MQ=255)
                                                                                        gCATAGCGGTGACAGACGGTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCt                                                               <  1:2786979/150‑1 (MQ=255)
                                                                                        gCATAGCGGTGACAGACGGTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCt                                                               <  2:2458726/150‑1 (MQ=255)
                                                                                        gCATAGCGGTGACAGACGGTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTAATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCt                                                               <  2:3151009/150‑1 (MQ=255)
                                                                                        gCATAGCGGTGACAGACGGTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTAATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCt                                                               <  2:2801352/150‑1 (MQ=255)
                                                                                         cATAGCGGTGACAGACGGTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTg                                                              <  2:1765072/150‑1 (MQ=255)
                                                                                            aGCGGTGACAGACGGTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTg                                                           <  1:4644299/150‑1 (MQ=255)
                                                                                               ggTGACAGACGGTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAAt                                                        >  2:1651251/1‑150 (MQ=255)
                                                                                                gTGACAGACGGTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATg                                                       >  1:682714/1‑150 (MQ=255)
                                                                                                gTGACAGACGGTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATg                                                       >  1:683419/1‑150 (MQ=255)
                                                                                                     aGACGGTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTaga                                                  >  2:4532501/1‑150 (MQ=255)
                                                                                                      gACGGTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTagag                                                 >  2:3910751/1‑150 (MQ=255)
                                                                                                       aCGGTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAgaga                                                >  2:240824/1‑150 (MQ=255)
                                                                                                        cGGTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGAt                                               <  2:1339397/150‑1 (MQ=255)
                                                                                                        cGGTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGAt                                               >  2:3194687/1‑150 (MQ=255)
                                                                                                         ggTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTAATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATa                                              <  2:1989841/150‑1 (MQ=255)
                                                                                                         ggTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTAATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATa                                              <  2:2894867/150‑1 (MQ=255)
                                                                                                           tCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACa                                            <  2:2116189/150‑1 (MQ=255)
                                                                                                             gAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAgg                                          >  1:2469549/1‑150 (MQ=255)
                                                                                                             gAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAgg                                          >  2:2610073/1‑150 (MQ=255)
                                                                                                             gAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAgg                                          >  1:3276260/1‑150 (MQ=255)
                                                                                                             gAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAgg                                          <  2:4567049/150‑1 (MQ=255)
                                                                                                               aTCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTg                                        <  1:4317500/150‑1 (MQ=255)
                                                                                                               aTCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTg                                        >  1:1577525/1‑150 (MQ=255)
                                                                                                               aTCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTg                                        <  1:375973/150‑1 (MQ=255)
                                                                                                               aTCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTg                                        <  2:3561144/150‑1 (MQ=255)
                                                                                                                  aaaTCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTGCGa                                     >  2:3010211/1‑150 (MQ=255)
                                                                                                                  aaaTCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTGCGa                                     <  1:4074005/150‑1 (MQ=255)
                                                                                                                  aaaTCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTGCGa                                     >  2:689397/1‑150 (MQ=255)
                                                                                                                    aTCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTGCGATc                                   <  1:4000823/150‑1 (MQ=255)
                                                                                                                         aaGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTGCGATCGGCCt                              >  1:3393507/1‑150 (MQ=255)
                                                                                                                         aaGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTGCGATCGGCCt                              >  1:3494426/1‑150 (MQ=255)
                                                                                                                          aGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTGCGATCGGCCTc                             >  1:4596935/1‑150 (MQ=255)
                                                                                                                            cTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTGCGATCGGCCTCaa                           >  1:4273625/1‑150 (MQ=255)
                                                                                                                              ggCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTGCGATCGGCCTCAACg                         <  1:4397960/150‑1 (MQ=255)
                                                                                                                              ggCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTGCGATCGGCCTCAACg                         <  1:1306322/150‑1 (MQ=255)
                                                                                                                                cATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTGCGATCGGCCTCAACggt                       >  2:524839/1‑150 (MQ=255)
                                                                                                                                      gCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTGCGATCGGCCTCAACGGTGGTCGc                 >  2:1214270/1‑150 (MQ=255)
                                                                                                                                      gCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTGCGATCGGCCTCAACGGTGGTCGc                 >  2:1502247/1‑150 (MQ=255)
                                                                                                                                       cctcctTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTGCGATCGGCCTCAACGGTGGTCGcc                >  1:4222989/1‑150 (MQ=255)
                                                                                                                                        ctcctTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTGCGATCGGCCTCAACGGTGGTCGCCa               <  1:1391727/150‑1 (MQ=255)
                                                                                                                                        ctcctTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTGCGATCGGCCTCAACGGTGGTCGCCa               >  1:3290816/1‑150 (MQ=255)
                                                                                                                                                cGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTGCGATCGGCCTCAACGGTGGTCGCCATACgcggc       <  2:2744065/150‑1 (MQ=255)
                                                                                                                                                 gggTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTGCGATCGGCCTCAACGGTGGTCGCCATACgcggcg      >  1:587214/1‑150 (MQ=255)
                                                                                                                                                 gggTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTGCGATCGGCCTCAACGGTGGTCGCCATACgcggcg      <  1:3890565/150‑1 (MQ=255)
                                                                                                                                                     cATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTGCGATCGGCCTCAACGGTGGTCGCCATACGCGGCGAGTa  >  1:1478440/1‑150 (MQ=255)
                                                                                                                                                     cATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTGCGATCGGCCTCAACGGTGGTCGCCATACGCGGCGAGTa  >  1:411446/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GAGTCGGCAGAAATCAGACCGAGGCTACGAATATGCGGCGGCAATTTCAGTTCACGCGCAAAATCGGCGTTAACGGAGGCAATCACCCGCATAGCGGTGACAGACGGTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTGCGATCGGCCTCAACGGTGGTCGCCATACGCGGCGAGTA  >  CP009273/2549629‑2549927

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: