Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 134,028 G→T 100% A153A (GCG→GCT cueO → multicopper oxidase (laccase)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273134,0280GT99.5% 648.7 / ‑4.1 207A153A (GCG→GCTcueOmulticopper oxidase (laccase)
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (99/107);  total (99/108)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 207 total aligned reads displayed.

TGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCGATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATAAGAAATTTAGCGCCGACGGGCAGATTGATTATCAACTGGA  >  CP009273/133883‑134171
                                                                                                                                                 |                                                                                                                                               
  gAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATgg                                                                                                                                             >  2:532863/1‑148 (MQ=255)
   aaGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATggg                                                                                                                                            >  1:2708972/1‑148 (MQ=255)
     gTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGggc                                                                                                                                          >  2:2658147/1‑148 (MQ=255)
     gTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGggc                                                                                                                                          >  1:2003212/1‑148 (MQ=255)
        ccGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGggctgg                                                                                                                                       <  1:1752740/148‑1 (MQ=255)
             tGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTggctggg                                                                                                                                  <  1:1229199/148‑1 (MQ=255)
             tGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTggctggg                                                                                                                                  <  1:784901/148‑1 (MQ=255)
             tGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTggctggg                                                                                                                                  <  1:781488/148‑1 (MQ=255)
             tGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTggctggg                                                                                                                                  <  1:3106710/148‑1 (MQ=255)
               aaGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTggctgggct                                                                                                                                <  2:963328/148‑1 (MQ=255)
               aaGTCGACGGCGGCCCGCAGGGAATTAGTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTggctgggct                                                                                                                                <  2:1860189/148‑1 (MQ=255)
                 gTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGgctgggctgg                                                                                                                              <  1:2730141/148‑1 (MQ=255)
                   cGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCGATGGGGCTGGCTGGGCtggtg                                                                                                                            <  1:2516616/148‑1 (MQ=255)
                     aCGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTggtggt                                                                                                                          >  2:2553417/1‑148 (MQ=255)
                     aCGGCGGCCCGCAGGGAATTATTCAGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTggtggt                                                                                                                          >  2:3268695/1‑148 (MQ=255)
                      cggcggCCCGCAGGGAATTAGTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGgtggtg                                                                                                                         <  2:3211455/148‑1 (MQ=255)
                        gcggcCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGAt                                                                                                                       <  2:3395876/148‑1 (MQ=255)
                        gcggcCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGAt                                                                                                                       >  2:175834/1‑148 (MQ=255)
                         cggcCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGAtt                                                                                                                      <  2:2683701/148‑1 (MQ=255)
                           gcCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGa                                                                                                                    <  2:404579/148‑1 (MQ=255)
                            cccGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGaa                                                                                                                   >  2:3464608/1‑148 (MQ=255)
                             ccGCAGGGAATGATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCGACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAg                                                                                                                  <  2:2509915/148‑1 (MQ=255)
                                 aGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGa                                                                                                              <  1:2809227/148‑1 (MQ=255)
                                 aGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGa                                                                                                              <  2:2957374/148‑1 (MQ=255)
                                  gggAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGAc                                                                                                             >  2:1999464/1‑148 (MQ=255)
                                  gggAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGAc                                                                                                             >  2:1405296/1‑148 (MQ=255)
                                  gggAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCAGCTGGTTCCATCCGCAACAGCACGGCAAAACCCGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGAc                                                                                                             >  1:1493050/1‑148 (MQ=255)
                                    gAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACga                                                                                                           >  2:1565034/1‑148 (MQ=255)
                                     aaTTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACgag                                                                                                          <  2:2201695/148‑1 (MQ=255)
                                     aaTGATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACgag                                                                                                          <  2:2009407/148‑1 (MQ=255)
                                      attattCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACgaga                                                                                                         >  1:2584987/1‑148 (MQ=255)
                                       ttattCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGAt                                                                                                        <  2:2833613/148‑1 (MQ=255)
                                        tattCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATc                                                                                                       <  1:3235337/148‑1 (MQ=255)
                                          ttCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCt                                                                                                     >  2:1521956/1‑148 (MQ=255)
                                             cGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGaa                                                                                                  <  1:109171/148‑1 (MQ=255)
                                                   gTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAAt                                                                                            >  1:2868334/1‑148 (MQ=255)
                                                   gTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAAt                                                                                            >  1:1645610/1‑148 (MQ=255)
                                                      gCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAAtgct                                                                                         >  1:593488/1‑148 (MQ=255)
                                                        aaGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAAtgctgc                                                                                       <  2:2259214/148‑1 (MQ=255)
                                                        aaGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAAtgctgc                                                                                       >  1:1543270/1‑148 (MQ=255)
                                                        aaGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAAtgatgc                                                                                       >  1:3284607/1‑148 (MQ=255)
                                                           cgcTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCaa                                                                                    >  1:3442406/1‑148 (MQ=255)
                                                           cgcTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCaa                                                                                    >  2:1046903/1‑148 (MQ=255)
                                                           cgcTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCaa                                                                                    >  1:377527/1‑148 (MQ=255)
                                                           cgcTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCGGGTGGTGATTGACGATGTCGAGATCCTGAAATTAATGCTGccca                                                                                    >  1:3348421/1‑146 (MQ=255)
                                                            gcTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCaaa                                                                                   <  1:2939899/148‑1 (MQ=255)
                                                              tCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAAc                                                                                 <  2:2729117/148‑1 (MQ=255)
                                                               cGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACa                                                                                <  1:1395839/148‑1 (MQ=255)
                                                                   gACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTgg                                                                            >  2:1228679/1‑148 (MQ=255)
                                                                    acgttgaacgttgaTCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTggg                                                                           <  2:2868334/148‑1 (MQ=255)
                                                                      gttgaacgttgaTCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGt                                                                         <  1:407914/148‑1 (MQ=255)
                                                                            cgttgaTCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCgat                                                                   >  2:2015206/1‑148 (MQ=255)
                                                                            cgttgaTCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCGGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCgat                                                                   >  2:2015369/1‑148 (MQ=255)
                                                                             gttgaTCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCgatg                                                                  <  1:3456700/148‑1 (MQ=255)
                                                                              ttgaTCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCgatga                                                                 >  2:1882819/1‑148 (MQ=255)
                                                                              ttgaTCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCgatga                                                                 >  2:1624827/1‑148 (MQ=255)
                                                                               tgaTCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCgatgat                                                                <  2:1970357/148‑1 (MQ=255)
                                                                               tgaTCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCgatgat                                                                <  2:2466207/148‑1 (MQ=255)
                                                                                gaTCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGatgatg                                                               <  2:56769/148‑1 (MQ=255)
                                                                                gaTCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGatgatg                                                               <  2:1834462/148‑1 (MQ=255)
                                                                                 aTCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGt                                                              >  1:469398/1‑148 (MQ=255)
                                                                                 aTCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGt                                                              <  2:3023091/148‑1 (MQ=255)
                                                                                  tCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGtt                                                             <  2:2233745/148‑1 (MQ=255)
                                                                                  tCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGtt                                                             <  1:321049/148‑1 (MQ=255)
                                                                                     aCCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCg                                                          >  1:1331343/1‑148 (MQ=255)
                                                                                      ccTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATTACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCgg                                                         >  1:865321/1‑148 (MQ=255)
                                                                                       cTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGt                                                        <  2:1861651/148‑1 (MQ=255)
                                                                                         gCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGa                                                      >  1:1676618/1‑148 (MQ=255)
                                                                                          ccGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGAt                                                     <  1:532863/148‑1 (MQ=255)
                                                                                              tACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGtt                                                 <  2:593488/148‑1 (MQ=255)
                                                                                              tACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGtt                                                 <  2:2708972/148‑1 (MQ=255)
                                                                                               aCCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTc                                                <  1:387571/148‑1 (MQ=255)
                                                                                                  tgctgGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAgg                                             <  1:2737125/148‑1 (MQ=255)
                                                                                                  tgctgGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAgg                                             <  1:191292/148‑1 (MQ=255)
                                                                                                  tgctgGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAgg                                             <  1:197385/148‑1 (MQ=255)
                                                                                                      ggTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGGGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATaa                                         >  2:803296/1‑148 (MQ=255)
                                                                                                       gTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATAAg                                        >  1:804519/1‑148 (MQ=255)
                                                                                                           cATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATAAGAAAt                                    <  2:377527/148‑1 (MQ=255)
                                                                                                              ccGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATAAGAAATTTa                                 <  1:3043780/148‑1 (MQ=255)
                                                                                                               cGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATAAGAAATTTAg                                <  1:2313295/148‑1 (MQ=255)
                                                                                                                gCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATAAGAAATTTAgc                               >  2:1641691/1‑148 (MQ=255)
                                                                                                                    cAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCGAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATAAGAAATTTAGCGCCg                           >  1:776805/1‑148 (MQ=255)
                                                                                                                    cAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATAAGAAATTTAGCGCCg                           >  1:2728141/1‑148 (MQ=255)
                                                                                                                      gCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATAAGAAATTTAGCGCCGAc                         >  1:2618295/1‑148 (MQ=255)
                                                                                                                       cACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATAAGAAATTTAGCGCCGACg                        <  2:2350051/148‑1 (MQ=255)
                                                                                                                            cAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATAAGAAATTTAGCGCCGACGGGCAg                   >  1:1859873/1‑148 (MQ=255)
                                                                                                                              aaaCCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATAAGAAATTTAGCGCCGACGGGCAgat                 >  2:2384493/1‑148 (MQ=255)
                                                                                                                               acccGGGCGACCGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATAAGAAATTTAGCGCCGACGGGCAgatt                <  1:131518/146‑1 (MQ=255)
                                                                                                                               aaCCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATAAGAAATTTAGCGCCGACGGGCAgatt                <  1:3090227/148‑1 (MQ=255)
                                                                                                                               aaCCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATAAGAAATTTAGCGCCGACGGGCAgatt                >  2:210166/1‑148 (MQ=255)
                                                                                                                                   gggCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATAAGAAATTTAGCGCCGACGGGCAgattgatt            <  1:2140728/148‑1 (MQ=255)
                                                                                                                                     gCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATAAGAAATTTAGCGCCGACGGGCAGATTGATTAt          <  2:2003767/148‑1 (MQ=255)
                                                                                                                                     gCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATAAGAAATTTAGCGCCGACGGGCAGATTGATTAt          <  1:1565034/148‑1 (MQ=255)
                                                                                                                                        aCAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATAAGAAATTTAGCGCCGACGGGCAGATTGATTATCaa       <  1:2553417/148‑1 (MQ=255)
                                                                                                                                        aCAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATAAGAAATTTAGCGCCGACGGGCAGATTGATTATCaa       >  1:961978/1‑148 (MQ=255)
                                                                                                                                         cAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATAAGAAATTTAGCGCCGACGGGCAGATTGATTATCAAc      >  2:73277/1‑148 (MQ=255)
                                                                                                                                         cAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATAAGAAATTTAGCGCCGACGGGCAGATTGATTATCAAc      >  2:1543077/1‑148 (MQ=255)
                                                                                                                                         cAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATAAGAAATTTAGCGCCGACGGGCAGATTGATTATCAAc      >  1:2901571/1‑148 (MQ=255)
                                                                                                                                           ggTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATTGATGATGTTCCGGTGATCGTTCAGGATAAGAAATTTAGCGCCGACGGGCAGATTGATTATCAACTg    >  1:561311/1‑148 (MQ=255)
                                                                                                                                             tGGCTATGGGGCTGGCTGGGCTGGTGGTGATTTAATATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCTTTCAGGATAAGAAATTTAGCGCCGACGGGCAGATTGATTATCAACTGGa  <  1:878789/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
TGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCGATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATAAGAAATTTAGCGCCGACGGGCAGATTGATTATCAACTGGA  >  CP009273/133883‑134171

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: