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.0% 788.8 / 5.5 251G256G (GGC→GGAsgcXputative endoglucanase with Zn‑dependent exopeptidase domain
Reads supporting (aligned to +/- strand):  ref base G (3/2);  new base T (100/146);  total (103/148)
Fisher's exact test for biased strand distribution p-value = 4.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 251 total aligned reads displayed.

TCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCCGCGGATATTAAATTCTTCCTGTACCGAAG  >  CP009273/4520555‑4520846
                                                                                                                                                  |                                                                                                                                                 
tCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTcc                                                                                                                                                 <  1:1572615/149‑1 (MQ=255)
 cTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTccc                                                                                                                                                >  2:1595611/1‑149 (MQ=255)
  tGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCt                                                                                                                                               >  2:3728139/1‑149 (MQ=255)
     tATCACGCCTGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCTAGAGCAGTCTGGTCCAACATCCGTATCAGGAGAGGAGGCGTGATTAATCCGGCCAGCGTTCACCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGAt                                                                                                                                            <  1:1490852/149‑1 (MQ=255)
     tATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGAt                                                                                                                                            <  2:341577/149‑1 (MQ=255)
     tATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGAt                                                                                                                                            <  1:3537451/149‑1 (MQ=255)
     tATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGAt                                                                                                                                            <  2:3031690/149‑1 (MQ=255)
       tCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTa                                                                                                                                          >  1:2501507/1‑149 (MQ=255)
        cACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTaa                                                                                                                                         >  2:3456694/1‑149 (MQ=255)
        cACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTaa                                                                                                                                         >  2:1800695/1‑149 (MQ=255)
        cACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTaa                                                                                                                                         <  1:1090737/149‑1 (MQ=255)
         aCGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAAt                                                                                                                                        <  2:2233413/149‑1 (MQ=255)
           gCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATcc                                                                                                                                      <  2:3753487/149‑1 (MQ=255)
             cGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGTCAGGTGATCCCAACTCCCTGATTAAACCTg                                                                                                                                    >  1:883360/1‑149 (MQ=255)
              gggAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGa                                                                                                                                   <  1:751184/149‑1 (MQ=255)
              gggAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGa                                                                                                                                   >  2:566083/1‑149 (MQ=255)
                 aGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCt                                                                                                                                <  1:2111862/149‑1 (MQ=255)
                  gCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTc                                                                                                                               >  1:2789383/1‑149 (MQ=255)
                     aCTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGa                                                                                                                            >  2:1408311/1‑149 (MQ=255)
                       ttCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGaat                                                                                                                          >  1:950134/1‑149 (MQ=255)
                         ctcGCTGCACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGaataa                                                                                                                        <  2:2925994/149‑1 (MQ=255)
                           cGCTGCACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATc                                                                                                                      <  2:3402941/149‑1 (MQ=255)
                             cTGCACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCgt                                                                                                                    <  1:1245854/149‑1 (MQ=255)
                             cTGCACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCgt                                                                                                                    <  2:1610561/149‑1 (MQ=255)
                             cTGCACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCgt                                                                                                                    <  1:3568673/149‑1 (MQ=255)
                              tGCACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCgtg                                                                                                                   <  2:2537663/149‑1 (MQ=255)
                                   cGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAAt                                                                                                              <  2:1080510/149‑1 (MQ=255)
                                    ggAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATc                                                                                                             <  2:2196834/149‑1 (MQ=255)
                                     gAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCa                                                                                                            <  1:3924911/149‑1 (MQ=255)
                                      aaTGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAg                                                                                                           <  1:174544/149‑1 (MQ=255)
                                      aaTGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAg                                                                                                           >  2:3067680/1‑149 (MQ=255)
                                        tgttgtGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAggg                                                                                                         >  2:1146337/1‑149 (MQ=255)
                                        tgttgtGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAggg                                                                                                         <  2:2613598/149‑1 (MQ=255)
                                         gttgtGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAgggg                                                                                                        <  1:646352/149‑1 (MQ=255)
                                           tgtgTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGgtg                                                                                                      <  1:4014740/149‑1 (MQ=255)
                                               ttCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTcac                                                                                                  <  1:362222/149‑1 (MQ=255)
                                                   agagCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGa                                                                                              <  1:1135755/149‑1 (MQ=255)
                                                       ccGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGa                                                                                          <  2:2828481/149‑1 (MQ=255)
                                                        cGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAg                                                                                         >  2:3841748/1‑149 (MQ=255)
                                                          tCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTg                                                                                       >  1:2132913/1‑149 (MQ=255)
                                                              ttCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGatat                                                                                   <  2:386943/149‑1 (MQ=255)
                                                                 cAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCaa                                                                                <  2:2266772/149‑1 (MQ=255)
                                                                 cAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCaa                                                                                <  2:3914073/149‑1 (MQ=255)
                                                                  aaCATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAAt                                                                               >  2:2824943/1‑149 (MQ=255)
                                                                   aCATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATa                                                                              >  1:2713066/1‑149 (MQ=255)
                                                                      tCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCACTACCa                                                                           >  2:2391943/1‑149 (MQ=255)
                                                                       ccGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCaa                                                                          <  2:3172061/149‑1 (MQ=255)
                                                                        cGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAAt                                                                         <  2:4056204/149‑1 (MQ=255)
                                                                        cGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAAt                                                                         <  2:1310780/149‑1 (MQ=255)
                                                                         gTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATc                                                                        <  2:848259/149‑1 (MQ=255)
                                                                         gTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATc                                                                        <  1:2969087/149‑1 (MQ=255)
                                                                         gTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATAAATACCAATc                                                                        <  1:2251695/149‑1 (MQ=255)
                                                                         gTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATc                                                                        <  2:2458565/149‑1 (MQ=255)
                                                                            tCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGcc                                                                     <  1:3004526/149‑1 (MQ=255)
                                                                             cAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCca                                                                    >  2:1982205/1‑149 (MQ=255)
                                                                              aGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCcag                                                                   >  2:1287165/1‑149 (MQ=255)
                                                                              aGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCcag                                                                   <  2:2131390/149‑1 (MQ=255)
                                                                              aGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCcag                                                                   <  2:876461/149‑1 (MQ=255)
                                                                                  gaggagGCGTGATTAATCCGGCCAGCGTTCCCCGACCTTGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCcaggtca                                                               <  2:1042648/149‑1 (MQ=255)
                                                                                   aggaggCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCcaggtcag                                                              >  1:1083613/1‑149 (MQ=255)
                                                                                      aggCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGtt                                                           <  1:271165/149‑1 (MQ=255)
                                                                                       ggCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGttt                                                          <  1:1800695/149‑1 (MQ=255)
                                                                                       ggCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGttt                                                          >  2:859653/1‑149 (MQ=255)
                                                                                       ggCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGttt                                                          <  2:935082/149‑1 (MQ=255)
                                                                                       ggCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGttt                                                          <  1:1517522/149‑1 (MQ=255)
                                                                                        gCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGtttt                                                         <  2:2274573/149‑1 (MQ=255)
                                                                                        gCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGtttt                                                         >  2:756574/1‑149 (MQ=255)
                                                                                         cGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTa                                                        <  1:2766463/149‑1 (MQ=255)
                                                                                           tGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGACAGGTCAGGTTTTaca                                                      >  2:3302136/1‑149 (MQ=255)
                                                                                            gATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTacac                                                     >  2:3633/1‑149 (MQ=255)
                                                                                            gATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTacac                                                     >  2:3602/1‑149 (MQ=255)
                                                                                            gATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTacac                                                     <  2:2724483/149‑1 (MQ=255)
                                                                                                    cGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAAt                                             >  1:1676378/1‑149 (MQ=255)
                                                                                                     ggCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATa                                            <  1:481957/149‑1 (MQ=255)
                                                                                                      gCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATAc                                           >  2:2094438/1‑149 (MQ=255)
                                                                                                      gCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATAc                                           >  2:662511/1‑149 (MQ=255)
                                                                                                        cAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACgg                                         <  1:3728139/149‑1 (MQ=255)
                                                                                                        cAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACgg                                         >  1:3065785/1‑149 (MQ=255)
                                                                                                          gCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGa                                       <  1:1144151/149‑1 (MQ=255)
                                                                                                          gCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGa                                       <  1:3843968/149‑1 (MQ=255)
                                                                                                             ttCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACa                                    >  1:1060741/1‑149 (MQ=255)
                                                                                                              tCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGACCaa                                   >  2:3163665/1‑149 (MQ=255)
                                                                                                                 ccGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGc                                >  2:1001680/1‑149 (MQ=255)
                                                                                                                 ccGACCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGc                                >  2:1991187/1‑149 (MQ=255)
                                                                                                                    aCCATGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCcgc                             <  1:1001533/149‑1 (MQ=255)
                                                                                                                       aTGGTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCCGCGGa                          <  2:4013883/149‑1 (MQ=255)
                                                                                                                         ggTAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCCGCGGata                        <  2:1852288/149‑1 (MQ=255)
                                                                                                                           tAGTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCCGCGGATAtt                      >  1:2194119/1‑149 (MQ=255)
                                                                                                                             gTTCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCCGCGGATATTaa                    <  2:1488198/149‑1 (MQ=255)
                                                                                                                              ttCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCCGCGGATATTaaa                   <  2:3327129/149‑1 (MQ=255)
                                                                                                                               tCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCCGCGGATATTAAAt                  >  2:2682660/1‑149 (MQ=255)
                                                                                                                               tCAGGCAGGTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCCGCGGATATTAAAt                  <  1:1375602/149‑1 (MQ=255)
                                                                                                                                 aggcaggTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCCGCGGATATTAAAttc                <  1:3470569/149‑1 (MQ=255)
                                                                                                                                     aggTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCCGCGGATATTAAATTCTTcc            <  2:3000747/149‑1 (MQ=255)
                                                                                                                                      ggTGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCCGCGGATATTAAATTCTTCCt           <  1:1410951/149‑1 (MQ=255)
                                                                                                                                        tGATCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCCGCGGATATTAAATTCTTCCTGt         >  1:1827041/1‑149 (MQ=255)
                                                                                                                                          aTCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCCGCGGATATTAAATTCTTCCTGTAc       <  2:2279005/149‑1 (MQ=255)
                                                                                                                                          aTCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCCGCGGATATTAAATTCTTCCTGTAc       <  1:3137538/149‑1 (MQ=255)
                                                                                                                                          aTCCCAACTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCCGCGGATATTAAATTCTTCCTGTAc       >  1:2045149/1‑149 (MQ=255)
                                                                                                                                               aaCTCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCCGCGGATATTAAATTCTTCCTGTACCGAAg  <  2:2264267/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
TCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCCGCGGATATTAAATTCTTCCTGTACCGAAG  >  CP009273/4520555‑4520846

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: