Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,835,094 C→A 37.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,0940CA37.5% 438.8 / inf 316P249Q (CCG→CAG) ascBcryptic 6‑phospho‑beta‑glucosidase
Reads supporting (aligned to +/- strand):  ref base C (104/93);  new base A (60/59);  total (164/152)
Fisher's exact test for biased strand distribution p-value = 7.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 318 total aligned reads displayed.

TGAAAATCAGGATCAGGTGAAATATCAGGCCGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCCGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCATCAACAAAGCACCGGGCGATGATGAAA  >  CP009273/2834945‑2835240
                                                                                                                                                     |                                                                                                                                                  
   aaaTCAGGATCAGGTGAAATATCAGGCCGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCCGGa                                                                                                                                                 >  1:1435526/1‑150 (MQ=255)
   aaaTCAGGATCAGGTGAAATATCAGGCCGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGa                                                                                                                                                 >  1:3209746/1‑150 (MQ=255)
     atcaggatcaggTGAAATATCAGGCCGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCCGGAAg                                                                                                                                               <  1:5581/150‑1 (MQ=255)
     atcaggatcaggTGAAATATCAGGCCGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCCGGAAg                                                                                                                                               <  2:4445994/150‑1 (MQ=255)
       caggatcaggTGAAATATCAGGCCGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCCGGAAGAt                                                                                                                                             <  2:3910371/150‑1 (MQ=255)
         ggatcaggTGAAATATCAGGCCGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCTGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGt                                                                                                                                           <  1:4610050/150‑1 (MQ=255)
          gatcaggTGAAATATCAGGCCGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCCGGAAGATGTc                                                                                                                                          <  2:3692005/150‑1 (MQ=255)
          gatcaggTGAAATATCAGGCCGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCCGGAAGATGTc                                                                                                                                          >  1:532069/1‑150 (MQ=255)
           atcaggTGAAATATCAGGCCGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCCGGAAGATGTCt                                                                                                                                         >  1:4054822/1‑150 (MQ=255)
            tcaggTGAAATATCAGGCCGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCACAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCCGGAAGATGTCTg                                                                                                                                        <  2:2691729/150‑1 (MQ=255)
            tcaggTGAAATATCAGGCCGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCCGGAAGATGTCTg                                                                                                                                        >  2:4375037/1‑150 (MQ=255)
               ggTGAAATATCAGGCCGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCCGGAAGATGTCTGggc                                                                                                                                     <  2:4403509/150‑1 (MQ=255)
               ggTGAAATATCAGGCCGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGggc                                                                                                                                     >  1:1947273/1‑150 (MQ=255)
                  gAAATATCAGGCCGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGggcggc                                                                                                                                  <  1:3087344/150‑1 (MQ=255)
                   aaaTATCAGGCCGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCCGGAAGATGTCTGGgcggcg                                                                                                                                 >  1:2362104/1‑150 (MQ=255)
                    aaTATCAGGCCGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGgcgc                                                                                                                                >  1:3766325/1‑150 (MQ=255)
                          aGGCCGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCCGGAAGATGTCTGGGCGGCGCTGgaga                                                                                                                          <  1:817162/150‑1 (MQ=255)
                           ggCCGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCCGGAAGATGTCTGGGCGGCGCTGGAGaa                                                                                                                         <  2:2294383/150‑1 (MQ=255)
                             ccGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAg                                                                                                                       <  2:673343/150‑1 (MQ=255)
                             ccGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAg                                                                                                                       <  2:1693489/150‑1 (MQ=255)
                             ccGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAg                                                                                                                       <  1:798578/150‑1 (MQ=255)
                               gcgcATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGAc                                                                                                                     >  2:4250243/1‑150 (MQ=255)
                                cgcATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGAcc                                                                                                                    >  2:2324169/1‑150 (MQ=255)
                                cgcATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGAGGTCTGGGCGGCGCTGGAGAAAGAcc                                                                                                                    >  2:3105133/1‑150 (MQ=255)
                                 gcATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCACGCCGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCg                                                                                                                   <  2:2111400/150‑1 (MQ=255)
                                 gcATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCCGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCg                                                                                                                   >  1:1359504/1‑150 (MQ=255)
                                 gcATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCg                                                                                                                   <  1:3384560/150‑1 (MQ=255)
                                 gcATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCg                                                                                                                   >  2:4709100/1‑150 (MQ=255)
                                 gcATCACCAGCTGGGTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCCGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCg                                                                                                                   <  2:4666361/150‑1 (MQ=255)
                                    tCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCCGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGa                                                                                                                <  1:4266506/150‑1 (MQ=255)
                                     cACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCCGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGaa                                                                                                               >  2:3485426/1‑150 (MQ=255)
                                     cACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGaa                                                                                                               >  1:438596/1‑150 (MQ=255)
                                       ccAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCCGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGaaaa                                                                                                             >  2:4832812/1‑150 (MQ=255)
                                       ccAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGaaaa                                                                                                             <  1:4992792/150‑1 (MQ=255)
                                         aGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCCGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAAcc                                                                                                           <  1:3884474/150‑1 (MQ=255)
                                              gTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCCGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGttt                                                                                                      <  1:1163640/150‑1 (MQ=255)
                                                 gCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGtttttt                                                                                                   >  1:3845692/1‑150 (MQ=255)
                                                    aGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCCGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATc                                                                                                >  1:2696675/1‑150 (MQ=255)
                                                     gTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCg                                                                                               <  2:3888022/150‑1 (MQ=255)
                                                            aGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCCGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCa                                                                                        <  1:3735275/150‑1 (MQ=255)
                                                                accaAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCCGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGgcg                                                                                    <  1:1272527/150‑1 (MQ=255)
                                                                accaAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGgcg                                                                                    >  1:1200928/1‑150 (MQ=255)
                                                                accaAAAGCGCCCATGAGGTTATCCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCCGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGgcn                                                                                    <  2:5033148/150‑2 (MQ=255)
                                                                  caAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCCGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGcgcg                                                                                  <  2:2403462/150‑1 (MQ=255)
                                                                     aaTAGCCCATGAGGTTAACCCCCAAAAAAAGGTGGGCTGTATGCTGGCGGGCCGTAACATCTACCCGTACAGTTTCAACCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCCGGAAAACCTGTTTTTTATAGATGTGCAGGCGCgggg                                                                               >  2:1057775/1‑150 (MQ=255)
                                                                         gCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACg                                                                           <  1:5094771/150‑1 (MQ=255)
                                                                           ccATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCCGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTa                                                                         <  2:1369271/150‑1 (MQ=255)
                                                                           ccATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCCGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTa                                                                         <  2:3933432/150‑1 (MQ=255)
                                                                           ccATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCCGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTa                                                                         >  2:3320262/1‑150 (MQ=255)
                                                                           ccATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCCGGAAGATGTCTGGGCCGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTa                                                                         >  1:4563098/1‑150 (MQ=255)
                                                                              tGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATcc                                                                      >  2:1154485/1‑150 (MQ=255)
                                                                               gAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCCGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCg                                                                     >  1:4812417/1‑150 (MQ=255)
                                                                                aGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCCGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCgg                                                                    >  1:1865628/1‑150 (MQ=255)
                                                                                aGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCCGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCgg                                                                    >  2:2798721/1‑150 (MQ=255)
                                                                                  gTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCCGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCt                                                                  >  2:4067407/1‑150 (MQ=255)
                                                                                  gTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCCGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCt                                                                  <  1:1126429/150‑1 (MQ=255)
                                                                                  gTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCt                                                                  <  1:4425416/150‑1 (MQ=255)
                                                                                    tAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCCGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTa                                                                >  1:616843/1‑150 (MQ=255)
                                                                                    tAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCCGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTa                                                                <  1:4945037/150‑1 (MQ=255)
                                                                                    tAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCCGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTa                                                                >  2:1531154/1‑150 (MQ=255)
                                                                                    tAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCCGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTa                                                                >  2:3570560/1‑150 (MQ=255)
                                                                                        ccGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTActct                                                            >  1:2290671/1‑150 (MQ=255)
                                                                                         cGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCCGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTg                                                           <  2:2751392/150‑1 (MQ=255)
                                                                                           cAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCCGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGcc                                                         <  1:2316728/150‑1 (MQ=255)
                                                                                            aaaaTCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGccc                                                        >  1:2945760/1‑150 (MQ=255)
                                                                                                 cAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCCGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTa                                                   <  1:2460928/150‑1 (MQ=255)
                                                                                                 cAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTa                                                   <  1:1745809/150‑1 (MQ=255)
                                                                                                     tGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCCGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTcc                                               >  2:4730734/1‑150 (MQ=255)
                                                                                                     tGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCCGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTcc                                               <  2:3528132/150‑1 (MQ=255)
                                                                                                      gggCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCCGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCcg                                              <  1:3010230/150‑1 (MQ=255)
                                                                                                       ggCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCcgc                                             <  1:326236/150‑1 (MQ=255)
                                                                                                          tGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCCGGAAGATGTCTGGGCGGCGCTGGGGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGaa                                          <  1:3516678/150‑1 (MQ=255)
                                                                                                             aTGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCCGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGaaaaa                                       <  1:4397351/150‑1 (MQ=255)
                                                                                                              tGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCCGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAg                                      >  1:103280/1‑150 (MQ=255)
                                                                                                              tGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCCGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAg                                      >  1:3858348/1‑150 (MQ=255)
                                                                                                              tGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAg                                      <  2:2905707/150‑1 (MQ=255)
                                                                                                              tGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAg                                      <  1:2597429/150‑1 (MQ=255)
                                                                                                                     gggcggTAACTTCTACCCTTACAGTTGCAAGCCGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAAc                               >  2:424858/1‑150 (MQ=255)
                                                                                                                      ggcggTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAAcc                              <  1:3365626/150‑1 (MQ=255)
                                                                                                                         ggTAACTTCTACCCTTACAGTTGCAAGCCGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCAta                           <  2:3868867/150‑2 (MQ=255)
                                                                                                                         ggTAACTTCTACCCTTACAGTTGCAAGCCGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCATc                           >  1:2018018/1‑150 (MQ=255)
                                                                                                                         ggTAACTTCTACCCTTACAGTTGCAAGCCGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCATc                           >  1:2040117/1‑150 (MQ=255)
                                                                                                                         ggTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCATc                           <  1:41171/150‑1 (MQ=255)
                                                                                                                           tAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCATcaa                         <  2:4780662/150‑1 (MQ=255)
                                                                                                                            aaCTTCTACCCTTACAGTTGCAAGCCGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCATcaac                        >  1:3104109/1‑150 (MQ=255)
                                                                                                                               ttCTACCCTTACAGTTGCAAGCCGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCATCAACaaa                     <  2:3009723/150‑1 (MQ=255)
                                                                                                                                tCTACCCTTACAGTTGCAAGCCGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAATGGGGTAACCATCAACAAAg                    >  2:4558544/1‑150 (MQ=255)
                                                                                                                                tCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCATCAACAAAg                    <  1:1753413/150‑1 (MQ=255)
                                                                                                                                  tACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCATCAACAAAGCa                  >  1:2186729/1‑150 (MQ=255)
                                                                                                                                      cTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCATCAACAAAGCACCgg              <  1:4905755/150‑1 (MQ=255)
                                                                                                                                      cTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCATCAACAAAGCACCgg              <  1:318584/150‑1 (MQ=255)
                                                                                                                                       ttACAGTTGCAAGCCGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCATCAACAAAGCACCggg             >  2:4583414/1‑150 (MQ=255)
                                                                                                                                        tACAGTTGCAAGCCGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCATCAACAAAGCACCGGGc            <  2:1865628/150‑1 (MQ=255)
                                                                                                                                        tACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCATCAACAAAGCACCGGGc            <  1:764149/150‑1 (MQ=255)
                                                                                                                                        tACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCATCAACAAAGCACCGGGc            >  2:758590/1‑150 (MQ=255)
                                                                                                                                         aCAGTTGCAAGCCGGAAGATGTCTGGGCGGTGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCATCAACAAAGCACCGGGCg           <  1:3071945/150‑1 (MQ=255)
                                                                                                                                         aCAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCATCAACAAAGCACCGGGCg           >  2:1956989/1‑150 (MQ=255)
                                                                                                                                          cAGTTGCAAGCCGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCATCAACAAAGCACCGGGCga          >  2:471215/1‑150 (MQ=255)
                                                                                                                                           aGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCATCAACAAAGCACCGGGCgat         >  1:2638470/1‑150 (MQ=255)
                                                                                                                                                 aaGCCGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCATCAACAAAGCACCGGGCGATGATGaa   <  2:1261034/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
TGAAAATCAGGATCAGGTGAAATATCAGGCCGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCCGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCATCAACAAAGCACCGGGCGATGATGAAA  >  CP009273/2834945‑2835240

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: