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,7010GT99.7% 937.7 / ‑4.4 285G256G (GGC→GGAsgcXputative endoglucanase with Zn‑dependent exopeptidase domain
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base T (133/151);  total (134/151)
Fisher's exact test for biased strand distribution p-value = 4.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-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 285 total aligned reads displayed.

TTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCCGCGGATATTAAATTCTTCCTGTACCGAAGC  >  CP009273/4520554‑4520847
                                                                                                                                                   |                                                                                                                                                  
  cTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCt                                                                                                                                                <  1:5137039/150‑1 (MQ=255)
  cTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCt                                                                                                                                                <  2:3584193/150‑1 (MQ=255)
     ttATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGAt                                                                                                                                             >  1:353231/1‑150 (MQ=255)
     ttATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGAt                                                                                                                                             <  2:1570455/150‑1 (MQ=255)
     ttATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTGGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGAt                                                                                                                                             >  2:3892852/1‑150 (MQ=255)
      tATCTCGCCGGGAGCCTCTTCTCGCTGCACCGTAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGAtt                                                                                                                                            <  2:421793/150‑1 (MQ=255)
      tATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGAtt                                                                                                                                            <  2:2260107/150‑1 (MQ=255)
           cGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATcc                                                                                                                                       <  1:4247211/150‑1 (MQ=255)
               gggAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGAc                                                                                                                                   >  1:3122357/1‑150 (MQ=255)
               gggAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGAc                                                                                                                                   <  2:1141494/150‑1 (MQ=255)
                 gAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCt                                                                                                                                 <  1:4904571/150‑1 (MQ=255)
                   gCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCg                                                                                                                               <  2:791399/150‑1 (MQ=255)
                     cACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGa                                                                                                                             >  2:1809471/1‑150 (MQ=255)
                     cACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGa                                                                                                                             >  1:269080/1‑150 (MQ=255)
                       cTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGaat                                                                                                                           <  1:2391436/150‑1 (MQ=255)
                       cTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGaat                                                                                                                           >  1:1638447/1‑150 (MQ=255)
                       cTTCTCGCTGCACCAGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGaat                                                                                                                           >  2:4930789/1‑150 (MQ=255)
                         tctcGCTGCACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGaataa                                                                                                                         <  2:1200736/150‑1 (MQ=255)
                          ctcGCTGCACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGaataat                                                                                                                        >  2:3070820/1‑150 (MQ=255)
                           tcGCTGCACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATc                                                                                                                       <  2:5295378/150‑1 (MQ=255)
                           tcGCTGCACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATc                                                                                                                       <  1:3651192/150‑1 (MQ=255)
                            cGCTGCACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCg                                                                                                                      <  1:3773260/150‑1 (MQ=255)
                               tGCACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCgtgt                                                                                                                   <  1:75798/150‑1 (MQ=255)
                               tGCACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATACGTATCAGGCGAGGAGGATTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCTACTCCCTGATTAATCCTGACCTCGTAATAATCgtgt                                                                                                                   >  2:2186074/1‑150 (MQ=255)
                                 cACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTaa                                                                                                                 <  2:1674000/150‑1 (MQ=255)
                                    cGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATc                                                                                                              >  1:1215539/1‑150 (MQ=255)
                                    cGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATc                                                                                                              >  1:2018667/1‑150 (MQ=255)
                                    cGGAAGGTTGTGTTCAAGAGCCGACTGTACCAACATCCGTATCAGGCGAGGAGGCGAGATTAAACCGGCCAGCGGTCCCCGACCCTCGTAGTTCAGGCAGGTGATCCCAACTCCCAGATTAATCCTGACCACGGACTAATCGTGTAAATc                                                                                                              >  2:4772492/1‑150 (MQ=255)
                                     ggAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCa                                                                                                             >  2:2649386/1‑150 (MQ=255)
                                       aaTGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAgg                                                                                                           >  1:3511811/1‑150 (MQ=255)
                                          gttgtGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGgt                                                                                                        >  2:2014425/1‑150 (MQ=255)
                                          gttgtGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGgt                                                                                                        >  2:2018172/1‑150 (MQ=255)
                                             gtgtTCAAGAGCCGTCTGGTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTc                                                                                                     <  2:4384483/150‑1 (MQ=255)
                                               gtTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTcac                                                                                                   <  2:1981666/150‑1 (MQ=255)
                                                  cAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACg                                                                                                <  2:2501543/150‑1 (MQ=255)
                                                       gCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGa                                                                                           <  1:2741822/150‑1 (MQ=255)
                                                          gTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTg                                                                                        >  1:47404/1‑150 (MQ=255)
                                                           tCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGa                                                                                       >  1:2534603/1‑150 (MQ=255)
                                                            cTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGat                                                                                      >  1:5299753/1‑150 (MQ=255)
                                                            cTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGat                                                                                      >  1:4056333/1‑150 (MQ=255)
                                                             tGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGata                                                                                     >  2:4486028/1‑150 (MQ=255)
                                                              gTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGTTGGAGTGatat                                                                                    <  1:1479767/150‑1 (MQ=255)
                                                              gTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGatat                                                                                    <  2:2335122/150‑1 (MQ=255)
                                                                  cAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATTCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAAt                                                                                <  1:1924909/150‑1 (MQ=255)
                                                                  cAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAAt                                                                                <  1:1657685/150‑1 (MQ=255)
                                                                     cATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATAcc                                                                             <  1:374094/150‑1 (MQ=255)
                                                                     cATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATAcc                                                                             >  1:1121067/1‑150 (MQ=255)
                                                                      aTCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCCGGGGTGTCACACGATGGAGTGATATCAATACCa                                                                            <  2:1552031/150‑1 (MQ=255)
                                                                          gTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCg                                                                        >  2:1070543/1‑150 (MQ=255)
                                                                              cAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCcag                                                                    <  2:1254992/150‑1 (MQ=255)
                                                                                ggCTAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATTATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCcaggt                                                                  <  2:2719571/150‑1 (MQ=255)
                                                                                  cGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCcaggtca                                                                >  1:5044203/1‑150 (MQ=255)
                                                                                        ggCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGtttt                                                          <  1:3304622/150‑1 (MQ=255)
                                                                                          cGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTac                                                        <  1:4165097/150‑1 (MQ=255)
                                                                                          cGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTac                                                        >  1:2258147/1‑150 (MQ=255)
                                                                                           gTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTaca                                                       >  2:2877216/1‑150 (MQ=255)
                                                                                           gTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTaca                                                       >  2:4990125/1‑150 (MQ=255)
                                                                                           gTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTaca                                                       >  2:4023907/1‑150 (MQ=255)
                                                                                           gTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTaca                                                       >  1:1275951/1‑150 (MQ=255)
                                                                                            ttaTTAATCCGGCAAGCGTTCCCCGACCATTGTAGTTCAGGCAGGTGAACCCAACTCCCTGATTAATCCTGACCTCGGATTATTCTTGTAAATCAGGGGTGTCACACGATGGAGTTATATCAATACCAATCGCCAGGTCAGGTTTTacac                                                      <  1:4487663/148‑1 (MQ=255)
                                                                                             gATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACg                                                     >  2:1299089/1‑150 (MQ=255)
                                                                                             gATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACg                                                     >  1:68220/1‑150 (MQ=255)
                                                                                              aTTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGACTAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACgg                                                    >  2:4963739/1‑150 (MQ=255)
                                                                                              aTTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACgg                                                    <  2:1904952/150‑1 (MQ=255)
                                                                                              aTTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACgg                                                    >  1:532069/1‑150 (MQ=255)
                                                                                               ttAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGc                                                   >  2:966652/1‑150 (MQ=255)
                                                                                                 aaTCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGt                                                 >  2:3202335/1‑150 (MQ=255)
                                                                                                  aTCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTa                                                <  2:4302166/150‑1 (MQ=255)
                                                                                                     cGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATa                                             <  2:3818147/150‑1 (MQ=255)
                                                                                                     cGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATa                                             >  2:3754846/1‑150 (MQ=255)
                                                                                                      ggCCAGCGGTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATAc                                            <  2:692396/150‑1 (MQ=255)
                                                                                                          aGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGa                                        <  1:2066908/150‑1 (MQ=255)
                                                                                                            cGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAAc                                      >  1:2652844/1‑150 (MQ=255)
                                                                                                               tCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAAGAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAAt                                   >  1:4616881/1‑150 (MQ=255)
                                                                                                                 cccGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGc                                 <  2:2071410/150‑1 (MQ=255)
                                                                                                                  ccGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAGCAATGcc                                >  2:3562221/1‑150 (MQ=255)
                                                                                                                    gACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCcgc                              >  2:1269659/1‑150 (MQ=255)
                                                                                                                         tGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGACCAATGCCGCGGata                         >  1:2189716/1‑150 (MQ=255)
                                                                                                                           gTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCCGCGGATAtt                       <  1:3005139/150‑1 (MQ=255)
                                                                                                                               ttCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCCGCGGATATTAAAt                   <  2:3122357/150‑1 (MQ=255)
                                                                                                                                tCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCCGCGGATATTAAAtt                  <  2:3567121/150‑1 (MQ=255)
                                                                                                                                tCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCCGCGGATATTAAAtt                  <  2:1195848/150‑1 (MQ=255)
                                                                                                                                tCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCCGCGGATATTAAAtt                  >  2:1451941/1‑150 (MQ=255)
                                                                                                                                tCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCCGCGGATATTAAAtt                  <  2:4099460/150‑1 (MQ=255)
                                                                                                                                  aggcaggTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGACCAATGCCGCGGATATTAAAttct                >  2:3031956/1‑150 (MQ=255)
                                                                                                                                  aggcaggTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCCGCGGATATTAAAttct                <  2:3483669/150‑1 (MQ=255)
                                                                                                                                   ggcaggTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCCGCGGATATTAAAttctt               >  1:1702037/1‑150 (MQ=255)
                                                                                                                                   ggcaggTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCCGCGGATATTAAAttctt               <  1:1446158/150‑1 (MQ=255)
                                                                                                                                    gcaggTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCCGCGGATATTAAAttcttc              <  2:2779156/150‑1 (MQ=255)
                                                                                                                                      aggTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCCGCGGATATTAAATTCTTCCt            <  1:3059058/150‑1 (MQ=255)
                                                                                                                                      aggTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCCGCGGATATTAAATTCTTCCt            <  1:1544835/150‑1 (MQ=255)
                                                                                                                                      aggTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCCGCGGATATTAAATTCTTCCt            >  2:1413480/1‑150 (MQ=255)
                                                                                                                                       ggTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCCGCGGATATTAAATTCTTCCTg           >  2:2434804/1‑150 (MQ=255)
                                                                                                                                        gTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCCGCGGATATTAAATTCTTCCTGt          <  1:4182095/150‑1 (MQ=255)
                                                                                                                                          gATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATTGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCCGCGGATATTAAATTCTTCCTGTAc        <  2:2394242/150‑1 (MQ=255)
                                                                                                                                          gATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCCGCGGATATTAAATTCTTCCTGTAc        >  2:3206342/1‑150 (MQ=255)
                                                                                                                                           aTCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCCGCGGATATTAAATTCTTCCTGTAcc       <  2:3647086/150‑1 (MQ=255)
                                                                                                                                           aTCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCCGCGGATATTAAATTCTTCCTGTAcc       <  1:4910405/150‑1 (MQ=255)
                                                                                                                                                aaCTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCCGCGGATATTAAATTCTTCCTGTACCGAAGc  >  2:2507285/1‑150 (MQ=255)
                                                                                                                                                aaCTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCCGCGGATATTAAATTCTTCCTGTACCGAAGc  >  1:1379511/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCCGCGGATATTAAATTCTTCCTGTACCGAAGC  >  CP009273/4520554‑4520847

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: