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,0940CA99.1% 737.0 / ‑1.9 230P249Q (CCG→CAG) ascBcryptic 6‑phospho‑beta‑glucosidase
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base A (128/100);  total (129/101)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-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 230 total aligned reads displayed.

AAAATCAGGATCAGGTGAAATATCAGGCCGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCCGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCATCAACAAAGCACCGGGCGATGATGAA  >  CP009273/2834947‑2835239
                                                                                                                                                   |                                                                                                                                                 
aaaaTCAGGATCAGGTGAAATATCAGGCCGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAgg                                                                                                                                                 <  1:2630027/150‑1 (MQ=255)
   atcaggatcaggTGAAATATCAGGCCGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAg                                                                                                                                              >  2:888648/1‑150 (MQ=255)
         atcaggTGAAAGATCAGGCCGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCt                                                                                                                                        <  2:2184138/150‑1 (MQ=255)
          tcaggTGAAATATCAGGCCGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACCGTTGCAAGCAGGAAGATGTCTg                                                                                                                                       <  2:2936948/150‑1 (MQ=255)
           caggTGAAATATCAGGCCGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTgg                                                                                                                                      >  2:4899903/1‑150 (MQ=255)
            aggTGAAATATCAGGCCGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTggg                                                                                                                                     <  2:3422080/150‑1 (MQ=255)
             ggTGAAATATCAGGCCGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGggc                                                                                                                                    <  1:4216209/150‑1 (MQ=255)
               tGAAATATCAGGCCGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCTGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGggcgg                                                                                                                                  >  1:1212079/1‑150 (MQ=255)
                 aaaTATCAGGCCGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGgcggcg                                                                                                                                >  1:1630684/1‑150 (MQ=255)
                  aaTATCAGGCCGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGgcgc                                                                                                                               >  1:3036509/1‑150 (MQ=255)
                   atatCAGGCCGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCt                                                                                                                              <  1:754856/150‑1 (MQ=255)
                   atatCAGGCCGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCt                                                                                                                              <  1:3979426/150‑1 (MQ=255)
                    tatCAGGCCGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTg                                                                                                                             >  1:3907705/1‑150 (MQ=255)
                     atCAGGCCGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTgg                                                                                                                            >  2:3295410/1‑150 (MQ=255)
                      tCAGGCCGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGga                                                                                                                           >  2:4430558/1‑150 (MQ=255)
                           ccGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAg                                                                                                                      <  2:5222095/150‑1 (MQ=255)
                           ccGCGCATCACCAGCTGGGTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAg                                                                                                                      <  2:821254/150‑1 (MQ=255)
                             gcgcATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGAc                                                                                                                    >  1:2381645/1‑150 (MQ=255)
                              cgcATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGAcc                                                                                                                   <  1:2941877/150‑1 (MQ=255)
                               gcATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCg                                                                                                                  >  1:4382618/1‑150 (MQ=255)
                                       aGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAAcc                                                                                                          <  1:1408509/150‑1 (MQ=255)
                                           ggTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGACAACCTGtt                                                                                                      >  1:1209804/1‑150 (MQ=255)
                                           ggTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGtt                                                                                                      <  2:803593/150‑1 (MQ=255)
                                           ggTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGtt                                                                                                      >  2:653180/1‑150 (MQ=255)
                                           ggTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGtt                                                                                                      >  2:2973586/1‑150 (MQ=255)
                                            gTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAACTCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGttt                                                                                                     <  2:1148938/150‑1 (MQ=255)
                                            gTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGttt                                                                                                     >  1:591578/1‑150 (MQ=255)
                                               gCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGtttttt                                                                                                  >  1:662993/1‑150 (MQ=255)
                                                   gTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCCGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCg                                                                                              <  1:312702/150‑1 (MQ=255)
                                                   gTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCg                                                                                              <  1:5237836/150‑1 (MQ=255)
                                                   gTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCg                                                                                              <  1:4787641/150‑1 (MQ=255)
                                                    tGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGa                                                                                             >  2:4092572/1‑150 (MQ=255)
                                                     gcgcTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGAt                                                                                            <  1:2869295/150‑1 (MQ=255)
                                                      cgcTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCGTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAGGACCGGGAAAACCTGTTTTTTATCGAtg                                                                                           >  2:4997034/1‑150 (MQ=255)
                                                       gcTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAACACCTGTTTTTTATCGAtgt                                                                                          >  2:4296368/1‑150 (MQ=255)
                                                          aGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCa                                                                                       <  1:41223/150‑1 (MQ=255)
                                                           gCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAg                                                                                      >  2:4849240/1‑150 (MQ=255)
                                                             caccaAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGgc                                                                                    >  1:2529204/1‑150 (MQ=255)
                                                             caccaAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGgc                                                                                    >  2:3254472/1‑150 (MQ=255)
                                                              accaAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGgcg                                                                                   <  2:2740802/150‑1 (MQ=255)
                                                                caAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGcgcg                                                                                 >  2:2800727/1‑150 (MQ=255)
                                                                  aaaTCGCCCATGAGGCTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCggg                                                                               <  1:4426129/150‑1 (MQ=255)
                                                                   aaTCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCgggg                                                                              <  2:4040699/150‑1 (MQ=255)
                                                                       gCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACg                                                                          >  2:3471771/1‑150 (MQ=255)
                                                                        cccATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGATTTTTATCGATGTGCAGGCGCGGGGCACGt                                                                         >  2:4565147/1‑150 (MQ=255)
                                                                         ccATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTa                                                                        >  1:101092/1‑150 (MQ=255)
                                                                           aTGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGtctc                                                                      >  1:1871859/1‑150 (MQ=255)
                                                                           aTGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATc                                                                      <  1:59940/150‑1 (MQ=255)
                                                                           aTGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATc                                                                      >  2:3807515/1‑150 (MQ=255)
                                                                            tgtggTTATCCCGCAATATGAGGTGGGCTGTATGCTGGCGGGCGGTTACTTGTACCCTTACTGGTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACGGGGATAACCTGTTGTTTATCGATGTGCAGGCGCGGGGCACGTTTcc                                                                     <  2:1999570/147‑1 (MQ=255)
                                                                             gAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCg                                                                    >  1:3215892/1‑150 (MQ=255)
                                                                                   aaCCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTAc                                                              <  1:2550841/150‑1 (MQ=255)
                                                                                    aCCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTAct                                                             <  2:23716/150‑1 (MQ=255)
                                                                                    aCCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTAct                                                             >  1:667619/1‑150 (MQ=255)
                                                                                     cccGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTActc                                                            <  1:4865131/150‑1 (MQ=255)
                                                                                      ccGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTActct                                                           >  2:1424159/1‑150 (MQ=255)
                                                                                      ccGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTActct                                                           >  1:4370356/1‑150 (MQ=255)
                                                                                          aaaaTCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGccc                                                       >  2:3217584/1‑150 (MQ=255)
                                                                                               cAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTa                                                  <  2:466452/150‑1 (MQ=255)
                                                                                               cAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTa                                                  >  2:2682780/1‑150 (MQ=255)
                                                                                                 ggTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTAtt                                                <  2:2129052/150‑1 (MQ=255)
                                                                                                 ggTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTAtt                                                >  2:5059653/1‑150 (MQ=255)
                                                                                                  gTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTc                                               <  1:73490/150‑1 (MQ=255)
                                                                                                   tGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTcc                                              <  2:4265180/150‑1 (MQ=255)
                                                                                                   tGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTcc                                              <  1:1724438/150‑1 (MQ=255)
                                                                                                   tGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTcc                                              <  1:2220438/150‑1 (MQ=255)
                                                                                                    gggCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCcg                                             >  1:3859275/1‑150 (MQ=255)
                                                                                                     ggCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCcgc                                            >  1:4101215/1‑150 (MQ=255)
                                                                                                       cTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGa                                          >  2:4344403/1‑150 (MQ=255)
                                                                                                            tGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAg                                     >  1:2103860/1‑150 (MQ=255)
                                                                                                            tGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCCTACTCTGCCCGCGTATTCCGCGAAAAAg                                     <  1:1960748/150‑1 (MQ=255)
                                                                                                             gCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAgg                                    <  1:4894149/150‑1 (MQ=255)
                                                                                                                ggcgggcgGTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGt                                 <  1:2445314/150‑1 (MQ=255)
                                                                                                                 gcgggcggTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTa                                <  1:329992/150‑1 (MQ=255)
                                                                                                                 gcgggcggTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTa                                <  1:4996652/150‑1 (MQ=255)
                                                                                                                  cgggcggTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTaa                               <  2:5245692/150‑1 (MQ=255)
                                                                                                                     gcggTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCa                            >  2:1897145/1‑150 (MQ=255)
                                                                                                                        gTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCATca                         <  2:1773560/150‑1 (MQ=255)
                                                                                                                        gTAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCATca                         >  1:956224/1‑150 (MQ=255)
                                                                                                                         tAACTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCATcaa                        >  2:2519051/1‑150 (MQ=255)
                                                                                                                         tAACTTCGACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCATcaa                        <  1:352418/150‑1 (MQ=255)
                                                                                                                          aaCTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCATcaac                       <  2:1842015/150‑1 (MQ=255)
                                                                                                                          aaCTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCATcaac                       >  2:1321363/1‑150 (MQ=255)
                                                                                                                          aaCTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCATcaac                       >  1:3706154/1‑150 (MQ=255)
                                                                                                                           aCTTCTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCATcaaca                      <  2:3182349/150‑1 (MQ=255)
                                                                                                                               cTACCCTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGCAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGACAAAGGGGTAACCATCAACAAAGc                  >  2:1462097/1‑150 (MQ=255)
                                                                                                                                 aCCCTTACAGTTGCAAGAAGGAAGATGTCTGGTCGGCGCTGGAGAAAGACCTTGAAAACCTGTTTTTTATCGATTTGCAGGCGCGTGGCACTTATCCGGCTTACTCTGCACGCGTATTCCGCGAAAAAGGGGTAACCATCAACAAAGCAc                >  1:2382525/1‑150 (MQ=255)
                                                                                                                                   ccTTACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCATCAACAAAGCACCg              >  1:5113346/1‑150 (MQ=255)
                                                                                                                                     ttACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCATCAACAAAGCACCggg            <  1:2798042/150‑1 (MQ=255)
                                                                                                                                     ttACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCATCAACAAAGCACCggg            <  1:4430558/150‑1 (MQ=255)
                                                                                                                                     ttACAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCATCAACAAAGCACCggg            <  1:619560/150‑1 (MQ=255)
                                                                                                                                       aCAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCATCAACAAAGCACCGGGCg          <  1:5317547/150‑1 (MQ=255)
                                                                                                                                        cAGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCATCAACAAAGCACCGGGCga         <  1:2439693/150‑1 (MQ=255)
                                                                                                                                         aGTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCATCAACAAAGCACCGGGCgat        <  2:1450592/150‑1 (MQ=255)
                                                                                                                                          gTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCATCAACAAAGCACCGGGCgatg       <  1:4415987/150‑1 (MQ=255)
                                                                                                                                          gTTGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCATCAACAAAGCACCGGGCgatg       <  1:1758459/150‑1 (MQ=255)
                                                                                                                                            tGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCATCAACAAAGCACCGGGCgatgat     <  2:373115/150‑1 (MQ=255)
                                                                                                                                            tGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCATCAACAAAGCACCGGGCgatgat     <  1:2048646/150‑1 (MQ=255)
                                                                                                                                            tGCAAGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCATCAACAAAGCACCGGGCgatgat     >  1:869668/1‑150 (MQ=255)
                                                                                                                                               aaGCAGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCATCAACAAAGCACCGGGCGATGATGaa  >  1:3787437/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
AAAATCAGGATCAGGTGAAATATCAGGCCGCGCATCACCAGCTGGTTGCCAGTGCGCTAGCCACCAAAATCGCCCATGAGGTTAACCCGCAAAATCAGGTGGGCTGTATGCTGGCGGGCGGTAACTTCTACCCTTACAGTTGCAAGCCGGAAGATGTCTGGGCGGCGCTGGAGAAAGACCGGGAAAACCTGTTTTTTATCGATGTGCAGGCGCGGGGCACGTATCCGGCTTACTCTGCCCGCGTATTCCGCGAAAAAGGGGTAACCATCAACAAAGCACCGGGCGATGATGAA  >  CP009273/2834947‑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: