Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,835,094 C→A 100% P249Q (CCG→CAG)  ascB → cryptic 6‑phospho‑beta‑glucosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,835,0940CA98.7% 494.9 / ‑1.6 156P249Q (CCG→CAG) ascBcryptic 6‑phospho‑beta‑glucosidase
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base A (89/65);  total (91/65)
Fisher's exact test for biased strand distribution p-value = 5.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-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 156 total aligned reads displayed.

AATCAGGATCAGGTGAAATATCAGGCCGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCCGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCATCAACAAAGCACCGGGCGATG  >  CP009273/2834949‑2835234
                                                                                                                                                 |                                                                                                                                            
aaTCAGGATCAGGTGAAATATCAGGCCGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAgg                                                                                                                                            >  2:1766819/1‑148 (MQ=255)
 atcaggatcaggTGAAATATCAGGCCGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGa                                                                                                                                           >  2:2061826/1‑148 (MQ=255)
   caggatcaggTGAAATATCAGGCCGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAg                                                                                                                                         >  2:2752795/1‑148 (MQ=255)
   caggatcaggTGAAATATCAGGCCGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAg                                                                                                                                         >  2:2752989/1‑148 (MQ=255)
    aggatcaggTGAAATATCAGGCCGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGa                                                                                                                                        >  1:741491/1‑148 (MQ=255)
       atcaggTGAAATATCAGGCCGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGt                                                                                                                                     <  2:275708/148‑1 (MQ=255)
           ggTGAAATATCAGGCCGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTgg                                                                                                                                 <  1:3227528/148‑1 (MQ=255)
                 atatCAGGCCGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGgcggcg                                                                                                                           >  2:602297/1‑148 (MQ=255)
                 atatCAGGCCGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGgcggcg                                                                                                                           >  2:2153032/1‑148 (MQ=255)
                 atatCAGGCCGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGATAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGgcggcg                                                                                                                           >  2:600522/1‑148 (MQ=255)
                  tatCAGGCCGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGgcgc                                                                                                                          >  2:3275428/1‑148 (MQ=255)
                  tatCAGGCCGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGgcgc                                                                                                                          <  2:2322798/148‑1 (MQ=255)
                   atCAGGCCGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCt                                                                                                                         <  1:3591106/148‑1 (MQ=255)
                     cAGGCCGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTgg                                                                                                                       >  1:3401353/1‑148 (MQ=255)
                        gCCGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGgaga                                                                                                                    >  2:2897785/1‑148 (MQ=255)
                           gcgcATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAg                                                                                                                 <  1:2385629/148‑1 (MQ=255)
                            cgcATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGa                                                                                                                <  2:2493335/148‑1 (MQ=255)
                               aTCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCg                                                                                                             <  1:2514884/148‑1 (MQ=255)
                               aTCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCg                                                                                                             <  1:3298942/148‑1 (MQ=255)
                                tCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCgg                                                                                                            >  1:1815769/1‑148 (MQ=255)
                                tCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCgg                                                                                                            <  1:1773781/148‑1 (MQ=255)
                                 cACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCggg                                                                                                           >  1:3543990/1‑148 (MQ=255)
                                 cACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCggg                                                                                                           >  1:50606/1‑148 (MQ=255)
                                 cACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCggg                                                                                                           >  1:3544392/1‑148 (MQ=255)
                                        tGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCt                                                                                                    >  2:3252089/1‑148 (MQ=255)
                                            tGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGttt                                                                                                <  1:1281225/148‑1 (MQ=255)
                                            tGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGttt                                                                                                <  2:3070448/148‑1 (MQ=255)
                                                aGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTa                                                                                            <  2:1200069/148‑1 (MQ=255)
                                                 gTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTAt                                                                                           <  2:2642816/148‑1 (MQ=255)
                                                  tGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATc                                                                                          >  1:37177/1‑148 (MQ=255)
                                                    cgcTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGACAGACCGGGAAAACCTGTTTTTTATCGa                                                                                        >  2:3413374/1‑148 (MQ=255)
                                                    cgcTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGACAACCTGTTTTTTATCGa                                                                                        >  2:3413707/1‑148 (MQ=255)
                                                      cTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGAtg                                                                                      >  1:1122529/1‑148 (MQ=255)
                                                        aGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGAtgtg                                                                                    <  1:443582/148‑1 (MQ=255)
                                                          ccaccaAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCa                                                                                  <  2:3384566/148‑1 (MQ=255)
                                                          ccaccaAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCa                                                                                  <  2:1828828/148‑1 (MQ=255)
                                                             ccaAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGgc                                                                               >  2:1406855/1‑148 (MQ=255)
                                                                 aaTCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCgg                                                                           >  2:1497228/1‑148 (MQ=255)
                                                                 aaTCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCgg                                                                           >  2:1498266/1‑148 (MQ=255)
                                                                 aaTCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCgg                                                                           >  2:2078661/1‑148 (MQ=255)
                                                                   tCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCgggg                                                                         <  1:262225/148‑1 (MQ=255)
                                                                   tCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCgggg                                                                         <  1:1615253/148‑1 (MQ=255)
                                                                        cATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGt                                                                    >  1:3536862/1‑148 (MQ=255)
                                                                        cATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGATGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACAGGGACAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGt                                                                    >  2:2276475/1‑148 (MQ=255)
                                                                          tGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTAt                                                                  >  1:1705644/1‑148 (MQ=255)
                                                                          tGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTAt                                                                  >  1:1945092/1‑148 (MQ=255)
                                                                          tGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTAt                                                                  >  1:1945074/1‑148 (MQ=255)
                                                                          tGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTAt                                                                  >  1:1945019/1‑148 (MQ=255)
                                                                          tGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTAt                                                                  >  1:1561769/1‑148 (MQ=255)
                                                                           gAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATc                                                                 >  2:1575765/1‑148 (MQ=255)
                                                                            aGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATcc                                                                <  1:2271367/148‑1 (MQ=255)
                                                                            aGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATcc                                                                <  1:345814/148‑1 (MQ=255)
                                                                             ggTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCg                                                               >  1:2117882/1‑148 (MQ=255)
                                                                              gTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCgg                                                              >  2:3137493/1‑148 (MQ=255)
                                                                               ttAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGc                                                             <  1:1678324/148‑1 (MQ=255)
                                                                               ttAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGc                                                             <  2:3401353/148‑1 (MQ=255)
                                                                                  ccccGCCAAATCACCTGTCCTGCACGCAGCCGGTCGGTAACGTCAACCCTTACACTTCCGAGCACGAAGATCCCTGGGCCGCGCTCGAGAAAGACCCGGAAAACCTGTTTTTTACCGATGTGAAGGCGCGGGGCACGTATCCGGCTTa                                                          <  2:2395993/147‑1 (MQ=255)
                                                                                  aCCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTa                                                          >  1:2379265/1‑148 (MQ=255)
                                                                                  aCCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTa                                                          <  2:115138/148‑1 (MQ=255)
                                                                                  aCCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTa                                                          <  2:3215762/148‑1 (MQ=255)
                                                                                  aCCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTa                                                          <  2:3215797/148‑1 (MQ=255)
                                                                                   cccGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTAc                                                         >  1:3281366/1‑148 (MQ=255)
                                                                                    ccGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTAct                                                        <  1:1952560/148‑1 (MQ=255)
                                                                                    ccGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTAct                                                        <  2:2879287/148‑1 (MQ=255)
                                                                                     cGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTActc                                                       >  1:2893763/1‑148 (MQ=255)
                                                                                        aaaaTCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGc                                                    <  2:1239583/148‑1 (MQ=255)
                                                                                        aaaaTCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGc                                                    <  1:2460555/148‑1 (MQ=255)
                                                                                          aaTCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGccc                                                  >  2:1216653/1‑148 (MQ=255)
                                                                                          aaTCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGccc                                                  >  1:2393730/1‑148 (MQ=255)
                                                                                            tCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCcgc                                                >  1:2151611/1‑148 (MQ=255)
                                                                                             cAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCcgcg                                               >  1:900865/1‑148 (MQ=255)
                                                                                                gTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTAt                                            >  1:1133794/1‑148 (MQ=255)
                                                                                                gTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTAt                                            <  1:1321187/148‑1 (MQ=255)
                                                                                                          tGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGaaaa                                  <  1:602297/148‑1 (MQ=255)
                                                                                                          tGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGaaaa                                  <  1:2153032/148‑1 (MQ=255)
                                                                                                           gCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGaaaaa                                 >  2:1917237/1‑148 (MQ=255)
                                                                                                           gCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGaaaaa                                 <  2:1571542/148‑1 (MQ=255)
                                                                                                            cTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAg                                <  1:1112412/148‑1 (MQ=255)
                                                                                                             tGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAgg                               <  2:1421212/148‑1 (MQ=255)
                                                                                                             tGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAgg                               <  2:741491/148‑1 (MQ=255)
                                                                                                             tGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAgg                               <  2:2472603/148‑1 (MQ=255)
                                                                                                              ggcgggcgGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAggg                              <  1:3226395/148‑1 (MQ=255)
                                                                                                               gcgggcggTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAgggg                             >  2:1066221/1‑148 (MQ=255)
                                                                                                               gcgggcggTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAgggg                             >  1:3651803/1‑148 (MQ=255)
                                                                                                                        aaCTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCATca                    >  1:544981/1‑148 (MQ=255)
                                                                                                                         aCTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCATcaa                   >  2:2682636/1‑148 (MQ=255)
                                                                                                                         aCTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCATcaa                   <  2:3467733/148‑1 (MQ=255)
                                                                                                                            tCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCATcaacaa                >  1:55937/1‑148 (MQ=255)
                                                                                                                            tCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCATcaacaa                <  1:1469775/148‑1 (MQ=255)
                                                                                                                                 ccTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCATCAACAAAGCAc           >  1:2542831/1‑148 (MQ=255)
                                                                                                                                  cTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCATCAACAAAGCAcc          >  1:3007386/1‑148 (MQ=255)
                                                                                                                                   ttACAGTTGCAAGCCGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCATCAACAAAGCACCg         >  2:1376648/1‑148 (MQ=255)
                                                                                                                                   ttACAGTTGCAAGCCGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCATCAACAAAGCACCg         >  2:1375885/1‑148 (MQ=255)
                                                                                                                                     aCAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCATCAACAAAGCACCggg       <  2:1068685/148‑1 (MQ=255)
                                                                                                                                      cAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCATCAACACAGCACCGGGc      >  2:2087813/1‑148 (MQ=255)
                                                                                                                                        gTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCATCAACAAAGCACCGGGCga    <  1:2296209/148‑1 (MQ=255)
                                                                                                                                         ttGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCATCAACAAAGCACCGGGCgat   >  1:2325902/1‑148 (MQ=255)
                                                                                                                                          tGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCATCAACAAAGCACCGGGCgatg  >  1:1204923/1‑148 (MQ=255)
                                                                                                                                          tGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCATCAACAAAGCACCGGGCgatg  <  2:1354180/148‑1 (MQ=255)
                                                                                                                                          tGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCATCAACAAAGCACCGGGCgatg  <  1:513913/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                            
AATCAGGATCAGGTGAAATATCAGGCCGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCCGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCATCAACAAAGCACCGGGCGATG  >  CP009273/2834949‑2835234

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: