Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,808,987 G→T 6.0% S271* (TCG→TAG)  gshA ← glutamate‑cysteine ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,808,9870GT 6.0% 1260.5 / 29.2 327S271* (TCG→TAG) gshAglutamate‑cysteine ligase
Reads supporting (aligned to +/- strand):  ref base G (127/179);  new base T (5/16);  total (132/195)
Fisher's exact test for biased strand distribution p-value = 1.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 332 total aligned reads displayed.

TAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCGCATACGGCAGGTAATACATACCGCACTC  >  CP009273/2808839‑2809135
                                                                                                                                                    |                                                                                                                                                    
  aCAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGAt                                                                                                                                                    <  1:3906066/149‑1 (MQ=255)
  aCAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGAt                                                                                                                                                    <  1:696481/149‑1 (MQ=255)
   cAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTCCCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATg                                                                                                                                                   <  2:955873/149‑1 (MQ=255)
      cgcATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCg                                                                                                                                                <  1:3445729/149‑1 (MQ=255)
       gcATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGt                                                                                                                                               <  2:1414701/149‑1 (MQ=255)
        cATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGtt                                                                                                                                              >  1:1540544/1‑149 (MQ=255)
         aTCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGttt                                                                                                                                             <  1:4479819/149‑1 (MQ=255)
          tCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGtttt                                                                                                                                            <  2:583740/149‑1 (MQ=255)
               aGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTg                                                                                                                                       >  1:579202/1‑149 (MQ=255)
               aGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTg                                                                                                                                       >  2:2406815/1‑149 (MQ=255)
                ggCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGc                                                                                                                                      <  1:4257442/149‑1 (MQ=255)
                   gACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTg                                                                                                                                   <  2:149518/149‑1 (MQ=255)
                   gACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTg                                                                                                                                   <  1:4004560/149‑1 (MQ=255)
                    aCTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGt                                                                                                                                  >  1:1268264/1‑149 (MQ=255)
                      tCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCATGttt                                                                                                                                >  2:3699427/1‑149 (MQ=255)
                          cgcTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAgg                                                                                                                            >  2:294322/1‑149 (MQ=255)
                          cgcTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAgg                                                                                                                            <  1:4663905/149‑1 (MQ=255)
                             tGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGccc                                                                                                                         <  1:1919658/149‑1 (MQ=255)
                             tGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGccc                                                                                                                         <  1:3057159/149‑1 (MQ=255)
                                gggTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCt                                                                                                                      >  2:2735297/1‑149 (MQ=255)
                                      cgcgTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTAc                                                                                                                <  2:2688368/149‑1 (MQ=255)
                                       gcgTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACt                                                                                                               >  2:36078/1‑149 (MQ=255)
                                         gTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCg                                                                                                             >  1:1004005/1‑149 (MQ=255)
                                          ttttGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGt                                                                                                            <  1:2628074/149‑1 (MQ=255)
                                           tttGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTa                                                                                                           <  1:3247203/149‑1 (MQ=255)
                                            ttGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTaa                                                                                                          <  1:2831534/149‑1 (MQ=255)
                                             tGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTaaa                                                                                                         <  1:94822/149‑1 (MQ=255)
                                             tGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTaaa                                                                                                         <  1:4571656/149‑1 (MQ=255)
                                               gACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGa                                                                                                       >  1:527351/1‑149 (MQ=255)
                                                aCGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGAt                                                                                                      <  2:6185/149‑1 (MQ=255)
                                                aCGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGAt                                                                                                      <  2:6041/149‑1 (MQ=255)
                                                aCGAATCGGCGCGTACAGTTAGTGTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTCCCGTCTTTCTCAATCCCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGAt                                                                                                      <  2:3878039/149‑1 (MQ=255)
                                                   aaTCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGt                                                                                                   >  2:3968521/1‑149 (MQ=255)
                                                      cGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGa                                                                                                <  1:772260/149‑1 (MQ=255)
                                                       ggCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGaa                                                                                               <  1:432039/149‑1 (MQ=255)
                                                          gcgTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGt                                                                                            <  1:3577958/149‑1 (MQ=255)
                                                            gTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTaa                                                                                          <  2:3503274/149‑1 (MQ=255)
                                                            gTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTaa                                                                                          <  2:3789635/149‑1 (MQ=255)
                                                             tACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAAt                                                                                         <  2:112735/149‑1 (MQ=255)
                                                             tACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAAt                                                                                         <  2:1802760/149‑1 (MQ=255)
                                                              aCAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATa                                                                                        <  1:1464337/149‑1 (MQ=255)
                                                              aCAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATa                                                                                        >  2:2158619/1‑149 (MQ=255)
                                                               catttagTTTTCACTCTTCAAAACGTTTCTGTTGATTAGCATCCTCTTACCGTATTTCTCAATAACAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCTTAAAGATCGTTGAAGGTAATAc                                                                                       <  2:4701979/143‑1 (MQ=255)
                                                                agttagttTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATAcc                                                                                      <  1:2142039/149‑1 (MQ=255)
                                                                   tagTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAg                                                                                   <  1:3269829/147‑1 (MQ=255)
                                                                    agTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGa                                                                                  >  2:1741107/2‑149 (MQ=255)
                                                                     gTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGAt                                                                                 <  2:2058749/149‑1 (MQ=255)
                                                                      ttttCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGAtt                                                                                <  1:2290354/149‑1 (MQ=255)
                                                                            aTCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTg                                                                          <  2:815155/149‑1 (MQ=255)
                                                                                gCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGAt                                                                      <  2:4406707/149‑1 (MQ=255)
                                                                                gCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATAACAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTAGTAAATATCGTTGAAGGTAATACCAAGAATGAttttagat                                                                      >  1:2057549/1‑144 (MQ=255)
                                                                                 cAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATAACAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGAtt                                                                     <  1:3196344/149‑1 (MQ=255)
                                                                                   acacGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTa                                                                   >  2:3593682/1‑149 (MQ=255)
                                                                                   acacGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTa                                                                   >  2:1053063/1‑149 (MQ=255)
                                                                                    cacGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTAt                                                                  <  2:1339470/149‑1 (MQ=255)
                                                                                    cacGTTGCTGTTGATTTGCAGCCTCTTACCGTCTGTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTAt                                                                  <  2:3420423/149‑1 (MQ=255)
                                                                                     acGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTAtt                                                                 <  2:2564900/149‑1 (MQ=255)
                                                                                       gTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTgg                                                               >  1:3944715/1‑149 (MQ=255)
                                                                                         tGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGta                                                             <  1:553461/149‑1 (MQ=255)
                                                                                          gCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGtat                                                            >  2:3893302/1‑149 (MQ=255)
                                                                                            tGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAg                                                          <  2:43350/149‑1 (MQ=255)
                                                                                             gTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGc                                                         <  2:3766371/149‑1 (MQ=255)
                                                                                                gATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCga                                                      <  1:1814495/149‑1 (MQ=255)
                                                                                                 aTTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCgag                                                     <  1:822745/149‑1 (MQ=255)
                                                                                                  tttGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCgaga                                                    <  1:4177955/149‑1 (MQ=255)
                                                                                                    tGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATc                                                  >  2:1987366/1‑149 (MQ=255)
                                                                                                       aGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCt                                               <  2:3506500/149‑1 (MQ=255)
                                                                                                       aGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCt                                               >  2:1944048/1‑149 (MQ=255)
                                                                                                       aGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCt                                               >  1:2228563/1‑149 (MQ=255)
                                                                                                          ctctTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCaa                                            >  2:4143576/1‑149 (MQ=255)
                                                                                                            ctgaCCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAAGACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAAc                                          <  1:4651229/146‑1 (MQ=255)
                                                                                                            ctTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAAc                                          <  1:2252801/149‑1 (MQ=255)
                                                                                                            ctTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAAc                                          <  1:2922101/149‑1 (MQ=255)
                                                                                                               aCCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGTGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGaa                                       <  1:3821061/149‑1 (MQ=255)
                                                                                                               aCCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGaa                                       >  1:2350873/1‑149 (MQ=255)
                                                                                                                 cGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAaga                                     >  2:3298165/1‑149 (MQ=255)
                                                                                                                  gTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAagag                                    <  1:3662315/149‑1 (MQ=255)
                                                                                                                   tCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAagaga                                   >  2:4739043/1‑149 (MQ=255)
                                                                                                                   tCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAagaga                                   >  2:3002561/1‑149 (MQ=255)
                                                                                                                   tCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAagaga                                   <  2:1697151/149‑1 (MQ=255)
                                                                                                                      ttCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGt                                <  2:766179/149‑1 (MQ=255)
                                                                                                                       tctcAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTc                               >  2:3919536/1‑149 (MQ=255)
                                                                                                                       tctcAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTc                               >  2:4294005/1‑149 (MQ=255)
                                                                                                                           aaTACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCGCAt                           >  2:2768128/1‑149 (MQ=255)
                                                                                                                            aTACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCGCATa                          >  2:3197083/1‑149 (MQ=255)
                                                                                                                             tACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCGCATAc                         >  1:4051082/1‑149 (MQ=255)
                                                                                                                              aCCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCGCATACg                        <  2:3330165/149‑1 (MQ=255)
                                                                                                                                cAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCGCATACGGc                      <  2:4501274/149‑1 (MQ=255)
                                                                                                                                 aaTCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCGCATACTGCa                     >  1:2496166/1‑149 (MQ=255)
                                                                                                                                   tCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCGCATACGGCAgg                   >  2:2248/1‑149 (MQ=255)
                                                                                                                                     ttCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCGCATACGGCAGGTa                 <  2:1177729/149‑1 (MQ=255)
                                                                                                                                     ttCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCGCATACGGCAGGTa                 >  2:4650085/1‑149 (MQ=255)
                                                                                                                                        gcgTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCGCATACGGCAGGTAata              <  1:1185169/149‑1 (MQ=255)
                                                                                                                                          gTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCGCATACGGCAGGTAataca            <  1:3405019/149‑1 (MQ=255)
                                                                                                                                          gTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCGCATACGGCAGGTAataca            <  1:1714372/149‑1 (MQ=255)
                                                                                                                                          gTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCGCATACGGCAGGTAataca            >  1:1076838/1‑149 (MQ=255)
                                                                                                                                           tACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCGCATACGGCAGGTAatacat           >  2:1951294/1‑149 (MQ=255)
                                                                                                                                             ctctTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCGCATACGGCAGGTAatacatac         <  1:1053063/149‑1 (MQ=255)
                                                                                                                                                ttCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCGCATACGGCAGGTAATACATACCGc      <  1:2369269/149‑1 (MQ=255)
                                                                                                                                                    gATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCGCATACGGCAGGTAATACATACCGCACTc  <  2:2364053/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCGCATACGGCAGGTAATACATACCGCACTC  >  CP009273/2808839‑2809135

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: