Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,835,094 C→A 89.5% 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,0940CA89.5% 457.4 / 51.4 181P249Q (CCG→CAG) ascBcryptic 6‑phospho‑beta‑glucosidase
Reads supporting (aligned to +/- strand):  ref base C (10/9);  new base A (76/86);  total (86/95)
Fisher's exact test for biased strand distribution p-value = 8.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

Only 100 of 181 total aligned reads displayed.

AATCAGGATCAGGTGAAATATCAGGCCGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCCGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCATCAACAAAGCACCGGGCGATGATGAA  >  CP009273/2834949‑2835239
                                                                                                                                                 |                                                                                                                                                 
aaTCAGGATCAGGTGAAATATCAGGCCGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAgg                                                                                                                                                 >  1:1248941/1‑148 (MQ=255)
 atcaggatcatgTGAAATATCAGGCCGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGa                                                                                                                                                <  2:589100/148‑1 (MQ=255)
  tcaggatcaggTGAAATATCAGGCCGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGaa                                                                                                                                               >  2:2458622/1‑148 (MQ=255)
   caggatcaggTGAAATATCAGGCCGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAg                                                                                                                                              <  2:678946/148‑1 (MQ=255)
    aggatcaggTGAAATATCAGGCCGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGa                                                                                                                                             <  2:874051/148‑1 (MQ=255)
     ggatcaggTGAAATATCAGGCCGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGAt                                                                                                                                            >  2:1428295/1‑148 (MQ=255)
     ggatcaggTGAAATATCAGGCCGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGAt                                                                                                                                            >  2:2969675/1‑148 (MQ=255)
     ggatcaggTGAAATATCAGGCCGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGAt                                                                                                                                            >  2:2969777/1‑148 (MQ=255)
     ggatcaggTGAAATATCAGGCCGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGAt                                                                                                                                            >  2:474930/1‑148 (MQ=255)
      gatcaggTGAAATATCAGGCCGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATg                                                                                                                                           <  1:2325537/148‑1 (MQ=255)
      gatcaggTGAAATATCAGGCCGCGCATCACCAGCTGGTGGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATg                                                                                                                                           <  1:2811454/148‑1 (MQ=255)
       atcaggTGAAATATCAGGCCGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGt                                                                                                                                          >  1:662225/1‑148 (MQ=255)
       atcaggTGAAATATCAGGCCGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGt                                                                                                                                          >  2:2348621/1‑148 (MQ=255)
        tcaggTGAAATATCAGGCCGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTc                                                                                                                                         <  1:2158452/148‑1 (MQ=255)
           ggTGAAATATCAGGCCGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTgg                                                                                                                                      <  1:1795283/148‑1 (MQ=255)
             tGAAATATCAGGCCGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCCGGAAGATGTCTGggc                                                                                                                                    >  1:3300691/1‑148 (MQ=255)
             tGAAATATCAGGCCGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGggc                                                                                                                                    >  1:1379692/1‑148 (MQ=255)
             tGAAATATCAGGCCGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGggc                                                                                                                                    >  1:321657/1‑148 (MQ=255)
               aaaTATCAGGCCGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGggcgg                                                                                                                                  <  1:1640082/148‑1 (MQ=255)
               aaaTATCAGGCCGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGggcgg                                                                                                                                  <  2:3272766/148‑1 (MQ=255)
               aaaTATCAGGCCGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGggcgg                                                                                                                                  <  1:298272/148‑1 (MQ=255)
                 atatCAGGCCGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGgcggcg                                                                                                                                <  1:119609/148‑1 (MQ=255)
                  tatCAGGCCGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGgcgc                                                                                                                               >  1:590515/1‑148 (MQ=255)
                   atCAGGCCGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCt                                                                                                                              >  2:40728/1‑148 (MQ=255)
                        gCCGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGgaga                                                                                                                         >  2:1368463/1‑148 (MQ=255)
                         ccGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCCGGAAGATGTCTGGGCGGCGCTGGAGaa                                                                                                                        >  1:3063090/1‑148 (MQ=255)
                         ccGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGaa                                                                                                                        <  2:3136919/148‑1 (MQ=255)
                           gcgcATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAg                                                                                                                      <  1:1309171/148‑1 (MQ=255)
                             gcATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGAc                                                                                                                    <  2:869673/148‑1 (MQ=255)
                              cATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGAcc                                                                                                                   <  2:3192809/148‑1 (MQ=255)
                              cATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGCTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGAcc                                                                                                                   <  2:402913/148‑1 (MQ=255)
                                  aCCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGa                                                                                                               >  2:944033/1‑148 (MQ=255)
                                   ccAGCTGGGTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGaa                                                                                                              <  2:1973921/148‑1 (MQ=255)
                                      gCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAAc                                                                                                           >  2:2163563/1‑148 (MQ=255)
                                      gCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAAc                                                                                                           <  1:1557865/148‑1 (MQ=255)
                                        tGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCt                                                                                                         <  2:1604611/148‑1 (MQ=255)
                                         ggTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTg                                                                                                        >  1:2915526/1‑148 (MQ=255)
                                           ttGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGtt                                                                                                      <  1:2681282/148‑1 (MQ=255)
                                           ttGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGtt                                                                                                      >  2:2132842/1‑148 (MQ=255)
                                            tGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGttt                                                                                                     >  1:1464229/1‑148 (MQ=255)
                                             gCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGtttt                                                                                                    <  1:2997124/148‑1 (MQ=255)
                                             gCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGtttt                                                                                                    <  1:2798254/148‑1 (MQ=255)
                                                aGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTa                                                                                                 <  1:2752497/148‑1 (MQ=255)
                                                 gTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTAt                                                                                                <  2:2493000/148‑1 (MQ=255)
                                                    cgcTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCCGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGa                                                                                             <  1:1065456/148‑1 (MQ=255)
                                                          ccaccaAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCa                                                                                       >  2:2105533/1‑148 (MQ=255)
                                                           caccaAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAg                                                                                      >  2:2241087/1‑148 (MQ=255)
                                                              caAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGgcg                                                                                   <  1:2809439/148‑1 (MQ=255)
                                                                  aTCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCggg                                                                               >  2:693948/1‑148 (MQ=255)
                                                                  aTCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCggg                                                                               <  1:2138043/148‑1 (MQ=255)
                                                                   tCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCgggg                                                                              <  2:1277781/148‑1 (MQ=255)
                                                                         aTGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTa                                                                        >  1:407555/1‑148 (MQ=255)
                                                                         aTGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGAAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTa                                                                        >  2:250920/1‑148 (MQ=255)
                                                                          tGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTAt                                                                       <  1:2216801/148‑1 (MQ=255)
                                                                           gAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATc                                                                      >  1:944878/1‑148 (MQ=255)
                                                                            aGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATcc                                                                     <  2:2312974/148‑1 (MQ=255)
                                                                                  aCCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTa                                                               <  1:1874041/148‑1 (MQ=255)
                                                                                    ccGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCCGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTAct                                                             <  2:2859672/148‑1 (MQ=255)
                                                                                       cAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCCGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTg                                                          >  1:2041089/1‑148 (MQ=255)
                                                                                       cAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTg                                                          <  2:1843887/148‑1 (MQ=255)
                                                                                         aaaTCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGcc                                                        >  2:742850/1‑148 (MQ=255)
                                                                                         aaaTCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGcc                                                        >  2:2402251/1‑148 (MQ=255)
                                                                                          aaTCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGccc                                                       >  1:2313052/1‑148 (MQ=255)
                                                                                           aTCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCCGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCcg                                                      >  2:1694295/1‑148 (MQ=255)
                                                                                               ggTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCCGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTa                                                  >  1:821084/1‑148 (MQ=255)
                                                                                               ggTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTa                                                  <  2:217824/148‑1 (MQ=255)
                                                                                                gTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCCTGAAGATGTCTGGGAGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCCGCTTACTCTGCCCGCGTAt                                                 >  1:2669365/1‑148 (MQ=255)
                                                                                                gTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTAt                                                 <  1:2540409/148‑1 (MQ=255)
                                                                                                    gCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCcg                                             >  1:838462/1‑148 (MQ=255)
                                                                                                    gCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCcg                                             <  1:102679/148‑1 (MQ=255)
                                                                                                           gCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGaaaaa                                      >  2:2740949/1‑148 (MQ=255)
                                                                                                           gCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGaaaaa                                      >  2:608843/1‑148 (MQ=255)
                                                                                                           gCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGAGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGaaaaa                                      <  1:777530/148‑1 (MQ=255)
                                                                                                            cTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAg                                     <  2:260168/148‑1 (MQ=255)
                                                                                                               gcgggcggTAACTTCTACCCTTACAGTTGCAAGCCGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAgggg                                  <  2:1079490/148‑1 (MQ=255)
                                                                                                               gcgggcggTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAgggg                                  <  2:2629834/148‑1 (MQ=255)
                                                                                                               gcgggcggTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAgggg                                  <  1:2583124/148‑1 (MQ=255)
                                                                                                                cgggcggTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGt                                 <  2:132595/148‑1 (MQ=255)
                                                                                                                  ggcggTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTaa                               <  2:1464660/148‑1 (MQ=255)
                                                                                                                  ggcggTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTaa                               <  2:661803/148‑1 (MQ=255)
                                                                                                                  ggcggTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTaa                               <  2:661740/148‑1 (MQ=255)
                                                                                                                   gcggTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAAc                              <  1:1312534/148‑1 (MQ=255)
                                                                                                                      gTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCAt                           >  2:830243/1‑148 (MQ=255)
                                                                                                                        aaCTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCATca                         <  1:1139887/148‑1 (MQ=255)
                                                                                                                         aCTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCATcaa                        >  2:2764728/1‑148 (MQ=255)
                                                                                                                              tACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCATCAACAAAg                   <  1:3078164/148‑1 (MQ=255)
                                                                                                                                  cTTACAGTTGCAAGCCGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCATCAACAAAGCAcc               >  2:1074232/1‑148 (MQ=255)
                                                                                                                                  cTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCATCAACAAAGCAcc               >  2:1361171/1‑148 (MQ=255)
                                                                                                                                     aCAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCATCAACAAAGCACCggg            >  2:609753/1‑148 (MQ=255)
                                                                                                                                     aCAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCATCAACAAAGCACCggg            <  2:2427096/148‑1 (MQ=255)
                                                                                                                                       aGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCATCAACAAAGCACCGGGCg          >  2:1471891/1‑148 (MQ=255)
                                                                                                                                       aGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCATCAACAAAGCACCGGGCg          >  2:1471487/1‑148 (MQ=255)
                                                                                                                                        gTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCATCAACAAAGCACCGGGCga         >  2:559021/1‑148 (MQ=255)
                                                                                                                                          tGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCATCAACAAAGCACCGGGCgatg       >  2:2220450/1‑148 (MQ=255)
                                                                                                                                          tGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCATCAACAAAGCACCGGGCgatg       >  2:3088991/1‑148 (MQ=255)
                                                                                                                                           gCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCATCAACAAAGCACCGGGCgatga      <  2:2915526/148‑1 (MQ=255)
                                                                                                                                             aaGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCATCAACAAAGCACCGGGCGatgatg    <  1:1573751/148‑1 (MQ=255)
                                                                                                                                              aGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCATCAACAAAGCACCGGGCGAtgatga   >  2:2762228/1‑148 (MQ=255)
                                                                                                                                               gCCGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCATCAACAAAGCACCGGGCGATGATGaa  <  2:1909563/148‑1 (MQ=255)
                                                                                                                                               gCCGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCATCAACAAAGCACCGGGCGATGATGaa  >  2:3095605/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
AATCAGGATCAGGTGAAATATCAGGCCGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCCGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCATCAACAAAGCACCGGGCGATGATGAA  >  CP009273/2834949‑2835239

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: