Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,520,701 G→T 100% G256G (GGC→GGA)  sgcX ← putative endoglucanase with Zn‑dependent exopeptidase domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,520,7010GT98.5% 890.8 / 2.0 272G256G (GGC→GGA) sgcXputative endoglucanase with Zn‑dependent exopeptidase domain
Reads supporting (aligned to +/- strand):  ref base G (3/1);  new base T (136/132);  total (139/133)
Fisher's exact test for biased strand distribution p-value = 6.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01
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 272 total aligned reads displayed.

TTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCCGCGGATATTAAATTCTTCCTGTACCGAAGC  >  CP009273/4520554‑4520847
                                                                                                                                                   |                                                                                                                                                  
ttCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTcc                                                                                                                                                  <  2:4385001/150‑1 (MQ=255)
  cTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCt                                                                                                                                                >  2:1787187/1‑150 (MQ=255)
     ttATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGAt                                                                                                                                             >  2:1570540/1‑150 (MQ=255)
       aTCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTa                                                                                                                                           <  1:2018347/150‑1 (MQ=255)
        tCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTaa                                                                                                                                          <  1:1695010/150‑1 (MQ=255)
          aCGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATc                                                                                                                                        >  1:2245833/1‑150 (MQ=255)
           cGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTACTcc                                                                                                                                       >  1:561789/1‑150 (MQ=255)
                ggAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGAcc                                                                                                                                  <  2:1021877/150‑1 (MQ=255)
                ggAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGAcc                                                                                                                                  <  2:3029334/150‑1 (MQ=255)
                  aGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGTCCTc                                                                                                                                >  1:2626082/1‑150 (MQ=255)
                  aGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTc                                                                                                                                <  1:3530565/150‑1 (MQ=255)
                     cACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGa                                                                                                                             >  1:4336202/1‑150 (MQ=255)
                     cACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGa                                                                                                                             <  1:1401793/150‑1 (MQ=255)
                      aCTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGaa                                                                                                                            <  1:3250924/150‑1 (MQ=255)
                        ttCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGaata                                                                                                                          <  2:1090247/150‑1 (MQ=255)
                        ttCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGaata                                                                                                                          >  1:3603286/1‑150 (MQ=255)
                         tctcGCTGCACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCTTGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGaataa                                                                                                                         <  2:5059626/150‑1 (MQ=255)
                           tcacTGCACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATc                                                                                                                       <  2:3852440/147‑1 (MQ=255)
                           tcGCTGCACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATc                                                                                                                       >  2:340410/1‑150 (MQ=255)
                           tcGCTGCACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATc                                                                                                                       >  2:1942406/1‑150 (MQ=255)
                             gCTGCACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCgt                                                                                                                     <  1:2895161/150‑1 (MQ=255)
                              cTGCACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGAATAATCgtg                                                                                                                    >  1:842975/1‑150 (MQ=255)
                              cTGCACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGAGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGAATATTCgtg                                                                                                                    >  2:2658456/1‑150 (MQ=255)
                                 cACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTaa                                                                                                                 >  1:2582497/1‑150 (MQ=255)
                                 cACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTaa                                                                                                                 >  1:3886728/1‑150 (MQ=255)
                                 cACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGAGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTaa                                                                                                                 >  2:1559550/1‑150 (MQ=255)
                                  aCCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTaaa                                                                                                                >  2:2835705/1‑150 (MQ=255)
                                  aCCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTaaa                                                                                                                >  1:3650991/1‑150 (MQ=255)
                                    cGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATc                                                                                                              >  1:1745719/1‑150 (MQ=255)
                                          gttgtGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGgt                                                                                                        <  2:2510613/150‑1 (MQ=255)
                                             gtgtTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTc                                                                                                     >  1:2167952/1‑150 (MQ=255)
                                               gtTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTcac                                                                                                   <  1:4080587/150‑1 (MQ=255)
                                                 tCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCacac                                                                                                 <  2:52095/150‑1 (MQ=255)
                                                  cAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCacaca                                                                                                <  1:5150535/150‑2 (MQ=255)
                                                  cAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCACAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTTAATCAGGGGTGTCACACg                                                                                                <  1:3458605/150‑1 (MQ=255)
                                                    agagCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGAt                                                                                              >  2:4528169/1‑150 (MQ=255)
                                                    agagCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGAt                                                                                              <  2:4585821/150‑1 (MQ=255)
                                                        ccGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAg                                                                                          >  2:1843123/1‑150 (MQ=255)
                                                        ccGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAg                                                                                          >  2:3653908/1‑150 (MQ=255)
                                                         cGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGt                                                                                         >  1:1077130/1‑150 (MQ=255)
                                                           tCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGa                                                                                       >  2:952223/1‑150 (MQ=255)
                                                             tGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGata                                                                                     >  2:3479218/1‑150 (MQ=255)
                                                             tGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGata                                                                                     <  1:3824107/150‑1 (MQ=255)
                                                              gTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGatat                                                                                    >  1:4266990/1‑150 (MQ=255)
                                                              gTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGatat                                                                                    >  2:2309599/1‑150 (MQ=255)
                                                              gTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGatat                                                                                    <  1:3718476/150‑1 (MQ=255)
                                                                tCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCa                                                                                  >  1:3059638/1‑150 (MQ=255)
                                                                tCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCa                                                                                  >  1:4086963/1‑150 (MQ=255)
                                                                  cAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAAt                                                                                <  2:2770833/150‑1 (MQ=255)
                                                                     cATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATAcc                                                                             <  2:116493/150‑1 (MQ=255)
                                                                      aTCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCa                                                                            <  1:5050801/150‑1 (MQ=255)
                                                                           tATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGc                                                                       >  2:2363681/1‑150 (MQ=255)
                                                                            aTCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGcc                                                                      >  1:468945/1‑150 (MQ=255)
                                                                               aGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCcagg                                                                   <  1:608347/150‑1 (MQ=255)
                                                                                  cGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCcaggtca                                                                <  2:3202134/150‑1 (MQ=255)
                                                                                   gaggagGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCcaggtcag                                                               >  2:62278/1‑150 (MQ=255)
                                                                                     ggaggCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCcaggtcaggt                                                             >  1:2281002/1‑150 (MQ=255)
                                                                                       aggCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGttt                                                           >  2:477786/1‑150 (MQ=255)
                                                                                       aggCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGttt                                                           <  2:3139545/150‑1 (MQ=255)
                                                                                        ggCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGtttt                                                          <  1:1289776/150‑1 (MQ=255)
                                                                                         gCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTa                                                         <  1:4145620/150‑1 (MQ=255)
                                                                                         gCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTa                                                         <  2:4977338/150‑1 (MQ=255)
                                                                                         gCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTa                                                         <  1:4507027/150‑1 (MQ=255)
                                                                                          cGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTac                                                        >  1:666558/1‑150 (MQ=255)
                                                                                           gTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTaca                                                       >  2:469688/1‑150 (MQ=255)
                                                                                             gATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACg                                                     <  1:4967441/150‑1 (MQ=255)
                                                                                             gATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACg                                                     >  1:1001648/1‑150 (MQ=255)
                                                                                                tAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCg                                                  <  1:3318853/150‑1 (MQ=255)
                                                                                                 aaTCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGt                                                 >  2:1908199/1‑150 (MQ=255)
                                                                                                 aaTCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGt                                                 >  2:1269632/1‑150 (MQ=255)
                                                                                                    ccGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAAt                                              <  2:445202/150‑1 (MQ=255)
                                                                                                      ggCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATAc                                            <  1:715439/150‑1 (MQ=255)
                                                                                                         cAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACggg                                         <  2:1705209/150‑1 (MQ=255)
                                                                                                          aGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGa                                        <  2:4984414/150‑1 (MQ=255)
                                                                                                            cGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAAc                                      >  1:2082592/1‑150 (MQ=255)
                                                                                                            cGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAAc                                      <  1:1787187/150‑1 (MQ=255)
                                                                                                              ttCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACaa                                    <  1:2779913/150‑1 (MQ=255)
                                                                                                              ttCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACaa                                    >  2:20897/1‑150 (MQ=255)
                                                                                                                ccccGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATg                                  >  2:1411087/1‑150 (MQ=255)
                                                                                                                 cccGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGc                                 >  2:4037346/1‑150 (MQ=255)
                                                                                                                      ccATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGACCACTGCCGCgg                            >  2:717853/1‑150 (MQ=255)
                                                                                                                       cATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCCGCGGa                           >  2:1808837/1‑150 (MQ=255)
                                                                                                                         tGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCCGCGGata                         <  2:5125030/150‑1 (MQ=255)
                                                                                                                         tGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCCGCGGata                         <  2:2934629/150‑1 (MQ=255)
                                                                                                                          ggTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCCGCGGatat                        <  1:2298534/150‑1 (MQ=255)
                                                                                                                           gTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCCGCGGATAtt                       >  1:2916763/1‑150 (MQ=255)
                                                                                                                            tAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCCGCGGATATTa                      >  1:4571699/1‑150 (MQ=255)
                                                                                                                            tAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCCGCGGATATTa                      <  2:3104971/150‑1 (MQ=255)
                                                                                                                             aGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCCGCGGATATTaa                     >  2:4841840/1‑150 (MQ=255)
                                                                                                                                tCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCCGCGGATATTAAAtt                  >  2:39585/1‑150 (MQ=255)
                                                                                                                                 caggcaggTGATCCCAACTCCCTGATTAATCCTGACATCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCCGCGGATATTAAAttc                 >  1:2740207/1‑150 (MQ=255)
                                                                                                                                  aggcaggTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCCGCGGATATTAAAttct                >  1:613523/1‑150 (MQ=255)
                                                                                                                                   ggcaggTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCCGCGGATATTAAAttctt               <  1:1190512/150‑1 (MQ=255)
                                                                                                                                    gcaggTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCCGCGGATATTAAAttcttc              <  1:3122916/150‑1 (MQ=255)
                                                                                                                                     caggTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCCGCGGATATTAAATTCTTcc             >  1:1483828/1‑150 (MQ=255)
                                                                                                                                      aggTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCCGCGGATATTAAATTCTTCCt            >  1:1675183/1‑150 (MQ=255)
                                                                                                                                        gTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCCGCGGATATTAAATTCTTCCTGt          >  2:4804206/1‑150 (MQ=255)
                                                                                                                                          gATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCCGCGGATATTAAATTCTTCCTGTAc        >  2:2574477/1‑150 (MQ=255)
                                                                                                                                            tCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCCGCGGATATTAAATTCTTCCTGTACCg      <  1:952484/150‑1 (MQ=255)
                                                                                                                                                aaCTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCCGCGGATATTAAATTCTTCCTGTACCGAAGc  >  1:30261/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCCGCGGATATTAAATTCTTCCTGTACCGAAGC  >  CP009273/4520554‑4520847

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: ‑