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,0280GT95.6% 660.3 / 20.4 228A153A (GCG→GCTcueOmulticopper oxidase (laccase)
Reads supporting (aligned to +/- strand):  ref base G (3/7);  new base T (96/122);  total (99/129)
Fisher's exact test for biased strand distribution p-value = 5.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.62e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 228 total aligned reads displayed.

AAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCGATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATAAGAAATTTAGCGCCGACGGGCAGATTGATTATCAACTGGA  >  CP009273/133886‑134171
                                                                                                                                              |                                                                                                                                               
aaGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATggg                                                                                                                                            >  1:149665/1‑148 (MQ=255)
 aGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATgggg                                                                                                                                           <  1:3529737/148‑1 (MQ=255)
 aGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATgggg                                                                                                                                           <  1:3454915/148‑1 (MQ=255)
   tACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGggct                                                                                                                                         >  1:551963/1‑148 (MQ=255)
   tACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGggct                                                                                                                                         <  2:3396888/148‑1 (MQ=255)
      cGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGggctggc                                                                                                                                      >  2:2871497/1‑148 (MQ=255)
           gAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTggctgggc                                                                                                                                 <  1:1582608/148‑1 (MQ=255)
            aaGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTggctgggct                                                                                                                                >  2:1120152/1‑148 (MQ=255)
             aGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTggctgggctg                                                                                                                               >  1:412416/1‑148 (MQ=255)
             aGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTggctgggctg                                                                                                                               <  1:2274101/148‑1 (MQ=255)
              gTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGgctgggctgg                                                                                                                              >  2:152544/1‑148 (MQ=255)
               tCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCtggt                                                                                                                             <  1:1140264/148‑1 (MQ=255)
               tCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCtggt                                                                                                                             >  2:2847969/1‑148 (MQ=255)
                  aCGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCGATGGGGCTGGCTGGGCTggtggt                                                                                                                          >  2:440112/1‑148 (MQ=255)
                   cggcggCCCGCAGGGAATTAGTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGgtggtg                                                                                                                         <  2:2435839/148‑1 (MQ=255)
                    ggcggcCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGa                                                                                                                        >  2:1493951/1‑148 (MQ=255)
                      cggcCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGAtt                                                                                                                      <  1:895376/148‑1 (MQ=255)
                        gcCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGa                                                                                                                    <  2:2013439/148‑1 (MQ=255)
                        gcCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGa                                                                                                                    >  1:2700851/1‑148 (MQ=255)
                          ccGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAg                                                                                                                  >  1:519381/1‑148 (MQ=255)
                           cGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGa                                                                                                                 <  1:2851044/148‑1 (MQ=255)
                             cAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATg                                                                                                               <  1:383776/148‑1 (MQ=255)
                             cAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATg                                                                                                               <  1:1671461/148‑1 (MQ=255)
                              aGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGa                                                                                                              >  1:1912301/1‑148 (MQ=255)
                                ggAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAACCCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACg                                                                                                            <  2:2082235/148‑1 (MQ=255)
                                ggAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACg                                                                                                            <  2:2581952/148‑1 (MQ=255)
                                ggAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACg                                                                                                            >  2:236863/1‑148 (MQ=255)
                                ggAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACg                                                                                                            <  1:2147058/148‑1 (MQ=255)
                                 gAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACga                                                                                                           >  2:3437549/1‑148 (MQ=255)
                                    ttattCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGAt                                                                                                        <  1:2765818/148‑1 (MQ=255)
                                     tattCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATc                                                                                                       <  1:2089907/148‑1 (MQ=255)
                                     tattCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATc                                                                                                       >  1:115611/1‑148 (MQ=255)
                                     tattCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATc                                                                                                       <  1:1090594/148‑1 (MQ=255)
                                       ttCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCt                                                                                                     <  2:268053/148‑1 (MQ=255)
                                       ttCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCt                                                                                                     <  1:216147/148‑1 (MQ=255)
                                       ttCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCt                                                                                                     >  2:1781966/1‑148 (MQ=255)
                                            ccAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAAt                                                                                                <  2:1407435/148‑1 (MQ=255)
                                             cAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAAtt                                                                                               >  2:2804577/1‑148 (MQ=255)
                                              aGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTa                                                                                              <  1:431933/148‑1 (MQ=255)
                                              aGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCGATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTa                                                                                              <  2:1763070/148‑1 (MQ=255)
                                              aGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCGATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTa                                                                                              <  2:2933160/148‑1 (MQ=255)
                                               ggTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTaa                                                                                             >  2:1384911/1‑148 (MQ=255)
                                               ggTGGCAAGCGCTAGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTaa                                                                                             >  2:754374/1‑148 (MQ=255)
                                                gTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAAt                                                                                            <  2:2741346/148‑1 (MQ=255)
                                                gTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAAt                                                                                            >  2:2853578/1‑148 (MQ=255)
                                                 tGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAAtg                                                                                           <  1:2617988/148‑1 (MQ=255)
                                                 tGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAAtg                                                                                           <  2:1727021/148‑1 (MQ=255)
                                                 tGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAAtg                                                                                           <  2:860631/148‑1 (MQ=255)
                                                       gcgcTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATcc                                                                                                      >  1:1750612/1‑131 (MQ=255)
                                                        cgcTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACCGCAAAAACGGGCGAAAGGTGGCTATGGGGCTGGCAGGGCTGTTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCaa                                                                                    >  1:1040347/1‑148 (MQ=255)
                                                         gcTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCaaa                                                                                   <  1:2829515/148‑1 (MQ=255)
                                                            cGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACa                                                                                <  1:2928203/148‑1 (MQ=255)
                                                            cGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACa                                                                                <  1:2927812/148‑1 (MQ=255)
                                                               tGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTg                                                                             <  1:771096/148‑1 (MQ=255)
                                                                gACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTgg                                                                            >  2:848877/1‑148 (MQ=255)
                                                                  cgttgaacgttgaTCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTgggg                                                                          >  2:2604694/1‑148 (MQ=255)
                                                                  cgttgaacgttgaTCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTgggg                                                                          >  1:2879386/1‑148 (MQ=255)
                                                                    ttgaacgttgaTCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTa                                                                        <  2:2560408/148‑1 (MQ=255)
                                                                    ttgaacgttgaTCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTa                                                                        <  2:3289587/148‑1 (MQ=255)
                                                                     tgaacgttgaTCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTAt                                                                       <  2:1238132/148‑1 (MQ=255)
                                                                     tgaacgttgaTCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCGATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTAt                                                                       <  1:197577/148‑1 (MQ=255)
                                                                      gaacgttgaTCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATTACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATc                                                                      <  2:3168615/148‑1 (MQ=255)
                                                                         cgttgaTCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCgat                                                                   >  1:76334/1‑148 (MQ=255)
                                                                         cgttgaTCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCGATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCgat                                                                   <  1:1303391/148‑1 (MQ=255)
                                                                          gttgaTCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCgatg                                                                  <  2:683966/148‑1 (MQ=255)
                                                                          gttgaTCAACCCGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCgatg                                                                  <  1:2858389/148‑1 (MQ=255)
                                                                           ttgaTCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCgatga                                                                 <  1:2976065/148‑1 (MQ=255)
                                                                               tCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGtt                                                             <  1:1471355/148‑1 (MQ=255)
                                                                                  aCCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCg                                                          >  2:658303/1‑148 (MQ=255)
                                                                                         gCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCg                                                   <  1:236863/148‑1 (MQ=255)
                                                                                            aCCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTc                                                <  1:1176844/148‑1 (MQ=255)
                                                                                             ccTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCa                                               >  2:491801/1‑148 (MQ=255)
                                                                                              ctgctgGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAg                                              <  2:3320594/148‑1 (MQ=255)
                                                                                                  tgGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATa                                          >  2:3222140/1‑148 (MQ=255)
                                                                                                   ggTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATaa                                         >  1:657843/1‑148 (MQ=255)
                                                                                                      tCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATAAGaa                                      <  2:3024537/148‑1 (MQ=255)
                                                                                                      tCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATAAGaa                                      >  1:1033194/1‑148 (MQ=255)
                                                                                                      tCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATAAGaa                                      >  1:577146/1‑148 (MQ=255)
                                                                                                         aTCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATAAGAAAtt                                   <  2:1865098/148‑1 (MQ=255)
                                                                                                         aTCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCGATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATAAGAAAtt                                   <  2:1181818/148‑1 (MQ=255)
                                                                                                           ccGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATAAGAAATTTa                                 >  1:905332/1‑148 (MQ=255)
                                                                                                           ccGCATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATAAGAAATTTa                                 >  2:791024/1‑148 (MQ=255)
                                                                                                              cATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATAAGAAATTTAgcg                              >  1:3185144/1‑148 (MQ=255)
                                                                                                              cATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATAAGAAATTTAgcg                              >  1:3202160/1‑148 (MQ=255)
                                                                                                              cATCAGCACGGCAAAACCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATAAGAAATTTAgcg                              <  2:3471783/148‑1 (MQ=255)
                                                                                                                            aaCCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCGGCCAAAACAGTGGGGTCTCGATGATGTTCCGGGGATCGTTCAGGATAAGAATTTTAGCGCCGACGGGCAgctg                >  2:2535996/1‑145 (MQ=255)
                                                                                                                            aaCCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCGGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATAAGAAATTTAGCGCCGACGGGCAgatt                >  2:3036024/1‑148 (MQ=255)
                                                                                                                             aCCGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATAAGAAATTTAGCGCCGACGGGCAgattg               <  1:1861496/148‑1 (MQ=255)
                                                                                                                              ccGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATAAGAAATTTAGCGCCGACGGGCAgattga              >  2:107818/1‑148 (MQ=255)
                                                                                                                              ccGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATAAGAAATTTAGCGCCGACGGGCAgattga              <  1:1120094/148‑1 (MQ=255)
                                                                                                                              ccGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATAAGAAATTTAGCGCCGACGGGCAgattga              <  1:848877/148‑1 (MQ=255)
                                                                                                                               cGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATAAGAAATTTAGCGCCGACGGGCAgattgat             >  2:44877/1‑148 (MQ=255)
                                                                                                                               cGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATAAGAAATTTAGCGCCGACGGGCAgattgat             <  2:3409976/148‑1 (MQ=255)
                                                                                                                               cGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATAAGAAATTTAGCGCCGACGGGCAgattgat             <  1:2937347/148‑1 (MQ=255)
                                                                                                                               cGGGCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATAAGAAATTTAGCGCCGACGGGCAgattgat             <  2:2700851/148‑1 (MQ=255)
                                                                                                                                gggCGACAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATAAGAAATTTAGCGCCGACGGGCAgattgatt            <  1:335077/148‑1 (MQ=255)
                                                                                                                                      cAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATAAGAAATTTAGCGCCGACGGGCAGATTGATTATCAAc      >  2:1922755/1‑148 (MQ=255)
                                                                                                                                      cAGGTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATAAGAAATTTAGCGCCGACGGGCAGATTGATTATCAAc      >  2:2255135/1‑148 (MQ=255)
                                                                                                                                         gTGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATAAGAAATTTAGCGCCGACGGGCAGATTGATTATCAACTgg   >  1:1430327/1‑148 (MQ=255)
                                                                                                                                          tGGCTATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATAAGAAATTTAGCGCCGACGGGCAGATTGATTATCAACTGGa  >  2:305163/1‑148 (MQ=255)
                                                                                                                                              |                                                                                                                                               
AAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCGATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATAAGAAATTTAGCGCCGACGGGCAGATTGATTATCAACTGGA  >  CP009273/133886‑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: