Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,808,646 G→T 17.3% R385S (CGC→AGC)  gshA ← glutamate‑cysteine ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,808,6460GT17.3% 794.4 / 103.9 265R385S (CGC→AGC) gshAglutamate‑cysteine ligase
Reads supporting (aligned to +/- strand):  ref base G (108/110);  new base T (26/21);  total (134/131)
Fisher's exact test for biased strand distribution p-value = 5.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 269 total aligned reads displayed.

CACACACTTTCTGATACGCTTCGCCGCCGTTAATACTATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGAATCGCACCTGCTGTTCATCTACACCAATCGGCGAGAACGGGT  >  CP009273/2808499‑2808795
                                                                                                                                                   |                                                                                                                                                     
   acacTTTCTGATACGCTTCGCCGCCGTTAATACTATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCtt                                                                                                                                                  <  1:1093294/150‑1 (MQ=255)
   acacTTTCTGATACGCTTCGCCGCCGTTAATACTATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCtt                                                                                                                                                  >  1:1468888/1‑150 (MQ=255)
     acTTTCTGATACGCTTCGCCGCCGTTAATACTATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCCGCGTCAGACCCGGTTTGCGACCTTCg                                                                                                                                                >  1:1080475/1‑150 (MQ=255)
       tttCTGATACGCTTCGCCGCCGTTAATACTATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAg                                                                                                                                              <  1:2783123/150‑1 (MQ=255)
       tgtCTGATACGCTTCGCCGCCGTTAATACTATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGAGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAg                                                                                                                                              <  1:4126037/148‑1 (MQ=255)
          cTGATACGCTTCGCCGCCGTTAATACTATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGAt                                                                                                                                           <  1:2020014/150‑1 (MQ=255)
           tGATACGCTTCGCCGCCGTTAATACTATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATc                                                                                                                                          >  1:2644770/1‑150 (MQ=255)
           tGATACGCTTCGCCGCCGTTAATACTATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATc                                                                                                                                          >  1:4343726/1‑150 (MQ=255)
            gATACGCTTCGCCGCCGTTAATACTATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCa                                                                                                                                         >  2:1100769/1‑150 (MQ=255)
            gATACGCTTCGCCGCCGTTAATACTATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCa                                                                                                                                         >  2:2189980/1‑150 (MQ=255)
                  cTTCGCCGCCGTTAATACTATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGt                                                                                                                                   <  2:3348026/150‑1 (MQ=255)
                  cTTCGCCGCCGTTAATACTATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGt                                                                                                                                   <  2:1153779/150‑1 (MQ=255)
                   ttCGCCGCCGTTAATACTATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACTGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGATGATCACCCGGtt                                                                                                                                  <  2:239569/150‑1 (MQ=255)
                    tCGCCGCCGTTAATACTATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTAGAGGATCACCCGGTTc                                                                                                                                 >  1:547726/1‑150 (MQ=255)
                     cgccgcCGTTAATACTATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTcc                                                                                                                                <  1:1657777/150‑1 (MQ=255)
                      gccgccGTTAATACTATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCa                                                                                                                               <  1:46771/150‑1 (MQ=255)
                      gccgccGTTAATACTATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCa                                                                                                                               <  2:433461/150‑1 (MQ=255)
                         gccgTTAATACTATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGtt                                                                                                                            <  2:1300927/150‑1 (MQ=255)
                         gccgTTAATACTATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGtt                                                                                                                            <  2:2755898/150‑1 (MQ=255)
                          ccgTTAATACTATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTa                                                                                                                           >  2:3417655/1‑150 (MQ=255)
                            gTTAATACTATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAAc                                                                                                                         >  1:141166/1‑150 (MQ=255)
                              tAATACTATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAAcgc                                                                                                                       <  1:68014/150‑1 (MQ=255)
                               aaTACTATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAAcgcg                                                                                                                      <  1:3940632/150‑1 (MQ=255)
                               aaTACTATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAAcgcg                                                                                                                      <  1:2389630/150‑1 (MQ=255)
                                 tACTATCCCGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCgtg                                                                                                                    <  2:2843881/150‑1 (MQ=255)
                                   cTATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTa                                                                                                                  <  1:4340446/150‑1 (MQ=255)
                                     aTCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACa                                                                                                                >  2:1691329/1‑150 (MQ=255)
                                          gCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCCCCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCaa                                                                                                           <  2:125498/150‑1 (MQ=255)
                                           cGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAg                                                                                                          >  2:4282215/1‑150 (MQ=255)
                                           cGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAg                                                                                                          <  2:2973681/150‑1 (MQ=255)
                                             tttGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCTACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGtt                                                                                                        <  2:2222194/150‑1 (MQ=255)
                                              ttGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTc                                                                                                       <  2:1895340/150‑1 (MQ=255)
                                              ttGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTc                                                                                                       <  1:2597987/150‑1 (MQ=255)
                                               tGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCg                                                                                                      >  1:2853564/1‑150 (MQ=255)
                                                 cgcgACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCt                                                                                                    <  2:1936947/150‑1 (MQ=255)
                                                 cgcgACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCt                                                                                                    >  1:2806240/1‑150 (MQ=255)
                                                 cgcgACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCt                                                                                                    <  2:3755682/150‑1 (MQ=255)
                                                  gcgACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTa                                                                                                   <  2:2926437/150‑1 (MQ=255)
                                                   cgACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCACGTTCGCTAc                                                                                                  >  2:2449532/1‑150 (MQ=255)
                                                       gcgTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCt                                                                                              <  2:1633210/150‑1 (MQ=255)
                                                             cAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCTACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTc                                                                                        >  1:2592681/1‑150 (MQ=255)
                                                             cAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTc                                                                                        >  2:1332527/1‑150 (MQ=255)
                                                                aTCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCgg                                                                                     >  2:2277374/1‑150 (MQ=255)
                                                                 tCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGt                                                                                    <  1:3482511/150‑1 (MQ=255)
                                                                     ggAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCAt                                                                                <  2:4128799/150‑1 (MQ=255)
                                                                       aaCAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATca                                                                              <  1:3292184/150‑1 (MQ=255)
                                                                        aCAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATcag                                                                             <  2:1051578/150‑1 (MQ=255)
                                                                        aCAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATcag                                                                             >  1:3716884/1‑150 (MQ=255)
                                                                        aCAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATcag                                                                             >  1:126430/1‑150 (MQ=255)
                                                                           gATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATcagcca                                                                          <  2:933865/150‑1 (MQ=255)
                                                                            aTCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATcagccag                                                                         >  1:3198504/1‑150 (MQ=255)
                                                                            aTCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATcagccag                                                                         >  1:3198677/1‑150 (MQ=255)
                                                                                 tACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGcac                                                                    >  2:481107/1‑150 (MQ=255)
                                                                                 tACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGcac                                                                    >  1:1848221/1‑150 (MQ=255)
                                                                                 tACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGcac                                                                    >  2:2653770/1‑150 (MQ=255)
                                                                                 tACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGcac                                                                    >  2:1838314/1‑150 (MQ=255)
                                                                                   cccACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATGTCCGGTGCATCAGCCAGCGCacac                                                                  >  1:2073033/1‑150 (MQ=255)
                                                                                    ccACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACAcc                                                                 >  1:2709645/1‑150 (MQ=255)
                                                                                      aCCTGCGGTAACGGGAACTGTGCTGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAg                                                               <  1:1834046/150‑1 (MQ=255)
                                                                                      aCCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAg                                                               >  2:4281426/1‑150 (MQ=255)
                                                                                       ccTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCTACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGa                                                              <  1:1497390/150‑1 (MQ=255)
                                                                                        cTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCTACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGAc                                                             <  2:1986289/150‑1 (MQ=255)
                                                                                          gCGGTAACGGGAACTGTGCGGTTTCTCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCa                                                           <  1:2533140/150‑1 (MQ=255)
                                                                                             gTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATaa                                                        >  2:1405965/1‑150 (MQ=255)
                                                                                             gTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATaa                                                        >  2:3599721/1‑150 (MQ=255)
                                                                                               aaCGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAAc                                                      <  1:1915201/150‑1 (MQ=255)
                                                                                                 cGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAg                                                    >  1:3798951/1‑150 (MQ=255)
                                                                                                   ggAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGt                                                  >  1:840539/1‑150 (MQ=255)
                                                                                                   ggAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGt                                                  <  1:3649248/150‑1 (MQ=255)
                                                                                                   ggAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGt                                                  >  1:3017809/1‑150 (MQ=255)
                                                                                                    gAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCTACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTc                                                 >  1:2170975/1‑150 (MQ=255)
                                                                                                    gAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTc                                                 <  1:763965/150‑1 (MQ=255)
                                                                                                     aaCTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCg                                                >  1:2329792/1‑150 (MQ=255)
                                                                                                      aCTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGa                                               >  2:620984/1‑150 (MQ=255)
                                                                                                      aCTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGa                                               >  2:1397105/1‑150 (MQ=255)
                                                                                                       cTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCTACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAg                                              <  1:3078737/150‑1 (MQ=255)
                                                                                                       cTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCTACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAg                                              <  1:1529408/150‑1 (MQ=255)
                                                                                                        tgtgCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAgg                                             <  1:2134802/150‑1 (MQ=255)
                                                                                                             ggTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGAATCg                                        >  1:1624850/1‑150 (MQ=255)
                                                                                                               tttCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGAATCGCa                                      >  1:288838/1‑150 (MQ=255)
                                                                                                                  cGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGAATCGCACct                                   <  2:4023187/150‑1 (MQ=255)
                                                                                                                     aGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGAATCGCACctgct                                <  1:2814822/150‑1 (MQ=255)
                                                                                                                     aGCCGATACCCAGCGTCAGACCCGGTGTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGAATCGCACctgct                                <  1:4044808/150‑1 (MQ=255)
                                                                                                                       ccGATACCCAGCGTCAGACCCGGTTTGCTACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGAATCGCACCTGCTGt                              >  2:2719076/1‑150 (MQ=255)
                                                                                                                       ccGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGACTCGCACCTGCTGt                              >  2:1488630/1‑150 (MQ=255)
                                                                                                                        cGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGAATCGCACCTGCTGtt                             <  2:320021/150‑1 (MQ=255)
                                                                                                                        cGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGAATCGCACCTGCTGtt                             >  2:1821347/1‑150 (MQ=255)
                                                                                                                        cGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGAATCGCACCTGCTGtt                             <  2:2494856/150‑1 (MQ=255)
                                                                                                                          aTACCCAGCGTCAGACCCGGTTTGCTACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGAATCGCACCTGCTGTTCa                           <  1:3439756/150‑1 (MQ=255)
                                                                                                                                 gCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGAATCGCACCTGCTGTTCATCTacac                    <  2:3300055/150‑1 (MQ=255)
                                                                                                                                  cGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGAATCGCACCTGCTGTTCATCTACAcc                   <  1:308582/150‑1 (MQ=255)
                                                                                                                                  cGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGAATCGCACCTGCTGTTCATCTACAcc                   <  2:3190680/150‑1 (MQ=255)
                                                                                                                                  cGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGAATCGCACCTGCTGTTCATCTACAcc                   <  1:922174/150‑1 (MQ=255)
                                                                                                                                     cAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGAATCGCACCTGCTGTTCATCTACACCAAt                >  1:3932722/1‑150 (MQ=255)
                                                                                                                                      aGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGAATCGCACCTGCTGTTCATCTACACCAATc               >  2:555805/1‑150 (MQ=255)
                                                                                                                                         cccGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGAATCGCACCTGCTGTTCATCTACACCAATCGGc            >  1:897533/1‑150 (MQ=255)
                                                                                                                                            ggTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGAATCGCACCTGCTGTTCATCTACACCAATCGGCgag         >  1:3733146/1‑150 (MQ=255)
                                                                                                                                              tttGCTACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGAATCGCACCTGCTGTTCATCTACACCAATCGGCGAGaa       <  2:3907196/150‑1 (MQ=255)
                                                                                                                                                 gCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGAATCGCACCTGCTGTTCATCTACACCAATCGGCGAGAACgg    <  1:2290041/150‑1 (MQ=255)
                                                                                                                                                   gACCTTCGTGGATCACCCGGTTCCAGTTCACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCCGCCAGCGCACACCAGACCATAACCCGGTCGAGGAATCGCACCTGCTGTGCATCGACACCAATCGGCGAGAACGGGt  >  2:238711/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CACACACTTTCTGATACGCTTCGCCGCCGTTAATACTATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGAATCGCACCTGCTGTTCATCTACACCAATCGGCGAGAACGGGT  >  CP009273/2808499‑2808795

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: