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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273134,0280GT88.4% 344.4 / 44.8 147A153A (GCG→GCTcueOmulticopper oxidase (laccase)
Reads supporting (aligned to +/- strand):  ref base G (5/12);  new base T (54/76);  total (59/88)
Fisher's exact test for biased strand distribution p-value = 4.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.95e-01

Only 100 of 148 total aligned reads displayed.

GCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCGATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATAAGAAATTTAGCGCCGACGGGCAGATTGATTATCAACTGGAT  >  CP009273/133881‑134172
                                                                                                                                                   |                                                                                                                                                
gCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCg                                                                                                                                                  <  2:664805/148‑1 (MQ=255)
  tGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTAt                                                                                                                                                <  2:2928706/148‑1 (MQ=255)
   ggAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATg                                                                                                                                               >  2:339248/1‑148 (MQ=255)
     aaGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATggg                                                                                                                                             >  2:487435/1‑148 (MQ=255)
      aGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTAGTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATgggg                                                                                                                                            <  1:1362134/148‑1 (MQ=255)
       gTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCGATGGggc                                                                                                                                           <  1:514474/148‑1 (MQ=255)
        tACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGggct                                                                                                                                          <  1:2180232/148‑1 (MQ=255)
        tACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCGATGGggct                                                                                                                                          <  2:955220/148‑1 (MQ=255)
          ccGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGggctgg                                                                                                                                        <  1:540972/148‑1 (MQ=255)
                  aGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTggctgggctg                                                                                                                                <  2:620332/148‑1 (MQ=255)
                   gTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGgctgggctgg                                                                                                                               <  2:2528823/148‑1 (MQ=255)
                      gACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCtggtgg                                                                                                                            <  1:1239551/148‑1 (MQ=255)
                      gACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCtggtgg                                                                                                                            <  1:731881/148‑1 (MQ=255)
                      gACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCtggtgg                                                                                                                            <  1:731778/148‑1 (MQ=255)
                          gcggcCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGAt                                                                                                                        <  2:3292646/148‑1 (MQ=255)
                          gcggcCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGAt                                                                                                                        >  2:3002637/1‑148 (MQ=255)
                          gcggcCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGAt                                                                                                                        >  2:2407349/1‑148 (MQ=255)
                             gcCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGa                                                                                                                     <  1:960285/148‑1 (MQ=255)
                             gcCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGa                                                                                                                     >  2:2700275/1‑148 (MQ=255)
                               ccGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAg                                                                                                                   >  2:1878745/1‑148 (MQ=255)
                               ccGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAg                                                                                                                   <  2:2750938/148‑1 (MQ=255)
                                cGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATagaaga                                                                                                                  <  1:1183965/148‑1 (MQ=255)
                                 gCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGAt                                                                                                                 <  2:1964164/148‑1 (MQ=255)
                                  cAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATg                                                                                                                <  1:2661863/148‑1 (MQ=255)
                                  cAGGGAATTAGTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATg                                                                                                                <  2:3282695/148‑1 (MQ=255)
                                  cAGGGAATTAGTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCGATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATg                                                                                                                <  1:2336510/148‑1 (MQ=255)
                                   aGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGa                                                                                                               <  1:1500707/148‑1 (MQ=255)
                                    gggAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTCTGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGAc                                                                                                              >  2:2303396/1‑148 (MQ=255)
                                       aaTTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACgag                                                                                                           <  2:3168130/148‑1 (MQ=255)
                                        attattCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACgaga                                                                                                          <  2:3159846/148‑1 (MQ=255)
                                         ttattCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGAt                                                                                                         <  1:1011272/148‑1 (MQ=255)
                                          tattCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATc                                                                                                        <  2:683551/148‑1 (MQ=255)
                                          tattCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATc                                                                                                        <  2:683620/148‑1 (MQ=255)
                                            ttCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCt                                                                                                      >  2:2664532/1‑148 (MQ=255)
                                            ttCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCt                                                                                                      >  2:2664938/1‑148 (MQ=255)
                                                gCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGaaa                                                                                                  <  1:1503736/148‑1 (MQ=255)
                                                 ccAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAAt                                                                                                 <  2:2973740/148‑1 (MQ=255)
                                                 ccAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAAt                                                                                                 <  2:2973619/148‑1 (MQ=255)
                                                  cAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAAtt                                                                                                <  1:1068696/148‑1 (MQ=255)
                                                   aGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTa                                                                                               <  1:2783248/148‑1 (MQ=255)
                                                   aGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTa                                                                                               <  2:2101298/148‑1 (MQ=255)
                                                   aGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCGGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTa                                                                                               >  2:3261738/1‑148 (MQ=255)
                                                    ggTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTaa                                                                                              <  2:595826/148‑1 (MQ=255)
                                                        gCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAAtgct                                                                                          >  1:1626545/1‑148 (MQ=255)
                                                         cAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAAtgctg                                                                                         <  1:1918056/148‑1 (MQ=255)
                                                         cAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAAtgctg                                                                                         <  2:1457075/148‑1 (MQ=255)
                                                            gcgcTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCNGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCa                                                                                      >  1:1677219/1‑148 (MQ=255)
                                                               cTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCaaaa                                                                                   >  1:1870518/1‑148 (MQ=255)
                                                               cTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCaaaa                                                                                   <  1:3170650/148‑1 (MQ=255)
                                                               cTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGGGGTGATTGAAGATGTCGAGATCCTGACAGTCAtgctgcgccag                                                                                   >  1:1870786/1‑143 (MQ=255)
                                                               cTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCGATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCaaaa                                                                                   <  1:1125617/148‑1 (MQ=255)
                                                                tCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAAc                                                                                  <  1:3298752/148‑1 (MQ=255)
                                                                tCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAAc                                                                                  <  1:3298748/148‑1 (MQ=255)
                                                                    tGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTg                                                                              <  2:3246999/148‑1 (MQ=255)
                                                                     gACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCGATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTgg                                                                             >  1:1447700/1‑148 (MQ=255)
                                                                       cgttgaacgttgaTCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTgggg                                                                           >  2:22375/1‑148 (MQ=255)
                                                                          tgaacgttgaTCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTAt                                                                        <  1:962173/148‑1 (MQ=255)
                                                                            aacgttgaTCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCg                                                                      >  1:2150019/1‑148 (MQ=255)
                                                                            aacgttgaTCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCGATGGGGCTGGCTGGGCGGGTGGGGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTCTCg                                                                      >  2:2491874/1‑148 (MQ=255)
                                                                                ttgaTCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCgatga                                                                  <  2:1156855/148‑1 (MQ=255)
                                                                                  gaTCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGatgatg                                                                <  2:2032843/148‑1 (MQ=255)
                                                                                   aTCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGt                                                               >  1:917933/1‑148 (MQ=255)
                                                                                    tCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGtt                                                              <  1:2384923/148‑1 (MQ=255)
                                                                                     cAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTc                                                             >  2:111450/1‑148 (MQ=255)
                                                                                       aCCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCg                                                           >  1:1209242/1‑148 (MQ=255)
                                                                                        ccTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCgg                                                          <  1:3002536/148‑1 (MQ=255)
                                                                                        ccTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCgg                                                          <  1:2407349/148‑1 (MQ=255)
                                                                                        ccTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCgg                                                          <  1:887765/148‑1 (MQ=255)
                                                                                           gCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGa                                                       <  1:2148852/148‑1 (MQ=255)
                                                                                              gCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCGATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCg                                                    <  2:2391805/148‑1 (MQ=255)
                                                                                               cTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGt                                                   <  2:2917498/148‑1 (MQ=255)
                                                                                                  ccTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCa                                                <  1:2201076/148‑1 (MQ=255)
                                                                                                    tgctgGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAgg                                              <  2:1995951/148‑1 (MQ=255)
                                                                                                       tgGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATa                                           <  2:2261473/148‑1 (MQ=255)
                                                                                                         gTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATAAg                                         >  1:2639933/1‑148 (MQ=255)
                                                                                                            ccATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATAAGaaa                                      <  1:2031377/148‑1 (MQ=255)
                                                                                                            ccATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATAAGaaa                                      <  1:2031376/148‑1 (MQ=255)
                                                                                                             cATCCGCATCAGCACGGCAAACCCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATAAGAAAt                                     <  2:945519/148‑1 (MQ=255)
                                                                                                             cATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATAAGAAAt                                     <  1:1878821/148‑1 (MQ=255)
                                                                                                               tCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCGGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATAAGAAAttt                                   >  2:2930147/1‑148 (MQ=255)
                                                                                                               tCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGATGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATAAGAAAttt                                   >  1:3171973/1‑148 (MQ=255)
                                                                                                                  gCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATAAGAAATTTAgc                                >  1:1537369/1‑148 (MQ=255)
                                                                                                                  gCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATAAGAAATTTAgc                                >  1:1537213/1‑148 (MQ=255)
                                                                                                                   cATCAGCACGGCAAAACCGGGCGACAGGTGGCGATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATAAGAAATTTAgcg                               <  1:2901746/148‑1 (MQ=255)
                                                                                                                     tCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATAAGAAATTTATCGcc                             >  2:1922071/1‑148 (MQ=255)
                                                                                                                     tCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATAAGAAATTTAGCGcc                             >  2:127055/1‑148 (MQ=255)
                                                                                                                     tCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATAAGAAATTTAGCGcc                             >  2:1922051/1‑148 (MQ=255)
                                                                                                                           cGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATAAGAAATTTAGCGCCGACggg                       <  2:3020167/148‑1 (MQ=255)
                                                                                                                           cGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATAAGAAATTTAGCGCCGACggg                       >  2:2228856/1‑148 (MQ=255)
                                                                                                                                aaaCCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATAAGAAATTTAGCGCCGACGGGCAgat                  >  2:2535993/1‑148 (MQ=255)
                                                                                                                                aaaCCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATAAGAAATTTAGCGCCGACGGGCAgat                  >  2:1106460/1‑148 (MQ=255)
                                                                                                                                aaaCCGGGCGACAGGTGGCGATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATAAGAAATTTAGCGCCGACGGGCAgat                  <  2:1207323/148‑1 (MQ=255)
                                                                                                                                 aaCCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATAAGAAATTTAGCGCCGACGGGCAgatt                 >  2:140406/1‑148 (MQ=255)
                                                                                                                                     gggCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATAAGAAATTTAGCGCCGACGGGCAgattgatt             <  2:861144/148‑1 (MQ=255)
                                                                                                                                         gACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATAAGAAATTTAGCGCCGACGGGCAGATTGATTATCa         <  2:3228925/148‑1 (MQ=255)
                                                                                                                                          aCAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATAAGAAATTTAGCGCCGACGGGCAGATTGATTATCaa        <  1:1047328/148‑1 (MQ=255)
                                                                                                                                             ggTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATAAGAAATTTAGCGCCGACGGGCAGATTGATTATCAACTg     >  2:1749186/1‑148 (MQ=255)
                                                                                                                                               tGGCGATGGGGCTTGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATAAGAAATTTAGCGCCGACGGGCAGATTGATTATCAACTGGa   <  2:249814/148‑1 (MQ=255)
                                                                                                                                               tGGCGATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATAAGAAATTTAGCGCCGACGGGCAGATTGATTATCAACTGGa   <  2:762086/148‑1 (MQ=255)
                                                                                                                                                ggCTATGGGGCTGGGTGGGCTGGTTGTGATTGAAGATTACGAGATCCTGAAATTAATGCTGCCAAAACTGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATAAGAAATTTAGCGCCGACGGGCAGATTGATTATCAACTGGAt  <  1:2986448/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
GCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCGATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATAAGAAATTTAGCGCCGACGGGCAGATTGATTATCAACTGGAT  >  CP009273/133881‑134172

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: