Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,808,987 G→T 52.7% 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,9870GT52.7% 48.1 / inf 216S271* (TCG→TAG) gshAglutamate‑cysteine ligase
Reads supporting (aligned to +/- strand):  ref base G (50/52);  new base T (51/63);  total (101/115)
Fisher's exact test for biased strand distribution p-value = 5.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.91e-01

Only 100 of 217 total aligned reads displayed.

ACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCGCATACGGCAGGTAATACATACCGCACTCG  >  CP009273/2808841‑2809136
                                                                                                                                                  |                                                                                                                                                     
  aGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGc                                                                                                                                                  >  2:1355821/1‑150 (MQ=255)
      cATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAAACTTCGCGTACTCTTCCGATGTCGttt                                                                                                                                              >  1:2134997/1‑150 (MQ=255)
       aTCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGtttt                                                                                                                                             >  1:2116684/1‑150 (MQ=255)
          aGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGAt                                                                                                                                          >  1:801952/1‑150 (MQ=255)
           gAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGAtt                                                                                                                                         <  1:159231/150‑1 (MQ=255)
           gAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGAtt                                                                                                                                         >  1:1940559/1‑150 (MQ=255)
            aaGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTg                                                                                                                                        >  2:2849625/1‑150 (MQ=255)
            aaGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTg                                                                                                                                        <  2:1223793/150‑1 (MQ=255)
              ggCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGcc                                                                                                                                      >  2:3147564/1‑150 (MQ=255)
              ggCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGcc                                                                                                                                      <  2:257014/150‑1 (MQ=255)
              ggCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGcc                                                                                                                                      <  2:256552/150‑1 (MQ=255)
              ggCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGcc                                                                                                                                      >  2:2543987/1‑150 (MQ=255)
                 gACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGt                                                                                                                                   <  1:809607/150‑1 (MQ=255)
                  aCTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGtt                                                                                                                                  >  1:3645579/1‑150 (MQ=255)
                  aCTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGtt                                                                                                                                  <  2:2209588/150‑1 (MQ=255)
                   cTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGttt                                                                                                                                 <  2:2173613/150‑1 (MQ=255)
                   cTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGttt                                                                                                                                 >  1:2862233/1‑150 (MQ=255)
                     cgccgcTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTaa                                                                                                                               <  2:2814679/150‑1 (MQ=255)
                        cgcTGCGGGTAACGCGTTTTGGTCGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACTCGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCACTCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGc                                                                                                                            <  1:2685545/150‑1 (MQ=255)
                        cgcTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGc                                                                                                                            >  1:1003507/1‑150 (MQ=255)
                             cGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCt                                                                                                                       <  2:1461317/150‑1 (MQ=255)
                              gtgTAACGCGTTTTGGACGAATCGGCGCGTAAAGTTAGTTTTCAATCTGAAACACGTTGCTGTTGATTTGCAGACTCTTACCGTATTTATAAATACCAATCTTCGCGTAATCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTa                                                                                                                      <  2:196230/148‑1 (MQ=255)
                              gggTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTa                                                                                                                      <  1:3291795/150‑1 (MQ=255)
                              gggTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTa                                                                                                                      <  1:2300143/150‑1 (MQ=255)
                               ggTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTAc                                                                                                                     >  2:247539/1‑150 (MQ=255)
                                gTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACg                                                                                                                    <  1:2212050/150‑1 (MQ=255)
                                gTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACg                                                                                                                    <  2:3866461/150‑1 (MQ=255)
                                      cgTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCg                                                                                                              <  1:2806531/150‑1 (MQ=255)
                                        ttttGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTa                                                                                                            <  2:3707651/150‑1 (MQ=255)
                                          ttGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTaaa                                                                                                          <  1:3751404/150‑1 (MQ=255)
                                            ggACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGa                                                                                                        >  2:288532/1‑150 (MQ=255)
                                            ggACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGa                                                                                                        >  2:46356/1‑150 (MQ=255)
                                            ggACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGa                                                                                                        >  2:2718372/1‑150 (MQ=255)
                                            ggACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGa                                                                                                        >  2:685100/1‑150 (MQ=255)
                                              aCGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATc                                                                                                      <  2:3104682/150‑1 (MQ=255)
                                              aCGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATc                                                                                                      <  1:236281/150‑1 (MQ=255)
                                              aCGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATc                                                                                                      <  1:3819044/150‑1 (MQ=255)
                                                    cGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGaa                                                                                                >  1:2839227/1‑150 (MQ=255)
                                                         cgTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTaa                                                                                           <  2:1603885/150‑1 (MQ=255)
                                                          gTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAAt                                                                                          <  2:1812774/150‑1 (MQ=255)
                                                           tACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATa                                                                                         >  1:139915/1‑150 (MQ=255)
                                                             cAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATAcc                                                                                       <  1:1733729/150‑1 (MQ=255)
                                                              aGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCa                                                                                      >  1:2094915/1‑150 (MQ=255)
                                                                ttCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAg                                                                                    <  2:943899/150‑1 (MQ=255)
                                                                  cGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGAt                                                                                  >  2:668365/1‑150 (MQ=255)
                                                                   gTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGAtt                                                                                 >  2:1728320/1‑150 (MQ=255)
                                                                   gTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTGTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGAtt                                                                                 >  1:2324438/1‑150 (MQ=255)
                                                                    ttttCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTg                                                                                <  1:725177/150‑1 (MQ=255)
                                                                    ttttCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTg                                                                                <  2:1754774/150‑1 (MQ=255)
                                                                      ttCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGTAGGTAATACCAAGATTGCt                                                                              >  2:2800680/1‑150 (MQ=255)
                                                                      ttCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCt                                                                              >  2:2388665/1‑150 (MQ=255)
                                                                      ttCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCt                                                                              <  1:3961320/150‑1 (MQ=255)
                                                                       tCAATCTGCAACACGTTGCGGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTGTTGATTGCCTGTTTAATGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGAATGAtt                                                                             >  2:1271374/1‑150 (MQ=255)
                                                                         aaTCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTg                                                                           >  2:3728289/1‑150 (MQ=255)
                                                                           tCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCg                                                                         <  1:591561/150‑1 (MQ=255)
                                                                            cTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGa                                                                        >  2:3988892/1‑150 (MQ=255)
                                                                               cAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGAttt                                                                     <  1:898486/150‑1 (MQ=255)
                                                                               cAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGAttt                                                                     <  2:1635064/150‑1 (MQ=255)
                                                                                aaCACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTa                                                                    >  1:1363723/1‑150 (MQ=255)
                                                                                aaCACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTa                                                                    <  1:2876425/150‑1 (MQ=255)
                                                                                aaCACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTa                                                                    <  1:1797915/150‑1 (MQ=255)
                                                                                  cacGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTAtt                                                                  <  1:3187694/150‑1 (MQ=255)
                                                                                   acGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTg                                                                 <  2:3841175/150‑1 (MQ=255)
                                                                                       tGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGtat                                                             <  2:3496912/150‑1 (MQ=255)
                                                                                        gCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGtata                                                            <  1:2911084/150‑1 (MQ=255)
                                                                                        gCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGtata                                                            >  1:4127075/1‑150 (MQ=255)
                                                                                          tGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGc                                                          <  1:516187/150‑1 (MQ=255)
                                                                                          tGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGc                                                          >  2:3266938/1‑150 (MQ=255)
                                                                                           gTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCACTCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGcc                                                         <  1:2191788/150‑1 (MQ=255)
                                                                                           gTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGcc                                                         <  1:3373916/150‑1 (MQ=255)
                                                                                           gTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGcc                                                         <  1:2118309/150‑1 (MQ=255)
                                                                                           gTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGcc                                                         >  2:3606661/1‑150 (MQ=255)
                                                                                            ttGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCg                                                        >  2:728051/1‑150 (MQ=255)
                                                                                                  tGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCg                                                  >  2:2943891/1‑150 (MQ=255)
                                                                                                    cAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCt                                                >  2:2306435/1‑150 (MQ=255)
                                                                                                    cAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCt                                                <  1:1847904/150‑1 (MQ=255)
                                                                                                       ccTCTTACCGTCTTTCTCAAGACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCACGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCaa                                             <  2:1371193/150‑1 (MQ=255)
                                                                                                             aCCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAag                                       <  1:2166627/150‑1 (MQ=255)
                                                                                                             aCCGTCTTTCTCAATACCAATCTTCGCGTAAACTTCATATGGCGTTTTGATTGCATGTTTAAGGACCGCTACGAACTCGTAAAGATAGTTGAAGGTAATAACAAGATTCCTTTGAGAATTATTGGTCTAGCCGAGATAGCTCAAACGAag                                       >  1:1107146/1‑150 (MQ=255)
                                                                                                                 tCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATATATTGGTATAGCCGAGATCGCTCAAACGAagagag                                   >  1:1465397/1‑150 (MQ=255)
                                                                                                                  cTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCACTCAAACGAAGAGAgg                                  <  2:3328295/150‑1 (MQ=255)
                                                                                                                     tctcAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCg                               <  2:234489/150‑1 (MQ=255)
                                                                                                                         aaTACCAATCTTCGCGTACTCTTCCGATGGCGTTGTGATTGCCTGTTTAATGCCCGCTCCGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCACACGAAGAGAGGTCGCATa                           >  2:721185/1‑150 (MQ=255)
                                                                                                                           tACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCGCATACg                         >  1:1268385/1‑150 (MQ=255)
                                                                                                                           tACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCGCATACg                         >  1:999644/1‑150 (MQ=255)
                                                                                                                            aCCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCGCATACgg                        >  2:3913712/1‑150 (MQ=255)
                                                                                                                             ccAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCGCATACGGc                       >  1:2256103/1‑150 (MQ=255)
                                                                                                                             ccAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCGCATACGGc                       >  1:3340887/1‑150 (MQ=255)
                                                                                                                                aTCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCGCATACGGCAgg                    <  2:4125057/150‑1 (MQ=255)
                                                                                                                                 tCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCGCATACGGCAGGt                   <  2:553401/150‑1 (MQ=255)
                                                                                                                                  cTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCGCATACGGCAGGTa                  >  1:670686/1‑150 (MQ=255)
                                                                                                                                  cTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCGCATACGGCAGGTa                  <  2:1103584/150‑1 (MQ=255)
                                                                                                                                   ttCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCGCATACGGCAGGTaa                 >  2:385660/1‑150 (MQ=255)
                                                                                                                                     cgcgTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCGCATACGGCAGGTAata               <  2:1204629/150‑1 (MQ=255)
                                                                                                                                     cgcgTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCGCATACGGCAGGTAata               <  2:3960936/150‑1 (MQ=255)
                                                                                                                                     cgcgTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCGCATACGGCAGGTAata               >  2:2426012/1‑150 (MQ=255)
                                                                                                                                      gcgTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCGCATACGGCAGGTAatac              <  2:2749406/150‑1 (MQ=255)
                                                                                                                                       cgTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCGCATACGGCAGGTAataca             >  2:2510944/1‑150 (MQ=255)
                                                                                                                                        gTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCGCATACGGCAGGTAatacat            >  1:1680472/1‑150 (MQ=255)
                                                                                                                                             ctTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCGCATACGGCAGGTAATACATACCGc       >  1:2148720/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
ACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCGCATACGGCAGGTAATACATACCGCACTCG  >  CP009273/2808841‑2809136

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: